A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Marine Macrophyte Wrack Inputs and Dissolved Nutrients in Beach Sands

We investigated the role of sandy beaches in nearshore nutrient cycling by quantifying macrophyte wrack inputs and examining relationships between wrack accumulation and pore water nutrients during the summer dry season. Macrophyte inputs, primarily giant kelp Macrocystis pyrifera, exceeded 2.3 kg m⁻¹ day⁻¹. Mean wrack biomass varied 100-fold among beaches (range = 0.41 to 46.43 kg m⁻¹). Mean concentrations of dissolved inorganic nitrogen (DIN), primarily NO_x⁻-N, and dissolved organic nitrogen (DON) in intertidal pore water varied significantly among beaches (ranges = 1 to 6,553 μM and 7 to 2,006 μM, respectively). Intertidal DIN and DON concentrations were significantly correlated with wrack biomass. Surf zone concentrations of DIN were also strongly correlated with wrack biomass and with intertidal DIN, suggesting export of nutrients from re-mineralized wrack. Our results suggest beach ecosystems can process and re-mineralize substantial organic inputs and accumulate dissolved nutrients, which are subsequently available to nearshore waters and primary producers.     

Long‐term sediment yield from a small catchment in southern Brazil affected by land use and soil management changes

Sediments produced from eroding cultivated land can cause on‐site and off‐site effects that cause considerable economic and social impacts. Despite the importance of soil conservation practices (SCP) for the control of soil erosion and improvements in soil hydrological functions, limited information is available regarding the effects of SCP on sediment yield (SY) at the catchment scale. This study aimed to investigate the long‐term relationships between SY and land use, soil management, and rainfall in a small catchment. To determine the effects of anthropogenic and climatic factors on SY, rainfall, streamflow, and suspended sediment concentration were monitored at 10‐min intervals for 14 years (2002–2016), and the land use and soil management changes were surveyed annually. Using a statistical procedure to separate the SY effects of climate, land use, and soil management, we observed pronounced temporal effects of land use and soil management changes on SY. During the first 2 years (2002–2004), the land was predominantly cultivated with tobacco under a traditional tillage system (no cover crops and ploughed soil) using animal traction. In that period, the SY reached approximately 400 t·km^?2·year^?1. From 2005 to 2009, a soil conservation programme introduced conservation tillage and winter cover crops in the catchment area, which lowered the SY to 50 t·km^?2·year^?1. In the final period (2010–2016), the SCP were partially abandoned by farmers, and reforested areas increased, resulting in an SY of 150 t·km^?2·year^?1. This study also discusses the factors associated with the failure to continue using SCP, including structural support and farmer attitudes.     

The information system of the French Peatland Observation Service: Service National d'Observation Tourbières – A valuable tool to assess the impact of global changes on the hydrology and biogeochemistry of temperate peatlands through long term monitoring

Mitigating and adapting to global changes requires a better understanding of the response of the Biosphere to these environmental variations. Human disturbances and their effects act in the long term (decades to centuries) and consequently, a similar time frame is needed to fully understand the hydrological and biogeochemical functioning of a natural system. To this end, the ‘Centre National de la Recherche Scientifique’ (CNRS) promotes and certifies long-term monitoring tools called national observation services or ‘Service National d'Observation’ (SNO) in a large range of hydrological and biogeochemical systems (e.g., cryosphere, catchments, aquifers). The SNO investigating peatlands, the SNO ‘Tourbières’, was certified in 2011 ( https://www.sno-tourbieres.cnrs.fr/ ). Peatlands are mostly found in the high latitudes of the northern hemisphere and French peatlands are located in the southern part of this area. Thus, they are located in environmental conditions that will occur in northern peatlands in coming decades or centuries and can be considered as sentinels. The SNO Tourbières is composed of four peatlands: La Guette (lowland central France), Landemarais (lowland oceanic western France), Frasne (upland continental eastern France) and Bernadouze (upland southern France). Thirty target variables are monitored to study the hydrological and biogeochemical functioning of the sites. They are grouped into four datasets: hydrology, fluvial export of organic matter, greenhouse gas fluxes and meteorology/soil physics. The data from all sites follow a common processing chain from the sensors to the public repository. The raw data are stored on an FTP server. After operator or automatic processing, data are stored in a database, from which a web application extracts the data to make them available ( https://data-snot.cnrs.fr/data-access/ ). Each year at least, an archive of each dataset is stored in Zenodo, with a digital object identifier (DOI) attribution ( https://zenodo.org/communities/sno_tourbieres_data/ ).     

Late Quaternary vegetation and climate dynamics in central-eastern Brazil: insights from a ~35k cal a bp peat record in the Cerrado biome

The late Quaternary evolution of central-eastern Brazil has been under-researched. Questions remain as to the origin of the Cerrado, a highly endangered biome, and other types of vegetation, such as the Capões – small vegetation islands of semi-deciduous and mountain forests. We investigated the factors that influenced the expansion and contraction of the Cerrado and Capões during the late Quaternary (last ~35 ka), using a multi-proxy approach: stable isotopes (δ¹³C, δ¹⁵N), geochemistry, pollen and multivariate statistics derived from a peat core (Pinheiro mire, Serra do Espinhaço Meridional). Five major shifts in precipitation, temperature, vegetation and landscape stability occurred at different timescales. Our study revealed that changes in the South Atlantic Convergence Zone (SACZ) seem to have been coeval with these shifts: from the Late Glacial Maximum to mid-Holocene the SACZ remained near (~29.6 to ~16.5k cal a bp ) and over (~16.5 to ~6.1 k cal a bp ) the study area, providing humidity to the region. This challenges previous research which suggested that climate was drier for this time period. At present, the Capões are likely to be a remnant of a more humid climate; meanwhile, the Cerrado biome seems to have stablished in the late Holocene, after ~3.1 k cal a bp .     

Defining the Potential of Nanoscale Re‐Os Isotope Systematics Using Atom Probe Microscopy

Atom probe microscopy (APM) is a relatively new in situ tool for measuring isotope fractions from nanoscale volumes (< 0.01 μm³). We calculate the theoretical detectable difference of an isotope ratio measurement result from APM using counting statistics of a hypothetical data set to be ± 4δ or 0.4% (2s). However, challenges associated with APM measurements (e.g., peak ranging, hydride formation and isobaric interferences), result in larger uncertainties if not properly accounted for. We evaluate these factors for Re‐Os isotope ratio measurements by comparing APM and negative thermal ionisation mass spectrometry (N‐TIMS) measurement results of pure Os, pure Re, and two synthetic Re‐Os‐bearing alloys from Schwander et al. (2015, Meteoritics and Planetary Science, 50, 893) [the original metal alloy (HSE) and alloys produced by heating HSE within silicate liquid (SYN)]. From this, we propose a current best practice for APM Re‐Os isotope ratio measurements. Using this refined approach, mean APM and N‐TIMS ¹⁸⁷Os/¹⁸⁹Os measurement results agree within 0.05% and 2s (pure Os), 0.6–2% and 2s (SYN) and 5–10% (HSE). The good agreement of N‐TIMS and APM ¹⁸⁷Os/¹⁸⁹Os measurements confirms that APM can extract robust isotope ratios. Therefore, this approach permits nanoscale isotope measurements of Os‐bearing alloys using the Re‐Os geochronometer that could not be measured by conventional measurement principles.     

Modelling in-situ upgrading of heavy oil using operator splitting method

The in-situ upgrading (ISU) of bitumen and oil shale is a very challenging process to model numerically because of the large number of components that need to be modelled using a system of equations that are both highly non-linear and strongly coupled. Operator splitting methods are one way of potentially improving computational performance. Each numerical operator in a process is modelled separately, allowing the best solution method to be used for the given numerical operator. A significant drawback to the approach is that decoupling the governing equations introduces an additional source of numerical error, known as the splitting error. The best splitting method for modelling a given process minimises the splitting error whilst improving computational performance compared to a fully implicit approach. Although operator splitting has been widely used for the modelling of reactive-transport problems, it has not yet been applied to the modelling of ISU. One reason is that it is not clear which operator splitting technique to use. Numerous such techniques are described in the literature and each leads to a different splitting error. While this error has been extensively analysed for linear operators for a wide range of methods, the results cannot be extended to general non-linear systems. It is therefore not clear which of these techniques is most appropriate for the modelling of ISU. In this paper, we investigate the application of various operator splitting techniques to the modelling of the ISU of bitumen and oil shale. The techniques were tested on a simplified model of the physical system in which a solid or heavy liquid component is decomposed by pyrolysis into lighter liquid and gas components. The operator splitting techniques examined include the sequential split operator (SSO), the Strang-Marchuk split operator (SMSO) and the iterative split operator (ISO). They were evaluated on various test cases by considering the evolution of the discretization error as a function of the time-step size compared with the results obtained from a fully implicit simulation. We observed that the error was least for a splitting scheme where the thermal conduction was performed first, followed by the chemical reaction step and finally the heat and mass convection operator (SSO-CKA). This method was then applied to a more realistic model of the ISU of bitumen with multiple components, and we were able to obtain a speed-up of between 3 and 5.