Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Spatial Statistics”

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 1,171 records · Page 65Linked to original sources

Geospatial association between adverse birth outcomes and arsenic in groundwater in New Hampshire, USA

There is increasing evidence of the role of arsenic in the etiology of adverse human reproductive outcomes. Because drinking water can be a major source of arsenic to pregnant women, the effect of arsenic exposure through drinking water on human birth may be revealed by a geospatial association between arsenic concentration in groundwater and birth problems, particularly in a region where private wells substantially account for water supply, like New Hampshire, USA. We calculated town-level rates of preterm birth and term low birth weight (term LBW) for New Hampshire, by using data for 1997&ndash;2009 stratified by maternal age. We smoothed the rates by using a locally weighted averaging method to increase the statistical stability. The town-level groundwater arsenic probability values are from three GIS data layers generated by the US Geological Survey: probability of local groundwater arsenic concentration >1 &micro;g/L, probability >5 &micro;g/L, and probability >10 &micro;g/L. We calculated Pearson&rsquo;s correlation coefficients ( r ) between the reproductive outcomes (preterm birth and term LBW) and the arsenic probability values, at both state and county levels. For preterm birth, younger mothers (maternal age <20) have a statewide r = 0.70 between the rates smoothed with a threshold = 2,000 births and the town mean arsenic level based on the data of probability >10 &micro;g/L; for older mothers, r = 0.19 when the smoothing threshold = 3,500; a majority of county level r values are positive based on the arsenic data of probability >10 &micro;g/L. For term LBW, younger mothers (maternal age <25) have a statewide r = 0.44 between the rates smoothed with a threshold = 3,500 and town minimum arsenic concentration based on the data of probability >1 &micro;g/L; for older mothers, r = 0.14 when the rates are smoothed with a threshold = 1,000 births and also adjusted by town median household income in 1999, and the arsenic values are the town minimum based on probability >10 &micro;g/L. At the county level for younger mothers, positive r values prevail, but for older mothers, it is a mix. For both birth problems, the several most populous counties&mdash;with 60&ndash;80% of the state&rsquo;s population and clustering at the southwest corner of the state&mdash;are largely consistent in having a positive r across different smoothing thresholds. We found evident spatial associations between the two adverse human reproductive outcomes and groundwater arsenic in New Hampshire, USA. However, the degree of associations and their sensitivity to different representations of arsenic level are variable. Generally, preterm birth has a stronger spatial association with groundwater arsenic than term LBW, suggesting an inconsistency in the impact of arsenic on the two reproductive outcomes. For both outcomes, younger maternal age has stronger spatial associations with groundwater arsenic.

New Hampshire↗

Candidate-penetrative-fracture mapping of the Grand Canyon area, Arizona, from spatial correlation of deep geophysical features and surficial lineaments

Some aquifers of the southwestern Colorado Plateaus Province are deeply buried and overlain by several impermeable shale layers, and so recharge to the aquifer probably is mainly by seepage down penetrative-fracture systems. The purpose of this 2-year study, sponsored by the U.S. National Park Service, was to map candidate deep penetrative fractures over a 120,000-km2 area, using gravity and aeromagnetic-anomaly data together with surficial-fracture data. The study area was on the Colorado Plateau south of the Grand Canyon and west of Black Mesa; mapping was carried out at a scale of 1:250,000. The resulting database constitutes a spatially registered estimate of deep-fracture locations. Candidate penetrative fractures were located by spatial correlation of horizontal- gradient and analytic-signal maximums of gravity and magnetic anomalies with major surficial lineaments obtained from geologic, topographic, side-looking-airborne-radar, and satellite imagery. The maps define a subset of candidate penetrative fractures because of limitations in the data coverage and the analytical technique. In particular, the data and analytical technique used cannot predict whether the fractures are open or closed. Correlations were carried out by using image-processing software, such that every pixel on the resulting images was coded to uniquely identify which datasets are correlated. The technique correctly identified known and many new deep fracture systems. The resulting penetrative-fracture-distribution maps constitute an objectively obtained, repeatable dataset and a benchmark from which additional studies can begin. The maps also define in detail the tectonic fabrics of the southwestern Colorado Plateaus Province. Overlaying the correlated lineaments on the normalized-density-of-vegetation-index image reveals that many of these lineaments correlate with the boundaries of vegetation zones in drainages and canyons and so may be controlling near-surface water availability in some places. Many derivative products can be produced from the database, such as fracture-density-estimate maps, and maps with the number of correlations color-coded to estimate the possible quality of correlation. The database contained in this report is designed to be used in a geographic information system and image-processing systems, and most data layers are in georeferenced tagged image format (Geotiff) or ARC grids. The report includes 163 map plates and various metadata, supporting, and statistical diagram files.

Data Series↗

Status and distribution of mangrove forests of the world using earth observation satellite data

Aim Our scientific understanding of the extent and distribution of mangrove forests of the world is inadequate. The available global mangrove databases, compiled using disparate geospatial data sources and national statistics, need to be improved. Here, we mapped the status and distributions of global mangroves using recently available Global Land Survey (GLS) data and the Landsat archive. Methods We interpreted approximately 1000 Landsat scenes using hybrid supervised and unsupervised digital image classification techniques. Each image was normalized for variation in solar angle and earth–sun distance by converting the digital number values to the top-of-the-atmosphere reflectance. Ground truth data and existing maps and databases were used to select training samples and also for iterative labelling. Results were validated using existing GIS data and the published literature to map ‘true mangroves’. Results The total area of mangroves in the year 2000 was 137,760 km2 in 118 countries and territories in the tropical and subtropical regions of the world. Approximately 75% of world's mangroves are found in just 15 countries, and only 6.9% are protected under the existing protected areas network (IUCN I-IV). Our study confirms earlier findings that the biogeographic distribution of mangroves is generally confined to the tropical and subtropical regions and the largest percentage of mangroves is found between 5° N and 5° S latitude. Main conclusions We report that the remaining area of mangrove forest in the world is less than previously thought. Our estimate is 12.3% smaller than the most recent estimate by the Food and Agriculture Organization (FAO) of the United Nations. We present the most comprehensive, globally consistent and highest resolution (30 m) global mangrove database ever created. We developed and used better mapping techniques and data sources and mapped mangroves with better spatial and thematic details than previous studies.

Global Ecology and Biogeography↗

Design of a shallow ground-water network to monitor agricultural chemicals, Lake Wales Ridge, central Florida

Extensive agricultural land use and dynamic hydraulic connections between the land surface and ground-water resources render many of Florida's aquifers vulnerable to chemical contamination. In these areas, there is a need to monitor shallow ground water for agricultural chemicals to evaluate potential migration of the chemicals to the subsurface and to deeper aquifers, and to assess the effects of agricultural practices on ground-water quality. Historically, efforts to monitor water-table aquifers have been minimal compared with monitoring deeper aquifers used for municipal drinking-water supply. The Florida Department of Agriculture and Consumer Services has proposed the establishment of long-term, water-quality networks to monitor shallow ground water in agricultural areas of Florida. These networks will bridge the current information gap between local-scale short-duration field studies, which are required for pesticide licensing, and statewide monitoring of deeper ground-water resources. Laboratory (experimental) and short-term field evaluations of the transport of agricultural chemicals to the subsurface contain some degree of uncertainty due to widespread variations in environmental conditions. The proposed regional-scale shallow ground-water networks will serve as "early warning" networks to avoid contamination of ground-water resources. The objectives of these shallow ground-water networks are to provide information for early detection of pesticides and nitrate in the subsurface, and for evaluating temporal trends in concentrations in relation to changes in land use and agricultural practices. Information on trends is important to the U.S. Environmental Protection Agency, State agencies, and agricultural managers, and is an integral part of Florida's Pesticide Management Plan, which evaluates the effectiveness of implementing "best management practices" to minimize environmental impacts. Lake Wales Ridge in Polk and Highlands Counties has been identified by the Florida Department of Agriculture and Consumer Services as a pilot study area for implementing a regional shallow ground-water network. Information gained from designing the Lake Wales Ridge Network will be used to develop plans as similar networks are considered for other areas of Florida. Lake Wales Ridge is heavily utilized for citrus production, but the well-drained, clean sands underlying the region, coupled with a hydraulic connection with deeper karst formations, render the area extremely vulnerable to the transport and migration of agricultural chemicals within the subsurface. The network design for Lake Wales Ridge was developed on the basis of the objectives of monitoring, factors affecting aquifer vulnerability, probability (statistical) sampling theory, and specifications regarding well construction, sampling, and laboratory analysis. The network will provide estimates of regional pesticide and nitrate concentrations in targeted ground water, monitor for trends over time, and will provide baseline, regional-scale information. The areas targeted for sampling were citrus groves located on soils classified as vulnerable to leaching of agricultural chemicals, and ground water in close proximity to the water table. The duration of monitoring is planned to be long term. Chemicals targeted for monitoring include selected nutrients, trace elements, and pesticides used for citrus management. Standardization of field sampling and laboratory methods, and construction of a computerized data-base repository will be important for long-term utility of the data. The proposed network design consists of existing and newly drilled wells. Existing wells were evaluated as candidates for the network based on location and well-construction specifications for water-quality monitoring. Using probability selection techniques, new well locations were selected in areas that lacked adequate existing wells. Implementation of the network will occur in phases. Quarterly sampling of 13 Phase I wells commenced in April 1999. An additional 11 Phase II wells were added to the network in April 2000. Proposed locations have been identified for 8 Phase III wells. Future topics for study include a regional water-quality analysis, a network design evaluation, and consideration of local-scale studies needed to support regional monitoring efforts. An analysis of historical water-quality data for Lake Wales Ridge, including Phase I wells, is needed for the purpose of examining trends, providing a preliminary description of spatial and short-term variability in water quality, and evaluating the proposed network design for cost effectively meeting program objectives. Local-scale, process-oriented studies are needed to provide important information for evaluating regional water-quality data. Coordination between individuals conducting such local-scale studies and those interpreting regional sampling results will be considered as a component of the network design, and will require coordination between multiple agencies.

Florida↗

A within-season approach for detecting early crop stage of corn and soybean using high temporal and spatial resolution imagery

Crop emergence is a critical stage for crop development and crop growth modeling. Mapping crop emergence using remote sensing data is challenging. Previous remote sensing phenology algorithms showed that crop stages could be detected around the V3-V4 (3 to 4 established leaves) vegetative stage. Traditional approaches have a strong assumption regarding the temporal evolution of plant growth and normally require a complete growth period of observations to define seasonal changes. Most approaches were not designed for the within-season mapping in the early growing season. In the current paper, we developed a new within-season emergence (WISE) approach to mapping crop green-up date using satellite observations during early growth stages. The approach was first optimized using high spatiotemporal resolution (10 m, 2 day revisit) imagery from the Vegetation and Environment monitoring New MicroSatellite (VENµS) research mission, and assessed using ground observations of early crop growth stages (emergence VE and one leaf V1 stages for corn, and emergence VE and unifoliolate VC stages for soybeans) collected over the Beltsville Agricultural Research Center (BARC) experimental fields in Beltsville, MD during the 2019 growing season. Results show that early crop growth stages can be reliably detected at sub-field scale about two weeks after crop emergence. The remote sensing green-up dates were about 4-5 days after crop emergence on average. Coefficients of determination (R2) between green-up dates and the mid-point dates of the early growth stages were above 0.90. The mean absolute differences, standard deviations, and root mean square errors comparing to the early growth stage mid-point dates were within six days. The maximum differences were within ±10 days across all fields. The WISE approach was assessed using operational Sentinel-2 data (10 m, 5 day revisit) in BARC. The detected green-up dates from Sentinel-2 were found close to VENµS results. Some fields were not detected due to the lack of observations during emergence dates. For independent evaluation, the WISE approach was applied over an agricultural watershed on the Maryland Eastern Shore using both VENµS and the harmonized Landsat and Sentinel-2 (HLS) data (30 m, 3-4 day revisit). The green-up dates were compared with crop progress reports of crop emergence dates from the National Agricultural Statistics Service (NASS) at the state-level. The WISE -detected green-up dates at the regional scale are within VE stage ranges but slightly earlier than NASS crop progress reports at the state-level. The WISE approach uses remote sensing observations during the early crop growth stages and has potential for operational application within the season using Sentinel-2 and HLS data.

Maryland↗

Ecological effects of pinyon-juniper removal in the Western United States—A synthesis of scientific research, January 2014–March 2021

Executive Summary Increasing density of pinyon ( Pinus spp. ) and juniper ( Juniperus spp. ) woodlands (hereinafter “pinyon-juniper”), as well as expansion of these woodlands into adjacent shrublands and grasslands, has altered ecosystem function and wildlife habitat across large areas of the interior western United States. Although there are many natural and human-caused drivers of woodland infilling and expansion, restoration of sagebrush ( Artemisia spp. ) habitat through removal of pinyon-juniper is considered an urgent management objective in many locations, particularly in support of sagebrush-dependent wildlife species of conservation concern. In December 2020, the Bureau of Land Management (BLM) established the Pinyon-Juniper Management Categorical Exclusion (PJCX) to expedite the regulatory process for pinyon-juniper removal projects on public lands, largely intended to benefit mule deer ( Odocoileus hemionus ) and greater sage-grouse ( Centrocercus urophasianus ) habitats. During final preparation of this report, the BLM discontinued use of the PJCX (as of November, 2022), but the pinyon-juniper tree removal techniques assessed in this report are commonly used and understanding their effects remains relevant to land use planning. To address areas of uncertainty relative to potential ecological effects of the PJCX, we conducted a review of the peer-reviewed science literature to better understand the likely responses of vegetation, environmental (for example, soils), and wildlife variables to specific tree removal techniques permitted by the PJCX. In brief, the PJCX permitted removal of trees by either manual cutting, mechanical cutting, or mastication; allowed certain methods to redistribute or remove resulting tree biomass after treatment; and prohibited broadcast burning, roadbuilding, removal of old-growth, and seeding of non-native species. Specifically, we conducted our review to address the following questions: How will PJCX removal techniques affect plant communities, soils, and abiotic resources? How do these pinyon-juniper removal techniques affect wildlife communities, including both woodland- and sagebrush-dependent species? What are the potential ecological implications of different pinyon-juniper removal treatment types and implementation strategies (for example, treatment sizes) over time? What are the most important gaps in our scientific understanding of how treatments might affect targeted ecosystems over space and time (for example, potential effects of climate change)? To answer these questions, we considered studies related to pinyon-juniper ecosystems, focusing on research that occurred over a large portion of the interior western United States that is the primary focus of the PJCX. We also focused on papers published from 2014 onward, to avoid excessive overlap with other recent reviews on pinyon-juniper management effects. Using strict criteria, including only considering research that tested responses for statistical significance, we identified 48 papers that primarily examined treatment effects on vegetation and other environmental variables (1,709 responses), and 11 papers that addressed effects on wildlife (132 responses). Responses to the PJCX-permitted treatments were summarized as either positive (that is, a significant increase), negative (that is, a significant decrease), or non-significant (that is, no significant difference). Responses were assigned to categories (for example, Native Annual Grass/Forb Abundance) and hierarchical treatment levels. We found that there were large proportions of non-significant responses among all categories combined, with roughly half or more of all responses non-significant (48 percent for wildlife, 60 percent for vegetation-environmental), comparable to other recent systematic reviews of pinyon-juniper treatment effects. However, we also found that when there were significant responses, some important trends potentially emerged. Important undesirable outcomes included far more positive than negative responses of exotic grass and forb abundance among nearly all treatment types. Cutting treatments were also more likely to decrease biocrust cover and microbial activity. Potentially beneficial outcomes included mostly positive responses among sagebrush obligate species, including more positive than negative responses for mule deer and sage-grouse. Some treatment types (for example, mastication) also resulted in more positive than negative responses for native grasses and forbs (although, non-significant responses were the majority). We also highlighted many limitations of this review, including how responses often come from few studies, and how some response-treatment category combinations lack adequate response data. Moreover, the existing research is often insufficient to address many key questions about treatment effects, largely owing to short time-scales and limited spatial extents of observations, which do not match the size of treatments being implemented by land managers, nor capture long-term, post-treatment ecological dynamics. We also identify a lack of research that addresses key interactions that could undermine restoration objectives, including potential effects of climate change and grazing on post-treatment environments. Thus, we emphasize the importance of integrating these factors into future pinyon-juniper treatment research, and we stress the need for use of monitoring programs and research studies that partake in data collection and analysis over long durations and broad spatial scales.

Open-File Report↗

Hydrology of the coastal springs ground-water basin and adjacent parts of Pasco, Hernando, and Citrus Counties, Florida

The coastal springs in Pasco, Hernando, and Citrus Counties, Florida consist of three first-order magnitude springs and numerous smaller springs, which are points of substantial ground-water discharge from the Upper Floridan aquifer. Spring flow is proportional to the water-level altitude in the aquifer and is affected primarily by the magnitude and timing of rainfall. Ground-water levels in 206 Upper Floridan aquifer wells, and surface-water stage, flow, and specific conductance of water from springs at 10 gaging stations were measured to define the hydrologic variability (temporally and spatially) in the Coastal Springs Ground-Water Basin and adjacent parts of Pasco, Hernando, and Citrus Counties. Rainfall at 46 stations and ground-water withdrawals for three counties, were used to calculate water budgets, to evaluate long-term changes in hydrologic conditions, and to evaluate relations among the hydrologic components. Predictive equations to estimate daily spring flow were developed for eight gaging stations using regression techniques. Regression techniques included ordinary least squares and multiple linear regression techniques. The predictive equations indicate that ground-water levels in the Upper Floridan aquifer are directly related to spring flow. At tidally affected gaging stations, spring flow is inversely related to spring-pool altitude. The springs have similar seasonal flow patterns throughout the area. Water-budget analysis provided insight into the relative importance of the hydrologic components expected to influence spring flow. Four water budgets were constructed for small ground-water basins that form the Coastal Springs Ground-Water Basin. Rainfall averaged 55 inches per year and was the only source of inflow to the Basin. The pathways for outflow were evapotranspiration (34 inches per year), runoff by spring flow (8 inches per year), ground-water outflow from upward leakage (11 inches per year), and ground-water withdrawal (2 inches per year). Recharge (rainfall minus evapotranspiration) to the Upper Floridan aquifer consists of vertical leakage through the surficial deposits. Discharge is primarily through springs and diffuse upward leakage that maintains the extensive swamps along the Gulf of Mexico. The ground-water basins had slightly different partitioning of hydrologic components, reflecting variation among the regions. Trends in hydrologic data were identified using nonparametric statistical techniques to infer long-term changes in hydrologic conditions, and yielded mixed results. No trend in rainfall was detected during the past century. No trend in spring flow was detected in 1931-98. Although monotonic trends were not detected, rainfall patterns are naturally variable from month to month and year to year; this variability is reflected in ground-water levels and spring flows. A decreasing trend in ground-water levels was detected in the Weeki Wachee well (1966-98), but the trend was statistically weak. At current ground-water withdrawal rates, there is no discernible affect on ground-water levels and spring flows. Sporadic data records, lack of continuous data, and inconsistent periods of record among the hydrologic components impeded analysis of long-term changes to the hydrologic system and interrelations among components. The ongoing collection of hydrologic data from index sites could provide much needed information to assess the hydrologic factors affecting the quantity and quality of spring flow in the Coastal Springs Ground-Water Basin.

Water-Resources Investigations Report↗

Model-based surveillance system design under practical constraints with application to white-nose syndrome

Infectious diseases are powerful ecological forces structuring ecosystems, causing devastating economic impacts and disrupting society. Successful prevention and control of pathogens requires knowledge of the current scope and severity of disease, as well as the ability to forecast future disease dynamics. Assessment of the current situation as well as prediction of the future conditions, rely on spatially referenced information regarding the presence or absence of a pathogen, and the prevalence of the pathogen in the population. In particular, knowledge about the location of the disease front is foundational for deploying disease countermeasures to prevent further disease spread and focusing control efforts to reduce disease intensity in affected areas. In this paper, we develop a model-based approach to designing sampling strategies for wildlife disease surveillance at the disease front. Specifically, we use a mechanistic spatio-temporal model based on an underlying partial differential equation to track the disease dynamics and predict the disease prevalence in the future. We also devise an optimal surveillance system design at the disease front that takes into account the practical constraints of sampling. We evaluate the effectiveness of our proposed design via a simulation study and demonstrate the application of the proposed approach by designing a surveillance strategy for the pathogen that causes white-nose syndrome.

Environmental and Ecological Statistics↗

Evaluation of radon occurrence in groundwater from 16 geologic units in Pennsylvania, 1986–2015, with application to potential radon exposure from groundwater and indoor air

Results from 1,041 groundwater samples collected during 1986‒2015 from 16 geologic units in Pennsylvania, associated with 25 or more groundwater samples with concentrations of radon-222, were evaluated in an effort to identify variations in radon-222 activities or concentrations and to classify potential radon-222 exposure from groundwater and indoor air. Radon-222 is hereafter referred to as “radon.” Radon concentrations in groundwater greater than or equal to the proposed U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for public-water supply systems of 300 picocuries per liter (pCi/L) were present in about 87 percent of the water samples, whereas concentrations greater than or equal to the proposed alternative MCL (AMCL) for public water-supply systems of 4,000 pCi/L were present in 14 percent. The highest radon concentrations were measured in groundwater from the schists, gneisses, and quartzites of the Piedmont Physiographic Province. In this study, conducted by the U.S. Geological Survey in cooperation with the Pennsylvania Department of Health and the Pennsylvania Department of Environmental Protection, groundwater samples were aggregated among 16 geologic units in Pennsylvania to identify units with high median radon concentrations in groundwater. Graphical plots and statistical tests were used to determine variations in radon concentrations in groundwater and indoor air. Median radon concentrations in groundwater samples and median radon concentrations in indoor air samples within the 16 geologic units were classified according to proposed and recommended regulatory limits to explore potential radon exposure from groundwater and indoor air. All of the geologic units, except for the Allegheny (Pa) and Glenshaw (Pcg) Formations in the Appalachian Plateaus Physiographic Province, had median radon concentrations greater than the proposed EPA MCL of 300 pCi/L, and the Peters Creek Schist (Xpc), which is in the Piedmont Physiographic Province, had a median radon concentration greater than the EPA proposed AMCL of 4,000 pCi/L. Median concentrations of radon in groundwater and indoor air were determined to differ significantly among the geologic units (Kruskal-Wallis test, significance probability, p<0.001), and Tukey’s test indicated that radon concentrations in groundwater and indoor air in the Peters Creek Schist (Xpc) were significantly higher than those in the other units. Also, the Peters Creek Schist (Xpc) was determined to be the area with highest potential of radon exposure from groundwater and indoor air and one of two units with the highest percentage of population assumed to be using domestic self-supplied water (81 percent), which puts the population at greater potential of exposure to radon from groundwater. Potential radon exposure determined from classification of geologic units by median radon concentrations in groundwater and indoor air according to proposed and recommended regulatory limits is useful for drawing general conclusions about the presence, variation, and potential radon exposure in specific geologic units, but the associated data and maps have limitations. The aggregated indoor air radon data have spatial accuracy limitations owing to imprecision of geocoded test locations. In addition, the associated data describing geologic units and the public water supplier’s service areas have spatial and interpretation accuracy limitations. As a result, data and maps associated with this report are not recommended for use in predicting individual concentrations at specific sites nor for use as a decision-making tool for property owners to decide whether to test for radon concentrations at specific locations. Instead, the data and maps are meant to promote awareness regarding potential radon exposure in Pennsylvania and to point out data gaps that exist throughout the State.

Pennsylvania↗

The United States National Climate Assessment - Alaska Technical Regional Report

The Alaskan landscape is changing, both in terms of effects of human activities as a consequence of increased population, social and economic development and their effects on the local and broad landscape; and those effects that accompany naturally occurring hazards such as volcanic eruptions, earthquakes, and tsunamis. Some of the most prevalent changes, however, are those resulting from a changing climate, with both near term and potential upcoming effects expected to continue into the future. Alaska's average annual statewide temperatures have increased by nearly 4&deg;F from 1949 to 2005, with significant spatial variability due to the large latitudinal and longitudinal expanse of the State. Increases in mean annual temperature have been greatest in the interior region, and smallest in the State's southwest coastal regions. In general, however, trends point toward increases in both minimum temperatures, and in fewer extreme cold days. Trends in precipitation are somewhat similar to those in temperature, but with more variability. On the whole, Alaska saw a 10-percent increase in precipitation from 1949 to 2005, with the greatest increases recorded in winter. The National Climate Assessment has designated two well-established scenarios developed by the Intergovernmental Panel on Climate Change (Nakicenovic and others, 2001) as a minimum set that technical and author teams considered as context in preparing portions of this assessment. These two scenarios are referred to as the Special Report on Emissions Scenarios A2 and B1 scenarios, which assume either a continuation of recent trends in fossil fuel use (A2) or a vigorous global effort to reduce fossil fuel use (B1). Temperature increases from 4 to 22&deg;F are predicted (to 2070-2099) depending on which emissions scenario (A2 or B1) is used with the least warming in southeast Alaska and the greatest in the northwest. Concomitant with temperature changes, by the end of the 21st century the growing season is expected to lengthen by 15-25 days in some areas of Alaska, with much of that corresponding with earlier spring snow melt. Future projections of precipitation (30-80 years) over Alaska show an increase across the State, with the largest changes in the northwest and smallest in the southeast. Because of increasing temperatures and growing season length, however, increased precipitation may not correspond with increased water availability, due to temperature related increased evapotranspiration. The extent of snow cover in the Northern Hemisphere has decreased by about 10 percent since the late 1960s, with stronger trends noted since the late 1980s. Alaska has experienced similar trends, with a strong decrease in snow cover extent occurring in May. When averaged across the State, the disappearance of snow in the spring has occurred from 4 to 6 days earlier per decade, and snow return in fall has occurred approximately 2 days later per decade. This change appears to be driven by climate warming rather than a decrease in winter precipitation, with average winter temperatures also increasing by about 2.5&deg;F. The extent of sea ice has been declining, as has been widely published in both national and scientific media outlets, and is projected to continue to decline during this century. The observed decline in annual sea ice minimum extent (September) has occurred more rapidly than was predicted by climate models and has been accompanied by decreases in ice thickness and in the presence of multi-year ice. This decrease was first documented by satellite imagery in the late 1970s for the Bering and Chukchi Seas, and is projected to continue, with the potential for the disappearance of summer sea ice by mid- to late century. A new phenomenon that was not reported in previous assessments is ocean acidification. Uptake of carbon dioxide (CO2) by oceans has a significant effect on marine biogeochemistry by reducing seawater pH. Ocean acidification is of particular concern in Alaska, because cold sea water absorbs CO2 more rapidly than warm water, and a decrease in sea ice extent has allowed increased sea surface exposure and more uptake of CO2 into these northern waters. Ocean acidification will likely affect the ability of organisms to produce and maintain shell material, such as aragonite or calcite (calcium carbonate minerals structured from carbonate ions), required by many shelled organism, from mollusks to corals to microscopic organisms at the base of the food chain. Direct biological effects in Alaska further along the food chain have yet to be studied and may vary among organisms. Some of the potentially most significant changes to Alaska that could result from a changing climate are the effects on the terrestrial cryosphere - particularly glaciers and permafrost. Alaskan glaciers are changing at a rapid rate, the primary driver appearing to be temperature. Statewide, glaciers lost 13 cubic miles of ice annually from the 1950s to the 1990s, and that rate doubled in the 2000s. However, like temperature and precipitation, glacier ice loss is not spatially uniform; most glaciers are losing mass, yet some are growing (for example Hubbard Glacier in southeast Alaska). Alaska glaciers with the most rapid loss are those terminating in sea water or lakes. With this increasing rate of melt, the contribution of surplus fresh water entering into the oceans from Alaska's glaciers, as well as those in neighboring British Columbia, Canada, is approximately 20 percent of that contributed by the Greenland Ice Sheet. Permafrost degradation (that is, the thawing of ice-rich soils) is currently (2012) impacting infrastructure and surface-water availability in areas of both discontinuous and continuous ground ice. Over most of the State, the permafrost is warming, with increasing temperatures broadly consistent with increasing air temperatures. On the Arctic coastal plain of Alaska, permafrost temperatures showed some cooling in the 1950s and 1960s but have been followed by a roughly 5&deg;F increase since the 1980s. Many areas in the continuous permafrost zone have seen increases in temperature in the seasonally active layer and a decrease in re-freezing rates. Changes in the discontinuous permafrost zone are initially much more observable due to the resulting thermokarst terrain (land surface formed as ice rich permafrost thaws), most notable in boreal forested areas. Climate warming in Alaska has potentially broad implications for human health and food security, especially in rural areas, as well as increased risk for injury with changing winter ice conditions. Additionally, such warming poses the potential for increasing damage to existing water and sanitation facilities and challenges for development of new facilities, especially in areas underlain by permafrost. Non-infectious and infectious diseases also are becoming an increasing concern. For example, from 1999 to 2006 there was a statistically significant increase in medical claims for insectbite reactions in five of six regions of Alaska, with the largest percentage increase occurring in the most northern areas. The availability and quality of subsistence foods, normally considered to be very healthy, may change due to changing access, changing habitats, and spoilage of meat in food storage cellars. These and other trends and potential outcomes resulting from a changing climate are further described in this report. In addition, we describe new science leadership activities that have been initiated to address and provide guidance toward conducting research aimed at making available information for policy makers and land management agencies to better understand, address, and plan for changes to the local and regional environment. This report cites data in both metric and standard units due to the contributions by numerous authors and the direct reference of their data.

Alaska↗

Nitrate in ground water and stream base flow in the lower Susquehanna River Basin, Pennsylvania and Maryland

High concentrations of nitrate in both ground and surface water have been identified as a significant water-quality issue in the Lower Susquehanna River Basin. This report uses data collected by the National Water Quality Assessment (NAWQA) Program in the basin and compares nitrate concentrations found in ground water and surface water on both a spatial and temporal basis and relates nitrate concentrations to land use. Nitrate concentrations in the Lower Susquehanna River Basin in Pennsylvania and Maryland were higher in ground water than in surface water in agricultural areas underlain by carbonate bedrock and agricultural areas underlain by crystalline bedrock. Nitrate concentrations were higher in surface water than in ground water in urban areas underlain by carbonate bedrock. Nitrate concentrations also were higher in surface water than ground water in both agricultural and forested areas underlain by sandstone and shale. Nitrate concentrations in ground water vary in areas with different land use and bedrock type. Ground-water nitrate concentrations were highest in agricultural areas underlain by carbonate bedrock, where 45 percent of the samples exceeded the U.S. Environmental Protection Agency (USEPA) Maximum Contaminant Level (MCL) of 10 mg/L (milligrams per liter as N). Waters from 36 percent of the wells in agricultural areas underlain by crystalline bedrock also had nitrate concentrations greater than 10 mg/L. Nitrate concentrations in water from wells in urban areas underlain by carbonate bedrock and in forested and agricultural areas underlain by sandstone and shale seldom exceeded the MCL. Nitrate concentrations were generally higher in surface water in areas underlain by carbonate bedrock than in areas underlain by noncarbonate bedrock; however, when an agricultural area underlain by carbonate bedrock and an agricultural area underlain by sandstone and shale with similar manure application rates were compared, nitrate concentrations in surface water were not significantly different. A comparison of three agricultural areas underlain by carbonate bedrock shows that the manure application rate is strongly correlated with nitrate concentration. Nitrate concentrations in stream base flow at seven sites where samples were collected throughout the year were commonly higher in the winter months than in the summer months. A statistically significant correlation between streamflow and nitrate concentration existed for six of the seven sites, indicating that seasonal variability in precipitation may be the cause of some of the seasonal variation in concentration. Other possible explanations for this variation include the seasonal cycle in plant uptake of nitrogen and seasonal fluctuations in uptake of nitrate by algae in streams. Because no information was available about the traveltime for ground water, interpretation of this temporal variation was not conclusive. Estimates of base-flow loads and yields of nitrate showed that agricultural areas underlain by carbonate bedrock provide the highest yield of nitrate when compared with the other areas studied. Agricultural areas underlain by sandstone and shale and crystalline bedrock also provide large amounts of nitrate to the river. The large amount of nitrate in the water from these areas cause a significant increase in nitrate loads transported by the Susquehanna River to the Chesapeake Bay. Urban areas underlain by carbonate bedrock had a high yield of nitrate but comprise such a small part of the basin that the nitrate load from these areas was small. In contrast, forested areas underlain by sandstone and shale bedrock had low base-flow nitrate yields, but these areas comprise a large percentage of the basin, making the overall nitrate load from these areas high.

Maryland, Pennsylvania↗

Quantification of geodetic strain rate uncertainties and implications for seismic hazard estimates

Geodetic velocity data provide first-order constraints on crustal surface strain rates, which in turn are linked to seismic hazard. Estimating the 2-D surface strain tensor everywhere requires knowledge of the surface velocity field everywhere, while geodetic data such as Global Navigation Satellite System (GNSS) only have spatially scattered measurements on the surface of the Earth. To use these data to estimate strain rates, some type of interpolation is required. In this study, we review methodologies for strain rate estimation and compare a suite of methods, including a new implementation based on the geostatistical method of kriging, to compare variation between methods with uncertainty based on one method. We estimate the velocity field and calculate strain rates in southern California using a GNSS velocity field and five different interpolation methods to understand the sources of variability in inferred strain rates. Uncertainty related to data noise and station spacing (aleatoric uncertainty) is minimal where station spacing is dense and maximum far from observations. Differences between methods, related to epistemic uncertainty, are usually highest in areas of high strain rate due to differences in how gradients in the velocity field are handled by different interpolation methods. Parameter choices, unsurprisingly, have a strong influence on strain rate field, and we propose the traditional L -curve approach as one method for quantifying the inherent trade-off between fit to the data and models that are reflective of tectonic strain rates. Doing so, we find total variability between five representative strain rate models to be roughly 40 per cent, a much lower value than roughly 100 per cent that was found in previous studies (Hearn et al .). Using multiple methods to tune parameters and calculate strain rates provides a better understanding of the range of acceptable models for a given velocity field. Finally, we present an open-source Python package (Materna et al .) for calculating strain rates, Strain_2D, which allows for the same data and model grid to be used in multiple strain rate methods, can be extended with other methods from the community, and provides an interface for comparing strain rate models, calculating statistics and estimating strain rate uncertainty for a given GNSS data set.

Geophysical Journal International↗

Assessing the role of climate and resource management on groundwater dependent ecosystem changes in arid environments with the Landsat archive

Groundwater dependent ecosystems (GDEs) rely on near-surface groundwater. These systems are receiving more attention with rising air temperature, prolonged drought, and where groundwater pumping captures natural groundwater discharge for anthropogenic use. Phreatophyte shrublands, meadows, and riparian areas are GDEs that provide critical habitat for many sensitive species, especially in arid and semi-arid environments. While GDEs are vital for ecosystem services and function, their long-term (i.e. ~ 30 years) spatial and temporal variability is poorly understood with respect to local and regional scale climate, groundwater, and rangeland management. In this work, we compute time series of NDVI derived from sensors of the Landsat TM, ETM +, and OLI lineage for assessing GDEs in a variety of land and water management contexts. Changes in vegetation vigor based on climate, groundwater availability, and land management in arid landscapes are detectable with Landsat. However, the effective quantification of these ecosystem changes can be undermined if changes in spectral bandwidths between different Landsat sensors introduce biases in derived vegetation indices, and if climate, and land and water management histories are not well understood. The objective of this work is to 1) use the Landsat 8 under-fly dataset to quantify differences in spectral reflectance and NDVI between Landsat 7 ETM + and Landsat 8 OLI for a range of vegetation communities in arid and semiarid regions of the southwestern United States, and 2) demonstrate the value of 30-year historical vegetation index and climate datasets for assessing GDEs. Specific study areas were chosen to represent a range of GDEs and environmental conditions important for three scenarios: baseline monitoring of vegetation and climate, riparian restoration, and groundwater level changes. Google's Earth Engine cloud computing and environmental monitoring platform is used to rapidly access and analyze the Landsat archive along with downscaled North American Land Data Assimilation System gridded meteorological data, which are used for both atmospheric correction and correlation analysis. Results from the cross-sensor comparison indicate a benefit from the application of a consistent atmospheric correction method, and that NDVI derived from Landsat 7 and 8 are very similar within the study area. Results from continuous Landsat time series analysis clearly illustrate that there are strong correlations between changes in vegetation vigor, precipitation, evaporative demand, depth to groundwater, and riparian restoration. Trends in summer NDVI associated with riparian restoration and groundwater level changes were found to be statistically significant, and interannual summer NDVI was found to be moderately correlated to interannual water-year precipitation for baseline study sites. Results clearly highlight the complementary relationship between water-year PPT, NDVI, and evaporative demand, and are consistent with regional vegetation index and complementary relationship studies. This work is supporting land and water managers for evaluation of GDEs with respect to climate, groundwater, and resource management.

Remote Sensing of Environment↗

Pesticides in surface water of the Mid-Atlantic region

Water-quality data from 463 surface-water sites were compiled and analyzed to document the occurrence and distribution of pesticides in surface water of the Mid-Atlantic region as part of the Mid-Atlantic Integrated Assessment program of the U.S. Environmental Protection Agency. Those data collected by the U.S. Geological Survey from October 1973 through March 1997 were used in the analyses. Data are available for a large part of the Mid-Atlantic region, but large spatial gaps in the data do exist. USGS data bases contained analyses of surface-water samples for 127 pesticide compounds, including 12 degradates, but only 16 of the compounds were commonly detected. Atrazine, metolachlor, simazine, prometon, alachlor, tebuthiuron, cyanazine, diazinon, carbaryl, chlorpyrifos, pendimethalin, 2,4-D, dieldrin, DCPA, metribuzin, and desethylatrazine (an atrazine degradate) were detected in more than 100 of the samples analyzed. At least one pesticide was detected in about 75 percent of the samples collected and at more than 90 percent of the sites sampled. Concentrations greater than the Federal Maximum Contaminant Level (MCL) for drinking water of 3 micrograms per liter (ug/L) for atrazine were found in 67 of 2,076 samples analyzed; concentrations greater than the MCL of 2ug/L for alachlor were found in 13 of 1,693 samples analyzed, and concentrations greater than the MCL of 4 ug/L for simazine were found in 17 of 1,995 samples analyzed. Concentrations of four pesticides were greater than Federal Health Advisory levels for drinking water, and concentrations of nine pesticides were greater than Federal Ambient Water-Quality Criteria for the Protection of Aquatic Organisms. Streams draining basins with different land uses tend to have different pesticide detection frequencies and median concentrations. Median concentrations of herbicides tend to be highest in streams draining basins in which the major land use is agriculture, whereas median concentrations of insecticides tend to be highest in streams draining extensively urbanized basins. Concentrations of both herbicides and insecticides are usually highest during the spring and summer, although many pesticides are present at low concentrations in surface water throughout the year. Pesticide concentrations vary greatly seasonally and over different hydrologic conditions, with overall variation sometimes exceeding four orders of magnitude. During periods of pesticide application (typically spring and summer), the occurrence of selected pesticides in some streams in the Mid-Atlantic region is related to streamflow. Correlations between concentrations of selected pesticides and streamflow are statistically significant during spring and summer for small (draining less than 55 square miles) streams. Concentrations of selected pesticides in small streams increase during high flows in the growing season, up to 30 times the concentrations present during low-flow conditions in the growing season. In small streams draining urban areas, concentrations of atrazine decrease during high-flow events but concentrations of the insecticides diazinon and chlorpyrifos increase. This may be due to the differences in the pesticides used in agricultural and urban areas and the amounts applied.

Water-Resources Investigations Report↗

Yuma Ridgway’s rail selenium exposure and occupancy within managed and unmanaged emergent marshes at the Salton Sea

Yuma Ridgway’s rail ( Rallus obsoletus yumanensis , hereafter, rail) is an endangered species for which patches of emergent marsh within the Salton Sea watershed comprise a substantial part of habitat for the species’ disjointed range in the southwestern United States. These areas of emergent marsh include (1) marshes managed by federal (particularly the U.S. Fish and Wildlife Service’s Sonny Bono Salton Sea National Wildlife Refuge), state (California Department of Fish and Wildlife), and local (Imperial Irrigation District) resource agencies that are sustained by direct deliveries of Colorado River water and (2) unmanaged marshes sustained by agricultural drainage water. Management of rail habitat in this arid environment is complicated by increasingly limited availability of unimpaired freshwater owing to water management decisions associated with the Quantification Settlement Agreement and risks posed by potentially harmful concentrations of selenium found in agricultural drainage water that can readily bioaccumulate in aquatic food webs. To provide timely science for managers, herein we report summary statistics for managed and unmanaged emergent marshes sampled at the Salton Sea during the rail breeding season of 2016 pertaining to (1) selenium concentrations in food webs representing dietary pathways of selenium exposure and (2) patterns of rail occupancy and inter-marsh movements, estimated abundance, and regional population size of rail. For selenium-specific objectives, we sampled unfiltered surface water, midge larvae (Chironomidae), water boatmen (Corixidae), mosquitofish ( Gambusia spp.), and crayfish (Astacidae). Selenium samples were collected from 15 fixed sampling points, each in managed and unmanaged marshes, during late February, April, and June 2016, which corresponded to rail pre-nesting, nesting, and fledgling reproductive life-stages, respectively. Two areas within the two treatment types (managed versus unmanaged marsh) were of particular interest to help assess risks associated with changing sea dynamics and different water-management strategies: (1) a large unmanaged marsh (Morton Bay) unintentionally created in approximately 2008 when it became separated from the Salton Sea as water inflows began to drop and a berm formed from accumulated sediment and (2) a restored marsh (HZ9A) managed by the Sonny Bono Salton Sea National Wildlife Refuge, which is currently supplied with Colorado River water but may be sustained in the future by a blend of clean (that is, low selenium) Colorado River and agricultural drainage water with higher selenium from the Alamo River. Hence, baseline data for these marshes are important for future management decisions. We also report selenium concentrations in rail blood, head feathers, and breast feathers from rails captured as part of the movement study. Results indicated relatively higher risks from dietary selenium exposure for rails occupying unmanaged marshes compared to managed marshes and similar risks among unmanaged marshes. However, risks also were potentially elevated for rails occupying some managed marshes (that is, the Hazard Marshes), where relatively high proportions of Chironomidae and mosquitofish exceeded dietary thresholds for selenium effects on avian reproduction. For rail-specific objectives, we quantified occupancy and spatial distribution using call count data analyzed with imperfect detection models. Imperfect detection models allowed us to jointly estimate detection probability and abundance of detected rails in association with habitats. We then used estimates of detection probability and abundance at the habitat level to extrapolate rail population abundance for the Salton Sea region. Inter- and intra-marsh movements were described from over 5,000 locations obtained from 15 radio-marked rails. Resultant space use patterns indicated that, in general, selenium risk to individuals is not equally shared because of high levels of territoriality and very limited movement throughout the landscape. Moreover, the largest contiguous blocks of habitat are associated with unmanaged marshlands located on the former southeastern shoreline and outside traditional management areas and authorities. Thus, a substantial proportion of the rail population that is using unmanaged marsh on the southeastern shoreline may have disproportionate risk of elevated selenium exposure, yet how that risk translates to population-level effects remains unknown.

California↗

Fena Valley Reservoir watershed and water-balance model updates and expansion of watershed modeling to southern Guam

In 2014, the U.S. Geological Survey, in cooperation with the U.S. Department of Defense’s Strategic Environmental Research and Development Program, initiated a project to evaluate the potential impacts of projected climate-change on Department of Defense installations that rely on Guam’s water resources. A major task of that project was to develop a watershed model of southern Guam and a water-balance model for the Fena Valley Reservoir. The southern Guam watershed model provides a physically based tool to estimate surface-water availability in southern Guam. The U.S. Geological Survey’s Precipitation Runoff Modeling System, PRMS-IV, was used to construct the watershed model. The PRMS-IV code simulates different parts of the hydrologic cycle based on a set of user-defined modules. The southern Guam watershed model was constructed by updating a watershed model for the Fena Valley watersheds, and expanding the modeled area to include all of southern Guam. The Fena Valley watershed model was combined with a previously developed, but recently updated and recalibrated Fena Valley Reservoir water-balance model. Two important surface-water resources for the U.S. Navy and the citizens of Guam were modeled in this study; the extended model now includes the Ugum River watershed and improves upon the previous model of the Fena Valley watersheds. Surface water from the Ugum River watershed is diverted and treated for drinking water, and the Fena Valley watersheds feed the largest surface-water reservoir on Guam. The southern Guam watershed model performed “very good,” according to the criteria of Moriasi and others (2007), in the Ugum River watershed above Talofofo Falls with monthly Nash-Sutcliffe efficiency statistic values of 0.97 for the calibration period and 0.93 for the verification period (a value of 1.0 represents perfect model fit). In the Fena Valley watershed, monthly simulated streamflow volumes from the watershed model compared reasonably well with the measured values for the gaging stations on the Almagosa, Maulap, and Imong Rivers—tributaries to the Fena Valley Reservoir—with Nash-Sutcliffe efficiency values of 0.87 or higher. The southern Guam watershed model simulated the total volume of the critical dry season (January to May) streamflow for the entire simulation period within –0.54 percent at the Almagosa River, within 6.39 percent at the Maulap River, and within 6.06 percent at the Imong River. The recalibrated water-balance model of the Fena Valley Reservoir generally simulated monthly reservoir storage volume with reasonable accuracy. For the calibration and verification periods, errors in end-of-month reservoir-storage volume ranged from 6.04 percent (284.6 acre-feet or 92.7 million gallons) to –5.70 percent (–240.8 acre-feet or –78.5 million gallons). Monthly simulation bias ranged from –0.48 percent for the calibration period to 0.87 percent for the verification period; relative error ranged from –0.60 to 0.88 percent for the calibration and verification periods, respectively. The small bias indicated that the model did not consistently overestimate or underestimate reservoir storage volume. In the entirety of southern Guam, the watershed model has a “satisfactory” to “very good” rating when simulating monthly mean streamflow for all but one of the gaged watersheds during the verification period. The southern Guam watershed model uses a more sophisticated climate-distribution scheme than the older model to make use of the sparse climate data, as well as includes updated land-cover parameters and the capability to simulate closed depression areas. The new Fena Valley Reservoir water-balance model is useful as an updated tool to forecast short-term changes in the surface-water resources of Guam. Furthermore, the now spatially complete southern Guam watershed model can be used to evaluate changes in streamflow and recharge owing to climate or land-cover changes. These are substantial improvements to the previous models of the Fena Valley watershed and Reservoir. Datasets associated with this report are available as a U.S. Geological Survey data release (Rosa and Hay, 2017; DOI:10.5066/F7HH6HV4).

Scientific Investigations Report↗

Effects of selected low-impact-development (LID) techniques on water quality and quantity in the Ipswich River Basin, Massachusetts: Field and modeling studies

During the months of August and September, flows in the Ipswich River, Massachusetts, dramatically decrease largely due to groundwater withdrawals needed to meet increased residential and commercial water demands. In the summer, rates of groundwater recharge are lower than during the rest of the year, and water demands are higher. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of low-impact-development (LID) enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of (1) replacing an impervious parking lot surface with a porous surface on groundwater quality, (2) installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and (3) installing a 3,000-square foot (ft2) green roof on the quantity and quality of stormwater runoff. In addition, the effects of broad-scale implementation of LID techniques, reduced water withdrawals, and water-conservation measures on streamflow in large areas of the basin were simulated using the U.S. Geological Survey’s Ipswich River Basin model. From June 2005 to 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, MA, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers. Changes in the concentrations of the water-quality constituents, phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons, were monitored. Increased infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking lot surface. In Wilmington, MA, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the above constituents to Silver Lake. Flow-proportional water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and Escherichia coli bacteria. In general, when all storms were considered, no substantial decreases were observed in runoff volume as a result of installing LID enhancements. However, the relation between rainfall and runoff did provide some insight into how the LID enhancements affected the effective impervious area for the neighborhood. A decrease in runoff was observed for storms of 0.2 inches (in.) or less of precipitation, which indicated a reduction in effective impervious area from approximately 10 percent to about 4.5 percent for the 3-acre area. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. Three factors were probably most important in minimizing differences: (1) the small decrease in effective impervious area, (2) the differences in the size of storms sampled for water-quality constituents before and after installation of the infiltration enhancing measures, and (3) small sample sizes. In a third field study, the characteristics of runoff from a vegetated “green” roof and a conventional, rubber-membrane roof were compared. The amount of precipitation and the length of the antecedent dry period were the two primary factors affecting the green roof’s water-storage capacity. The green roof retained more than 50 percent of the precipitation from storms with 0.04 to 1.0 in. of rain. Approximately 95 percent of the precipitation from one storm of nearly 2 in. was retained by the green roof. On the rubber-membrane roof, only a small, shallow puddle of insubstantial volume ever remained after a storm. Bulk precipitation from 10 storms was monitored for the same constituents (nutrients, metals, and total petroleum hydrocarbons) as the roof runoff, and the results were compared with those for roof-runoff samples. The use of fertilizers to help establish the vegetation during the study probably distorted any effect the plants and growing medium may have had on the retention of target analytes. As a result of the fertilizer and growing medium chemistry, median concentrations of total nitrogen, total phosphorus, cadmium, copper, and nickel in runoff from the green roof were greater than in the runoff from the conventional roof or in bulk precipitation. Concentrations of lead and zinc were greater in runoff from the conventional roof, probably a result of passage through the old, metal drainpipes. Simulations of the effects of LID on streamflow in the Ipswich River Basin were conducted with a previously calibrated Hydrological Simulation Program-FORTRAN (HSPF) precipitation-runoff model. Simulations were conducted at multiple spatial scales to evaluate the effects of (1) updated water withdrawals for the towns of Reading and Wilmington; (2) potential land-use changes at buildout (potential future development); (3) effective impervious area reductions upstream from the South Middleton streamgage to represent the effects of widespread implementation of LID retrofit techniques; (4) basin-scale water withdrawal reductions scaled up (expanded to the town level) from water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation; and (5) land-use change and LID techniques at a local scale, which is smaller than the HSPF subbasin. Effects on streamflow generally were evaluated by comparing results of two or more related simulations for selected reaches in the basin; thus, relative rather than absolute changes in simulated flow were the focus of the assessment. Simulations indicated that reduced withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage. Simulations of water-conservation measures resulted in negligible effects on streamflow. Overall, simulations indicated that spatial scale is an important factor in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest (percent differences of less than 20 percent) effects on streamflow in most subbasins because relatively little land in the basin was available for development (about 17 percent); moreover, most of the available open land is zoned for low-density residential development, and this land-use category was simulated to contain relatively little effective impervious area and to be similar hydrologically to the forested land in place prior to development. Results of the simulations conducted to evaluate widespread effective impervious area reductions upstream from the South Middleton streamgage indicated that the percentage of urban land use and associated effective impervious area was too small for a 50-percent reduction of effective impervious area to appreciably affect streamflow (percent differences of less than 20 percent) in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, where the percentage of urban land use and associated effective impervious area in the drainage area may be substantially higher, land-use change, development patterns, and LID practices potentially have much greater effects on streamflow. Modeling results also indicated that LID was potentially most beneficial for minimizing streamflow alteration when applied to dense urban development, largely because larger tracts of effective impervious area were available for reduction than were available for other land-use categories. For example, commercial-industrial-transportation land use is composed of 37 percent pervious area and 63 percent effective impervious area in the HSPF model, whereas low-density residential area is composed of 97.5 percent pervious area and only 2.5 percent effective impervious area. Field and modeling studies concurred in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas. A measurable effect for small rainfall events (less than 0.25 inch) was determined in the small, highly pervious area that was monitored in this study, but the volume difference was not great.

Massachusetts↗

Lahar hazard zones for eruption-generated lahars in the Lassen Volcanic Center, California

Lahar deposits are found in drainages that head on or near Lassen Peak in northern California, demonstrating that these valleys are susceptible to future lahars. In general, lahars are uncommon in the Lassen region. Lassen Peak's lack of large perennial snowfields and glaciers limits its potential for lahar development, with the winter snowpack being the largest source of water for lahar generation. The most extensive lahar deposits are related to the May 1915 eruption of Lassen Peak, and evidence for pre-1915 lahars is sparse and spatially limited. The May 1915 eruption of Lassen Peak was a small-volume eruption that generated a snow and hot-rock avalanche, a pyroclastic flow, and two large and four smaller lahars. The two large lahars were generated on May 19 and 22 and inundated sections of Lost and Hat Creeks. We use 80 years of snow depth measurements from Lassen Peak to calculate average and maximum liquid water depths, 2.02 meters (m) and 3.90 m respectively, for the month of May as estimates of the 1915 lahars. These depths are multiplied by the areal extents of the eruptive deposits to calculate a water volume range, 7.05-13.6x10 6 cubic meters (m 3 ). We assume the lahars were a 50/50 mix of water and sediment and double the water volumes to provide an estimate of the 1915 lahars, 13.2-19.8x10 6 m 3 . We use a representative volume of 15x106 m 3 in the software program LAHARZ to calculate cross-sectional and planimetric areas for the 1915 lahars. The resultant lahar inundation zone reasonably portrays both of the May 1915 lahars. We use this same technique to calculate the potential for future lahars in basins that head on or near Lassen Peak. LAHARZ assumes that the total lahar volume does not change after leaving the potential energy, H/L , cone (the height of the edifice, H , down to the approximate break in slope at its base, L ); therefore, all water available to initiate a lahar is contained inside this cone. Because snow is the primary source of water for lahar generation, we assume that the maximum historical water equivalent, 3.90 m, covers the entire basin area inside the H/L cone. The product of planimetric area of each basin inside the H/L and the maximum historical water equivalent yields the maximum water volume available to generate a lahar. We then double the water volumes to approximate maximum lahar volumes. The maximum lahar volumes and an understanding of the statistical uncertainties inherent to the LAHARZ calculations guided our selection of six hypothetical volumes, 1, 3, 10, 30, 60, and 90x10 6 m 3 , to delineate concentric lahar inundation zones. The lahar inundation zones extend, in general, tens of kilometers away from Lassen Peak. The small, more-frequent lahar inundation zones (1 and 3x10 6 m 3 ) are, on average, 10 km long. The exceptions are the zones in Warner Creek and Mill Creek, which extend much further. All but one of the small, more-frequent lahar inundation zones reach outside of the Lassen Volcanic National Park boundary, and the zone in Mill Creek extends well past the park boundary. All of the medium, moderately frequent lahar inundation zones (10 and 30x10 6 m 3 ) extend past the park boundary and could potentially impact the communities of Viola and Old Station and State Highways 36 and 44, both north and west of Lassen Peak. The approximately 27-km-long on average, large, less-frequent lahar inundation zones (60 and 90x10 6 m 3 ) represent worst-case lahar scenarios that are unlikely to occur. Flood hazards continue downstream from the toes of the lahars, potentially affecting communities in the Sacramento River Valley.

California↗