USGS ScienceSearch

Geology topics

Kelly O. Maloney

Publications and source records attributed to Kelly O. Maloney.

At least 37 records · Page 2Linked to original sources

Techniques to improve ecological interpretability of black box machine learning models

Statistical modeling of ecological data is often faced with a large number of variables as well as possible nonlinear relationships and higher-order interaction effects. Gradient boosted trees (GBT) have been successful in addressing these issues and have shown a good predictive performance in modeling nonlinear relationships, in particular in classification settings with a categorical response variable. They also tend to be robust against outliers. However, their black-box nature makes it difficult to interpret these models. We introduce several recently developed statistical tools to the environmental research community in order to advance interpretation of these black-box models. To analyze the properties of the tools, we applied gradient boosted trees to investigate biological health of streams within the contiguous USA, as measured by a benthic macroinvertebrate biotic index. Based on these data and a simulation study, we demonstrate the advantages and limitations of partial dependence plots (PDP), individual conditional expectation (ICE) curves and accumulated local effects (ALE) in their ability to identify covariate–response relationships. Additionally, interaction effects were quantified according to interaction strength (IAS) and Friedman’s H 2 "> H 2 statistic. Interpretable machine learning techniques are useful tools to open the black-box of gradient boosted trees in the environmental sciences. This finding is supported by our case study on the effect of impervious surface on the benthic condition, which agrees with previous results in the literature. Overall, the most important variables were ecoregion, bed stability, watershed area, riparian vegetation and catchment slope. These variables were also present in most identified interaction effects. In conclusion, graphical tools (PDP, ICE, ALE) enable visualization and easier interpretation of GBT but should be supported by analytical statistical measures. Future methodological research is needed to investigate the properties of interaction tests. Supplementary materials accompanying this paper appear on-line.

Journal of Agricultural, Biological, and Environme

Using expert knowledge to support Endangered Species Act decision‐making for data‐deficient species

Many questions relevant to conservation decision making are characterized by extreme uncertainty due to lack of empirical data and complexity of the underlying ecological processes, leading to a rapid increase in the use of structured protocols to elicit expert knowledge. Published ecological applications often employ a modified Delphi method, where experts provide judgments anonymously and mathematical aggregation techniques are used to combine judgments. The Sheffield Elicitation Framework (SHELF) differs in its behavioral approach to synthesizing individual judgments into a fully specified probability distribution for an unknown quantity. This study demonstrates the remote use of the SHELF protocol for an extinction risk assessment of three subterranean aquatic species petitioned for listing under the US Endangered Species Act. Experts were provided an empirical threat assessment for each known locality using video conferencing and asked for judgments on the probability of population persistence over four generations using online submission forms and R‐shiny apps available through the SHELF package. Despite large uncertainty for all populations, results reveal key differences between species’ risk of extirpation based on spatial variation in dominant threats, local land use and management practices, and microhabitat use. The resulting probability distributions provide decision makers with a full picture of uncertainty that is consistent with the probabilistic nature of risk assessments, and discussions during the behavioral aggregation stage clearly document dominant threats (e.g., development, timber harvest, animal agriculture, and cave visitation) and their interactions with local cave geology and species’ habitat preferences. Our virtual implementation of the SHELF protocol demonstrates the flexibility of this approach for conservation applications operating on budgets and timelines that can limit in‐person meetings of geographically dispersed experts.

Conservation Biology

Linking altered flow regimes to biological condition: An example using benthic macroinvertebrates in small streams of the Chesapeake Bay watershed

Regionally scaled assessments of hydrologic alteration for small streams and its effects on freshwater taxa are often inhibited by a low number of stream gages. To overcome this limitation, we paired modeled estimates of hydrologic alteration to a benthic macroinvertebrate index of biotic integrity data for 4522 stream reaches across the Chesapeake Bay watershed. Using separate random-forest models, we predicted flow status (inflated, diminished, or indeterminant) for 12 published hydrologic metrics (HMs) that characterize the main components of flow regimes. We used these models to predict each HM status for each stream reach in the watershed, and linked predictions to macroinvertebrate condition samples collected from streams with drainage areas less than 200 km 2 . Flow alteration was calculated as the number of HMs with inflated or diminished status and ranged from 0 (no HM inflated or diminished) to 12 (all 12 HMs inflated or diminished). When focused solely on the stream condition and flow-alteration relationship, degraded macroinvertebrate condition was, depending on the number of HMs used, 3.8–4.7 times more likely in a flow-altered site; this likelihood was over twofold higher in the urban-focused dataset (8.7–10.8), and was never significant in the agriculture-focused dataset. Logistic regression analysis using the entire dataset showed for every unit increase in flow-alteration intensity, the odds of a degraded condition increased 3.7%. Our results provide an indication of whether altered streamflow is a possible driver of degraded biological conditions, information that could help managers prioritize management actions and lead to more effective restoration efforts.

Chesapeake Bay watershed

Time marches on, but do the causal pathways driving instream habitat and biology remain consistent?

Stream ecosystems are complex networks of interacting terrestrial and aquatic drivers. To untangle these ecological networks, efforts evaluating the direct and indirect effects of landscape, climate, and instream predictors on biological condition through time are needed. We used structural equation modeling and leveraged a stream survey program to identify and compare important predictors driving condition of benthic macroinvertebrate and fish assemblages. We used data resampled 14 years apart at 252 locations across Maryland, USA. Sample locations covered a wide range of conditions that varied spatiotemporally. Overall, the relationship directions were consistent between sample periods, but their relative strength varied temporally. For benthic macroinvertebrates, we found that the total effect of natural landscape ( e.g ., elevation, longitude, latitude, geology) and land use ( i.e ., forest, development, agriculture) predictors was 1.4 and 1.5 times greater in the late 2010s compared to the 2000s. Moreover, the total effect of water quality ( e.g ., total nitrogen and conductivity) and habitat ( e.g ., embeddedness, riffle quality) was 1.2 and 4.8 times lower in the 2010s, respectively. For fish assemblage condition, the total effect of land use-land cover predictors was 2.3 times greater in the 2010s compared to the 2000s, while the total effect of local habitat was 1.4 times lower in the 2010s, respectively. As expected, we found biological assemblages in catchments with more agriculture and urban development were generally comprised of tolerant, generalist species, while assemblages in catchments with greater forest cover had more-specialized, less-tolerant species ( e.g., Ephemeroptera, Plecoptera, and Trichoptera taxa, clingers, benthic and lithophilic spawning fishes). Changes in the relative importance of landscape and land-use predictors suggest other correlated, yet unmeasured, proximal factors became more important over time. By untangling these ecological networks, stakeholders can gain a better understanding of the spatiotemporal relationships driving biological condition to implement management practices aimed at improving stream condition.

Maryland

Sediment dynamics and implications for management: State of the science from long‐term research in the Chesapeake Bay watershed, USA

This review aims to synthesize the current knowledge of sediment dynamics using insights from long‐term research conducted in the watershed draining to the Chesapeake Bay, the largest estuary in the U.S., to inform management actions to restore the estuary and its watershed. The sediment dynamics of the Chesapeake are typical of many impaired watersheds and estuaries around the world, and this synthesis is intended to be relevant and transferable to other sediment‐impaired systems. The watershed's sediment sources, transport, delivery, and impacts are discussed with implications for effectively implementing best management practices (BMPs) to mitigate sediment issues. This synthesis revealed three key issues to consider when planning actions to reduce sediment loading: Scale, time, and land use. Geology and historical land use generated a template that current land use and climate, in addition to management, are acting upon to control sediment delivery. Important sediment sources in the Chesapeake include the Piedmont physiographic region, urban, and agricultural land use, and streambank erosion of headwater streams, whereas floodplain trapping is important along larger streams and rivers. Implementation of BMPs is widespread and is predicted to lead to decreased sediment loading; however, reworking of legacy sediment stored in stream valleys, with potentially long residence times in storage, can delay and complicate detection of the effects of BMPs on sediment loads. In conclusion, the improved understanding of sediment sources, storage areas, and transport lag times reviewed here can help target choices of BMP types and locations to better manage sediment problems—for both local streams and receiving waters.

New York, Pennsylvania, Maryland, Delaware, Virgin

A random forest approach for bounded outcome variables

Random forests have become an established tool for classication and regres- sion, in particular in high-dimensional settings and in the presence of non-additive predictor-response relationships. For bounded outcome variables restricted to the unit interval, however, classical modeling approaches based on mean squared error loss may severely suer as they do not account for heteroscedasticity in the data. To address this issue, we propose a random forest approach for relating a beta dis- tributed outcome to a set of explanatory variables. Our approach explicitly makes use of the likelihood function of the beta distribution for the selection of splits dur- ing the tree-building procedure. In each iteration of the tree-building algorithm it chooses one explanatory variable in combination with a split point that maximizes the log-likelihood function of the beta distribution with the parameter estimates de- rived from the nodes of the currently built tree. Results of several simulation studies and an application using data from the U.S.A. National Lakes Assessment Survey demonstrate the properties and usefulness of the method, in particular when com- pared to random forest approaches based on mean squared error loss and parametric regression models.

Journal of Computational and Graphical Statistics

Shale gas development has limited effects on stream biology and geochemistry in a gradient-based, multiparameter study in Pennsylvania

The number of horizontally drilled shale oil and gas wells in the United States has increased from nearly 28,000 in 2007 to nearly 127,000 in 2017, and research has suggested the potential for the development of shale resources to affect nearby stream ecosystems. However, the ability to generalize current studies is limited by the small geographic scope as well as limited breadth and integration of measured chemical and biological indicators parameters. This study tested the hypothesis that a quantifiable, significant relationship exists between the density of oil and gas (OG) development, increasing streamwater concentrations of known geochemical tracers of OG extraction, and the composition of benthic macroinvertebrate and microbial communities. 25 headwater streams that drain lands across a gradient of shale gas development intensity were sampled. Our strategy included comprehensive measurements across multiple seasons of sampling to account for temporal variability of geochemical parameters, including known shale OG geochemical tracers, and microbial and benthic macroinvertebrate communities. No significant relationships were found between the intensity of OG development, shale OG geochemical tracers, or benthic macroinvertebrate or microbial community composition, whereas significant seasonal differences in stream chemistry were observed. These results highlight the importance of considering spatial and temporal variability in stream chemistry and biota and not only the presence of anthropogenic activities in a watershed. This comprehensive, integrated study of geochemical and biological variability of headwater streams in watersheds undergoing OG development provides a robust framework for examining the effects of energy development at a regional scale.

Pennsylvania

Disentangling the potential effects of land-use and climate change on stream conditions

Land‐use and climate change are significantly affecting stream ecosystems, yet understanding of their long‐term impacts is hindered by the few studies that have simultaneously investigated their interaction and high variability among future projections. We modeled possible effects of a suite of 2030, 2060, and 2090 land‐use and climate scenarios on the condition of 70,772 small streams in the Chesapeake Bay watershed, United States. The Chesapeake Basin‐wide Index of Biotic Integrity, a benthic macroinvertebrate multimetric index, was used to represent stream condition. Land‐use scenarios included four Special Report on Emissions Scenarios (A1B, A2, B1, and B2) representing a range of potential landscape futures. Future climate scenarios included quartiles of future climate changes from downscaled Coupled Model Intercomparison Project ‐ Phase 5 (CMIP5) and a watershed‐wide uniform scenario (Lynch2016). We employed random forests analysis to model individual and combined effects of land‐use and climate change on stream conditions. Individual scenarios suggest that by 2090, watershed‐wide conditions may exhibit anywhere from large degradations (e.g., scenarios A1B, A2, and the CMIP5 25th percentile) to small degradations (e.g., scenarios B1, B2, and Lynch2016). Combined land‐use and climate change scenarios highlighted their interaction and predicted, by 2090, watershed‐wide degradation in 16.2% (A2 CMIP5 25th percentile) to 1.0% (B2 Lynch2016) of stream kilometers. A goal for the Chesapeake Bay watershed is to restore 10% of stream kilometers over a 2008 baseline; our results suggest meeting and sustaining this goal until 2090 may require improvement in 11.0%–26.2% of stream kilometers, dependent on land‐use and climate scenario. These results highlight inherent variability among scenarios and the resultant uncertainty of predicted conditions, which reinforces the need to incorporate multiple scenarios of both land‐use (e.g., development, agriculture, etc.) and climate change in future studies to encapsulate the range of potential future conditions.

Delaware, Maryland, New York, Pennsylvania, Virgin

A Bayesian framework for assessing extinction risk based on ordinal categories of population condition and projected landscape change

Many at-risk species lack standardized surveys across their range or quantitative data capable of detecting demographic trends. As a result, extinction risk assessments often rely on ordinal categories of risk based on explicit criteria or expert elicitation. This study demonstrates a Bayesian approach to assessing extinction risk based on this common data structure, using three freshwater mussel species being considered for listing under the US Endangered Species Act. The probability that a population is classified under each risk category was modeled as a function of projected landscape change using ordered probit regression, assuming observed categories reflect a latent, continuous probability of persistence. All three species were more likely than not (mean probability >0.5) to be classified as extirpated or low condition throughout their range based on effects of urban development and hydrologic alteration. Spatial variation in estimates revealed strongholds and high-risk areas relevant to conservation decision making. Projected change in probabilities of each risk category based on multiple land-use and climate models was generally small relative to high baseline risk resulting from past landscape changes. Assessing extinction risk based on probabilities of ordinal condition as a function of landscape patterns may provide a flexible and robust approach for many at-risk taxa by adjusting species' demographic criteria to match relative risk categories, following standardized criteria, or using expert elicitation for data-deficient species. This approach provides decision makers with a useful measure of uncertainty around ordinal classifications and provides a framework for estimating future risk based on projections of anthropogenic stressors.

Alabama, Georgia, Kentucky, North Carolina, South

A preliminary assessment of hyperspectral remote sensing technology for mapping submerged aquatic vegetation in the Upper Delaware River National Parks

Hyperspectral remote sensing of submerged aquatic vegetation is a complex and difficult process that is affected by unique constraints on the energy flow profile near and below the water surface. In addition, shallow, winding, lotic systems, such as the Upper Delaware River, present additional remote sensing problems in the form of specular reflectance, variable depth and constituents in the water column and sometimes extremely weak signal strength due to absorption and scattering in the water column that can be statistically overwhelmed by the reflectance from upland vegetation in any individual image scene. Here we test hyperspectral imagery from the Civil Air Patrol’s (CAP), Airborne Real-time Cueing Hyperspectral Enhanced Recon (ARCHER) system in the scenic waters of two National Parks on the Upper Delaware River. A number of unique image processing problems were encountered, including specular reflectance from winding lotic systems, variable depth and flow dynamics of the riverine environment, and disproportionate signal strength from surface reflectance in this riverine environment. This were solved by applying a specular reflectance removal algorithm, applying field data collections to classification results and masking upland vegetation so as to not statistically overwhelm the weak reflectance signal from surface and near-surface water. Much was learned about conducting imaging spectroscopy in such difficult conditions. Significant results include successful mapping of SAV presenece/absence, advantages of upland masking of the reflectance signal, and a numkber of processing approaches that are unique to this environment. In this paper we summarize our results and identify unique issues that must be addressed in this environment

New Jersey, New York, Pennsylvania

Water stress from high-volume hydraulic fracturing potentially threatens aquatic biodiversity and ecosystem services in Arkansas, United States

Demand for high-volume, short duration water withdrawals could create water stress to aquatic organisms in Fayetteville Shale streams sourced for hydraulic fracturing fluids. We estimated potential water stress using permitted water withdrawal volumes and actual water withdrawals compared to monthly median, low, and high streamflows. Risk for biological stress was considered at 20% of long-term median and 10% of high- and low-flow thresholds. Future well build-out projections estimated potential for continued stress. Most water was permitted from small, free-flowing streams and “frack” ponds (dammed streams). Permitted 12-h pumping volumes exceeded median streamflow at 50% of withdrawal sites in June, when flows were low. Daily water usage, from operator disclosures, compared to median streamflow showed possible water stress in 7–51% of catchments from June–November, respectively. If 100% of produced water was recycled, per-well water use declined by 25%, reducing threshold exceedance by 10%. Future water stress was predicted to occur in fewer catchments important for drinking water and species of conservation concern due to the decline in new well installations and increased use of recycled water. Accessible and precise withdrawal and streamflow data are critical moving forward to assess and mitigate water stress in streams that experience high-volume withdrawals.

Arkansas

Canopy volume removal from oil and gas development activity in the upper Susquehanna River basin in Pennsylvania and New York (USA): An assessment using lidar data

Oil and gas development is changing the landscape in many regions of the United States and globally. However, the nature, extent, and magnitude of landscape change and development, and precisely how this development compares to other ongoing land conversion (e.g. urban/sub-urban development, timber harvest) is not well understood. In this study, we examine land conversion from oil and gas infrastructure development in the upper Susquehanna River basin in Pennsylvania and New York, an area that has experienced much oil and gas development over the past 10 years. We quantified land conversion in terms of forest canopy geometric volume loss in contrast to previous studies that considered only areal impacts. For the first time in a study of this type, we use fine-scale lidar forest canopy geometric models to assess the volumetric change due to forest clearing from oil and gas development and contrast this land change to clear cut forest harvesting, and urban and suburban development. Results show that oil and gas infrastructure development removed a large volume of forest canopy from 2006 to 2013, and this removal spread over a large portion of the study area. Timber operations (clear cutting) on Pennsylvania State Forest lands removed a larger total volume of forest canopy during the same time period, but this canopy removal was concentrated in a smaller area. Results of our study point to the need to consider volumetric impacts of oil and gas development on ecosystems, and to place potential impacts in context with other ongoing land conversions.

New York, Pennsylvania

A detailed risk assessment of shale gas development on headwater streams in the Pennsylvania portion of the Upper Susquehanna River Basin, U.S.A.

The development of unconventional oil and gas (UOG) involves infrastructure development (well pads, roads and pipelines), well drilling and stimulation (hydraulic fracturing), and production; all of which have the potential to affect stream ecosystems. Here, we developed a fine-scaled (1:24,000) catchment-level disturbance intensity index (DII) that included 17 measures of UOG capturing all steps in the development process (infrastructure, water withdrawals, probabilistic spills) that could affect headwater streams (< 200 km 2 in upstream catchment) in the Upper Susquehanna River Basin in Pennsylvania, U.S.A. The DII ranged from 0 (no UOG disturbance) to 100 (the catchment with the highest UOG disturbance in the study area) and it was most sensitive to removal of pipeline cover, road cover and well pad cover metrics. We related this DII to three measures of high quality streams: Pennsylvania State Exceptional Value (EV) streams, Class A brook trout streams and Eastern Brook Trout Joint Venture brook trout patches. Overall only 3.8% of all catchments and 2.7% of EV stream length, 1.9% of Class A streams and 1.2% of patches were classified as having medium to high level DII scores (> 50). Well density, often used as a proxy for development, only correlated strongly with well pad coverage and produced materials, and therefore may miss potential effects associated with roads and pipelines, water withdrawals and spills. When analyzed with a future development scenario, 91.1% of EV stream length, 68.7% of Class A streams and 80.0% of patches were in catchments with a moderate to high probability of development. Our method incorporated the cumulative effects of UOG on streams and can be used to identify catchments and reaches at risk to existing stressors or future development.

Pennsylvania

Brook trout distributional response to unconventional oil and gas development: Landscape context matters

We conducted a large-scale assessment of unconventional oil and gas (UOG) development effects on brook trout ( Salvelinus fontinalis ) distribution. We compiled 2231 brook trout collection records from the Upper Susquehanna River Watershed, USA. We used boosted regression tree (BRT) analysis to predict occurrence probability at the 1:24,000 stream-segment scale as a function of natural and anthropogenic landscape and climatic attributes. We then evaluated the importance of landscape context (i.e., pre-existing natural habitat quality and anthropogenic degradation) in modulating the effects of UOG on brook trout distribution under UOG development scenarios. BRT made use of 5 anthropogenic (28% relative influence) and 7 natural (72% relative influence) variables to model occurrence with a high degree of accuracy [Area Under the Receiver Operating Curve (AUC) = 0.85 and cross-validated AUC = 0.81]. UOG development impacted 11% ( n = 2784) of streams and resulted in a loss of predicted occurrence in 126 (4%). Most streams impacted by UOG had unsuitable underlying natural habitat quality ( n = 1220; 44%). Brook trout were predicted to be absent from an additional 26% ( n = 733) of streams due to pre-existing non-UOG land uses (i.e., agriculture, residential and commercial development, or historic mining). Streams with a predicted and observed (via existing pre- and post-disturbance fish sampling records) loss of occurrence due to UOG tended to have intermediate natural habitat quality and/or intermediate levels of non-UOG stress. Simulated development of permitted but undeveloped UOG wells ( n = 943) resulted in a loss of predicted occurrence in 27 additional streams. Loss of occurrence was strongly dependent upon landscape context, suggesting effects of current and future UOG development are likely most relevant in streams near the probability threshold due to pre-existing habitat degradation.

New York, Pennsylvania

Predicting biological conditions for small headwater streams in the Chesapeake Bay watershed

A primary goal for Chesapeake Bay watershed restoration is to improve stream health and function in 10% of stream miles by 2025. Predictive spatial modeling of stream conditions, when accurate, is one method to fill gaps in monitoring coverage and estimate baseline conditions for restoration goals. Predictive modeling can also monitor progress as additional data become available. We developed a random forests model to predict biological condition of small streams (<200 km 2 in drainage) in the Chesapeake Bay watershed. Biological condition was measured with the Chesapeake Bay Basin-wide Index of Biotic Integrity (Chessie BIBI), a stream macroinvertebrate index. Our goal was to predict biological condition in all unsurveyed small streams present in a 1:24,000 scale catchment layer as a 2004–2008 baseline. We reclassified the 5-category Chessie BIBI ratings into two categories, poor and fair/good, to align with management goals of the Chesapeake Bay Program. The model included 12 geospatial predictor variables including measures on spatial location, bioregion, land cover, soil, precipitation, and number of dams in local catchments. We trained the model with a random 75% subset of Chessie BIBI data ( n = 1449), and used the remaining 25% of Chessie BIBI data ( n = 484) as test data. The model performed well, correctly predicting 72% of samples in training data and 73% of samples in test data, but model accuracy varied among bioregions. We performed uncertainty analyses by adding bands of either ±0.05 or ±0.10 BIBI units to the cutoff between poor and fair/good. These uncertainty analyses resulted in 14.5% (±0.05 band) and 24.8% (±0.10 band) of samples in test data being classified as in uncertain condition. For 95,877 small stream reaches in the Chesapeake Bay watershed, the model predicted 64% in fair/good condition, the ±0.05 uncertainty analyses predicted 57% in fair/good condition, and the ±0.10 uncertainty analysis predicted 50% in fair/good condition. These reported values have different implications for the number of improved stream miles required to meet the goal of improving 10%. Incorporating uncertainty provides an assessment of model strength as well as confidence in predictions. We, therefore, suggest increased reporting of uncertainty in studies that spatially predict stream conditions.

Chesapeake Bay watershed

Unconventional oil and gas spills: Risks, mitigation priorities, and state reporting requirements

Rapid growth in unconventional oil and gas (UOG) has produced jobs, revenue, and energy, but also concerns over spills and environmental risks. We assessed spill data from 2005 to 2014 at 31 481 UOG wells in Colorado, New Mexico, North Dakota, and Pennsylvania. We found 2–16% of wells reported a spill each year. Median spill volumes ranged from 0.5 m 3 in Pennsylvania to 4.9 m 3 in New Mexico; the largest spills exceeded 100 m 3 . Seventy-five to 94% of spills occurred within the first three years of well life when wells were drilled, completed, and had their largest production volumes. Across all four states, 50% of spills were related to storage and moving fluids via flowlines. Reporting rates varied by state, affecting spill rates and requiring extensive time and effort getting data into a usable format. Enhanced and standardized regulatory requirements for reporting spills could improve the accuracy and speed of analyses to identify and prevent spill risks and mitigate potential environmental damage. Transparency for data sharing and analysis will be increasingly important as UOG development expands. We designed an interactive spills data visualization tool ( http://snappartnership.net/groups/hydraulic-fracturing/webapp/spills.html ) to illustrate the value of having standardized, public data.

Environmental Science & Technology

Production and disposal of waste materials from gas and oil extraction from the Marcellus Shale Play in Pennsylvania

The increasing world demand for energy has led to an increase in the exploration and extraction of natural gas, condensate, and oil from unconventional organic-rich shale plays. However, little is known about the quantity, transport, and disposal method of wastes produced during the extraction process. We examined the quantity of waste produced by gas extraction activities from the Marcellus Shale play in Pennsylvania for 2011. The main types of wastes included drilling cuttings and fluids from vertical and horizontal drilling and fluids generated from hydraulic fracturing [i.e., flowback and brine (formation) water]. Most reported drill cuttings (98.4%) were disposed of in landfills, and there was a high amount of interstate (49.2%) and interbasin (36.7%) transport. Drilling fluids were largely reused (70.7%), with little interstate (8.5%) and interbasin (5.8%) transport. Reported flowback water was mostly reused (89.8%) or disposed of in brine or industrial waste treatment plants (8.0%) and largely remained within Pennsylvania (interstate transport was 3.1%) with little interbasin transport (2.9%). Brine water was most often reused (55.7%), followed by disposal in injection wells (26.6%), and then disposed of in brine or industrial waste treatment plants (13.8%). Of the major types of fluid waste, brine water was most often transported to other states (28.2%) and to other basins (9.8%). In 2011, 71.5% of the reported brine water, drilling fluids, and flowback was recycled: 73.1% in the first half and 69.7% in the second half of 2011. Disposal of waste to municipal sewage treatment plants decreased nearly 100% from the first half to second half of 2011. When standardized against the total amount of gas produced, all reported wastes, except flowback sands, were less in the second half than the first half of 2011. Disposal of wastes into injection disposal wells increased 129.2% from the first half to the second half of 2011; other disposal methods decreased. Some issues with data were uncovered during the analytical process (e.g., correct geospatial location of disposal sites and the proper reporting of end use of waste) that obfuscated the analyses; correcting these issues will help future analyses.

Pennsylvania

Unconventional oil and gas spills: Materials, volumes, and risks to surface waters in four states of the U.S.

Extraction of oil and gas from unconventional sources, such as shale, has dramatically increased over the past ten years, raising the potential for spills or releases of chemicals, waste materials, and oil and gas. We analyzed spill data associated with unconventional wells from Colorado, New Mexico, North Dakota and Pennsylvania from 2005 to 2014, where we defined unconventional wells as horizontally drilled into an unconventional formation. We identified materials spilled by state and for each material we summarized frequency, volumes and spill rates. We evaluated the environmental risk of spills by calculating distance to the nearest stream and compared these distances to existing setback regulations. Finally, we summarized relative importance to drinking water in watersheds where spills occurred. Across all four states, we identified 21,300 unconventional wells and 6622 reported spills. The number of horizontal well bores increased sharply beginning in the late 2000s; spill rates also increased for all states except PA where the rate initially increased, reached a maximum in 2009 and then decreased. Wastewater, crude oil, drilling waste, and hydraulic fracturing fluid were the materials most often spilled; spilled volumes of these materials largely ranged from 100 to 10,000 L. Across all states, the average distance of spills to a stream was highest in New Mexico (1379 m), followed by Colorado (747 m), North Dakota (598 m) and then Pennsylvania (268 m), and 7.0, 13.3, and 20.4% of spills occurred within existing surface water setback regulations of 30.5, 61.0, and 91.4 m, respectively. Pennsylvania spills occurred in watersheds with a higher relative importance to drinking water than the other three states. Results from this study can inform risk assessments by providing improved input parameters on volume and rates of materials spilled, and guide regulations and the management policy of spills.

Colorado, New Mexico, North Dakota, Pennsylvania