Machine Learning Predictions of pH in the Glacial Aquifer System,Northern USA

A boosted regression tree model was developed to predict pH conditions in three dimensions throughout the glacial aquifer system of the contiguous United States using pH measurements in samples from 18,386 wells and predictor variables that represent aspects of the hydrogeologic setting. Model results indicate that the carbonate content of soils and aquifer materials strongly controls pH and, when coupled with long flowpaths, results in the most alkaline conditions. Conversely, in areas where glacial sediments are thin and carbonate-poor, pH conditions remain acidic. At depths typical of drinking-water supplies, predicted pH >7.5—which is associated with arsenic mobilization—occurs more frequently than predicted pH <6—which is associated with water corrosivity and the mobilization of other trace elements. A novel aspect of this model was the inclusion of numerically based estimates of groundwater flow characteristics (age and flowpath length) as predictor variables. The sensitivity of pH predictions to these variables was consistent with hydrologic understanding of groundwater flow systems and the geochemical evolution of groundwater quality. The model was not developed to provide precise estimates of pH at any given location. Rather, it can be used to more generally identify areas where contaminants may be mobilized into groundwater and where corrosivity issues may be of concern to prioritize areas for future groundwater monitoring.     

Comparison of particle-tracking and lumped-parameter age-distribution models for evaluating vulnerability of production wells to contamination

Environmental age tracers have been used in various ways to help assess vulnerability of drinking-water production wells to contamination. The most appropriate approach will depend on the information that is available and that which is desired. To understand how the well will respond to changing nonpoint-source contaminant inputs at the water table, some representation of the distribution of groundwater ages in the well is needed. Such information for production wells is sparse and difficult to obtain, especially in areas lacking detailed field studies. In this study, age distributions derived from detailed groundwater-flow models with advective particle tracking were compared with those generated from lumped-parameter models to examine conditions in which estimates from simpler, less resource-intensive lumped-parameter models could be used in place of estimates from particle-tracking models. In each of four contrasting hydrogeologic settings in the USA, particle-tracking and lumped-parameter models yielded roughly similar age distributions and largely indistinguishable contaminant trends when based on similar conceptual models and calibrated to similar tracer data. Although model calibrations and predictions were variably affected by tracer limitations and conceptual ambiguities, results illustrated the importance of full age distributions, rather than apparent tracer ages or model mean ages, for trend analysis and forecasting.     

Esterification of vertebrate-like steroids in the eastern oyster (Crassostrea virginica)

The esterification of two model vertebrate steroid hormones - estradiol (E2) and dehidroepiandrosterone (DHEA) - was studied in the oyster Crassostrea virginica. The activity of acyl-CoA:steroid acyltransferase was characterized in microsomal fractions isolated from oyster digestive glands. The apparent Km and Vmax values changed with the fatty acid acyl-CoA used (C20:4, C18:2, C18:1, C16:1, C18:0 or C16:0), and were in the range of 9-17 microM, and 35-74 pmol/min/mg protein for E2, and in the range of 45-120 microM, and 30-182 pmol/min/mg protein for DHEA. Kinetic parameters were also assessed in gonadal tissue. The enzyme saturated at similar concentrations, although conjugation rates were lower than in digestive gland. Preliminary data shows that tributyltin (TBT) in the low microM range (1-50) strongly inhibits E2 and DHEA esterification, the esterification of E2 being more sensitive to inhibition than that of DHEA. Overall, results indicate that apolar conjugation occurs in oysters, in both digestive gland and gonads, at a very similar rate to mammals, suggesting that this is a well conserved conjugation pathway during evolution. Esterification, together with other mechanisms, can modulate endogenous steroid levels in C. virginica, and might be a target for endocrine disrupters, such as TBT.     

Comparison of age distributions estimated from environmental tracers by using binary-dilution and numerical models of fractured and folded karst: Shenandoah Valley of Virginia and West Virginia, USA

Measured concentrations of environmental tracers in spring discharge from a karst aquifer in the Shenandoah Valley, USA, were used to refine a numerical groundwater flow model. The karst aquifer is folded and faulted carbonate bedrock dominated by diffuse flow along fractures. The numerical model represented bedrock structure and discrete features (fault zones and springs). Concentrations of ³H, ³He, ⁴He, and CFC-113 in spring discharge were interpreted as binary dilutions of young (0–8  years) water and old (tracer-free) water. Simulated mixtures of groundwater are derived from young water flowing along shallow paths, with the addition of old water flowing along deeper paths through the model domain that discharge to springs along fault zones. The simulated median age of young water discharged from springs (5.7  years) is slightly older than the median age estimated from ³H/³He data (4.4  years). The numerical model predicted a fraction of old water in spring discharge (0.07) that was half that determined by the binary-dilution model using the ³H/³He apparent age and ³H and CFC-113 data (0.14). This difference suggests that faults and lineaments are more numerous or extensive than those mapped and included in the numerical model.     

On Modeling Weak Sinks in MODPATH

Regional groundwater flow systems often contain both strong sinks and weak sinks. A strong sink extracts water from the entire aquifer depth, while a weak sink lets some water pass underneath or over the actual sink. The numerical groundwater flow model MODFLOW may allow a sink cell to act as a strong or weak sink, hence extracting all water that enters the cell or allowing some of that water to pass. A physical strong sink can be modeled by either a strong sink cell or a weak sink cell, with the latter generally occurring in low‐resolution models. Likewise, a physical weak sink may also be represented by either type of sink cell. The representation of weak sinks in the particle tracing code MODPATH is more equivocal than in MODFLOW. With the appropriate parameterization of MODPATH, particle traces and their associated travel times to weak sink streams can be modeled with adequate accuracy, even in single layer models. Weak sink well cells, on the other hand, require special measures as proposed in the literature to generate correct particle traces and individual travel times and hence capture zones. We found that the transit time distributions for well water generally do not require special measures provided aquifer properties are locally homogeneous and the well draws water from the entire aquifer depth, an important observation for determining the response of a well to non‐point contaminant inputs.     

Meteoroid Engineering Model (MEM): A Meteoroid Model For The Inner Solar System

In an attempt to overcome some of the deficiencies of existing meteoroid models, NASA’s Space Environments and Effects (SEE) Program sponsored a 3 year research effort at the University of Western Ontario. The resulting understanding of the sporadic meteoroid environment – particularly the nature and distribution of the sporadic sources – were then incorporated into a new Meteoroid Engineering Model (MEM) by members of the Space Environments Team at NASA’s Marshall Space Flight Center. This paper discusses some of the revolutionary aspects of MEM which include (a) identification of the sporadic radiants with real sources of meteoroids, such as comets, (b) a physics-based approach which yields accurate fluxes and directionality for interplanetary spacecraft anywhere from 0.2 to 2.0 astronomical units (AU), and (c) velocity distributions obtained from theory and validated against observation. Use of the model, which gives penetrating fluxes and average impact speeds on the surfaces of a cube-like structure, is also described along with its current limitations and plans for future improvements.     

An appraisal of Indonesia’s immense peat carbon stock using national peatland maps: uncertainties and potential losses from conversion

Background

A large proportion of the world’s tropical peatlands occur in Indonesia where rapid conversion and associated losses of carbon, biodiversity and ecosystem services have brought peatland management to the forefront of Indonesia’s climate mitigation efforts. We evaluated peat volume from two commonly referenced maps of peat distribution and depth published by Wetlands International (WI) and the Indonesian Ministry of Agriculture (MoA), and used regionally specific values of carbon density to calculate carbon stocks.

Results

Peatland extent and volume published in the MoA maps are lower than those in the WI maps, resulting in lower estimates of carbon storage. We estimate Indonesia’s total peat carbon store to be within 13.6 GtC (the low MoA map estimate) and 40.5 GtC (the high WI map estimate) with a best estimate of 28.1 GtC: the midpoint of medium carbon stock estimates derived from WI (30.8 GtC) and MoA (25.3 GtC) maps. This estimate is about half of previous assessments which used an assumed average value of peat thickness for all Indonesian peatlands, and revises the current global tropical peat carbon pool to 75 GtC. Yet, these results do not diminish the significance of Indonesia’s peatlands, which store an estimated 30% more carbon than the biomass of all Indonesian forests. The largest discrepancy between maps is for the Papua province, which accounts for 62–71% of the overall differences in peat area, volume and carbon storage. According to the MoA map, 80% of Indonesian peatlands are <300 cm thick and thus vulnerable to conversion outside of protected areas according to environmental regulations. The carbon contained in these shallower peatlands is conservatively estimated to be 10.6 GtC, equivalent to 42% of Indonesia’s total peat carbon and about 12 years of global emissions from land use change at current rates.

Conclusions

Considering the high uncertainties in peatland extent, volume and carbon storage revealed in this assessment of current maps, a systematic revision of Indonesia’s peat maps to produce a single geospatial reference that is universally accepted would improve national peat carbon storage estimates and greatly benefit carbon cycle research, land use management and spatial planning.

     

Esterification of vertebrate-like steroids in the eastern oyster (Crassostrea virginica)

The esterification of two model vertebrate steroid hormones – estradiol (E2) and dehidroepiandrosterone (DHEA) – was studied in the oyster Crassostrea virginica. The activity of acyl-CoA:steroid acyltransferase was characterized in microsomal fractions isolated from oyster digestive glands. The apparent K_m and V_max values changed with the fatty acid acyl-CoA used (C20:4, C18:2, C18:1, C16:1, C18:0 or C16:0), and were in the range of 9–17 μM, and 35–74 pmol/min/mg protein for E2, and in the range of 45–120 μM, and 30–182 pmol/min/mg protein for DHEA. Kinetic parameters were also assessed in gonadal tissue. The enzyme saturated at similar concentrations, although conjugation rates were lower than in digestive gland. Preliminary data shows that tributyltin (TBT) in the low μM range (1–50) strongly inhibits E2 and DHEA esterification, the esterification of E2 being more sensitive to inhibition than that of DHEA. Overall, results indicate that apolar conjugation occurs in oysters, in both digestive gland and gonads, at a very similar rate to mammals, suggesting that this is a well conserved conjugation pathway during evolution. Esterification, together with other mechanisms, can modulate endogenous steroid levels in C. virginica, and might be a target for endocrine disrupters, such as TBT.