Magmatic zircon oxygen isotopes of 1.88–1.87 Ga orogenic and 1.65–1.54 Ga anorogenic magmatism in Finland

Summary Oxygen isotope ratios of igneous zircon from magmatic rocks in Finland provide insights into the evolution and growth of the Precambrian crust during the Svecofennian orogeny. These data preserve magmatic δ¹⁸O values and correlate with major discontinuities in the lower crust. Oxygen isotope ratios of zircon across the 1.88–1.87 Ga Central Finland granitoid complex (CFGC) range from 5.50‰ to 6.84‰, except for three plutons in contact with the adjacent greenstone and metasedimentary belts (δ¹⁸O(Zrc) = 7.60‰–7.78‰). There is a systematic variation in δ¹⁸O(Zrc) with respect to geographic location in the CFGC, ranging from 6.60±0.23‰ (σ) in the northeast to 5.90±0.40‰ in the west-southwest. These values correlate with a change in crustal thickness and shift in geochemical composition. The oxygen isotope composition of the 1.65–1.54 Ga rapakivi granites and related rocks in southern Finland show a decreasing trend from north to south, independent of their emplacement age. The southern anorogenic granite group has an average δ¹⁸O in zircon of 6.14±0.07‰ and the northern anorogenic group has an average δ¹⁸O in zircon of 8.14±0.59‰. This difference reflects the boundary between island arc terrains accreted during the Paleoproterozoic. Deceased     

Concentration, isotopic composition, and sources of lead in Southern Ocean air during 1999/2000, measured at the Cape Grim Baseline Air Pollution Station, Tasmania

Southern Ocean aerosols were collected at the Cape Grim Baseline Air Pollution Station from onshore air under baseline conditions between February 1999 and April 2000. Thermal ionization techniques (TIMS) and isotope dilution mass spectrometry (IDMS) were used to measure the isotopic composition and concentration of lead in the air giving concentrations as low as 0.6 ± 0.1 pg · m⁻³. Air collected under baseline conditions for 12 months (May 1999-April 2000) yielded an overall lead concentration of 11.0 ± 0.2 pg · m⁻³ and isotopic composition of ²⁰⁶Pb/²⁰⁷Pb = 1.154, ²⁰⁸Pb/²⁰⁷Pb = 2.387 and ²⁰⁶Pb/²⁰⁴Pb = 17.93. The range in isotopic ratios was consistent with the mixing of lead from major population centers in the Southern Hemisphere in the mid to high latitudes, except for the presence of highly radiogenic lead in some samples. Contributions from radiogenic lead of up to ∼0.8% were observed. Three periods with the highest percentage contribution of radiogenic lead (>0.5%) were investigated in more detail, and 4-d back-trajectories and radon concentrations were used to help identify the sources. The sources are probably associated with the mining and processing of uranium rich ores in southern Africa and possibly South Australia.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

Luminescence dating of ice‐marginal deposits in northern Germany: evidence for repeated glaciations during the Middle Pleistocene (MIS 12 to MIS 6)

The exact number, extent and chronology of the Middle Pleistocene Elsterian and Saalian glaciations in northern Central Europe are still controversial. This study presents new luminescence data from Middle Pleistocene ice‐marginal deposits in northern Germany, giving evidence for repeated glaciations during the Middle Pleistocene (MIS 12 to MIS 6). The study area is located in the Leine valley south of the North German Lowlands. The data set includes digital elevation models, high‐resolution shear wave seismic profiles, outcrop and borehole data integrated into a 3D subsurface model to reconstruct the bedrock relief surface. For numerical age determination, we performed luminescence dating on 12 ice‐marginal and two fluvial samples. Luminescence ages of ice‐marginal deposits point to at least two ice advances during MIS 12 and MIS 10 with ages ranging from 461±34 to 421±25 ka and from 376±27 to 337±21 ka. The bedrock relief model and different generations of striations indicate that the older ice advance came from the north and the younger one from the northeast. During rapid ice‐margin retreat, subglacial overdeepenings were filled with glaciolacustrine deposits, partly rich in re‐worked Tertiary lignite and amber. During MIS 8 and MIS 6, the study area may have been affected by two ice advances. Luminescence ages of glaciolacustrine delta deposits point to a deposition during MIS 8 or early MIS 6, and late MIS 6 (250±20 to 161±10 ka). The maximum extent of both the Elsterian (MIS 12 and MIS 10) and Saalian glaciations (MIS 8? and MIS 6) approximately reached the same position in the Leine valley and was probably controlled by the formation of deep proglacial lakes in front of the ice sheets, preventing a further southward advance.     

Energy partitioning over a semi‐arid shrubland in northern China

During the last decade, the widely distributed shrublands in northern China have shown significant signs of recovery from desertification, the result of widespread conservation practices. However, to support the current efforts in conservation, more knowledge is needed on surface energy partitioning and its biophysical controls. Using eddy‐covariance measurements made over a semi‐arid shrubland in northwest China in 2012, we examined how surface energy‐balance components vary on diurnal and seasonal scales, and how biophysical factors control bulk surface parameters and energy exchange. Sensible heat flux (H) exceeded latent heat flux (λE) during most of the year, resulting in an annual Bowen ratio (β, i.e. H/λE) of 2.0. λE exceeded H only in mid‐summer when frequent rainfall co‐occurred with the seasonal peak in leaf area index (LAI). Evapotranspiration reached a daily maximum of 3.3 mm day^?1, and summed to 283 mm yr^?1. The evaporative fraction (EF, i.e. λE/R_n), Priestley–Taylor coefficient (α), surface conductance (g_s) and decoupling coefficient (Ω) were all positively correlated with soil water content (SWC) and LAI. The direct enhancement of λE by high vapour pressure deficit (VPD) was buffered by a concurrent suppression of g_s. The g_s played a direct role in controlling EF and α by mediating the effects of LAI, SWC and VPD. Our results highlight the importance of adaptive plant responses to water scarcity in regulating ecosystem energy partitioning, and suggest an important role for revegetation in the reversal of desertification in semi‐arid areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.