Characterizing sources of groundwater to a tropical coastal lagoon in a karstic area using radium isotopes and water chemistry

Radium isotopes (²²³Ra, ²²⁴Ra, ²²⁶Ra, and ²²⁸Ra) and water chemistry were used to identify two chemically distinct sources of submarine groundwater discharge (SGD) in Celestún Lagoon, Yucatán, Mexico. Low salinity groundwater discharging from springs within the lagoon has previously been identified and extensively sampled for nutrient concentrations. However, a second type of groundwater discharging into the lagoon was detected during this study using radium isotope activity measurements. This second type of groundwater is characterized by moderate salinities (within the range of lagoon salinities) and very elevated radium activities in comparison to the low salinity groundwater, mixed lagoon water, and seawater. Further analysis showed that the two types of groundwater also have distinct chloride, strontium, and sulfate ratios, along with slightly different nutrient concentrations. Groundwater discharge occurs through large and small springs scattered throughout the lagoon, and both types of groundwater were detected discharging from one of the larger springs. The relative proportions of low salinity groundwater and brackish high radium groundwater varied over the tidal cycle. In order to better understand the relative contributions of each type of groundwater to the lagoon, a three end-member mixing model based on the distinct chemical and isotopic compositions of both types of groundwater and of seawater was used to estimate the distribution of each water type throughout the lagoon in different seasons. This study suggests that substantial groundwater discharge to the lagoon can occur during both dry and rainy seasons. The presence of two groundwater sources has implications for monitoring and protection of the Celestún Lagoon Biosphere Reserve, since the two sources may have different susceptibilities to anthropogenic contamination depending on their respective recharge area and recharge rates.     

4DVAR data assimilation in the Intra-Americas Sea with the Regional Ocean Modeling System (ROMS)

We present the background, development, and preparation of a state-of-the-art 4D variational (4DVAR) data assimilation system in the Regional Ocean Modeling System (ROMS) with an application in the Intra-Americas Sea (IAS). This initial application with a coarse model shows the efficacy of the 4DVAR methodology for use within complex ocean environments, and serves as preparation for deploying an operational, real-time assimilation system onboard the Royal Caribbean Cruise Lines ship Explorer of the Seas. Assimilating satellite sea surface height and temperature observations with in situ data from the ship in 14 day cycles over 2 years from January 2005 through March 2007, reduces the observation-model misfit by over 75%. Using measures of the Loop Current dynamics, we show that the assimilated solution is consistent with observed statistics.     

Transport and retention of a conservative tracer in an isolated creek–marsh system

A study of tracer transport and retention in a small tidal creek and marsh system located in the southeastern US was conducted using a three-dimensional hydrodynamic model and data from a dye tracer release. The model simulated tidally driven flow, including inundation and drying out of the marshes and the dispersal of the dye tracer. Flow measurements in the tidal creek showed that the simulations appeared to generally duplicate the tidally driven flow into and out of the tidal creeks and marshes. The dye tracer experiment was conducted to test the hydrodynamic model's ability to simulate dispersal of a point release of pollutants into the creek during an incoming tide. The simulations reproduced much of the experimental measurements, but bathymetric errors and lack of resolution of the smallest arms of the tidal creeks affected the ability to faithfully reproduce the initial peak in measured dye tracer concentrations at a sampling location far from the boundary. The lack of the smallest tidal creeks led to some trapping of water and dye in the marshes. Two independent estimates of flushing time yielded values between 1.6 and 2.4 days. An uncertainty analysis indicates that model simulations are sensitive to variations in parameters such as water depth and marsh grass density. On the other hand, omission of the smallest tidal creeks in simulations may be partially offset by decreasing marsh grass density. Further improvements must rely on more accurate and detailed bathymetric data that resolves the smaller arms of the tidal creeks and includes quantitative information about the density and distribution of marsh grass.     

Hydrocarbon‐Rich Groundwater above Shale‐Gas Formations: A Karoo Basin Case Study

Horizontal drilling and hydraulic fracturing have enhanced unconventional hydrocarbon recovery but raised environmental concerns related to water quality. Because most basins targeted for shale‐gas development in the USA have histories of both active and legacy petroleum extraction, confusion about the hydrogeological context of naturally occurring methane in shallow aquifers overlying shales remains. The Karoo Basin, located in South Africa, provides a near‐pristine setting to evaluate these processes, without a history of conventional or unconventional energy extraction. We conducted a comprehensive pre‐industrial evaluation of water quality and gas geochemistry in 22 groundwater samples across the Karoo Basin, including dissolved ions, water isotopes, hydrocarbon molecular and isotopic composition, and noble gases. Methane‐rich samples were associated with high‐salinity, NaCl‐type groundwater and elevated levels of ethane, ⁴He, and other noble gases produced by radioactive decay. This endmember displayed less negative δ¹³C‐CH₄ and evidence of mixing between thermogenic natural gases and hydrogenotrophic methane. Atmospheric noble gases in the methane‐rich samples record a history of fractionation during gas‐phase migration from source rocks to shallow aquifers. Conversely, methane‐poor samples have a paucity of ethane and ⁴He, near saturation levels of atmospheric noble gases, and more negative δ¹³C‐CH₄; methane in these samples is biogenic and produced by a mixture of hydrogenotrophic and acetoclastic sources. These geochemical observations are consistent with other basins targeted for unconventional energy extraction in the USA and contribute to a growing data base of naturally occurring methane in shallow aquifers globally, which provide a framework for evaluating environmental concerns related to unconventional energy development (e.g., stray gas).     

Tidal pumping and biogeochemical processes: Dissolution within the tidal capillary fringe of eogenetic coastal carbonates

Growing evidence suggests microbial respiration of dissolved organic carbon (DOC) may be a principal driver of subsurface dissolution and cave formation in eogenetic carbonate rock. Analyses of samples of vadose zone gasses, and geochemical and hydrological data collected from shallow, uncased wells on San Salvador Island, Bahamas, suggest tidally varying water tables may help fuel microbial respiration and dissolution through oxygenation. Respiration of soil organic carbon transported to water tables generates dysaerobic to anaerobic groundwater, limiting aerobic microbial processes. Positive correlations of carbon dioxide (CO₂), radon-222 (²²²Rn) and water table elevation indicate, however, that tidal pumping of water tables pulls atmospheric air that is rich in oxygen, and low in CO₂ and ²²²Rn, into contact with the tidal capillary fringe during falling tides. Ratios of CO₂ and O₂ in vadose gas relative to the atmosphere indicate this atmospheric oxygen fuels respiration within newly-exposed, wetted bedrock. Deficits of expected CO₂ relative to O₂ concentrations indicate some respired CO₂ is likely removed by carbonate mineral dissolution. Tidal pumping also appears capable of transferring oxygen to the freshwater lens, where it could also contribute to respiration and dissolution; dissolved oxygen concentrations at the water table are at least 5% saturated and decline to anaerobic conditions 1–2 m below. Our results demonstrate how tidal pumping of air to vadose zones can drive mineral dissolution reactions that are focused near water tables and may contribute to the formation of laterally continuous vuggy horizons and potentially caves. © 2020 John Wiley & Sons, Ltd.     

Distribution of diffuse flow megafauna in two sites on the Eastern Lau Spreading Center, Tonga

Hydrothermal vent environments are characterized by large gradients of toxic chemicals and high temperatures, which play a significant role in defining species’ distributions. We used high-resolution imagery and spatially explicit in-situ physico-chemical measurements analyzed within a Geographic Information System (GIS) in order to characterize the spatial relations among different groups of megafauna, temperature, and chemistry within two discrete vent communities (40 and 50 m²) on the Eastern Lau Spreading Center (ELSC). Chemical (sulfide and O₂ concentrations) and temperature data were obtained from approximately 75 different locations within each community using in-situ instruments. All data were integrated into a GIS, which served as a visualization tool and enabled the data to be analyzed in a spatial context. Our results confirm the importance of abiotic variables in defining the distributions of some fauna and elucidate several biological associations that are consistent between the two communities. The provannid snail, Alviniconcha spp., appears to actively avoid temperatures above 32–46 °C and/or sulfide concentrations exceeding approximately 260 μM. Slightly higher average sulfide concentrations and temperatures were measured among aggregations of Ifremeria nautilei compared to aggregations of the mussel Bathymodiolus brevior; however, the presence of mixed aggregations of the two species indicates an overlap in requirements. The brachyuran crab, Austinograea spp., was consistently observed directly on symbiont-containing species, particularly Alviniconcha spp. The solitary snail, Eosipho desbruyeresi, was rarely observed on biological substrata, but was often (60% of its population at the most active site) within 5 cm of symbiont-containing fauna, indicating a tolerance and preference for proximity to areas of high productivity. Densities and coverage of species differed substantially between the two communities despite high species overlap. Symbiont-containing species covered much larger areas at the more hydrothermally active site, ABE1, while shrimp and anemones occurred in relatively higher densities within the less-active site, TM1. This is the first study to thoroughly characterize realized distributions of megafauna at vent sites along the ELSC.     

Speleothem and glacier records of latest Pleistocene–early Holocene climate change in the western North American interior

Rapid and high-magnitude North Atlantic climate oscillations following the Last Glacial Maximum have been correlated to climate change events in western North America. However, the strength of teleconnections between the North Atlantic and the interior of western North America remains poorly understood. We present a U-series calibrated speleothem record from Timpanogos Cave National Monument, located at 2040 m asl in the Wasatch Mountains of Utah, spanning 13.5–10.6 ka. Additionally, we carried out a climate reconstruction for a coeval glacier advance in the Wind River Range of Wyoming. Our results indicate that between 13.5 and 12.8 ka, the Wasatch was probably first cool and dry and then warmed. After 12.8 ka, our record suggests cool and/or wetter conditions followed by reduced moisture until 11.8 ka, followed by an early Holocene wet period. The Timpanogos record exhibits few similarities with those from the North Atlantic. Climate reconstructions of the Titcomb Basin glacier suggest modest temperature depressions relative to modern (<−3 °C) were necessary to sustain the glacier with a moderate increase in precipitation (>150%). The high-altitude speleothem record presented here provides a valuable basis for understanding latest Pleistocene–early Holocene glacial episodes in western North America.     

Scientific briefing: quantifying streambed heat advection associated with groundwater–surface water interactions

Stream thermal regimes are controlled by the interactions of external and internal energy fluxes with the water in the channel. Solar radiation is typically the dominant driver of stream water temperature, but streambed heat fluxes can be important in forested headwater streams. Past studies have presented seemingly disparate formulae for quantifying streambed heat advection from upwelling groundwater. This note details the sources of the differences in these alternative formulations. The equations illustrate the difficulties of attempting to isolate the thermal influence of groundwater–surface water interactions and highlight future research opportunities. Copyright © 2015 John Wiley & Sons, Ltd.     

Suitability of North American Regional Reanalysis (NARR) output for hydrologic modelling and analysis in mountainous terrain

Meteorological observations at high elevations in mountainous regions are often lacking. One opportunity to fill this data gap is through the use of downscaled output from weather reanalysis models. In this study, we tested the accuracy of downscaled output from the North American Regional Reanalysis (NARR) against high‐elevation surface observations at four ridgetop locations in the southern Coast Mountains of British Columbia, Canada. NARR model output was downscaled to the surface observation locations through three‐dimensional interpolation for air temperature, vapour pressure and wind speed and two‐dimensional interpolation for radiation variables. Accuracy was tested at both the 3‐hourly and daily time scales. Air temperature displayed a high level of agreement, especially at the daily scale, with root mean square error (RMSE) values ranging from 0.98 to 1.21 °C across all sites. Vapour pressure downscaling accuracy was also quite high (RMSE of 0.06 to 0.11 hPa) but displayed some site specific bias. Although NARR overestimated wind speed, there were moderate to strong linear relations (r² from 0.38 to 0.84 for daily means), suggesting that the NARR output could be used as an index and bias‐corrected. NARR output reproduced the seasonal cycle for incoming short‐wave radiation, with Nash–Sutcliffe model efficiencies ranging from 0.78 to 0.87, but accuracy suffered on days with cloud cover, resulting in a positive bias and RMSE ranged from 42 to 46 Wm^{? 2}. Although fewer data were available, incoming long‐wave radiation from NARR had an RMSE of 19 Wm^{? 2} and outperformed common methods for estimating incoming long‐wave radiation. NARR air temperature showed potential to assist in hydrologic analysis and modelling during an atmospheric river storm event, which are characterized by warm and wet air masses with atypical vertical temperature gradients. The incorporation of a synthetic NARR air temperature station to better represent the higher freezing levels resulted in increased predicted peak flows, which better match the observed run‐off during the event. Copyright © 2016 John Wiley & Sons, Ltd.