Constraints on Upward Migration of Hydraulic Fracturing Fluid and Brine

Recent increases in the use of hydraulic fracturing (HF) to aid extraction of oil and gas from black shales have raised concerns regarding potential environmental effects associated with predictions of upward migration of HF fluid and brine. Some recent studies have suggested that such upward migration can be large and that timescales for migration can be as short as a few years. In this article, we discuss the physical constraints on upward fluid migration from black shales (e.g., the Marcellus, Bakken, and Eagle Ford) to shallow aquifers, taking into account the potential changes to the subsurface brought about by HF. Our review of the literature indicates that HF affects a very limited portion of the entire thickness of the overlying bedrock and therefore, is unable to create direct hydraulic communication between black shales and shallow aquifers via induced fractures. As a result, upward migration of HF fluid and brine is controlled by preexisting hydraulic gradients and bedrock permeability. We show that in cases where there is an upward gradient, permeability is low, upward flow rates are low, and mean travel times are long (often >10⁶ years). Consequently, the recently proposed rapid upward migration of brine and HF fluid, predicted to occur as a result of increased HF activity, does not appear to be physically plausible. Unrealistically high estimates of upward flow are the result of invalid assumptions about HF and the hydrogeology of sedimentary basins.     

Sm–Nd isotopic mapping of lithospheric growth and stabilization in the eastern Kaapvaal craton

Whole-rock Sm–Nd isotope systematics of 79 Archean granitoids from the eastern Kaapvaal craton, southern Africa, are used to delineate lithospheric boundaries and to constrain the timescale of crustal growth, assembly and geochemical differentiation c. 3.66–2.70 Ga. Offsets in ε_Nd values for 3.2–3.3 Ga granitoids across the Barberton greenstone belt (BGB) are consistent with existing models for c. 3.23 Ga accretion of newly formed lithosphere north of the BGB onto pre-existing c. 3.66 Ga lithosphere south of the BGB along a doubly verging subduction margin. The Nd isotopic signature of c. 3.3–3.2 Ga magmatic rocks show that significant crustal growth occurred during subduction–accretion. After c. 3.2 Ga, however, the Nd signature of intrusive rocks c. 3.1 and 2.7 Ga is dominated by intracrustal recycling rather than by new additions from the mantle, signalling cratonic stability.     

A Comparison of Hyporheic Transport at a Cross‐Vane Structure and Natural Riffle

While restoring hyporheic flowpaths has been cited as a benefit to stream restoration structures, little documentation exists confirming that constructed restoration structures induce comparable hyporheic exchange to natural stream features. This study compares a stream restoration structure (cross‐vane) to a natural feature (riffle) concurrently in the same stream reach using time‐lapsed electrical resistivity (ER) tomography. Using this hydrogeophysical approach, we were able to quantify hyporheic extent and transport beneath the cross‐vane structure and the riffle. We interpret from the geophysical data that the cross‐vane and the natural riffle induced spatially and temporally unique hyporheic extent and transport, and the cross‐vane created both spatially larger and temporally longer hyporheic flowpaths than the natural riffle. Tracer from the 4.67‐h injection was detected along flowpaths for 4.6 h at the cross‐vane and 4.2 h at the riffle. The spatial extent of the hyporheic zone at the cross‐vane was 12% larger than that at the riffle. We compare ER results of this study to vertical fluxes calculated from temperature profiles and conclude significant differences in the interpretation of hyporheic transport from these different field techniques. Results of this study demonstrate a high degree of heterogeneity in transport metrics at both the cross‐vane and the riffle and differences between the hyporheic flowpath networks at the two different features. Our results suggest that restoration structures may be capable of creating sufficient exchange flux and timescales of transport to achieve the same ecological functions as natural features, but engineering of the physical and biogeochemical environment may be necessary to realize these benefits.     

Mangrove‐forest evolution in a sediment‐rich estuarine system: opportunists or agents of geomorphic change?

The majority of the world's mangrove forests occur on mostly mineral sediments of fluvial origin. Two perspectives exist on the biogeomorphic development of these forests, i.e. that mangroves are opportunistic, with forest development primarily driven by physical processes, or alternatively that biophysical feedbacks strongly influence sedimentation and resulting geomorphology. On the Firth of Thames coast, New Zealand, we evaluate these two possible scenarios for sediment accumulation and forest development using high‐resolution sedimentary records and a detailed chronology of mangrove‐forest (Avicennia marina) development since the 1950s. Cores were collected along a shore‐normal transect of known elevation relative to mean sea level (MSL). Activities for lead‐210 (²¹⁰Pb), caesium‐137 (¹³⁷Cs) and beryllium‐7 (⁷Be), and sediment properties were analysed, with ²¹⁰Pb sediment accumulation rates (SARs), compensated for deep subsidence (~8 mm yr^?1) used as a proxy for elevation gain. At least four phases of forest development since the 1950s are recognized. An old‐growth forest developed by the late‐1970s with more recent seaward forest expansion thereafter. Excess ²¹⁰Pb profiles from the old‐growth forest exhibit relatively low SARs near the top (7–12 mm yr^?1) and bottom (10–22 mm yr^?1) of cores, separated by an interval of higher SARs (33–100 mm yr^?1). A general trend of increasing SAR over time characterizes the recent forest. Biogeomorphic evolution of the system is more complex than simple mudflat accretion/progradation and mangrove‐forest expansion. Surface‐elevation gain in the old‐growth forest displays an asymptotic trajectory, with a secondary depocentre developing on the seaward mudflat from the mid‐1970s. Two‐ to ten‐fold increases in ²¹⁰Pb SARs are unambiguously large and occurred years to decades before seedling recruitment, demonstrating that mangroves do not measurably enhance sedimentation over annual to decadal timescales. This suggests that mangrove‐forest development is largely dependent on physical processes, with forests occupying mudflats once they reach a suitable elevation in the intertidal. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of Sill Matrices in the Linear Model of Coregionalization

When using least squares to fit the linear model of coregionalization to multivariate geostatistical data, the sill matrices for the different regions must be estimated, subject to the constraint that they be non-negative definite. In 1992, Goulard and Voltz proposed and empirically examined an iterative algorithm for doing this. Although no proof was given for its convergence or for the uniqueness of the solution to the problem, the algorithm has subsequently been extensively and successfully used. In this paper, we prove that the minimization problem, in fact, has a unique solution and that the algorithm is guaranteed to converge to it from any starting point. We also discuss the effect of the starting point on the speed of convergence.     

An application of Spartan spatial random fields in environmental mapping: focus on automatic mapping capabilities

This paper investigates the potential of Spartan spatial random fields (SSRFs) in real-time mapping applications. The data set that we study focuses on the distribution of daily gamma dose rates over part of Germany. Our goal is to determine a Spartan spatial model from the data, and then use it to generate “predictive” maps of the radioactivity. In the SSRF framework, the spatial dependence is determined from sample functions that focus on short-range correlations. A recently formulated SSRF predictor is used to derive isolevel contour maps of the dose rates. The SSRF predictor is explicit. Moreover, the adjustments that it requires by the user are reduced compared to classical geostatistical methods. These features present clear advantages for an automatic mapping system. The performance of the SSRF predictor is evaluated by means of various cross-validation measures. The values of the performance measures are similar to those obtained by classical geostatistical methods. Application of the SSRF method to data that simulate a radioactivity release scenario is also discussed. Hot spots are detected and removed using a heuristic method. The extreme values that appear in the path of the simulated plume are not captured by the currently used Spartan spatial model. Modeling of the processes leading to extreme values can enhance the predictive capabilities of the spatial model, by incorporating physical information.     

Dense water formation in the Gulf of Lions shelf: Impact of atmospheric interannual variability and climate change

Dense water formed over the continental shelf and cascading down the slope is responsible for shelf-slope exchanges in many parts of the world ocean, and transports large amounts of sediment and organic matter into the deep ocean. Here we perform numerical modeling experiments to investigate the impact of atmospheric interannual variability and climate change on dense water formation over the Gulf of Lions shelf, in the Northwestern Mediterranean Sea. Results obtained for a 140 years eddy-permitting simulation (1960–2100) performed over the whole Mediterranean Sea under IPCC A2 scenario forcings are used to force a regional eddy-resolving model of the Northwestern Mediterranean Sea.