Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling and Using Context”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39Linked to original sources

Documentation and mapping of flooding from the January and March 2018 nor’easters in coastal New England

In January and March 2018, coastal Massachusetts experienced flooding from two separate nor’easters. To put the January and March floods into historical context, the USGS computed statistical stillwater elevations. Stillwater elevations recorded in January 2018 in Boston (9.66 feet relative to the North American Vertical Datum of 1988) have an annual exceedance probability of between 2 and 1 percent (between a 50- and 100-year recurrence interval). Stillwater elevations recorded in March 2018 in Boston (9.17 feet relative to the North American Vertical Datum of 1988) have an annual exceedance probability of between 4 and 2 percent (between a 25- and 50-year recurrence interval). Flood maps show that the area inundated by the January storm is slightly more extensive than that of the March storm, reflecting the respective profiles of the two storms. On the basis of a limited dataset, the attenuation of peak water levels was estimated as a function of the hydraulic distance inland and the starting stillwater elevation computed for the flood within 0.6 foot of what was measured in the field. A simple one-dimensional model was calibrated using flood elevation data collected after the January flood, and the results of the model were validated using flood elevation data collected after the March flood to model the attenuation of the flood elevations as the storms move inland.

Connecticut, Massachusetts, Maine, New Hampshire, ↗

Multiple dimensions define thresholds for population resilience of the eastern oyster, Crassostrea virginica

A species' distribution depends on its tolerance to environmental conditions. These conditions are defined by a minimum, maximum, and optimal ranges of single and combined factors. Forays into environmental conditions outside the minimum or maximum tolerance of a species (i.e., thresholds) are predicted to have large effects on a species' population and may help predict population resilience in the face of changing conditions. Here, we explore ecological thresholds for an important fisheries species and ecosystem engineer, Crassostrea virginica (eastern oyster). In coastal Louisiana, extreme freshwater inputs from rivers and precipitation events impact estuarine salinity, which is a key driver of oyster population dynamics. Using daily salinity and monthly oyster abundance monitoring data across Louisiana estuaries, we explore low salinity exposure threshold levels for oysters. Two statistical approaches were applied, with each model highlighting a different operational definition of a threshold: random forest models identified a threshold as an abrupt change in the oyster abundance- salinity relationship, while Bayesian models identified an increased probability of oyster abundance dropping below a critical threshold, defined here as less than 50% of the 5-year mean. All model results indicate oysters in coastal Louisiana experience low salinity exposure thresholds, defined as the number of consecutive summer days of salinity levels less than 5. However, actual number of days and salinity threshold differed by statistical approach, oyster life stage, and estuary highlighting the multiple dimensions defining ecological thresholds. While thresholds are considered important benchmarks to inform management and assess population or ecosystem vulnerability, our results reveal the need to carefully relate threshold definition to management goals and to acknowledge that thresholds may be highly context dependent.

Louisiana↗

Trends in midwinter counts of bald eagles in the contiguous United States, 1986-2000

We estimated statewide, regional, and national trends in counts of Bald Eagles ( Haliaeetus leucocephalus ) along selected routes in the contiguous United States during midwinter, 1986-2000. Each January, several hundred observers collected data as part of a survey initiated by the National Wildlife Federation in 1979. To analyze these data, we used only those routes surveyed consistently in at least four years and on which at least four eagles were counted in a single year. We included surveys conducted during fog or precipitaion after determining that changes in weather conditions probably did not affect trend estimates. Our final analysis, using a hierarchical mixed model, was based on 101,777 eagle sightings during 5,180 surveys of 563 routes in 42 states. In the model, fixed effects were year, region, and route-length category; the random effect was the route itself. Model-based estimates of Bald Eagle counts throughout the U.S. increased 1.9% yr -1 , but trend estimates varied by region. Estimated trends were statistically significant, and positive, in the northeastern U.S. (6.1% yr -1 ), but were not significant in other regions. The proportion of increasing counts was higher north of 40° N and east of 100° W. Trends in numbers of adults and immatures showed similar geographic patterns, but counts of adults increased at a higher rate. Overall, trends were more similar to those identified by the Christmas Bird Count than the Breeding Bird Survey. In spite of limitations, the survey is a cost-effective way to monitor wintering eagles in the lower 48 states. We discuss estimated trends in the context of increased urbanization, changed weather, and recovery from pesticide pollution.

Bird Populations↗

Movement behavior preceding autumn mortality for white-tailed deer in central New York

A common yet largely untested assumption in the theory of animal movements is that increased rates and a wider range of movements, such as occurs during breeding, make animals more vulnerable to mortality. We examined mortality among 34 white-tailed deer (Odocoileus virginianus) wearing GPS collars during the autumn breeding season of 2006 and 2007 in a heavily hunted, forest-agricultural landscape of central New York state. We evaluated whether individuals having higher rates of movement incurred higher rates of mortality and whether mortality risk was higher when deer were in less familiar areas. We used a Cox proportional hazards model to analyze how mortality risk changes with movement rates measured over 3 time periods: < 1 day, up to 2 weeks prior to death, and 3–4 weeks prior to death. Overall, deer increased their movement rates as autumn progressed, males more so than females. However, deer that died moved at a slower rate relative to surviving deer up to 2 weeks prior to death (ß = -2.22 ± 0.81; 95% confidence interval [CI] = -3.91 to -0.51) and a slower rate on their day of death compared to deer that survived (ß = -1.77 ± 0.73; 95% CI = -3.19 to -0.33). Site familiarity was not significantly related to mortality risk. Deer were equally likely to die within their 50% core use area as elsewhere within their autumn home range. We hypothesize that increased sociality associated with breeding may make animals more vulnerable to harvest mortality. Our findings contradict general assumptions about the influences of movement behavior on mortality risk, suggesting that patterns may be sensitive to the spatiotemporal context of the movement analysis.

New York↗

Implications of community concordance for assessing stream integrity at three nested spatial scales in Minnesota, U.S.A.

1. Fish and invertebrate assemblage data collected from 670 stream sites in Minnesota (U.S.A.) were used to calculate concordance across three nested spatial scales (statewide, ecoregion and catchment). Predictive taxa richness models, calibrated using the same data, were used to evaluate whether concordant communities exhibited similar trends in human‐induced taxa loss across all three scales. Finally, we evaluated the strength of the relationship between selected environmental variables and the composition of both assemblages at all three spatial scales. 2. Significant concordance between fish and invertebrate communities occurred at the statewide scale as well as in six of seven ecoregions and 17 of the 21 major catchments. However, concordance was not consistently indicative of significant relationships between rates of fish and invertebrate taxa loss at those same scales. 3. Fish and invertebrate communities were largely associated with different environmental variables, although the composition of both communities was strongly correlated with stream size across all three scales. 4. Predictive taxa‐loss models for fish assemblages were less sensitive and precise than models for invertebrate assemblages, likely because of the relatively low number of common fish taxa in our data set. Both models, however, distinguished reference from non‐reference sites. 5. The importance of concordance, geographic context and scale are discussed in relation to the design and interpretation of stream integrity indicators. In particular, our findings suggest that community concordance should not be viewed as a substitute for an evaluation of how assemblages respond to environmental stressors.

Minnesota↗

Intraspecific variation in surface water uptake in a perennial desert shrub

Despite broad recognition that water is a major limiting factor in arid ecosystems, we lack an empirical understanding of how this resource is shared and distributed among neighbouring plants. Intraspecific variability can further contribute to this variation via divergent life‐history traits, including root architecture. We investigated these questions in the shrub Artemisia tridentata and hypothesized that the ability to access and utilize surface water varies among subspecies and cytotypes. We used an isotope tracer to quantify below‐ground zone of influence in A. tridentata , and tested whether spatial neighbourhood characteristics can alter plant water uptake. We introduced deuterium‐enriched water to the soil in plant interspaces in a common garden experiment and measured deuterium composition of plant stems. We then applied spatially explicit models to test for differential water uptake by A. tridentata , including intermingled populations of three subspecies and two ploidy levels. The results suggest that lateral root functioning in A. tridentata is associated with intraspecific identity and ploidy level. Subspecies adapted to habitats with deep soils generally had a smaller horizontal reach, and polyploid cytotypes were associated with greater water uptake compared to their diploid variants. We also found that plant crown volume was a weak predictor of water uptake, and that neighbourhood crowding had no discernable effect on water uptake. Intraspecific variation in lateral root functioning can lead to differential patterns of resource acquisition, an essential process in arid ecosystems in the contexts of changing climate and seasonal patterns of precipitation. Altogether, we found that lateral root development and activity are more strongly related to genetic variability within A. tridentata than to plant size. Our study highlights how intraspecific variation in life strategies is linked to mechanisms of resource acquisition.

Functional Ecology↗

Hydroclimatology of the Mississippi River Basin

Model estimated monthly water balance (WB) components (i.e., potential evapotranspiration, actual evapotranspiration, and runoff [ R ]) for 848 United States (U.S.) Geological Survey 8-digit hydrologic units located in the Mississippi River Basin (MRB) are used to examine the temporal and spatial variability of the MRB WB for water years 1901 through 2014. Results indicate the MRB can be divided into nine subregions with similar temporal variability in R . The WB analyses indicated ~79% of total water-year MRB runoff is generated by four of the nine subregions and most of the R in the basin is derived from surplus ( S ) water during the months of December through May. Furthermore, the analyses showed temporal variability in S is largely controlled by the occurrence of negative atmospheric pressure anomalies over the western U.S. and positive atmospheric pressure anomalies over the eastern U.S. coast. This combination of atmospheric pressure anomalies results in an anomalous flow of moist air from the Gulf of Mexico into the MRB. In the context of paleo-climate reconstructions of the Palmer Drought Severity Index, since about 1900 the MRB has experienced wetter conditions than were experienced during the previous 500 years.

Mississippi River Basin↗

Two-dimensional streamflow simulations of the Jordan River, Midvale and West Jordan, Utah

The Jordan River in Midvale and West Jordan, Utah, flows adjacent to two U.S. Environmental Protection Agency Superfund sites: Midvale Slag and Sharon Steel. At both sites, geotechnical caps extend to the east bank of the river. The final remediation tasks for these sites included the replacement of a historic sheet-pile dam and the stabilization of the river banks adjacent to the Superfund sites. To assist with these tasks, two hydraulic modeling codes contained in the U.S. Geological Survey (USGS) Multi-Dimensional Surface-Water Modeling System (MD_SWMS), System for Transport and River Modeling (SToRM) and Flow and Sediment Transport and Morphological Evolution of Channels (FaSTMECH), were used to provide predicted water-surface elevations, velocities, and boundary shear-stress values throughout the study reach of the Jordan River. A SToRM model of a 0.7 mile subreach containing the sheet-pile dam was used to compare water-surface elevations and velocities associated with the sheet-pile dam and a proposed replacement structure. Maps showing water-surface elevation and velocity differences computed from simulations of the historic sheet-pile dam and the proposed replacement structure topographies for streamflows of 500 and 1,000 cubic feet per second (ft 3 /s) were created. These difference maps indicated that the velocities associated with the proposed replacement structure topographies were less than or equal to those associated with the historic sheet-pile dam. Similarly, water-surface elevations associated with the proposed replacement structure topographies were all either greater than or equal to water-surface elevations associated with the sheet-pile dam. A FaSTMECH model was developed for the 2.5-mile study reach to aid engineers in bank stabilization designs. Predicted water-surface elevations, velocities and shear-stress values were mapped on an aerial photograph of the study reach to place these parameters in a spatial context. Profile plots of predicted cross-stream average water-surface elevations and cross-stream maximum and average velocities showed how these parameters change along the study reach for two simulated discharges of 1,040 ft 3 /s and 2,790 ft 3 /s. The profile plots for the simulated streamflow of 1,040 ft 3 /s show that the highest velocities are associated with the constructed sheet-pile replacement structure. Results for the simulated streamflow of 2,790 ft 3 /s indicate that the geometry of the 7800 South Bridge causes more backwater and higher velocities than the constructed sheet-pile replacement structure.

Utah↗

Estimating forest canopy cover dynamics in Valles Caldera National Preserve, New Mexico, using LiDAR and Landsat data

Increasing tree canopy cover has led to increasing wildfire activity in conifer dominated areas of the southwestern United States. Estimating historical changes in the spatial distribution of tree canopy cover can provide further insights into the dynamics of forest and fuel conditions in these landscapes and help prioritize areas for restoration to mitigate wildfire risks and restore biological functioning. In this study, we explored the relationship between LiDAR derived canopy cover data and Landsat reflectance values, and derived a model to estimate percent canopy cover (PCC) on historical Landsat data from 1987 to 2015 for the Valles Caldera National Preserve (VCNP), located in the southwest Jemez Mountains of New Mexico. We developed a regression model between LiDAR generated canopy cover collected in June 2010 and Landsat Thematic Mapper (TM) reflectance values (bands 1–7 except band 6) and vegetation indices collected for the same date. About 5% (17,000) of the total LiDAR points (329,102) were used as training points and a separate, non-overlapping set of 17,000 points as test points to validate the regression model. A simple linear model with the red band (band 3; R 2 = 0.70) was selected as the best model to predict PCC in the rest of the images for 1987–2015. In general, we found a strong consistency between the spatial dynamics of modelled tree canopy cover based on historical Landsat data, wildfire events and forest management practicesthat occurred during the same period. Results showed that about 11% of the study area experienced an increase in PCC for the period of 1987–2015 while 41% of the study area experienced a reduction in PCC during the same time period, mostly in the areas which were affected by stand replacing wildfires in 2011 and 2013. The results indicate an overall increase in medium and high canopy cover classes in specific regions of the study area, which could lead to hazardous wildfires such as those in 2011 and 2013. In the context of ongoing ecological restoration of these montane forests, predicted PCC of contemporary forests could help local managers to identify the areas in the need of immediate restoration efforts by focusing management practices on the areas with closed canopy.

New Mexico↗

Identifying factors that affect mountain lake sensitivity to atmospheric nitrogen deposition across multiple scales

Increased nitrogen (N) deposition rates over the past century have affected both North American and European mountain lake ecosystems. Ecological sensitivity of mountain lakes to N deposition varies, however, because chemical and biological responses are modulated by local watershed and lake properties. We evaluated predictors of mountain lake sensitivity to atmospheric N deposition across North American and European mountain ranges and included as response variables dissolved inorganic N (DIN = N NH 4 + + N NO 3 – ) concentrations and phytoplankton biomass. Predictors of these responses were evaluated at three different spatial scales (hemispheric, regional, subregional) using regression tree, random forest, and generalized additive model (GAM) analysis. Analyses agreed that Northern Hemisphere mountain lake DIN was related to N deposition rates and smaller scale spatial variability (e.g., regional variability between North American and European lakes, and subregional variability between mountain ranges). Analyses suggested that DIN, N deposition, and subregional variability were important for Northern Hemisphere mountain lake phytoplankton biomass. Together, these findings highlight the need for finer-scale, subregional analyses (by mountain range) of lake sensitivity to N deposition. Subregional analyses revealed differences in predictor variables of lake sensitivity. In addition to N deposition rates, lake and watershed features such as land cover, bedrock geology, maximum lake depth (Z max ), and elevation were common modulators of lake DIN. Subregional phytoplankton biomass was consistently positively related with total phosphorus (TP) in Europe, while North American locations showed variable relationships with N or P. This study reveals scale-dependent watershed and lake characteristics modulate mountain lake ecological responses to atmospheric N deposition and provides important context to inform empirically based management strategies.

Water Research↗

Impacts of climate change and renewable energy development on habitat of an endemic squirrel, Xerospermophilus mohavensis , in the Mojave Desert, USA

Predicting changes in species distributions under a changing climate is becoming widespread with the use of species distribution models (SDMs). The resulting predictions of future potential habitat can be cast in light of planned land use changes, such as urban expansion and energy development to identify areas with potential conflict. However, SDMs rarely incorporate an understanding of dispersal capacity, and therefore assume unlimited dispersal in potential range shifts under uncertain climate futures. We use SDMs to predict future distributions of the Mojave ground squirrel, Xerospermophilus mohavensis Merriam, and incorporate partial dispersal models informed by field data on juvenile dispersal to assess projected impact of climate change and energy development on future distributions of X. mohavensis . Our models predict loss of extant habitat, but also concurrent gains of new habitat under two scenarios of future climate change. Under the B1 emissions scenario- a storyline describing a convergent world with emphasis on curbing greenhouse gas emissions- our models predicted losses of up to 64% of extant habitat by 2080, while under the increased greenhouse gas emissions of the A2 scenario, we suggest losses of 56%. New potential habitat may become available to X. mohavensis , thereby offsetting as much as 6330 km 2 (50%) of the current habitat lost. Habitat lost due to planned energy development was marginal compared to habitat lost from changing climates, but disproportionately affected current habitat. Future areas of overlap in potential habitat between the two climate change scenarios are identified and discussed in context of proposed energy development.

California↗

Post-drought groundwater storage recovery in California’s Central Valley

Groundwater depletion is a major threat to agricultural and municipal water supply in California's Central Valley. Recent droughts during 2007–2009 and 2012–2016 exacerbated chronic groundwater depletion. However, it is unclear how much groundwater storage recovered from drought-related overdrafts during post-drought years, and how climatic conditions and water management affected recovery times. We estimated groundwater storage change in the Central Valley for April 2002 through September 2019 using four methods: GRACE satellite data, a water balance approach, a hydrologic simulation model, and monitoring wells. We also evaluated the sensitivity of drought recovery to different climate scenarios (recent climate ± droughts and future climate change scenarios: 20 GCMs and 2 RCPs) using water balance method and statistical sampling of historical climate data. Estimated Central Valley groundwater loss during the two droughts ranged from 19 km 3 (2007–2009) to 28 km 3 (2012–2016) (median of four methods). Median aquifer storage recovery was 34% and 19% of the overdraft during the 2010–2011 and 2017–2019 post-drought years, respectively. Numerical experiments show that recovery times are sensitive to climate forcing, with longer recovery times for a future climate scenario that replicate historical climatology relative to historical forcing with no droughts. Overdraft recovery times decrease by ∼2× with implementation of pumping restrictions (30th to 50th percentiles of historical groundwater depletion) to constrain groundwater depletion relative to no restrictions with a no-drought future climatology. This study highlights the importance of considering water management implications for future drought recoveries within the context of climate change scenarios.

California↗

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York↗

Coastal habitats of the Elwha River, Washington- Biological and physical patterns and processes prior to dam removal

This report includes chapters that summarize the results of multidisciplinary studies to quantify and characterize the current (2011) status and baseline conditions of the lower Elwha River, its estuary, and the adjacent nearshore ecosystems prior to the historic removal of two long-standing dams that have strongly influenced river, estuary, and nearshore conditions. The studies were conducted as part of the U.S. Geological Survey Multi-disciplinary Coastal Habitats in Puget Sound (MD-CHIPS) project. Chapter 1 is the introductory chapter that provides background and a historical context for the Elwha River dam removal and ecosystem restoration project. In chapter 2, the volume and timing of sediment delivery to the estuary and nearshore are discussed, providing an overview of the sediment stored in the two reservoirs and the expected erosion mechanics of the reservoir sediment deposits after removal of the dams. Chapter 3 describes the geological background of the Olympic Peninsula and the geomorphology of the Elwha River and nearshore. Chapter 4 details a series of hydrological data collected by the MD-CHIPS Elwha project. These include groundwater monitoring, surface water-groundwater interactions in the estuary, an estimated surface-water budget to the estuary, and a series of temperature and salinity measurements. Chapter 5 details the work that has been completed in the nearshore, including the measurement of waves, tides, and currents; the development of a numerical hydrodynamic model; and a description of the freshwater plume entering the Strait of Juan de Fuca. Chapter 6 includes a characterization of the nearshore benthic substrate developed using sonar, which formed a habitat template used to design scuba surveys of the benthic biological communities. Chapter 7 describes the ecological studies conducted in the lower river and estuary and includes characterization of juvenile salmon diets and seasonal estuary utilization patterns using otolith analysis to determine habitat specific and hatchery compared with wild patterns in juvenile Chinook salmon, assessment of benthic and terrestrial macroinvertebrate communities, and seasonal patterns of water nutrients. In Chapter 8, the vegetation communities of the eastern estuary are characterized by mapped vegetation cover types and samples collected for vegetation composition and diversity. Chapter 9 summarizes the existing conditions of the study area as detailed in this report and describes some of the possible outcomes of river restoration on the coastal ecosystems of the Elwha River. Together, these different scientific perspectives form a basis for understanding the Elwha River ecosystem, an environment that has and will undergo substantial change. A century of change began with the start of dam construction in 1910; additional major change will result from dam removal scheduled to begin in September 2011. This report provides a scientific snapshot of the lower Elwha River, its estuary, and adjacent nearshore ecosystems prior to dam removal that can be used to evaluate the responses and dynamics of various system components following dam removal.

Scientific Investigations Report↗

Predicting fine-scale distributions of peripheral aquatic species in headwater streams

Headwater species and peripheral populations that occupy habitat at the edge of a species range may hold an increased conservation value to managers due to their potential to maximize intraspecies diversity and species' adaptive capabilities in the context of rapid environmental change. The southern Appalachian Mountains are the southern extent of the geographic range of native Salvelinus fontinalis and naturalized Oncorhynchus mykiss and Salmo trutta in eastern North America. We predicted distributions of these peripheral, headwater wild trout populations at a fine scale to serve as a planning and management tool for resource managers to maximize resistance and resilience of these populations in the face of anthropogenic stressors. We developed correlative logistic regression models to predict occurrence of brook trout, rainbow trout, and brown trout for every interconfluence stream reach in the study area. A stream network was generated to capture a more consistent representation of headwater streams. Each of the final models had four significant metrics in common: stream order, fragmentation, precipitation, and land cover. Strahler stream order was found to be the most influential variable in two of the three final models and the second most influential variable in the other model. Greater than 70% presence accuracy was achieved for all three models. The underrepresentation of headwater streams in commonly used hydrography datasets is an important consideration that warrants close examination when forecasting headwater species distributions and range estimates. Additionally, it appears that a relative watershed position metric (e.g., stream order) is an important surrogate variable (even when elevation is included) for biotic interactions across the landscape in areas where headwater species distributions are influenced by topographical gradients.

Ecology and Evolution↗

Understanding the spatiotemporal distribution of snow refugia in the rain-snow transition zone of north-central Idaho

Knowledge of snow cover distribution and disappearance dates over a wide range of scales is imperative for understanding hydrological dynamics and for habitat management of wildlife species that rely on snow cover. Identification of snow refugia, or places with relatively late snow disappearance dates (SDDs) compared to surrounding areas, is especially important as climate change alters snow cover timing and duration. The purpose of this study was to increase understanding of snow refugia in complex terrain spanning the rain-snow transition zone at fine spatial and temporal scales. To accomplish this objective, we used remote cameras to provide relatively high temporal and spatial resolution measurements on snowpack conditions. We built linear models to relate SDDs at the monitoring sites to topoclimatic and canopy cover metrics. One model to quantify SDDs included elevation, aspect, and an interaction between canopy cover and cold-air pooling potential. High-elevation, north-facing sites in cold-air pools (CAPs) had the latest SDDs, but isolated lower-elevation points also exhibited relatively late potential SDDs. Importantly, canopy cover had a much stronger effect on SDDs in CAPs than in non-CAPs, indicating that best practices in forest management for snow refugia could vary across microtopography. A second model that included in situ hydroclimate observations (December – February (DJF) temperature and March 1 snow depth) indicated that March 1 snow depth had little impact on SDD at the coldest winter temperatures, and that DJF temperatures had a stronger effect on SDD at lower snow depths, implying that the relative importance of snowfall and temperature could vary across hydroclimatic contexts in their impact on snow refugia. This new understanding of factors influencing snow refugia can guide forest management actions to increase snow retention and inform management of snow-dependent wildlife species in complex terrain.

Environmental Research Letters↗

Hydrogeology and hydrologic conditions of the Ozark Plateaus aquifer system

The hydrogeology and hydrologic characteristics of the Ozark Plateaus aquifer system were characterized as part of ongoing U.S. Geological Survey efforts to assess groundwater availability across the Nation. The need for such a study in the Ozark Plateaus physiographic province (Ozark Plateaus) is highlighted by increasing demand on groundwater resources by the 5.3 million people of the Ozark Plateaus, water-level declines in some areas, and potential impacts of climate change on groundwater availability. The subject study integrates knowledge gained through local investigation within a regional perspective to develop a regional conceptual model of groundwater flow in the Ozark Plateaus aquifer system (Ozark system), a key phase of groundwater availability assessment. The Ozark system extends across much of southern Missouri and northwestern and north-central Arkansas and smaller areas of southeastern Kansas and northeastern Oklahoma. The region is one of the major karst landscapes in the United States, and karst aquifers are predominant in the Ozark system. Groundwater flow is ultimately controlled by aquifer and confining unit lithologies and stratigraphic relations, geologic structure, karst development, and the character of surficial lithologies and regolith mantle. The regolith mantle is a defining element of Ozark Plateaus karst, affecting recharge, karst development, and vulnerability to surface-derived contaminants. Karst development is more advanced—as evidenced by larger springs, hydraulic characteristics, and higher well yields—in the Salem Plateau and in the northern part of the Springfield Plateau (generally north of the Arkansas-Missouri border) as compared with the southern part of the Springfield Plateau in Arkansas, largely due to thinner, less extensive regolith and purer carbonate lithology. Precipitation is the ultimate source of all water to the Ozark system, and the hydrologic budget for the Ozark system includes inputs from recharge, losing-stream sections, and groundwater inflows and losses of water to gaining-stream sections, groundwater withdrawals, and surface-water and groundwater outflows to neighboring systems. Groundwater recharge, estimated by a soil-water-balance model, represents about 24 percent, or 11 inches, of 43.9 inches annual precipitation. Recharge is spatially variable, being greater in the northern Springfield Plateau and Salem Plateau than in the southern Springfield Plateau (generally south of the Arkansas border) because of differences in regolith mantle extent and thickness and carbonate lithology and hydraulic properties. Increased precipitation and decreased agricultural land use during the period 1951 through 2011 increased recharge by approximately 5 percent. Although all Ozark streams have losing, neutral, and gaining sections, they are dominantly gaining and are a net sink for groundwater with nearly 90 percent of groundwater recharge returned to springs and streams. Groundwater pumping is a small but important loss of water in the Ozark system hydrologic budget; water-level declines and local cones of depression have been observed around pumping centers and strong concerns exist over potential effects on stream and spring flow. Data indicate that societal needs for freshwater resources in the Ozark Plateaus will continue to increase and will do so in the context of changing climate and hydrology. Groundwater will continue to be an important part of supporting these societal needs and also local ecosystems. The unique character and hydrogeologic variability across the Ozark system will control how the system responds to future stress. Groundwater of the Ozark system in the northern study area is more dynamic, has greater storage and larger flux, and has greater potential for further development than in the part of the study area south of the Arkansas-Missouri border. Further south in Arkansas, a line exists, roughly defined as 5 miles south of the Springfield Plateau-Boston Mountains boundary, beyond which further extensive municipal or commercial development appears unlikely under current economic and resource-need conditions. A small part of the Ozark system groundwater budget is currently drafted for use, leaving an apparently large component available for further development and use—particularly in the northern Springfield Plateau and Salem Plateau; however, the effects of increased pumping on groundwater’s role in maintaining ecosystems and ecosystem services are not quantitatively well understood, and the close relation between groundwater and surface water highlights the importance of further quantitative assessment.

Arkansas, Kansas, Missouri, Oklahoma↗

The changing shapes of active volcanoes: History, evolution, and future challenges for volcano geodesy

At the very heart of volcanology lies the search for the 'plumbing systems' that form the inner workings of Earth’s active volcanoes. By their very nature, however, the magmatic reservoirs and conduits that underlie these active volcanic systems are elusive; mostly they are observable only through circumstantial evidence, using indirect, and often ambiguous, surficial measurements. Of course, we can infer much about these systems from geologic investigation of materials brought to the surface by eruptions and of the exposed roots of ancient volcanoes. But how can we study the magmatic processes that are occurring beneath Earth’s active volcanoes? What are the geometry, scale, physical, and chemical characteristics of magma reservoirs? Can we infer the dynamics of magma transport? Can we use this information to better forecast the future behavior of volcanoes? These questions comprise some of the most fundamental, recurring themes of modern research in volcanology. The field of volcano geodesy is uniquely situated to provide critical observational constraints on these problems. For the past decade, armed with a new array of technological innovations, equipped with powerful computers, and prepared with new analytical tools, volcano geodesists have been poised to make significant advances in our fundamental understanding of the behavior of active volcanic systems. The purpose of this volume is to highlight some of these recent advances, particularly in the collection and interpretation of geodetic data from actively deforming volcanoes. The 18 papers that follow report on new geodetic data that offer valuable insights into eruptive activity and magma transport; they present new models and modeling strategies that have the potential to greatly increase understanding of magmatic, hydrothermal, and volcano-tectonic processes; and they describe innovative techniques for collecting geodetic measurements from remote, poorly accessible, or hazardous volcanoes. To provide a proper context for these studies, we offer a short review of the evolution of volcano geodesy, as well as a case study that highlights recent advances in the field by comparing the geodetic response to recent eruptive episodes at Mount St. Helens. Finally, we point out a few areas that continue to challenge the volcano geodesy community, some of which are addressed by the papers that follow and which undoubtedly will be the focus of future research for years to come.

Journal of Volcanology and Geothermal Research↗