USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

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 937 records · Page 52Linked to original sources

Elk monitoring in Mount Rainier and Olympic national parks: 2008-2011 synthesis report

In 2008, the USGS collaborated with the NPS, the Muckleshoot and Puyallup Indian Tribes, and WDFW to develop a protocol tor monitor changes in abundance, population composition, and spatial distribution of elk on summer ranges in MORA and OLYM and winter ranges in OLYM. We developed double-observer sightability (DO-S) models that adjusted raw counts of elk as a function of factors influencing detection probabilities from the air, e.g. vegetation, elk group size, light, elk activity, and pilot experience. We plan to develop DO-S models for both MORA and OLYM summer ranges, but due to radiotelemetry collar failures in OLYM, we do not yet have enough data to model detection probabilities in OLYM. We analyzed results of the first 4 years of elk monitoring conducted under the new protocol from 2008-2011. Objectives of this first synthetic analysis are to: • update the DO-S model for MORA aerial survey results • examine abundance, composition, and distribution of elk trends in MORA summer ranges • establish a baseline of counts, population composition, and distribution of elk in OLYM summer ranges • examine trends in counts and distribution of elk in OLYM low-elevation winter ranges during early spring • determine environmental factors influencing abundance and composition of elk in selected MORA summer ranges and unadjusted counts of elk on selected OLYM winter ranges • review progress in developing a DO-S model for OLYM elk surveys • examine aerial survey operations and provide suggestions for future surveys. There was no trend in elk numbers in the N. Rainier TCA from 2008-2011; the trend in the S. Rainier TCA was not statistically significant but increased 3.3%/year. Maximum counts increased in the N. Rainier TCA ~6%/year. Maximum counts in the S. Rainier TCA increased at a rate of 17% annually. Due to failed radiocollars, weather, and other problems, we completed surveys in 2 of 5 OLYM summer range TCAs; no trend data are available. In OLYM winter ranges, we surveyed the Hoh TCA during early spring 2008-10 and of the S. Fork Hoh and Queets in 2008 and 2010. No surveys were done for early-spring counts in 2011 and 2012 due to high snowfall and lack of funding. Legacy early-spring surveys in OLYM since 1985 allowed us to assess trends in counts in the early-spring TCAs from 1985-2010. Counts of elk in the early-spring TCAs declined: 63% in the S. Fork Hoh, 18% in the Hoh, and 22% in the Queets Valley. We continue to develop and improve the DO-S model for application to OLYM summer surveys. In the next synthesis report, we will update findings with additional data following the 2015 field season, based on 8 years of survey results; it will be a more complete analysis of elk population trends.

Washington↗

Modeling flow and transport in unsaturated fractured rock: An evaluation of the continuum approach

Because the continuum approach is relatively simple and straightforward to implement, it has been commonly used in modeling flow and transport in unsaturated fractured rock. However, the usefulness of this approach can be questioned in terms of its adequacy for representing fingering flow and transport in unsaturated fractured rock. The continuum approach thus needs to be evaluated carefully by comparing simulation results with field observations directly related to unsaturated flow and transport processes. This paper reports on such an evaluation, based on a combination of model calibration and prediction, using data from an infiltration test carried out in a densely fractured rock within the unsaturated zone of Yucca Mountain, Nevada. Comparisons between experimental and modeling results show that the continuum approach may be able to capture important features of flow and transport processes observed from the test. The modeling results also show that matrix diffusion may have a significant effect on the overall transport behavior in unsaturated fractured rocks, which can be used to estimate effective fracture-matrix interface areas based on tracer transport data. While more theoretical, numerical, and experimental studies are needed to provide a conclusive evaluation, this study suggests that the continuum approach is useful for modeling flow and transport in unsaturated, densely fractured rock. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Contaminant Hydrology↗

Development, calibration, and analysis of a hydrologic and water-quality model of the Delaware Inland Bays watershed

Excessive nutrients and sediment are among the most significant environmental stressors in the Delaware Inland Bays (Rehoboth, Indian River, and Little Assawoman Bays). Sources of nutrients, sediment, and other contaminants within the Inland Bays watershed include point-source discharges from industries and wastewater-treatment plants, runoff and infiltration to ground water from agricultural fields and poultry operations, effluent from on-site wastewater disposal systems, and atmospheric deposition. To determine the most effective restoration methods for the Inland Bays, it is necessary to understand the relative distribution and contribution of each of the possible sources of nutrients, sediment, and other contaminants. A cooperative study involving the Delaware Department of Natural Resources and Environmental Control, the Delaware Geological Survey, and the U.S. Geological Survey was initiated in 2000 to develop a hydrologic and water-quality model of the Delaware Inland Bays watershed that can be used as a water-resources planning and management tool. The model code Hydrological Simulation Program - FORTRAN (HSPF) was used. The 719-square-kilometer watershed was divided into 45 model segments, and the model was calibrated using streamflow and water-quality data for January 1999 through April 2000 from six U.S. Geological Survey stream-gaging stations within the watershed. Calibration for some parameters was accomplished using PEST, a model-independent parameter estimator. Model parameters were adjusted systematically so that the discrepancies between the simulated values and the corresponding observations were minimized. Modeling results indicate that soil and aquifer permeability, ditching, dominant land-use class, and land-use practices affect the amount of runoff, the mechanism or flow path (surface flow, interflow, or base flow), and the loads of sediment and nutrients. In general, the edge-of-stream total suspended solids yields in the Inland Bays watershed are low in comparison to yields reported for the Eastern Shore from the Chesapeake Bay watershed model. The flatness of the terrain and the low annual surface runoff are important factors in determining the amount of detached sediment from the land that is delivered to streams. The highest total suspended solids yields were found in the southern part of the watershed, associated with high total streamflow and a high surface runoff component, and related to soil and aquifer permeability and land use. Nutrient yields from watershed model segments in the southern part of the Inland Bays watershed were the highest of all calibrated segments, due to high runoff and the substantial amount of available organic fertilizer (animal waste), which results in over-application of organic fertilizer to crops. Time series of simulated hourly total nitrogen concentrations and observed instantaneous values indicate a seasonal pattern, with the lowest values occurring during the summer and the highest during the winter months. Total phosphorus and total suspended solids concentrations are somewhat less seasonal. During storm events, total nitrogen concentrations tend to be diluted and total phosphorus concentrations tend to rise sharply. Nitrogen is transported mainly in the aqueous phase and primarily through ground water, whereas phosphorus is strongly associated with sediment, which washes off during precipitation events.

Water-Resources Investigations Report↗

Evaluation of models and data for assessing whooping crane habitat in the central Platte River, Nebraska

The primary objectives of this evaluation were to improve the performance of the Whooping Crane Habitat Suitability model (C4R) used by the U.S. Fish and Wildlife Service (Service) for defining the relationship between river discharge and habitat availability, and to assist the Service in implementing improved model(s) with existing hydraulic files. The C4R habitat model is applied at the scale of individual river cross-sections, but the model outputs are scaledup to larger reaches of the river using a decision support “model” comprised of other data and procedures. Hence, the validity of the habitat model depends at least partially on how its outputs are incorporated into this larger context. For that reason, we also evaluated other procedures including the PHABSIM data files, the FORTRAN computer programs used to implement the model, and other parameters used to simulate the relationship between river flows and the availability of Whooping Crane roosting habitat along more than 100 miles of heterogeneous river channels. An equally important objective of this report was to fully document these related procedures as well as the model and evaluation results so that interested parties could readily understand the technical basis for the Service’s recommendations.

Nebraska↗

Development of simulated groundwater-contributing areas to selected streams, ponds, coastal water bodies, and production wells in the Plymouth-Carver region and Cape Cod, Massachusetts

Introduction The U.S. Geological Survey (USGS), in support of the Massachusetts Estuaries Project (MEP), delineated groundwater-contributing areas to various hydrologic receptors including ponds, streams, and coastal water bodies throughout southeastern Massachusetts, including portions of the Plymouth-Carver aquifer system and all of Cape Cod. These contributing areas were delineated over a 6-year period from 2003 through 2008 by using previously published regional USGS groundwater-flow models for the Plymouth-Carver region (Masterson and others, 2009), the Sagamore (western) and Monomoy (eastern) flow lenses of Cape Cod (Walter and Whealan, 2005), and lower Cape Cod (Masterson, 2004). The original USGS groundwater-contributing areas were subsequently revised in some locations by the MEP to remove modeling artifacts or to make the contributing areas more consistent with site-specific hydrologic conditions without further USGS review. This report describes the process used to create the USGS groundwater-contributing areas and provides these model results in their original format in a single, publicly accessible publication.

Massachusetts↗

Implementation of the Next Generation Attenuation (NGA) ground-motion prediction equations in Fortran and R

This report presents two methods for implementing the earthquake ground-motion prediction equations released in 2008 as part of the Next Generation Attenuation of Ground Motions (NGA-West, or NGA) project coordinated by the Pacific Earthquake Engineering Research Center (PEER). These models were developed for predicting ground-motion parameters for shallow crustal earthquakes in active tectonic regions (such as California). Of the five ground-motion prediction equations (GMPEs) developed during the NGA project, four models are implemented: the GMPEs of Abrahamson and Silva (2008), Boore and Atkinson (2008), Campbell and Bozorgnia (2008), and Chiou and Youngs (2008a); these models are abbreviated as AS08, BA08, CB08, and CY08, respectively. Since site response is widely recognized as an important influence of ground motions, engineering applications typically require that such effects be modeled. The model of Idriss (2008) is not implemented in our programs because it does not explicitly include site response, whereas the other four models include site response and use the same variable to describe the site condition (VS30). We do not intend to discourage the use of the Idriss (2008) model, but we have chosen to implement the other four NGA models in our programs for those users who require ground-motion estimates for various site conditions. We have implemented the NGA models by using two separate programming languages: Fortran and R (R Development Core Team, 2010). Fortran, a compiled programming language, has been used in the scientific community for decades. R is an object-oriented language and environment for statistical computing that is gaining popularity in the statistical and scientific community. Derived from the S language and environment developed at Bell Laboratories, R is an open-source language that is freely available at http://www.r-project.org/ (last accessed 11 January 2011). In R, the functions for computing the NGA equations can be loaded as an add-on user-contributed code, which is referred to as a ?package? in R. The details of the nga package (Kaklamanos and Thompson, 2010) are presented in this report. In addition, differences between the R and Fortran implementations will be discussed later in this report. The NGA models have established a new baseline for seismic hazard assessments, and they have been incorporated into the most recent national seismic hazard maps published by the U.S. Geological Survey (Petersen and others, 2008). However, many of the new models are significantly more complicated than previous GMPEs and, therefore, require a substantial investment of time to implement and validate. We hope that the scientific and engineering communities find our implementations to be useful in research and practice. Our implementations may be considered as an alternate to the Microsoft Excel spreadsheet implementation available on the PEER NGA project Web site (http://peer.berkeley.edu/ngawest/, last accessed 11 January 2011). The implementations in Fortran and R are more appropriate for performing batch calculations than the implementation in Microsoft Excel. Spreadsheets and Fortran code for some of the individual models also are available on the PEER NGA project Web site; our programs implement the four GMPEs simultaneously. Our programs give the same results as the programs on the PEER NGA Web site, but we offer some additional flexibility of input, additional methods of estimating unknown input parameters, and additional options for output. Although these programs have been used by the U.S. Geological Survey (USGS), Tufts University, and others, no warranty, expressed or implied, is made by Tufts or the USGS as to the accuracy or functioning of the programs and related material, nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by Tufts or the USGS in connection therewith.

Open-File Report↗

A regression model for computing index flows describing the median flow for the summer month of lowest flow in Michigan

In 2006, Michigan enacted laws to prevent new large capacity withdrawals from decreasing flows to the extent that they would functionally impair a stream's ability to support characteristic fish populations. The median streamflow for the summer month of lowest flow was specified by state decision makers as the index flow on which likely impacts of withdrawals would be assessed. At sites near long-term streamflow-gaging stations, analysis of streamflow records during July, August, and September was used to determine the index flow. At ungaged sites, an alternate method for computing the index flow was needed. This report documents the development of a method for computing index flows at ungaged stream sites in Michigan. The method is based on a regression model that computes the index water yield, which is the index flow divided by the drainage area. To develop the regression model, index flows were determined on the basis of daily flows measured during July, August, and September at 147 streamflow-gaging stations having 10 or more years of record (considered long-term stations) in Michigan. The corresponding index water yields were statistically related to climatic and basin characteristics upstream from the stations in the regression model. Climatic and basin characteristics selected as explanatory variables in the regression model include two aquifer-transmissivity and hydrologic-soil groups, forest land cover, and normal annual precipitation. Regression model estimates of water yield explain about 70.8 percent of the variability in index water yields indicated by streamflow-gaging station records. Index flows computed on the basis of regression-model estimates of water yield and corresponding drainage areas explain about 94.0 percent of the variability in index flows indicated by streamflow-gaging station records. No regional bias was detected in the regression-based estimates of water yield within seven hydrologic subregions spanning Michigan. Thus, the single regression model developed in this report can be used to produce unbiased estimates of index water yield and flow statewide. In addition, a technique is presented for computing prediction intervals about the index flow estimates.

Michigan↗

Simulated water budgets and ground-water/surface-water interactions in Bushkill and parts of Monocacy Creek watersheds, Northampton County, Pennsylvania: A preliminary study with identification of data needs

This report, prepared in cooperation with the Department of Environmental Protection, Office of Mineral Resources Management, provides a preliminary analysis of water budgets and generalized ground-water/surface-water interactions for Bushkill and parts of Monocacy Creek watersheds in Northampton County, Pa., by use of a ground-water flow model. Bushkill Creek watershed was selected for study because it has areas of rapid growth, ground-water withdrawals from a quarry, and proposed stream-channel modifications, all of which have the potential for altering ground-water budgets and the interaction between ground water and streams. Preliminary 2-dimensional, steady-state simulations of ground-water flow by the use of MODFLOW are presented to show the status of work through September 2005 and help guide ongoing data collection in Bushkill Creek watershed. Simulations were conducted for (1) predevelopment conditions, (2) a water table lowered for quarry operations, and (3) anthropogenic changes in hydraulic conductivity of the streambed and aquifer. Preliminary results indicated under predevelopment conditions, the divide between the Bushkill and Monocacy Creek ground-water basins may not have been coincident with the topographic divide and as much as 14 percent of the ground-water discharge to Bushkill Creek may have originated from recharge in the Monocacy Creek watershed. For simulated predevelopment conditions, Schoeneck Creek and parts of Monocacy Creek were dry, but Bushkill Creek was gaining throughout all reaches. Simulated lowering of the deepest quarry sump to an altitude of 147 feet for quarry operations caused ground-water recharge and streamflow leakage to be diverted to the quarry throughout about 14 square miles and caused reaches of Bushkill and Little Bushkill Creeks to change from gaining to losing streams. Lowering the deepest quarry sump to an altitude of 100 feet caused simulated ground-water discharge to the quarry to increase about 4 cubic feet per second. Raising the deepest sump to an altitude of 200 feet caused the simulated discharge to the quarry to decrease about 14 cubic feet per second.Decreasing the hydraulic conductivity of the streambed of Bushkill Creek in the reach of large losses of flow caused simulated ground-water levels to decline and ground-water discharge to a quarry to decrease from 74 to 45 cubic feet per second. Decreasing the hydraulic conductivity of a hypothesized highly transmissive zone with a plug of relatively impermeable material caused ground-water levels to increase east of the plug and decline west of the plug, and decreased the discharge to a quarry from 74 to 53 cubic feet per second. Preliminary results of the study have significant limitations, which need to be recognized by the user. The results demonstrated the usefulness of ground-water modeling with available data sets, but as more data become available through field studies, a more complete evaluation could be conducted of the preliminary assumptions in the conceptual model, model sensitivity, and effects of boundary conditions. Additional streamflow and ground-water-level measurements would be needed to better quantify recharge and aquifer properties, particularly the anisotropy of carbonate rocks. Measurements of streamflow losses at average, steady-state hydrologic conditions could provide a more accurate estimate of ground-water recharge from this source, which directly affects water budgets and contributing areas simulated by the model.

Pennsylvania↗

Incorporating induced seismicity in the 2014 United States National Seismic Hazard Model: results of the 2014 workshop and sensitivity studies

The U.S. Geological Survey National Seismic Hazard Model for the conterminous United States was updated in 2014 to account for new methods, input models, and data necessary for assessing the seismic ground shaking hazard from natural (tectonic) earthquakes. The U.S. Geological Survey National Seismic Hazard Model project uses probabilistic seismic hazard analysis to quantify the rate of exceedance for earthquake ground shaking (ground motion). For the 2014 National Seismic Hazard Model assessment, the seismic hazard from potentially induced earthquakes was intentionally not considered because we had not determined how to properly treat these earthquakes for the seismic hazard analysis. The phrases “potentially induced” and “induced” are used interchangeably in this report, however it is acknowledged that this classification is based on circumstantial evidence and scientific judgment. For the 2014 National Seismic Hazard Model update, the potentially induced earthquakes were removed from the NSHM’s earthquake catalog, and the documentation states that we would consider alternative models for including induced seismicity in a future version of the National Seismic Hazard Model. As part of the process of incorporating induced seismicity into the seismic hazard model, we evaluate the sensitivity of the seismic hazard from induced seismicity to five parts of the hazard model: (1) the earthquake catalog, (2) earthquake rates, (3) earthquake locations, (4) earthquake Mmax (maximum magnitude), and (5) earthquake ground motions. We describe alternative input models for each of the five parts that represent differences in scientific opinions on induced seismicity characteristics. In this report, however, we do not weight these input models to come up with a preferred final model. Instead, we present a sensitivity study showing uniform seismic hazard maps obtained by applying the alternative input models for induced seismicity. The final model will be released after further consideration of the reliability and scientific acceptability of each alternative input model. Forecasting the seismic hazard from induced earthquakes is fundamentally different from forecasting the seismic hazard for natural, tectonic earthquakes. This is because the spatio-temporal patterns of induced earthquakes are reliant on economic forces and public policy decisions regarding extraction and injection of fluids. As such, the rates of induced earthquakes are inherently variable and nonstationary. Therefore, we only make maps based on an annual rate of exceedance rather than the 50-year rates calculated for previous U.S. Geological Survey hazard maps.

Open-File Report↗

Modeling sandhill crane population dynamics

The impact of sport hunting on the Central Flyway population of sandhill cranes (Grus canadensis) has been a subject of controversy for several years. A recent study (Buller 1979) presented new and important information on sandhill crane population dynamics. The present report is intended to incorporate that and other information into a mathematical model for the purpose of assessing the long-range impact of hunting on the population of sandhill cranes.The model is a simple deterministic system that embodies density-dependent rates of survival and recruitment. The model employs four kinds of data: (1) spring population size of sandhill cranes, estimated from aerial surveys to be between 250,000 and 400,000 birds; (2) age composition in fall, estimated for 1974-76 to be 11.3% young; (3) annual harvest of cranes, estimated from a variety of sources to be about 5 to 7% of the spring population; and (4) age composition of harvested cranes, which was difficult to estimate but suggests that immatures were 2 to 4 times as vulnerable to hunting as adults.Because the true nature of sandhill crane population dynamics remains so poorly understood, it was necessary to try numerous (768 in all) combinations of survival and recruitment functions, and focus on the relatively few (37) that yielded population sizes and age structures comparable to those extant in the real population. Hunting was then applied to those simulated populations. In all combinations, hunting resulted in a lower asymptotic crane population, the decline ranging from 5 to 54%. The median decline was 22%, which suggests that a hunted sandhill crane population might be about three-fourths as large as it would be if left unhunted. Results apply to the aggregate of the three subspecies in the Central Flyway; individual subspecies or populations could be affected to a greater or lesser degree.

Special Scientific Report - Wildlife↗

Glacier modeling in support of field observations of mass balance at South Cascade Glacier, Washington, USA

The long-term USGS measurement and reporting of mass balance at South Cascade Glacier was assisted in balance years 2006 and 2007 by a new mass balance model. The model incorporates a temperature-index melt computation and accumulation is modeled from glacier air temperature and gaged precipitation at a remote site. Mass balance modeling was used with glaciological measurements to estimate dates and magnitudes of critical mass balance phenomena. In support of the modeling, a detailed analysis was made of the "glacier cooling effect" that reduces summer air temperature near the ice surface as compared to that predicted on the basis of a spatially uniform temperature lapse rate. The analysis was based on several years of data from measurements of near-surface air temperature on the glacier. The 2006 and 2007 winter balances of South Cascade Glacier, computed with this new, model-augmented methodology, were 2.61 and 3.41 mWE, respectively. The 2006 and 2007 summer balances were -4.20 and -3.63 mWE, respectively, and the 2006 and 2007 net balances were -1.59 and -0.22 mWE. PDF version of a presentation on the mass balance of South Cascade Glacier in Washington state. Presented at the American Geophysical Union Fall Meeting 2010.

Washington↗

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗

A Bayesian network to predict vulnerability to sea-level rise: data report

During the 21st century, sea-level rise is projected to have a wide range of effects on coastal environments, development, and infrastructure. Consequently, there has been an increased focus on developing modeling or other analytical approaches to evaluate potential impacts to inform coastal management. This report provides the data that were used to develop and evaluate the performance of a Bayesian network designed to predict long-term shoreline change due to sea-level rise. The data include local rates of relative sea-level rise, wave height, tide range, geomorphic classification, coastal slope, and shoreline-change rate compiled as part of the U.S. Geological Survey Coastal Vulnerability Index for the U.S. Atlantic coast. In this project, the Bayesian network is used to define relationships among driving forces, geologic constraints, and coastal responses. Using this information, the Bayesian network is used to make probabilistic predictions of shoreline change in response to different future sea-level-rise scenarios.

Data Series↗

Verification of water-quality model to simulate effects of discharging treated wastewater during ice-cover conditions to the Red River of the North at Fargo, North Dakota, and Moorhead, Minnesota

The Red River at Fargo Water-Quality (RRatFGO QW) Model, which used the Enhanced Stream Water Quality Model (QUAL2E) computer program, was calibrated and verified for icefree conditions. The purpose of this study was to verify the model for ice-cover conditions using the same Red River of the North study reach that was used for ice-free conditions. The study reach begins about 0.1 mile downstream of the 12th Avenue North bridge in Fargo, North Dakota, and extends 30.8 miles downstream to a site 0.8 mile upstream of the confluence of the Buffalo River and the Red River of the North. The study reach receives treated wastewater outflow from municipal wastewater-treatment plants at Fargo, North Dakota, and Moorhead, Minnesota, and inflow from the Sheyenne River. For simulations conducted for ice-cover conditions, the RRatFGO QW Model will be referred to as the Red River at Fargo Ice-Cover Water-Quality (RRatFGOIC QW) Model. Streamflow measurements were made at 10 sites during February 21-24, 1995, and water-quality samples were collected and field properties were measured at 12 sites during February 23-24, 1995. Properties and constituents analyzed for include specific conductance, water temperature, dissolved oxygen, 5-day carbonaceous biochemical oxygen demand, total nitrite (reported as nitrogen), total nitrite plus nitrate (reported as nitrogen), total ammonia (reported as nitrogen), total organic nitrogen (reported as nitrogen), total phosphorus (reported as phosphorus), chlorophyll a, and algal biomass. The RRatFGOIC QW Model simulated streamflow, specific conductance, total organic nitrogen, total ammonia, total nitrite, total nitrite plus nitrate, 5-day carbonaceous biochemical oxygen demand, and dissolved oxygen. The model was considered verified for ice-cover conditions for all of the values or concentrations simulated except for the total organic nitrogen concentrations. Based on the results of this study, the QUAL2E Model computer program that was calibrated for ice-free conditions is capable of simulating water quality for both ice-free and ice-cover conditions.

Water-Resources Investigations Report↗

Global circulation as the main source of cloud activity on Titan

Clouds on Titan result from the condensation of methane and ethane and, as on other planets, are primarily structured by circulation of the atmosphere. At present, cloud activity mainly occurs in the southern (summer) hemisphere, arising near the pole and at mid-latitudes from cumulus updrafts triggered by surface heating and/or local methane sources, and at the north (winter) pole, resulting from the subsidence and condensation of ethane-rich air into the colder troposphere. General circulation models predict that this distribution should change with the seasons on a 15-year timescale, and that clouds should develop under certain circumstances at temperate latitudes (40??) in the winter hemisphere. The models, however, have hitherto been poorly constrained and their long-term predictions have not yet been observationally verified. Here we report that the global spatial cloud coverage on Titan is in general agreement with the models, confirming that cloud activity is mainly controlled by the global circulation. The non-detection of clouds at latitude 40??N and the persistence of the southern clouds while the southern summer is ending are, however, both contrary to predictions. This suggests that Titans equator-to-pole thermal contrast is overestimated in the models and that its atmosphere responds to the seasonal forcing with a greater inertia than expected. ?? 2009 Macmillan Publishers Limited. All rights reserved.

Nature↗

Hydrogeologic data for the Blaine aquifer and associated units in southwestern Oklahoma and northwestern Texas

This report is a compilation of hydrogeologic data collected for an areal ground-water investigation of the Blaine aquifer and associated units in southwestern Oklahoma and northwestern Texas. The study area includes parts of Greer, Harmon, and Jackson counties in Oklahoma and parts of Childress, Collingsworth, Hall, Hardeman, and Wilbarger counties in Texas. The Blaine aquifer consists of cavernous gypsum and dolomite beds. Water from the Blaine aquifer supports a local agriculture based mainly on irrigated cotton and wheat. The purpose of the study was to determine the availability, quantity, and quality of ground water from the Blaine aquifer and associated units. This report provides a reference for some of the data that was used as input into a computer ground-water flow model that simulates ground-water flow in the Blaine aquifer. The data in this report consists of: (1) Monthly or periodic water-level measurements in 134 wells; (2) daily mean water-level measurements for 11 wells equipped with water-level recorders; (3) daily total precipitation measurements from five precipitation gages; (4) low-flow stream-discharge measurements for 89 stream sites; (5) miscellaneous stream-discharge measurements at seven stream sites; (6) chemical analyses of surface water from 78 stream sites during low-flow periods; (7) chemical analyses of ground water from 41 wells; and (8) chemical analyses of runoff water collected at five sites.

Open-File Report↗

Executive summary

Fossil fuels from the Appalachian basin region have been major contributors to the Nation’s energy needs over much of the last three centuries. Early records indicate that Appalachian coal was first mined in the middle 1700s (Virginia and Pennsylvania) and was used sparingly to fuel colonial settlements and, later, a fledgling industrial-based economy along the eastern seaboard of the United States (de Witt and Milici, 1989). In 2011, central Appalachian basin coal production accounted for approximately 77 percent of all U.S. metallurgical (or coking) coal and 29 percent of total U.S. production (U.S. Energy Information Administration, 2013). Following initial discoveries and commercial use in western New York (1821) and Ohio and West Virginia (mid-1830s), the Appalachian petroleum (oil and gas) industry began in earnest in 1859 with the discovery of oil at the Drake well in northwestern Pennsylvania. Between 1860 and 1989, the Appalachian basin produced more than 2.5 billion barrels of oil (BBO) and more than 30 trillion cubic feet of gas (TCFG) from more than 500,000 wells (de Witt and Milici, 1989). Although both oil and gas continue to be produced in the Appalachian basin, most new wells in the region are drilled in shale reservoirs to produce natural gas. Appalachian coal and petroleum resources are still available in sufficient quantities to contribute significantly to the Nation’s energy needs. For example, the U.S. Energy Information Administration (2010) estimated that there are 6,484 million short tons of recoverable coal reserves in the Appalachian basin. Similarly, about 14.7 billion barrels of oil equivalent (BBOE) (1.2 BBO+81 TCFG [or 13.5 BBOE]) of recoverable Appalachian basin oil and gas remain available of an estimated ultimate endowment of approximately 25.5 billion BBOE (cumulative production + reserves + estimated recoverable undiscovered resources) (this volume, chap. C.1). U.S. Geological Survey (USGS) Professional Paper 1708 is a modern, indepth collection of reports, cross sections, and maps that describe the geology of the Appalachian basin and its fossil fuel resources. Several of the chapters have been published in outside journals or as other USGS publications. Although this volume is not a comprehensive regional treatment of all notable geologic and fossil fuel localities in the Appalachian basin, the selected study areas and topics presented in the chapters cover large segments of the basin and a wide range of stratigraphic intervals. As the title implies, this volume addresses topics that refer to the locations of coal and petroleum accumulations, the stratigraphic and structural framework, and the geochemical characteristics of the coal beds and petroleum in the basin, as well as the results and documentation of recent USGS assessments of coal, oil, and gas resources in the basin. Many of the maps and accompanying data supporting the reports in this volume are available as downloadable geographic information system (GIS) data files (such as selected coal beds, selected oil and gas fields, locations of oil and gas wells, coal production, coal chemistry, total petroleum system (TPS) boundaries, and bedrock geology). Log ASCII Standard (LAS) files for geophysical (gamma ray) wireline well logs also are included. This publication supplements and updates older USGS regional studies of Appalachian basin coal and petroleum resources such as those by Arndt and others (1968) and the numerous contributors to USGS Miscellaneous Map Series I−917 (for example, Harris and others, 1978), respectively. USGS Professional Paper 1708 is intended primarily for geoscientists in academia, industry, and government who are interested in Appalachian basin geology and its coal and petroleum resources. Other users, however, may find the wide variety of topics, papers, and digital images of value for landuse and policy planning issues. Among the anticipated benefits of the report are improvements in (1) resource assessment estimates and methodology, (2) exploration strategies, (3) basin models, and (4) energy use policies.

Appalachian Basin↗

Climate change and tree-line ecosystems in the Sierra Nevada: Habitat suitability modelling to inform high-elevation forest dynamics monitoring

Whitebark pine and foxtail pine serve foundational roles in the subalpine zone of the Sierra Nevada. They provide the dominant structure in tree-line forests and regulate key ecosystem processes and community dynamics. Climate change models suggest that there will be changes in temperature regimes and in the timing and magnitude of precipitation within the current distribution of these species, and these changes may alter the species’ distributional limits. Other stressors include the non-native pathogen white pine blister rust and mountain pine beetle, which have played a role in the decline of whitebark pine throughout much of its range. The National Park Service is monitoring status and trends of these species. This report provides complementary information in the form of habitat suitability models to predict climate change impacts on the future distribution of these species within Sierra Nevada national parks. We used maximum entropy modeling to build habitat suitability models by relating species occurrence to environmental variables. Species occurrence was available from 328 locations for whitebark pine and 244 for foxtail pine across the species’ distributions within the parks. We constructed current climate surfaces for modeling by interpolating data from weather stations. Climate surfaces included mean, minimum, and maximum temperature and total precipitation for January, April, July, and October. We downscaled five general circulation models for the 2050s and the 2090s from ~125 km2 to 1 km2 under both an optimistic and an extreme climate scenario to bracket potential climatic change and its influence on projected suitable habitat. To describe anticipated changes in the distribution of suitable habitat, we compared, for each species, climate scenario, and time period, the current models with future models in terms of proportional change in habitat size, elevation distribution, model center points, and where habitat is predicted to expand or contract. Overall, models indicated that suitable habitats for whitebark and foxtail pine are more likely to shift geographically within the parks by 2100 rather than decline precipitously. This implies park managers might focus conservation efforts on stressors other than climate change, working toward species resilience in the face of threats from introduced disease and elevated native insect damage. More specifically, further understanding of the incidence and severity of white pine blister rust and other stressors in high elevation white pines would help assess vulnerability from threats other than climate change.

California↗