Modelling nitrogen and oxygen isotope fractionation during denitrification in a lacustrine redox-transition zone

The stable isotope composition (δ¹⁵N and δ¹⁸O) of nitrate was measured during Summer 1999 in the anaerobic hypolimnion of eutrophic Lake Lugano (Switzerland). Denitrification was demonstrated by a progressive nitrate depletion coupled to increasing δ¹⁵N and δ¹⁸O values for residual nitrate. Maximum δ¹⁵N and δ¹⁸O values amounted to 27.2 and 15.7‰, respectively.¹⁵N and ¹⁸O enrichment factors for denitrification (ε) were estimated using a closed-system model and a dynamic diffusion-reaction model. Using the Rayleigh equation (closed-system approach), we obtained ε values of −11.2 and −6.6‰ for nitrogen and oxygen, respectively. The average ε values derived using the diffusion-reaction model were determined to be −20.7 ± 3.8 for nitrogen and −11.0 ± 1.7 for oxygen. Both N and O isotope fractionation appeared to be lower when denitrification rates where high, possibly in association with high organic carbon availability. In addition, variations in the isotope effects may be attributed to the variable importance of sedimentary denitrification having only a small isotope effect on the water column. The combined measurement of N and O isotope ratios in nitrate revealed that coupled nitrification-denitrification in the open-water was of minor importance. This is the first study of nitrogen and oxygen isotope effects associated with microbial denitrification in a natural lake. Moreover, this study confirms the high potential of δ¹⁸O of nitrate as a valuable biogeochemical tracer in aquatic systems, complementing nitrate δ¹⁵N.     

Cation ordering of orthopyroxenes from the Skaergaard Intrusion: implications for the subsolidus cooling rates and permeabilities

 We have determined the quenched cation ordering states of five orthopyroxene crystals collected from the marginal border group and the lower zone a and b of the Skaergaard intrusion, and modeled these data to retrieve their closure temperatures (T _c) of Fe–Mg ordering and cooling rates. According to existing thermal models for the Skaergaard pluton, conductive cooling dominated the high and low temperature regimes, which were separated by an intermediate temperature regime in which the cooling was controlled primarily by convective fluid circulation. The cooling rates retrieved from the quenched cation ordering states of the orthopyroxene crystals strictly apply to temperatures around the closure temperatures of the ordering states, ∼340–400^° C, which fall at the transition from convective to the lower temperature conductive cooling. The cooling rates obtained from the cation ordering states of orthopyroxene vary from ∼1 to 270 K/ka. These results agree well with a thermal model calculated using an assumed average permeability of 10^-12 cm² for the pluton, but not completely with a model calculated on the basis of an average permeability of 10^-13 cm², although both values produced shifts of δ¹⁸O that are comparable to those observed in the pluton. Received: 27 February 1995/Accepted: 20 July 1995     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

Exhumation patterns along shallow low‐angle normal faults: an example from the Altotiberina active fault system (Northern Apennines,Italy)

A multi‐method approach (palaeothermal and thermochronological analyses; thermal modelling) is applied to reconstruct the exhumation history of the Altotiberina Fault (ATF), a representative example of crustal‐scale active low‐angle normal faulting in the Northern Apennines (Italy). Thermal maturity and thermochronological data yield similar burial histories but different exhumation patterns for the sedimentary successions in the hangingwall and the footwall of the ATF. Since 3.8 Ma, the ATF footwall has exhumed at rates of 0.90 mm a^?1. Exhumation led to bending and deactivation of the ATF uppermost portion as a result of tectonic unloading and isostatic adjustment, followed by migration of extension and the development of a set of domino‐like, east‐dipping normal faults, rooting on the buried portion of the ATF. ATF activity and isostatic rebound exhumed Triassic rock units from depths of about 4 km. We suggest that isostatic instability is accommodated at shallow crustal levels, in a similar way to what is observed on larger structures at mid‐low crustal levels.     

Ressi experimental catchment: Ecohydrological research in the Italian pre-Alps

Ressi is a small (2.4 ha) forested catchment located in the Italian pre-Alps. The site became an experimental catchment to investigate the water fluxes in the soil–plant–atmosphere continuum and the impact of vegetation on runoff generation in 2012. The elevation of the catchment ranges from 598 to 721 m a.s.l. and the climate is humid temperate. The bedrock consists of rhyolites and dacites; the soil is a Cambisol. The catchment is covered by a dense forest, dominated by beech, chestnut, maple, and hazel trees. The field set up includes measurements of the rainfall in an open area, streamflow at the outlet, soil moisture at various depths and locations, and depth to water table in six piezometers at a 5- or 10-min interval. Samples of precipitation, stream water, shallow groundwater and soil water are collected monthly for tracer analysis (stable isotopes (²H and ¹⁸O), electrical conductivity and major ions), and during selected rainfall–runoff events to determine the contribution of the various sources to runoff. Since 2017, soil and plant water samples have been collected to determine the sources of tree transpiration. Data collected in the period 2012–2016 are publicly available. Data collection is ongoing, and the data set is expected to be updated on an annual basis to include the most recent measurements.