Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.

Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

The Value of Adding Optics to Ecosystem Models: A Case Study

Many ecosystem models have been developed to study the ocean's biogeochemical properties, but most of these models use simple formulations to describe light penetration and spectral quality. Here, an optical model is coupled with a previously published ecosystem model that explicitly represents two phytoplankton (picoplankton and diatoms) and two zooplankton functional groups, as well as multiple nutrients and detritus. Surface ocean color fields and subsurface light fields are calculated by coupling the ecosystem model with an optical model that relates biogeochemical standing stocks with inherent optical properties (absorption, scattering); this provides input to a commercially available radiative transfer model (Ecolight). We apply this bio-optical model to the equatorial Pacific upwelling region, and find the model to be capable of reproducing many measured optical properties and key biogeochemical processes in this region. Our model results suggest that non-algal particles largely contribute to the total scattering or attenuation (> 50% at 660 nm) but have a much smaller contribution to particulate absorption (< 20% at 440 nm), while picoplankton dominate the total phytoplankton absorption (> 95% at 440 nm). These results are consistent with the field observations. In order to achieve such good agreement between data and model results, however, key model parameters, for which no field data are available, have to be constrained. Sensitivity analysis of the model results to optical parameters reveals a significant role played by colored dissolved organic matter through its influence on the quantity and quality of the ambient light. Coupling explicit optics to an ecosystem model provides advantages in generating: (1) a more accurate subsurface light-field, which is important for light sensitive biogeochemical processes such as photosynthesis and photo-oxidation, (2) additional constraints on model parameters that help to reduce uncertainties in ecosystem model simulations, and (3) model output which is comparable to basic remotely-sensed properties. In addition, the coupling of biogeochemical models and optics paves the road for future assimilation of ocean color and in-situ measured optical properties into the models.

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.

Ocean (De)oxygenation Across the Last Deglaciation: Insights for the Future

Anthropogenic warming is expected to drive oxygen out of the ocean as the water temperature rises and the rate of exchange between subsurface waters and the atmosphere slows due to enhanced upper ocean density stratification. Observations from recent decades are tantalizingly consistent with this prediction, though these changes remain subtle in the face of natural variability. Earth system model projections unanimously predict a long-term decrease in the global ocean oxygen inventory, but show regional discrepancies, particularly in the most oxygen-depleted waters, owing to the complex interplay between oxygen supply pathways and oxygen consumption. The geological record provides an orthogonal perspective, showing how the oceanic oxygen content varied in response to prior episodes of climate change. These past changes were much slower than the current, anthropogenic change, but can help to appraise sensitivities, and point toward potentially dominant mechanisms of change. Consistent with the model projections, marine sediments recorded an overall expansion of low-oxygen waters in the upper ocean as it warmed at the end of the last ice age. This expansion was not linearly related with temperature, though, but reached a deoxygenation extreme midway through the warming. Meanwhile, the deep ocean became better oxygenated, opposite the general expectation. These observations require that significant changes in apparent oxygen utilization occurred, suggesting that they will also be important in the future.

Erosive tooth wear: a multifactorial condition of growing concern and increasing knowledge.

Dental erosion is often described solely as a surface phenomenon, unlike caries where it has been established that the destructive effects involve both the surface and the subsurface region. However, besides removal of the surface, erosion shows dissolution of mineral within the softened layer - beneath the surface. In order to distinguish this process from the carious process it is now called 'near surface demineralization'. Erosion occurs in low pH, but there is no fixed critical pH value concerning dental erosion. The critical pH value for enamel concerning caries (pH 5.5-5.7) has to be calculated from calcium and phosphate concentrations of plaque fluid. In the context of dental erosion, the critical pH value is calculated from the calcium and phosphate concentrations in the erosive solution itself. Thus, critical pH for enamel with regard to erosion will vary according to the erosive solution. Erosive tooth wear is becoming increasingly significant in the management of the long-term health of the dentition. What is considered as an acceptable amount of wear is dependent on the anticipated lifespan of the dentition and is, therefore, different for deciduous compared to permanent teeth. However, erosive damage to the teeth may compromise the patient's dentition for their entire lifetime and may require repeated and increasingly complex and expensive restorations. Therefore, it is important that diagnosis of the tooth wear process in children and adults is made early and that adequate preventive measures are undertaken. These measures can only be initiated when the risk factors are known and interactions between them are present.

Sea surface temperature and sea ice variability in the subpolar North Atlantic from explosive volcanism of the late thirteenth century

In this study, we use IP and alkenone biomarker proxies to document the subdecadal variations of sea ice and sea surface temperature in the subpolar North Atlantic induced by the decadally paced explosive tropical volcanic eruptions of the second half of the thirteenth century. The short- and long-term evolutions of both variables were investigated by cross analysis with a simulation of the IPSL-CM5A LR model. Our results show short-term ocean cooling and sea ice expansion in response to each volcanic eruption. They also highlight that the long response time of the ocean leads to cumulative surface cooling and subsurface heat buildup due to sea ice capping. As volcanic forcing relaxes, the surface ocean rapidly warms, likely amplified by subsurface heat, and remains almost ice free for several decades

Tracking water pathways in steep hillslopes by d18O depth profiles of soil water

Assessing temporal variations in soil water flow is important, especially at the hillslope scale, to identify mechanisms of runoff and flood generation and pathways for nutrients and pollutants in soils. While surface processes are well considered and parameterized, the assessment of subsurface processes at the hillslope scale is still challenging since measurement of hydrological pathways is connected to high efforts in time, money and personnel work. The latter might not even be possible in alpine environments with harsh winter processes. Soil water stable isotope profiles may offer a time-integrating fingerprint of subsurface water pathways. In this study, we investigated the suitability of soil water stable isotope (d18O) depth profiles to identify water flow paths along two transects of steep subalpine hillslopes in the Swiss Alps. We applied a one-dimensional advection–dispersion model using d18O values of precipitation (ranging from _24.7 to _2.9‰) as input data to simulate the d18O profiles of soil water. The variability of d18O values with depth within each soil profile and a comparison of the simulated and measured d18O profiles were used to infer information about subsurface hydrological pathways. The temporal pattern of d18O in precipitation was found in several profiles, ranging from _14.5 to _4.0‰. This suggests that vertical percolation plays an important role even at slope angles of up to 46_. Lateral subsurface flow and/or mixing of soil water at lower slope angles might occur in deeper soil layers and at sites near a small stream. The difference between several observed and simulated d18O profiles revealed spatially highly variable infiltration patterns during the snowmelt periods: The d18O value of snow (_17.7 ± 1.9‰) was absent in several measured d18O profiles but present in the respective simulated d18O profiles. This indicated overland flow and/or preferential flow through the soil profile during the melt period. The applied methods proved to be a fast and promising tool to obtain time-integrated information on soil water flow paths at the hillslope scale in steep subalpine slopes.

From In Situ towards In Operando Conditions: Scanning Tunneling Microscopy Study of Hydrogen Intercalation in Cu(111) during Hydrogen Evolution

We used electrochemical scanning tunneling microscopy to study the intercalation of hydrogen into a Cu(111) model electrode under reactive (in operando) conditions. Hydrogen evolution causes hydrogen intermediates to migrate into the copper lattice as function of the applied potential and the resulting current density. This H-inclusion is demonstrated to be reversible. The presence of subsurface hydrogen leads to a significant surface relaxation/reconstruction affecting both the geometric and electronic structure of the electrode surface.

Radar investigations of Apollinaris Mons on Mars: Exploring the origin of the fan deposits

Apollinaris Mons is an isolated volcano on Mars straddling the boundary between the southern highlands and the northern plains. One of its most distinctive features is its massive fan-shaped deposit that extends from a breach on its summit to distances of more than 150 km and drapes its entire southern flank. The composition and formation mechanism of these deposits remains controversial. We investigate the radar properties of the fan deposits (FD) of Apollinaris Mons using low-frequency sounding radar data in combination with high-resolution images and crater-size frequency analysis to constrain their inner shape and bulk composition. Our analysis indicates that the FD attains an irregular thickness and is gradually thinner towards their lateral margins. The crater-size frequency analysis shows that they may have undergone repeated resurfacing, which is suggestive of long-term evolution. Our analysis of Shallow Radar (SHARAD) radargrams traversing different sections of the FD reveals multiple and different subsurface interfaces among the radargrams crossing the thinnest part, which suggests a layered and complex inner shape. Our estimates for the bulk real part of the dielectric constant of the FD ranges from 3 to 5, which is consistent with an icy-silicate mixture or pyroclastic composition. Therefore, we conclude that lahars or pyroclastic flows are the most likely mechanism that created the FD, yet we cannot rule out additional contributions from lava flows. A combination of multiple processes is also possible since the deposits appear to have been modified by fluvial processes at a later stage of their formation.

The future of fluorides and other protective agents in erosion prevention

The effectiveness of fluoride in caries prevention has been convincingly proven. In recent years, researchers have investigated the preventive effects of different fluoride formulations on erosive tooth wear with positive results, but their action on caries and erosion prevention must be based on different requirements, because there is no sheltered area in the erosive process as there is in the subsurface carious lesions. Thus, any protective mechanism from fluoride concerning erosion is limited to the surface or the near surface layer of enamel. However, reports on other protective agents show superior preventive results. The mechanism of action of tin-containing products is related to tin deposition onto the tooth surface, as well as the incorporation of tin into the near-surface layer of enamel. These tin-rich deposits are less susceptible to dissolution and may result in enhanced protection of the underlying tooth. Titanium tetrafluoride forms a protective layer on the tooth surface. It is believed that this layer is made up of hydrated hydrogen titanium phosphate. Products containing phosphates and/or proteins may adsorb either to the pellicle, rendering it more protective against demineralization, or directly to the dental hard tissue, probably competing with H(+) at specific sites on the tooth surface. Other substances may further enhance precipitation of calcium phosphates on the enamel surface, protecting it from additional acid impacts. Hence, the future of fluoride alone in erosion prevention looks grim, but the combination of fluoride with protective agents, such as polyvalent metal ions and some polymers, has much brighter prospects.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

On the seasonal variations of salinity of the tropical Atlantic mixed layer

The physical processes controlling the mixed layer salinity (MLS) seasonal budget in the tropical Atlantic Ocean are investigated using a regional configuration of an ocean general circulation model. The analysis reveals that the MLS cycle is generally weak in comparison of individual physical processes entering in the budget because of strong compensation. In evaporative regions, around the surface salinity maxima, the ocean acts to freshen the mixed layer against the action of evaporation. Poleward of the southern SSS maxima, the freshening is ensured by geostrophic advection, the vertical salinity diffusion and, during winter, a dominant contribution of the convective entrainment. On the equatorward flanks of the SSS maxima, Ekman transport mainly contributes to supply freshwater from ITCZ regions while vertical salinity diffusion adds on the effect of evaporation. All these terms are phase locked through the effect of the wind. Under the seasonal march of the ITCZ and in coastal areas affected by river (7°S:15°N), the upper ocean freshening by precipitations and/or runoff is attenuated by vertical salinity diffusion. In the eastern equatorial regions, seasonal cycle of wind forced surface currents advect freshwaters, which are mixed with subsurface saline water because of the strong vertical turbulent diffusion. In all these regions, the vertical diffusion presents an important contribution to the MLS budget by providing, in general, an upwelling flux of salinity. It is generally due to vertical salinity gradient and mixing due to winds. Furthermore, in the equator where the vertical shear, associated to surface horizontal currents, is developed, the diffusion depends also on the sheared flow stability.

A deep crust–mantle boundary in the asteroid 4 Vesta

The asteroid 4 Vesta was recently found to have two large impact craters near its south pole, exposing subsurface material. Modelling suggested that surface material in the northern hemisphere of Vesta came from a depth of about 20 kilometres, whereas the exposed southern material comes from a depth of 60 to 100 kilometres. Large amounts of olivine from the mantle were not seen, suggesting that the outer 100 kilometres or so is mainly igneous crust. Here we analyse the data on Vesta and conclude that the crust–mantle boundary (or Moho) is deeper than 80 kilometres.