Geochemical Implications of Brine Leakage into Freshwater Aquifers

CO₂ injection into deep saline formations as a way to mitigate climate change raises concerns that leakage of saline waters from the injection formations will impact water quality of overlying aquifers, especially underground sources of drinking water (USDWs). This paper aims to characterize the geochemical composition of deep brines, with a focus on constituents that pose a human health risk and are regulated by the U.S. Environmental Protection Agency (USEPA). A statistical analysis of the NATCARB brine database, combined with simple mixing model calculations, show total dissolved solids and concentrations of chloride, boron, arsenic, sulfate, nitrate, iron and manganese may exceed plant tolerance or regulatory levels. Twelve agricultural crops evaluated for decreased productivity in the event of brine leakage would experience some yield reduction due to increased TDS at brine‐USDW ratios of < 0.1, and a 50% yield reduction at < 0.2 brine‐USDW ratio. A brine‐USDW ratio as low as 0.004 may result in yield reduction in the most sensitive crops. The USEPA TDS secondary standard is exceeded at a brine fraction of approximately 0.002. To our knowledge, this is the first study to consider agricultural impacts of brine leakage, even though agricultural withdrawals of groundwater in the United States are almost three times higher than public and domestic withdrawals.     

Multistage growth of Fe–Mg–carpholite and Fe–Mg–chloritoid,from field evidence to thermodynamic modelling

We provide new insights into the prograde evolution of HP/LT metasedimentary rocks on the basis of detailed petrologic examination, element-partitioning analysis, and thermodynamic modelling of well-preserved Fe–Mg–carpholite- and Fe–Mg–chloritoid-bearing rocks from the Afyon Zone (Anatolia). We document continuous and discontinuous compositional (ferromagnesian substitution) zoning of carpholite (overall X _Mg = 0.27–0.73) and chloritoid (overall X _Mg = 0.07–0.30), as well as clear equilibrium and disequilibrium (i.e., reaction-related) textures involving carpholite and chloritoid, which consistently account for the consistent enrichment in Mg of both minerals through time, and the progressive replacement of carpholite by chloritoid. Mg/Fe distribution coefficients calculated between carpholite and chloritoid vary widely within samples (2.2–20.0). Among this range, only values of 7–11 correlate with equilibrium textures, in agreement with data from the literature. Equilibrium phase diagrams for metapelitic compositions are calculated using a newly modified thermodynamic dataset, including most recent data for carpholite, chloritoid, chlorite, and white mica, as well as further refinements for Fe–carpholite, and both chloritoid end-members, as required to reproduce accurately petrologic observations (phase relations, experimental constraints, Mg/Fe partitioning). Modelling reveals that Mg/Fe partitioning between carpholite and chloritoid is greatly sensitive to temperature and calls for a future evaluation of possible use as a thermometer. In addition, calculations show significant effective bulk composition changes during prograde metamorphism due to the fractionation of chloritoid formed at the expense of carpholite. We retrieve P–T conditions for several carpholite and chloritoid growth stages (1) during prograde stages using unfractionated, bulk-rock XRF analyses, and (2) at peak conditions using compositions fractionated for chloritoid. The P–T paths reconstructed for the Kütahya and Afyon areas shed light on contrasting temperature conditions for these areas during prograde and peak stages.     

Simplifying the interpretation of continuous time models for spatio-temporal networks

Autoregressive and moving average models for temporally dynamic networks treat time as a series of discrete steps which assumes even intervals between data measurements and can introduce bias if this assumption is not met. Using real and simulated data from the London Underground network, this paper illustrates the use of continuous time multilevel models to capture temporal trajectories of edge properties without the need for simultaneous measurements, along with two methods for producing interpretable summaries of model results. These including extracting ‘features’ of temporal patterns (e.g. maxima, time of maxima) which have utility in understanding the network properties of each connection and summarising whole-network properties as a continuous function of time which allows estimation of network properties at any time without temporal aggregation of non-simultaneous measurements. Results for temporal pattern features in the response variable were captured with reasonable accuracy. Variation in the temporal pattern features for the exposure variable was underestimated by the models. The models showed some lack of precision. Both model summaries provided clear ‘real-world’ interpretations and could be applied to data from a range of spatio-temporal network structures (e.g. rivers, social networks). These models should be tested more extensively in a range of scenarios, with potential improvements such as random effects in the exposure variable dimension.

     

An estimate of the influence of loading effects on tectonic velocities in the Pyrenees

Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.     

Tracking emission sources of sulfur and elemental carbon in Hong Kong/Pearl River Delta region

The Tagged Species Source Apportionment (TSSA) algorithm is applied to study contributions to pollutants PM sulfate, SO2 and elemental carbon (EC) in Hong Kong by emitting sectors as well as non-emitting sources within and beyond the Pearl River Delta (PRD) region. We look at three areas of Hong Kong: western new towns, central downtown, and eastern countryside. Pollutant transport from beyond the PRD influenced all three areas during January and October 2004 but the source sectors impacting the three areas were different. Specifically, power plant SO2 from Hong Kong and Shenzhen, and vehicle EC from Shenzhen contribute to 11?~?66?% of SO2 concentration, and 33?~?75?% of EC concentration in the western new towns, respectively. Ships in and around Hong Kong contribute 8?~?56?% to the sulfate concentration in the downtown area, while local moving vehicles take up 28?~?81?% of the EC concentration there. The study also shows that diurnal variation of planetary boundary layer (PBL) causes day-night difference in SO2 by as much as 50?%. In addition, 13?~?38?% of all SO2 emissions have been converted to PM sulfate.     

Quantification of Storm-Induced Bathymetric Change in a Back-Barrier Estuary

Geomorphology is a fundamental control on ecological and economic function of estuaries. However, relative to open coasts, there has been little quantification of storm-induced bathymetric change in back-barrier estuaries. Vessel-based and airborne bathymetric mapping can cover large areas quickly, but change detection is difficult because measurement errors can be larger than the actual changes over the storm timescale. We quantified storm-induced bathymetric changes at several locations in Chincoteague Bay, Maryland/Virginia, over the August 2014 to July 2015 period using fixed, downward-looking altimeters and numerical modeling. At sand-dominated shoal sites, measurements showed storm-induced changes on the order of 5 cm, with variability related to stress magnitude and wind direction. Numerical modeling indicates that the predominantly northeasterly wind direction in the fall and winter promotes southwest-directed sediment transport, causing erosion of the northern face of sandy shoals; southwesterly winds in the spring and summer lead to the opposite trend. Our results suggest that storm-induced estuarine bathymetric change magnitudes are often smaller than those detectable with methods such as LiDAR. More precise fixed-sensor methods have the ability to elucidate the geomorphic processes responsible for modulating estuarine bathymetry on the event and seasonal timescale, but are limited spatially. Numerical modeling enables interpretation of broad-scale geomorphic processes and can be used to infer the long-term trajectory of estuarine bathymetric change due to episodic events, when informed by fixed-sensor methods.     

Ontogenetic shifts and temporal changes in the trophic patterns of the deep‐sea red shrimp,Aristaeomorpha foliacea (Decapods: Aristeidae), in the Eastern Ionian Sea (Eastern Mediterranean)

The purpose of this paper is to provide the first detailed data concerning the diet and feeding activity of the giant red shrimp, Aristaeomorpha foliacea, in the Eastern Ionian Sea (Eastern Mediterranean), in relation to season, size class and sex. Feeding activity in A. foliacea was intense, based on its low vacuity index and high prey diversity, with a diet dominated by mesopelagic prey and less frequent occurrence of benthic taxa. Giant red shrimp displayed a highly diversified diet that exhibited slight seasonal fluctuations. The diets of both sexes consisted of 60 different prey categories belonging chiefly to three groups: crustaceans (e.g. decapods, such as Plesionika spp. and Pasiphaeidae, amphipods), cephalopods (mainly Enoploteuthidae) and fishes (Myctophidae, Macrouridae). These three prey categories accounted for 72–82% of the relative abundance and total occurrence for males and 70–88% for females, respectively. Variation in food availability, as well as increased energy demands related to gonad development and breeding activity, appear to be critical factors driving temporal changes in feeding strategy. Feeding activity increased during spring and summer, which coincides with reproductive activities (mating, gonad maturation, egg‐laying). Females seem to be more active predators than males, consuming prey with greater swimming ability. However, ontogenetic shifts in diet were also apparent, despite high dietary overlap among small, medium and large females. Large individuals, which are more efficient predators, selected highly mobile prey (e.g. fishes), whereas small individuals consumed low‐mobility prey (e.g. copepods, ostracods, tanaids and sipunculans).     

Streamflow response to the rapid retreat of a lake‐calving glacier

There has been increasing attention over the last decade to the potential effects of glacier retreat on downstream discharge and aquatic habitat. This study focused on streamflow variability downstream of Bridge Glacier in the southern Coast Mountains of BC between 1979 and 2014, prior to and during a period in which the glacier experienced enhanced calving and rapid retreat across a lake‐filled basin. Here we combined empirical trend detection and a conceptual‐parametric hydrological model to address the following hypotheses: (1) streamflow trends in late summer and early autumn should reflect the opposing influences of climatic warming (which would tend to increase unit‐area meltwater production) and the reduction in glacier area (which would tend to reduce the total volume of meltwater generated), and (2) winter streamflow should increase because of displacement of lake water as ice flows past the grounding line and calves into the lake basin. In relation to the first hypothesis, we found no significant trends in monthly discharge during summer. However, applying regression analysis to account for air temperature and precipitation variations, weak but statistically significant negative trends were detected for August and melt season discharge. The HBV‐EC model was applied using time‐varying glacier cover, as derived from Landsat imagery. Relative to simulations based on constant glacier extent, model results indicated that glacier recession caused a decline in mean monthly streamflow of 9% in August and 11% in September. These declines in late‐summer streamflow are consistent with the results from our empirical analysis. The second hypothesis is supported by the finding of positive trends for December, January, and February discharge. Despite the modelled declines in late‐summer mean monthly streamflow, recorded discharge data exhibited neither positive nor negative trends during the melt season, suggesting that Bridge Glacier may currently be at or close to the point of peak water. Further analysis of the impact of lake‐terminating glaciers on downstream discharge is needed to refine the peak water model. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of the seismicity in the central part of the Pyrenees (France), and tectonic implications

Two temporary seismological networks have been set up in 2000 and 2002 in the central part of the Pyrenees, in a region, which appears as a transition between two domains where both the seismic activity and the tectonic regime are different. Together with the permanent networks, they allowed us to obtain precise hypocenter locations for more than 400 events with local magnitudes ranging from 1.5 to 4.6, as well as 30 new focal mechanisms. The seismicity is distributed in several clusters, which are not located along the North Pyrenean Fault, considered as the major tectonic accident resulting from the suture of the Iberian and Eurasian convergent plates when the range formed. Several small fault segments dipping to the north are identified. The maximum focal depth varies from 10 to 20 km, with variations which are roughly parallel to those of the Moho, indicating a thickening of the seismogenic layer to the east of the studied area. The obtained focal solutions reveal a predominance of normal faulting to the West and reverse faulting to the East, with strike-slip motions in between. The largest fault segment to the East, with a length of about 20 km, could possibly be related to a large historical event which occurred in 1660, with intensity IX, close to cities which have become since then important touristic centres.     

Gen*: a generic toolkit to generate spatially explicit synthetic populations

Agent-based models tend to integrate more and more data that can deeply impact their outcomes. Among these data, the ones that deal with agent attributes and localization are particularly important, but are very difficult to collect. In order to tackle this issue, we propose a complete generic toolkit called Gen* dedicated to generating spatially explicit synthetic populations from global (census and GIS) data. This article focuses on the localization methods provided by Gen* that are based on regression, geometrical constraints and spatial distributions. The toolkit is applied for a case study concerning the generation of the population of Rouen (France) and shows the capabilities of Gen* regarding population spatialization.