What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

Groundwater recharge sources in semiarid irrigated mountain fronts

High-elevation mountains often constitute for basins important groundwater recharge sources through mountain-front recharge processes. These processes include streamflow losses and subsurface inflow from the mountain block. However, another key recharge process is from irrigation practices, where mountain streamflow is distributed across the irrigated piedmont. In this study, coupled groundwater fluctuation measurements and environmental tracers (¹⁸O, ²H, and major ions) were used to identify and compare the natural mountain-front recharge to the anthropogenically induced irrigation recharge. Within the High Atlas mountain front of the Ourika Basin, Central Morocco, the groundwater fluctuation mapping from the dry to wet season showed that recharge beneath the irrigated area was higher than the recharge along the streambed. Irrigation practices in the region divert more than 65% of the stream water, thereby reducing the potential for in-stream groundwater recharge. In addition, the irrigation areas close to the mountain front had greater water table increases (up to 3.5 m) compared with the downstream irrigation areas (<1 m increase). Upstream crops have priority to irrigation with stream water over downstream areas. The latter are only irrigated via stream water during large flood events and are otherwise supplemented by groundwater resources. These changes in water resources used for irrigation practices between upstream and downstream areas are reflected in the spatiotemporal evolution of the stable isotopes of groundwater. In the upstream irrigation area, the groundwater stable isotope values (δ¹⁸O: −8.4‰ to −7.4‰) reflect recharge by the diverted stream water. In the downstream irrigation area, the groundwater isotope values are lower (δ¹⁸O: −8.1‰ to −8.4‰) due to recharge via the flood water. In the nonirrigation area, the groundwater has the highest stable isotope values (δ¹⁸O: −6.8‰ to −4.8‰). This might be due to recharge via subsurface inflow from the mountain block to the mountain front and/or recharge via local low altitude rainfall. These findings highlight that irrigation practices can result in the dominant mountain-front recharge process for groundwater.     

Lateral extrusion along the Altyn Tagh Fault,Qilian Shan (NE Tibet): insight from a 3D crustal budget

The lithospheric strike‐slip Altyn Tagh Fault has accommodated hundreds of kilometres of displacement between the Qaidam and Tarim blocks since its Eocene reactivation. However, the way the deformation is accommodated in the Qilian Shan and further east remains uncertain. Based on 360 km of north‐eastward migration of the relatively rigid Qaidam block along the Altyn Tagh Fault and 3D isovolumetric balancing of the crustal deformation within the Altyn Tagh Fault–Qilian Shan system, we demonstrate that 250 ± 28 km (43.8–49.4%) of N20E directed crustal shortening and an additional ~250–370 km of eastward motion of the Qilian Shan crust must be accounted for by strike‐slip faulting in the Qilian Shan and crustal thickening in the Qinling area, as well as by extension in the adjoining North China block graben systems.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.     

Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.