Impact of multi-annual drought on streamflow and habitat in coastal California salmonid streams

The 2012–2015 drought in north-central coastal California ranks among the three most prolonged periods of below-median annual rainfall in the past 65 years. In three critical coho salmon streams, summer baseflow was less each additional dry year; streams with summer flow early in the drought had no flow for more than two months in latter years. By the third dry year, summer discharge was 1–5% of recent wet-type years, and 10–20% of the first dry year. Multiannual drought also caused increased dry channel conditions: the percentage of flowing channel reduced from 28 to 55% from the first to the third dry years among three study streams. In the first year following drought, dry-season streamflow resembled early to-mid-drought conditions, while in the second, it approached pre-drought discharge. This multiannual drought foreshadows how multi-annual drought predicted under future climate scenarios may affect critical salmonid streams later this century.     

A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

The Effect of Undetected Barriers on Groundwater Drawdown and Recovery

In large-scale pumping projects, such as mine dewatering, predictions are often made about the rate of groundwater level recovery after pumping has ceased. However, these predictions may be impacted by geological uncertainty—including the presence of undetected impermeable barriers. During pumping, an impermeable barrier may be undetected if it is located beyond the maximum extent of the cone of depression; yet it may still control drawdown during the recovery phase. This has implications for regional-scale modeling and monitoring of groundwater level recovery. In this article, non-dimensional solutions are developed to show the conditions under which a barrier may be undetected during pumping but still significantly impact groundwater level recovery. The magnitude of the impact from an undetected barrier will increase as the ratio of pumping rate to aquifer transmissivity increases. The results are exemplified for a hypothetical aquifer with an unknown barrier 3 km from a pumping well. The difference in drawdown between a model with and without a barrier may be <1 m in the 10 years while pumping is occurring, but up to 50 m after pumping has ceased.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

A Simple Method for Simulating Groundwater Interactions with Fens to Forecast Development Effects

Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.     

Zooplankton community structure in a cyclonic and mode-water eddy in the Sargasso Sea

Mesoscale eddies are important suppliers of nutrients to the surface waters of oligotrophic gyres, but little is known about the biological response, particularly that of higher trophic levels, to these physical perturbations. During the summers of 2004 and 2005, we followed the development of a cyclonic eddy and an anti-cyclonic mode-water eddy in the Sargasso Sea. Zooplankton (>150 μm) were collected across both eddies in 9 discrete depth intervals between 0 and 700 m. Comparison of the abundance of major taxa of mesozooplankton in the upper 150 m at eddy center and outside the eddies (day and night) indicated that the cyclone and mode-water eddy supported similar mesozooplankton communities, with several taxa significantly higher in abundance inside than outside the eddies, when compared with the Bermuda Atlantic Time-series Study site as representative of mean conditions. In both eddies copepod peak abundance occurred in the 50-100 m depth interval, coincident with the chlorophyll a maximum, suggesting elevated food concentration in the eddies may have influenced zooplankton vertical distribution. The two eddies differed in the strength of diel vertical migration of zooplankton, as indicated by the ratio of night:day abundance in the epipelagic zone, which was higher at the center of the mode-water eddy for most taxa. Over the sampling interval of 1-2 months, abundance of the three most common taxa (copepods, chaetognaths, and ostracods) decreased in the cyclone and increased in the mode-water eddy. This further supports previous findings that over the sampling period the cyclone was in a decay phase, while the mode-water eddy was sustaining nutrient fluxes and high phytoplankton concentrations. A more detailed analysis of community structure in the mode-water eddy indicated the 0-700 m integrated abundance of doliolids was significantly higher inside the mode-water eddy than outside. The presence of a mesopelagic (200-700 m) layer of lepadid barnacle cyprids in this eddy highlights the potential of eddies to transport and disperse biota. We conclude that when compared with average ambient conditions (as measured at BATS), mesoscale eddies can influence zooplankton behavior and alter zooplankton community structure which can affect food-web interactions and biogeochemical cycling in the open ocean.     

Tidal decay and stable Roche-lobe overflow of short-period gaseous exoplanets

Many gaseous exoplanets in short-period orbits are on the verge or are in the process of Roche-lobe overflow (RLO). Moreover, orbital stability analysis shows tides can drive many hot Jupiters to spiral inevitably toward their host stars. Thus, the coupled processes of orbital evolution and RLO likely shape the observed distribution of close-in exoplanets and may even be responsible for producing some of the short-period rocky planets. However, the exact outcome for an overflowing planet depends on its internal response to mass loss, and the accompanying orbital evolution can act to enhance or inhibit RLO. In this study, we apply the fully-featured and robust Modules for Experiments in Stellar Astrophysics suite to model RLO of short-period gaseous planets. We show that, although the detailed evolution may depend on several properties of the planetary system, it is largely determined by the core mass of the overflowing gas giant. In particular, we find that the orbital expansion that accompanies RLO often stops and reverses at a specific maximum period that depends on the core mass. We suggest that RLO may often strand the remnant of a gas giant near this orbital period, which provides an observational prediction that can corroborate the hypothesis that short-period gas giants undergo RLO. We conduct a preliminary comparison of this prediction to the observed population of small, short-period planets and find some planets in orbits that may be consistent with this picture. To the extent that we can establish some short-period planets are indeed the remnants of gas giants, that population can elucidate the properties of gas giant cores, the properties of which remain largely unconstrained.