Integrating lysimeter drainage and eddy covariance flux measurements in a groundwater recharge model

Abstract

Field-scale water balance is difficult to characterize because controls exerted by soils and vegetation are mostly inferred from local-scale measurements with relatively small support volumes. Eddy covariance flux and lysimeters have been used to infer and evaluate field-scale water balances because they have larger footprint areas than local soil moisture measurements. This study quantifies heterogeneity of soil deep drainage (D) in four 12.5-m² repacked lysimeters, compares evapotranspiration from eddy covariance (ET_EC) and mass balance residuals of lysimeters (ETwb_Lys), and models D to estimate groundwater recharge. Variation in measured D was attributed to redirection of snowmelt infiltration and differences in lysimeter hydraulic properties caused by surface soil treatment. During the growing seasons of 2010, 2011 and 2012, ETwb_Lys (278, 289 and 269 mm, respectively) was in good agreement with ET_EC (298, 301 and 335 mm). Annual recharge estimated from modelled D was 486, 624 and 613 mm for three calendar years 2010, 2011 and 2012, respectively. In summary, lysimeter D and ET_EC can be integrated to estimate and model groundwater recharge.

Editor D. Koutsoyiannis     

Soil water residence times and solute uptake on a dolomite bedrock—preliminary results

Calcium and magnesium levels have been monitored in slope foot drainage waters on a dolomite bedrock. Both calcium and magnesium rich pulses occur. Short term dissolution experiments demonstrate high calcium levels in solution while other authors have suggested that long residence time groundwater has relatively high levels of magnesium due to calcite precipitation. Patterns of field fluctuations in Ca: Mg ratios can thus be tentatively interpreted in terms of short residence time water of high calcium content mixing with long residence time water of high magnesium content. Fluorometric dye tracing has been used to indicate the orders of magnitude of soil water residence times, suggesting that quickflow components are resident in the system for a few hours to a few days. Further work is in progress.     

Carbon fluxes from eroding peatlands – the carbon benefit of revegetation following wildfire

Peatlands are among the largest long‐term soil carbon stores, but their degradation can lead to significant carbon losses. This study considers the carbon budget of peat‐covered sites after restoration, following degradation by past wildfires. The study measured the carbon budget of eight sites: four restored‐revegetated sites, two unrestored bare soil control sites, and two intact vegetated controls over two years (2006–2008). The study considered the following flux pathways: dissolved organic carbon (DOC); particulate organic carbon (POC); dissolved carbon dioxide (CO₂); primary productivity; net ecosystem respiration, and methane (CH₄). The study shows that unrestored, bare peat sites can have significant carbon losses as high as 522 ± 3 tonnes C/km²/yr. Most sites showed improved carbon budgets (decreased source and/or increased sink of carbon) after restoration; this improvement was mainly in the form of a reduction in the size of the net carbon source, but for one restored site the measured carbon budget after four years of restoration was greater than observed for vegetated controls. The carbon sequestration benefit of peatland restoration would range between 122 and 833 tonnes C/km²/yr. Copyright © 2011 John Wiley & Sons, Ltd.     

The age as a diagnostic of the dynamics of marine ecosystem models

The constituent-oriented age theory (CAT) worked out by Delhez et al. (1999) is a flexible tool that can be applied to diagnose complex models. It is shown here how this can be used to quantify the pace at which an ecosystem model works. At the cost of the introduction of one additional evolution equation for each compartment of the ecosystem model, the mean age of the biological material forming these compartments can be computed. The information obtained in this way complements the information provided by the concentration data; while the latter measures the standing stocks, the former provides an integrated assessment of the interaction rates and matter fluxes. The benefits of the method are demonstrated with a simple Lotka–Volterra system and a one-dimensional vertical model of the nitrogen cycle in the Ligurian Sea. The theory can be used to study the biological compartments individually or the ecosystem as a whole. In particular, the age is a valuable tool to quantify the overall cycling rate of nitrogen in the food web.Responsible Editor: Phil Dyke     

ANISOTROPIC TRAVELTIME TOMOGRAPHY

Velocity estimation technique using seismic data is often based on time/distance equations which are independent of direction, and even though we now know that many rocks are quite anisotropic, useful results have been obtained over the years from these isotropic estimates. Nevertheless, if velocities are significantly direction-dependent, then the isotropic assumption may lead to serious structural interpretation errors. Additionally, information on angle-dependence may lead to a better understanding of the lithology of the rocks under measurement. VSP and cross-well data may each lack the necessary aperture to estimate more than two velocity parameters for each wave type, and if the data straddle a symmetry axis, then these may be usefully chosen to be the direct velocities (from time-and-distance measurements along the axis) and NMO velocities (from differential time-offset measurements). These sets of two parameters define ellipses, and provide intermediate models for the variation of velocity with angle which can later be assembled and translated into estimates of the elastic moduli of the rocks under scrutiny. If the aperture of the measurements is large enough e.g. we have access to both VSP and cross-well data, we divide the procedure into two independent steps, first fitting best ellipses around one symmetry axis and then fitting another set around the orthogonal axis. These sets of four elliptical parameters are then combined into a new, double elliptical approximation. This approximation keeps the useful properties of elliptical anisotropy, in particular the simple relation between group and phase velocities which simplifies the route from the traveltimes measurements to the elastic constants of the medium. The inversion proposed in this paper is a simple extension of well-known isotropic schemes and it is conceptually identical for all wave types. Examples are shown to illustrate the application of the technique to cross-well synthetic and field P-wave data. The examples demonstrate three important points: (a) When velocity anisotropy is estimated by iterative techniques such as conjugate gradients, early termination of the iterations may produce artificial anisotropy. (b) Different components of the velocity are subject to different type of artifacts because of differences in ray coverage, (c) Even though most rocks do not have elliptical dispersion relations, our elliptical schemes represent a useful intermediate step in the full characterization of the elastic properties.     

Hydrothermal alteration of Tertiary igneous rocks from the Isle of Skye,northwest Scotland

Hydrothermal alteration of Tertiary gabbros from Skye involved the reaction of igneous olivine, augite, hypersthene, plagioclase, magnetite, and ilmenite with aqueous fluid primarily to combinations of talc, chlorite, montmorillonite, calcic amphibole, biotite, and secondary magnetite. Lesser amounts of calcite, epidote, quartz, sphene, prehnite, and garnet also developed. During mineralogical alteration of gabbro there was a net addition to rock of K, Na, Sr, and H₂O and a net loss of Mg. Gabbro was oxidized early in the hydrothermal event and later reduced. Iron and silicon were probably initially lost and later added. There is no evidence for significant change in the Al or Ca content of the gabbros. Hydrothermal alteration of Skye gabbro involved not only large-scale migration of ¹⁸O, ¹⁶O, D and H but also of K, Na, Sr, Mg, and probably Fe and Si.Mineral thermometry indicates that pyroxenes in the gabbros crystallized at 1000° C–1150° C and were very resistent chemically as well as isotopically to later hydrothermal alteration. Hypothetical equilibrium between primary and secondary mafic silicates suggests that mineralogical alteration of gabbro occurred at 450°–550° C. The lack of correlation between mineralogical and isotopic alteration of gabbro requires that much isotopic alteration occurred at temepratures above those at which the secondary minerals developed, 550°–1000° C. The chemical alteration of gabbro is correlated with its mineralogical alteration and therefore occurred at 450°–550° C.Measured progress of the mineral-fluid reactions was used to estimate the amount of H₂O fluid that infiltrated the gabbro as primary olivine was converted to talc+magnetite at 525°–550° C. Calculated fluid-rock ratios are in the range 0.2–6 (volume basis) and are smaller than values estimated from isotopic data (fluid/rock 1–10, volume basis). Both isotopic and petrologic data point to pervasive flow of fluid through crystalline rock at elevated temperatures of 500°–1000° C. Isotopic fluid-rock ratios are larger than petrologic fluid-rock ratios because isotopic alteration of cooling gabbro began earlier and at higher temperatures than did the mineralogical alteration.     

Interaction of the Atmosphere and Lithosphere in the Minute Range of Periods

Doklady Earth Sciences - The interaction between the atmosphere and the lithosphere in the minute range of periods is considered. The technique of assessment of the influence of atmospheric...     

Olivine in ultramafic lamprophyres: chemistry,crystallisation, and melt sources of Siberian Pre- and post-trap aillikites

We studied olivines from the Devonian pre-trap (the Ilbokich occurrence) and the Triassic post-trap (the Chadobets occurrence) carbonate-rich ultramafic lamprophyres (UMLs) in the southwestern portion of the Siberian craton. On the basis of detailed investigations of major, minor, and trace-element distributions, we have reconstructed the main processes that control the origins of these olivines. These include fractional crystallisation from melt, assimilation, and fractional crystallisation processes with orthopyroxene assimilation, melt-reaction diffusive re-equilibration, alkali enrichment, and CO₂ degassing of the melt. Furthermore, we inferred the composition of the sources of the primary UML melt and their possible correlations with proto-kimberlitic melts, as well as the influence of the Triassic Siberian plume on the composition of the lithospheric mantle. The main differences between olivines from the Ilbokich and the Chadobets aillikites were that the olivines from the former had more magnesium-rich cores (Mg# = 89.2?±?0.2), had Mg- and Cr-rich transition zones (Mg# = 89.7?±?0.2 and 300–500 ppm Cr), had lower Ni (up to 3100 ppm) and Li (1.4–1.5 ppm), and had higher B (0.8–2.6 ppm) contents, all at higher Fo values (90–86), relative to the olivines from the latter (Mg# = 88–75; 200–300 ppm Cr; up to 3400 ppm Ni; 1.4–2.4 ppm Li; 0.4–2.2 ppm B). The Siberian aillikite sources contained a significant amount of metasomatic material. Phlogopite-rich MARID-type veins provided the likely metasomatic component in the pre-trap Devonian Ilbokich aillikite source, whereas the Triassic Chadobets aillikitic post-trap melts were derived from a source with a significant carbonate component. A comparison of UML olivines with olivines from the pre-trap and post-trap Siberian kimberlites shows a striking similarity. This suggests that the carbonate component in the aillikitic source could have been produced by evolved kimberlite melts. The differences in the lithospheric metasomatic component that contributed to pre-trap and post-trap aillikitic melts can be interpreted as reflections of the thermal impact of the Siberian Traps, which reduced phlogopite-bearing metasomes within the southwestern Siberian sub-continental lithospheric mantle.