Towards the construction of a validated numerical system to study the mesoscale dynamics of the North East Atlantic (2003-2006)

Máster en Oceanografía

External drivers of biogeochemical cycles in a tropical montane forest in Ecuador

Successful conservation of tropical montane forest, one of the most threatened ecosystems on earth, requires detailed knowledge of its biogeochemistry. Of particular interest is the response of the biogeochemical element cycles to external influences such as element deposition or climate change. Therefore the overall objective of my study was to contribute to improved understanding of role and functioning of the Andean tropical montane forest. In detail, my objectives were to determine (1) the role of long-range transported aerosols and their transport mechanisms, and (2) the role of short-term extreme climatic events for the element budget of Andean tropical forest. In a whole-catchment approach including three 8-13 ha microcatchments under tropical montane forest on the east-exposed slope of the eastern cordillera in the south Ecuadorian Andes at 1850-2200 m above sea level I monitored at least in weekly resolution the concentrations and fluxes of Ca, Mg, Na, K, NO3-N, NH4-N, DON, P, S, TOC, Mn, and Al in bulk deposition, throughfall, litter leachate, soil solution at the 0.15 and 0.3 m depths, and runoff between May 1998 and April 2003. I also used meteorological data from my study area collected by cooperating researchers and the Brazilian meteorological service (INPE), as well as remote sensing products of the North American and European space agencies NASA and ESA. My results show that (1) there was a strong interannual variation in deposition of Ca [4.4-29 kg ha-1 a-1], Mg [1.6-12], and K [9.8-30]) between 1998 and 2003. High deposition changed the Ca and Mg budgets of the catchments from loss to retention, suggesting that the additionally available Ca and Mg was used by the ecosystem. Increased base metal deposition was related to dust outbursts of the Sahara and an Amazonian precipitation pattern with trans-regional dry spells allowing for dust transport to the Andes. The increased base metal deposition coincided with a strong La Niña event in 1999/2000. There were also significantly elevated H+, N, and Mn depositions during the annual biomass burning period in the Amazon basin. Elevated H+ deposition during the biomass burning period caused elevated base metal loss from the canopy and the organic horizon and deteriorated already low base metal supply of the vegetation. Nitrogen was only retained during biomass burning but not during non-fire conditions when deposition was much smaller. Therefore biomass burning-related aerosol emissions in Amazonia seem large enough to substantially increase element deposition at the western rim of Amazonia. Particularly the related increase of acid deposition impoverishes already base-metal scarce ecosystems. As biomass burning is most intense during El Niño situations, a shortened ENSO cycle because of global warming likely enhances the acid deposition at my study forest. (2) Storm events causing near-surface water flow through C- and nutrient-rich topsoil during rainstorms were the major export pathway for C, N, Al, and Mn (contributing >50% to the total export of these elements). Near-surface flow also accounted for one third of total base metal export. This demonstrates that storm-event related near-surface flow markedly affects the cycling of many nutrients in steep tropical montane forests. Changes in the rainfall regime possibly associated with global climate change will therefore also change element export from the study forest. Element budgets of Andean tropical montane rain forest proved to be markedly affected by long-range transport of Saharan dust, biomass burning-related aerosols, or strong rainfalls during storm events. Thus, increased acid and nutrient deposition and the global climate change probably drive the tropical montane forest to another state with unknown consequences for its functions and biological diversity.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Development of an In Situ Measurement Device for Obtaining Material Thermal Properties

This thesis presents a methodology for measuring thermal properties in situ, with a special focus on obtaining properties of layered stack-ups commonly used in armored vehicle components. The technique involves attaching a thermal source to the surface of a component, measuring the heat flux transferred between the source and the component, and measuring the surface temperature response. The material properties of the component can subsequently be determined from measurement of the transient heat flux and temperature response at the surface alone. Experiments involving multilayered specimens show that the surface temperature response to a sinusoidal heat flux forcing function is also sinusoidal. A frequency domain analysis shows that sinusoidal thermal excitation produces a gain and phase shift behavior typical of linear systems. Additionally, this analysis shows that the material properties of sub-surface layers affect the frequency response function at the surface of a particular stack-up. The methodology involves coupling a thermal simulation tool with an optimization algorithm to determine the material properties from temperature and heat flux measurement data. Use of a sinusoidal forcing function not only provides a mechanism to perform the frequency domain analysis described above, but sinusoids also have the practical benefit of reducing the need for instrumentation of the backside of the component. Heat losses can be minimized by alternately injecting and extracting heat on the front surface, as long as sufficiently high frequencies are used.

A Preliminary Study of the Micro-Fauna of Certain Paleozoic Formations of Montana

The study of the micro-fauna of Montana formations has been almost entirely neglected. Because the petroleum industry of this state has not felt the necessity for using micro-paleontology in its sub-surface correlations, the science has been but little used. The Montana Power Company has had an examination made of some of its well cuttings by a competent micro-paleontologist who found some foraminifera in Mesozoic sediments. However, no investiga­tions have been made to determine the presence and character of the micro-fauna of the Paleozoic formations of Montana.

Insoluble Residues of the Lower Mississippian Limestones of the Madison Group

The correlation of non-fossiliferous drill samples is one of the difficult problems that is en­countered in sub-surface stratigraphy. In order to truly correlate a formation, it must have some dis tinctive features and have an areal persistence of these features. These requirements are probably met best by limestone.

Seasonal Climatology of Hydrographic Conditions in the Upwelling Region Off Northern Chile

Over 30 years of hydrographic data from the northern Chile (18 degreesS-24 degreesS) upwelling region are used to calculate the surface and subsurface seasonal climatology extending 400 km offshore. The data are interpolated to a grid with sufficient spatial resolution to preserve cross-shelf gradients and then presented as means within four seasons: austral winter (July-September), spring (October-December), summer (January-March), and fall (April-June). Climatological monthly wind forcing, surface temperature, and sea level from three coastal stations indicate equatorward (upwelling favorable) winds throughout the year, weakest in the north. Seasonal maximum alongshore wind stress is in late spring and summer (December-March). Major water masses of the region are identified in climatological T-S plots and their sources and implied circulation discussed. Surface fields and vertical transects of temperature and salinity confirm that upwelling occurs year-round, strongest in summer and weakest in winter, bringing relatively fresh water to the surface nearshore. Surface geostrophic flow nearshore is equatorward throughout the year. During summer, an anticyclonic circulation feature in the north which extends to at least 200 m depth is evident in geopotential anomaly and in both temperature and geopotential variance fields. Subsurface fields indicate generally poleward flow throughout the year, strongest in an undercurrent near the coast. This undercurrent is strongest in summer and most persistent and organized in the south (south of 21 degreesS), A subsurface oxygen minimum, centered at similar to 250 m, is strongest at lower latitudes. Low-salinity subsurface water intrudes into the study area near 100 m, predominantly in offshore regions, strongest during summer and fall and in the southernmost portion of the region. The climatological fields are compared to features off Baja within the somewhat analogous California Current and to measurements from higher latitudes within the Chile-Peru Current system.

Plutonium and americium concentrations in sea water and in sinking particles

Vertical fluxes of 239+240Pu and 241Am and temporal changes in their inventories in the northwestern Mediterranean Sea have been examined through high-resolution water column sampling coupled with direct measurements of the vertical flux of particle-bound transuranics using time-series sediment traps. Water column profiles of both radionuclides showed well-defined sub-surface maxima (2391240Pu between 100-400 m; 241Am at 100-200 m and 800 m), the depths of which are a result of the different biogeochemical scavenging behavior of the two radionuclides. Comparison of deep water column (0-2,000 m) transuranic inventories with those derived from earlier measurements demonstrate that the total 2391240Pu inventory had not substantially changed between 1976-1990 whereas 241Am had decreased by approximately 24%. Enhanced scavenging of 241Am and a resultant, more rapid removal from the water column relative to 239+240Pu was also supported by the observation of elevated Am/Pu activity ratios in sinking particles collected in sediment traps at depth. Direct measurements of the downward flux of particulate 239+240Pu and 241Am compared with transuranic removal rates derived from observed total water column inventory differences over time, show that particles sinking out of deep waters (1,000-2,000 m) could account for 26-72% of the computed total annual 239+240Pu loss and virtually all of the 241Am removal from the water column. Upper water column (0-200 m) residence times based on direct flux measurements ranged from 20-30 yr for 239+240Pu and 5-10 yr for 241Am. The observation that 241Am/239+240Pu activity ratios in unfiltered Mediterranean seawater are six times lower than those in the north Pacific suggests the existence of a specific mechanism for enhanced scavenging and removal of 241Am from the generally oligotrophic waters of the open Mediterranean. It is proposed that atmospheric inputs of aluminosilicate particles transported by Saharan dust events which frequently occur in the Mediterranean region could enhance the geochemical scavenging and resultant removal of 241Am to the sediments.

(Table 1.4) Ratios of petrogenic components in glassy and crystalline basalts from the north part of the Bougault anomaly, Mid-Atlantic Ridge 15°N

Structure and composition of sub-surface bottom sediments from the southwest Barents Sea have been under study. The study has revealed heterogeneity of sediment structure resulted from temporal irregularity and variability of sedimentation processes. The study of the heavy minerals from 0.1-0.01 mm grain size fraction has shown prevalence of green hornblende, epidote, garnet, and ilmenite in all types of sediments; these minerals are the basis of terrigenous-mineralogical province. At the same time in different areas local terrigenous-mineralogical associations have been identified. Clay mineral composition of in the sediments was quite uniform: biotite, chlorite, hydromica, smectite. Despite this, a number of features indicating initial stages of clay mineral transformation has been identified. Differences in material composition and structure of the studied sediments are associated with rapid change in paleogeographic situation on the land - ice cover melting on the Kola Peninsula and subsequent Holocene climatic situation.

(Table 1.5) Contents of some volatile components in glassy MORB from the north part of the Bougault anomaly, Mid-Atlantic Ridge 15°N

Structure and composition of sub-surface bottom sediments from the southwest Barents Sea have been under study. The study has revealed heterogeneity of sediment structure resulted from temporal irregularity and variability of sedimentation processes. The study of the heavy minerals from 0.1-0.01 mm grain size fraction has shown prevalence of green hornblende, epidote, garnet, and ilmenite in all types of sediments; these minerals are the basis of terrigenous-mineralogical province. At the same time in different areas local terrigenous-mineralogical associations have been identified. Clay mineral composition of in the sediments was quite uniform: biotite, chlorite, hydromica, smectite. Despite this, a number of features indicating initial stages of clay mineral transformation has been identified. Differences in material composition and structure of the studied sediments are associated with rapid change in paleogeographic situation on the land - ice cover melting on the Kola Peninsula and subsequent Holocene climatic situation.

Late Miocene development of the western Pacific warm pool: Planktonic foraminifer and oxygen isotopic evidence

The disappearance at ~10 Ma of the deep dwelling planktonic foraminifer Globoquadrina dehiscens from the western Pacific including the South China Sea was about 3 Myr earlier than its final extinction elsewhere. Accompanying this event at ~10 Ma was a series of faunal turnover characterized by increase in mixed layer, warm-water species and decrease to a minimum in deepwater species. Paleobiological and isotopic evidence indicates sea surface warming and a deepened local thermocline that we interpret as related to the development of an early western Pacific warm pool. The stepwise decline of G. dehiscens and other deep dwelling species from the NW and SW Pacific suggests more intensive warm water pileup than equatorial localities where surface bypass flow through the narrowing Indonesia seaway appears to remain efficient during the late Miocene. Planktonic delta18O values from the South China Sea consistently lighter than the tropical western Pacific during the Miocene also suggest, similar to today, more variable hydrologic conditions along the periphery than in the core of the warm pool. Stronger hydrologic variability affected mainly by monsoons and increased thermal gradient along the western margin of the late Miocene warm pool may have contributed to the decline of deep dwelling planktonic species including the early extinction of G. dehiscens from the South China Sea region. The late Miocene warm pool became influential and paleobiologically detectable from ~10 Ma, but the modern warm pool did not appear until about 4 Ma, in the middle Pliocene.

Past circulation along the western Iberian margin

Fifteen Iberian margin sediment cores, distributed between 43°12'N and 35°53'N, have been used to reconstruct spatial and temporal (sub)surface circulation along the Iberian margin since the Last Glacial period. Time-slice maps of planktonic foraminiferal derived summer sea surface temperature (SST) and export productivity (Pexp) were established for specific time intervals within the last 35 ky: the Holocene (Recent and last 8 ky), Younger Dryas (YD), Heinrich Stadials (HS) 1, 2a, 2b, 3, and the Last Glacial Maximum (LGM). The SST during the Holocene shows the same latitudinal gradient along the western Iberian margin as present-day with cold but productive areas that reflect the influence of coastal upwelling centers. The LGM appears as a slightly less warm, but more productive period relative to the Holocene and present-day conditions, suggesting that sea-level minima forced a westward displacement of the coastal upwelling centers possibly accompanied by a strengthening of northward winds. During the YD, a longitudinal thermal front is depicted at 10°W, with cold polar waters offshore and warmer subtropical waters inshore, suggesting that the subtropical Paleo-Iberian Poleward Current more likely flowed at a more inshore location masking the local SST signal and amplitude of variation. A substantial cooling and drop in productivity is observed during all HS, in particular HS1 and HS3, reflecting the penetration of icebergs-derived meltwater. These most extreme southward extensions of very cold waters define a strong SST gradient that marks a possible Paleo-Azores Front. Higher production south of this front was likely fed by frontal nutrient advection.

Uranium concentrations and isotope ratios Bahamas seawater and ODP holes pore-water

The geometry, timing, and rate of fluid-flow through carbonate margins and platforms is not well constrained. In this study, we use U concentrations and isotope ratios measured on small volumes of pore-water from Bahamas slope sediment, coupled with existing chlorinity data, to place constraints on the fluid-flow in this region and, by implication, other carbonate platforms. These data also allow an assessment of the behaviour of U isotopes in an unusually well constrained water-rock system. We report pore-water U concentrations which are controlled by dissolution of high-U organic material at shallow depths in the sediment and by reduction of U to its insoluble 4+ state at greater depths. The dominant process influencing pore-water (234U/238U) is alpha recoil. In Holocene sediments, the increase of pore-water (234U/238U) due to recoil provides an estimate of the horizontal flow rate of 11 cm/year, but with considerable uncertainty. At depths in the sediment where conditions are reducing, features in the U concentration and (234U/238U) profiles are offset from one another which constrains the effective diffusivity for U in these sediments to be c. 1-2 * 10**-8 cm**2/s. At depths between the Holocene and these reducing sediments, pore-water (234U/238U) values are unusually low due to a recent increase in the dissolution rate of grain surfaces. This suggests a strengthening of fluid flow, probably due to the flooding of the banks at the last deglaciation and the re-initiation of thermally-driven venting of fluid on the bank top and accompanying recharge on the slopes. Interpretation of existing chlorinity data, in the light of this change in flow rate, constrain the recent horizontal flow rate to be 10.6 ( 3.4) cm/year. Estimates of flow rate from (234U/238U) and Cl[-] are therefore in agreement and suggest flow rates close to those predicted by thermally-driven models of fluid flow. This agreement supports the idea that flow within the Bahamas Banks is mostly thermally driven and suggests that flow rates on the order of 10 cm/year are typical for carbonate platforms where such flow occurs.

Strontium isotopes, ages and carbonate content of sediments from the Agulhas Current

The Agulhas Leakage represents a significant portion of the warm, surface return flow of the global overturning circulation and thus may be an important feedback in the ocean climate system. Models indicate that reduced leakage could be caused by a stronger Agulhas Current and/or a more upstream (eastward) Agulhas Retroflection, while a weaker Agulhas Current would result in a more westward retroflection and increased leakage. However, data for the Last Glacial Maximum support both a weaker Agulhas Current and less leakage, implying a possible displacement of the retroflection. We present new 87Sr/86Sr results for modern sediments within this region, confirming that the modern pathway of the Agulhas Current, Retroflection, and Leakage can be traced by terrigenous sediment provenance using Sr isotopes. New 87Sr/86Sr data from sediments deposited during the Last Glacial Maximum suggest that the glacial Agulhas Current and Retroflection followed nearly their modern trajectory. The provenance data appear to rule out both a stronger Agulhas Current and a more upstream Agulhas Retroflection. We conclude that the reduced glacial leakage was caused by the weakened Agulhas Current, with no significant change in the retroflection position. This is inconsistent with the model predictions and thus emphasizes the need for further work in this region.

Objective Mapping of Argo data in the Weddell Gyre: a gridded dataset of upper ocean water properties, link to data files in NetCDF format

The Weddell Gyre plays a crucial role in the regulation of climate by transferring heat into the deep ocean through deep and bottom water mass formation. However, our understanding of Weddell Gyre water mass properties is limited to regions of data availability, primarily along the Prime Meridian. The aim is to provide a dataset of the upper water column properties of the entire Weddell Gyre. Objective mapping was applied to Argo float data in order to produce spatially gridded, time composite maps of temperature and salinity for fixed pressure levels ranging from 50 to 2000 dbar, as well as temperature, salinity and pressure at the level of the sub-surface temperature maximum. While the data are currently too limited to incorporate time into the gridded structure, the data are extensive enough to produce maps of the entire region across three time composite periods (2002-2005, 2006-2009 and 2010-2013), which can be used to determine how representative conclusions drawn from data collected along general RV transect lines are on a gyre scale perspective. The time composite data sets are provided as netCDF files; one for each time period. Mapped fields of conservative temperature, absolute salinity and potential density are provided for 41 vertical pressure levels. The above variables as well as pressure are provided at the level of the sub-surface temperature maximum. Corresponding mapping errors are also included in the netCDF files. Further details are provided in the global attributes, such as the unit variables and structure of the corresponding data array (i.e. latitude x longitude x vertical pressure level). In addition, all files ending in "_potTpSal" provide mapped fields of potential temperature and practical salinity.