USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geographical Analysis”

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,279 records · Page 71Linked to original sources

Cubic map algebra functions for spatio-temporal analysis

We propose an extension of map algebra to three dimensions for spatio-temporal data handling. This approach yields a new class of map algebra functions that we call "cube functions." Whereas conventional map algebra functions operate on data layers representing two-dimensional space, cube functions operate on data cubes representing two-dimensional space over a third-dimensional period of time. We describe the prototype implementation of a spatio-temporal data structure and selected cube function versions of conventional local, focal, and zonal map algebra functions. The utility of cube functions is demonstrated through a case study analyzing the spatio-temporal variability of remotely sensed, southeastern U.S. vegetation character over various land covers and during different El Nin??o/Southern Oscillation (ENSO) phases. Like conventional map algebra, the application of cube functions may demand significant data preprocessing when integrating diverse data sets, and are subject to limitations related to data storage and algorithm performance. Solutions to these issues include extending data compression and computing strategies for calculations on very large data volumes to spatio-temporal data handling.

Cartography and Geographic Information Science↗

The Badain Jaran Desert: Remote sensing investigations.

Approximately half the Badain Jaran Desert in the north-western Alashan Plain of northern China is a sand sea. The remainder is gravel or bedrock. The north-western border of the desert is a playa. The desert has been imaged by both Landsat and the Shuttle Imaging Radar (SIR-A). Landsat analysis indicates there are two dune patterns in the sand sea, north-east oriented crescentic dunes along the northern and western borders, and complex star dunes in the central and eastern desert. Although the orientation and morphology of the dunes are easily visible on Landsat, they cannot be determined with the radar image obtained from the aspect angle used during the SIR-A mission. An abrupt change in wavelength of the dune pattern near the Badain Jaran Playa is mappable on Landsat, but not seen on the radar image. The playa appears to be considerably larger on radar than on the Landsat, and we may be seeing subsurface penetration of dry surficial sands with the radar. Archaeological evidence suggests the playa was the location of prehistoric and historic human activity. SIR-A data indicate the playa was formerly a considerably larger inland lake.

Badain Jaran Desert↗

Strategic conservation planning for the Eastern North Carolina/Southeastern Virginia Strategic Habitat Conservation Team

The Eastern North Carolina/Southeastern Virginia Strategic Habitat Conservation Team (ENCSEVA) is a partnership among local federal agencies and programs with a mission to apply Strategic Habitat Conservation to accomplish priority landscape-level conservation within its geographic region. ENCSEVA seeks to further landscape-scale conservation through collaboration with local partners. To accomplish this mission, ENCSEVA is developing a comprehensive Strategic Habitat Conservation Plan (Plan) to provide guidance for its members, partners, and collaborators by establishing mutual conservation goals, objectives, strategies, and metrics to gauge the success of conservation efforts. Identifying common goals allows the ENCSEVA team to develop strategies that leverage joint resources and are more likely to achieve desired impacts across the landscape. The Plan will also provide an approach for ENCSEVA to meet applied research needs (identify knowledge gaps), foster adaptive management principles, identify conservation priorities, prioritize threats (including potential impacts of climate change), and identify the required capacity to implement strategies to create more resilient landscapes. ENCSEVA seeks to support the overarching goals of the South Atlantic Landscape Conservation Cooperative (SALCC) and to provide scientific and technical support for conservation at landscape scales as well as inform the management of natural resources in response to shifts in climate, habitat fragmentation and loss, and other landscape-level challenges (South Atlantic LCC 2012). The ENCSEVA ecoregion encompasses the northern third of the SALCC geography and offers a unique opportunity to apply landscape conservation at multiple scales through the guidance of local conservation and natural resource management efforts and by reporting metrics that reflect the effectiveness of those efforts (Figure 1). The Environmental Decision Analysis Team, housed within the North Carolina Cooperative Fish and Wildlife Research Unit at North Carolina State University, is assisting the ENCSEVA team in developing a scientifically sound basis for the Plan though the elicitation of expert knowledge and the organization of that knowledge using the Open Standards for the Practice of Conservation. The Open Standards for the Practice of Conservation is a framework that is well suited to incorporating decision-making tools such as Structured Decision Making and provides a multi-step process to conceptually organize conservation projects in a manner that enhances the rigor and transparency of expert and knowledge-based plans. It helps define explicit pathways from 2 planned conservation activities and ultimate impact, as well as indicators to measure success (Stem et al. 2005). Specifically, the framework identifies conservation targets, key ecological attributes, threats, and associated indicators to monitor responses given the implementation of a conservation action (Conservation Measures Partnership 2007). This report serves to provide a scientific foundation for the Plan by summarizing the expert opinion of wildlife biologists, ecologists, hydrologists, researchers, natural resource managers, and conservation practitioners regarding five environments (wetlands, riverine systems, estuaries, uplands, and barrier islands) within the ENCSEVA geography. Specifically, this report describes (1) the approach to elicit expert knowledge meant to support the strategic plan, (2) how this knowledge can inform collaborative conservation planning, and (3) a summary of opportunities available for the ENCSEVA team to address threats and impacts associated with climate change within the ecoregion.

Technical Bulletin↗

Influence of habitat and intrinsic characteristics on survival of neonatal pronghorn

Increased understanding of the influence of habitat (e.g., composition, patch size) and intrinsic (e.g., age, birth mass) factors on survival of neonatal pronghorn ( Antilocapra americana ) is a prerequisite to successful management programs, particularly as they relate to population dynamics and the role of population models in adaptive species management. Nevertheless, few studies have presented empirical data quantifying the influence of habitat variables on survival of neonatal pronghorn. During 2002–2005, we captured and radiocollared 116 neonates across two sites in western South Dakota. We documented 31 deaths during our study, of which coyote ( Canis latrans ) predation ( n = 15) was the leading cause of mortality. We used known fate analysis in Program MARK to investigate the influence of intrinsic and habitat variables on neonatal survival. We generated a priori models that we grouped into habitat and intrinsic effects. The highest-ranking model indicated that neonate mortality was best explained by site, percent grassland, and open water habitat; 90-day survival (0.80; 90% CI = 0.71–0.88) declined 23% when grassland and water increased from 80.1 to 92.3% and 0.36 to 0.40%, respectively, across 50% natal home ranges. Further, our results indicated that grassland patch size and shrub density were important predictors of neonate survival; neonate survival declined 17% when shrub density declined from 5.0 to 2.5 patches per 100 ha. Excluding the site covariates, intrinsic factors (i.e., sex, age, birth mass, year, parturition date) were not important predictors of survival of neonatal pronghorns. Further, neonatal survival may depend on available land cover and interspersion of habitats. We have demonstrated that maintaining minimum and maximum thresholds for habitat factors (e.g., percentages of grassland and open water patches, density of shrub patches) throughout natal home ranges will in turn, ensure relatively high (>0.50) neonatal survival rates, especially as they relate to coyote predation. Thus, landscape level variables (particularly percentages of open water, grassland habitats, and shrub density) should be incorporated into the development or implementation of pronghorn management plans across sagebrush steppe communities of the western Dakotas, and potentially elsewhere within the geographic range of pronghorn.

South Dakota↗

The Louisiana Amphibian Monitoring Program from 1997 to 2017: Results, analyses, and lessons learned

To determine trends in either frog distribution or abundance in the State of Louisiana, we reviewed and analyzed frog call data from the Louisiana Amphibian Monitoring Program (LAMP). The data were collected between 1997 and 2017 using North American Amphibian Monitoring Program protocols. Louisiana was divided into three survey regions for administration and analysis: the Florida Parishes, and 2 areas west of the Florida parishes called North and South. Fifty-four routes were surveyed with over 12,792 stops and 1,066 hours of observation. Observers heard 26 species of the 31 species reported to be in Louisiana. Three of the species not heard were natives with ranges that did not overlap with survey routes. The other two species were introduced species, the Rio Grande Chirping Frog ( Eleutherodactylus cystignathoides) and the Cuban Treefrog ( Osteopilus septentrionalis) . Both seem to be limited to urban areas with little to no route coverage. The 15 most commonly occurring species were examined in detail using the percentage of stops at which they observed along a given survey and their call indices. Most species exhibited a multimodal, concave, or convex pattern of abundance over a 15-year period. Among LAMP survey regions, none of the species had synchronous population trends. Only one group of species, winter callers, regularly co-occur. Based on the species lists, the North region could be seen as a subset of the South. However, based on relative abundance, the North was more similar to Florida parishes for both the winter and summer survey runs. Our analyses demonstrate that long-term monitoring (10 years or more) may be necessary to determine population and occupancy trends, and that frog species may have different local demographic patterns across large geographic areas.

Louisiana↗

Estimated 2012 groundwater potentiometric surface and drawdown from predevelopment to 2012 in the Santa Fe Group aquifer system in the Albuquerque metropolitan area, central New Mexico

Historically, the water-supply requirements of the Albuquerque metropolitan area of central New Mexico were met almost exclusively by groundwater withdrawal from the Santa Fe Group aquifer system. In response to water-level declines, the Albuquerque Bernalillo County Water Utility Authority (ABCWUA) began diverting water from the San Juan-Chama Drinking Water Project in December 2008 to reduce the use of groundwater to meet municipal demand. Modifications in the demand for water and the source of the supply of water for the Albuquerque metropolitan area have resulted in a variable response in the potentiometric surface of the production zone (the interval of the aquifer, from within about 200 feet below the water table to 900 feet or more, in which supply wells generally are screened) of the Santa Fe Group aquifer system. Analysis of the magnitude and spatial distribution of water-level change can help improve the understanding of how the groundwater system responds to withdrawals and variations in the management of the water supply and can support water-management agencies’ efforts to minimize future water-level declines and improve sustainability. The U.S. Geological Survey (USGS), in cooperation with the ABCWUA, has developed an estimate of the 2012 potentiometric surface of the production zone of the Santa Fe Group aquifer system in the Albuquerque metropolitan area. This potentiometric surface is the latest in a series of reports depicting the potentiometric surface of the area. This report presents the estimated potentiometric surface during winter (from December to March) of water year 2012 and the estimated changes in potentiometric surface between predevelopment (pre-1961) and water year 2012 for the production zone of the Santa Fe Group aquifer system in the Albuquerque metropolitan area. Hydrographs from selected piezometers are included to provide details of historical water-level changes. In general, water-level measurements used for this report were collected in small-diameter observation wells screened over short intervals near the middle of the production zone and were considered to best represent the potentiometric head in the production zone. The water-level measurements were collected by various local and Federal agencies. The water year 2012 potentiometric surface map was created in a geographic information system, and the change in water-level altitude from predevelopment to water year 2012 was calculated. The 2012 potentiometric surface indicates that the general direction of groundwater flow is from the Rio Grande towards clusters of supply wells in the east, north, and west. Water-level changes from predevelopment to 2012 were variable across the Albuquerque metropolitan area. Estimated drawdown from 2008 was spatially variable across the Albuquerque metropolitan area. Hydrographs from piezometers on the east side of the river indicate an increase in the annual highest water-level measurement from 2008 to 2012. Hydrographs from piezometers in the northwest part of the study area indicate either steady decline of the water-level altitude over the period of record or recently variable trends in which water-level altitudes increased for a number of years but have declined since water year 2012.

New Mexico↗

A spatially referenced regression model (SPARROW) for suspended sediment in streams of the Conterminous U.S.

Suspended sediment has long been recognized as an important contaminant affecting water resources. Besides its direct role in determining water clarity, bridge scour and reservoir storage, sediment serves as a vehicle for the transport of many binding contaminants, including nutrients, trace metals, semi-volatile organic compounds, a nd numerous pesticides (U.S. Environmental Protection Agency, 2000a). Recent efforts to addr ess water-quality concerns through the Total Maximum Daily Load (TMDL) process have iden tified sediment as the single most prevalent cause of impairment in the Nation’s streams a nd rivers (U.S. Environmental Protection Agency, 2000b). Moreover, sediment has been identified as a medium for the tran sport and sequestration of organic carbon, playing a potentia lly important role in understa nding sources and sinks in the global carbon budget (Stallard, 1998). A comprehensive understanding of sediment fate a nd transport is considered essential to the design and implementation of effective plans for sediment management (Osterkamp and others, 1998, U.S. General Accounting Office, 1990). An exte nsive literature addr essing the problem of quantifying sediment transport has produced a nu mber of methods for estimating its flux (see Cohn, 1995, and Robertson and Roerish, 1999, for us eful surveys). The accuracy of these methods is compromised by uncertainty in the concentration measurements and by the highly episodic nature of sediment movement, particul arly when the methods are applied to smaller basins. However, for annual or decadal flux es timates, the methods are generally reliable if calibrated with extended periods of data (Robertson and Roerish, 1999). A substantial literature also supports the Universal Soil Loss Equation (U SLE) (Soil Conservation Service, 1983), an engineering method for estimating sheet and rill erosion, although the empirical credentials of the USLE have recently been questioned (Tri mble and Crosson, 2000). Conversely, relatively little direct evidence is available concerning the fate of sediment. The common practice of quantifying sediment fate with a sediment deliv ery ratio, estimated from a simple empirical relation with upstream basin area, does not artic ulate the relative importance of individual storage sites within a basin (Wolman, 1977). Rates of sediment deposition in reservoirs and flood plains can be determined from empirical measurement s , but only a limited number of sites have been monitored, and net rates of deposition or loss from other potential sinks and sources is largely unknown (Stallard, 1998). In particular, little is known about how much sediment loss from fields ultimately makes its way to stream channels, and how much sediment is subsequently stored in or lost from th e streambed (Meade and Parker, 1985, Trimble and Crosson, 2000). This paper reports on recent progress made to a ddress empirically the question of sediment fate and transport on a national scale. The model pres ented here is based on the SPAtially Referenced Regression On Watershed attr ibutes (SPARROW) methodology, fi rst used to estimate the distribution of nutrients in str eams and rivers of the United Stat es, and subsequently shown to describe land and stream processes affecting the delivery of nutrients (Smith and others, 1997, Alexander and others, 2000, Preston and Brakeb ill, 1999). The model makes use of numerous spatial datasets, available at the national level, to explain long-term sediment water-quality conditions in major streams and rivers throughou t the United States. Sediment sources are identified using sediment erosion rates from the National Resources I nventory (NRI) (Natural Resources Conservation Service, 2000) and apportioned over the landscape according to 30- meter resolution land-use information from th e National Land Cover Data set (NLCD) (U.S. Geological Survey, 2000a). More than 76,000 reservoirs from the National Inventory of Dams (NID) (U.S. Army Corps of Engin eers, 1996) are identified as pot ential sediment sinks. Other, non-anthropogenic sources and sinks are identified using soil in formation from the State Soil Survey Geographic (STATSGO) data base (Schwarz and Alexander, 1995) and spatial coverages representing surficial rock t ype and vegetative cover. The SPA RROW model empirically relates these diverse spatial datasets to estimates of long-term, mean annual sediment flux computed from concentration and flow measurements co llected over the period 1985 -95 from more than 400 monitoring stations maintained by the Na tional Stream Quality Accounting Network (Alexander and others, 1998), the National Wa ter Quality Assessment Program, and U.S. Geological Survey District offices (Turcios and Gray, in press). Th e calibrated model is used to estimate sediment flux for over 60,000 stream segments included in the River Reach File 1 (RF1) stream network (Alexander and others, 1999). SPARROW uses statis tical methods to calibrate a simple, structural model of riverine water quality, one that imposes mass ba lance in accounting for changes in contaminant flux. As applied here, the mass-balance approach facilitates the interpretation of model results in terms of physical processes affecting sediment transport, and makes possible the estimation of various rates of sediment generation and loss associated with stream channels and features of the landscape. The statistical approach provides a basi s for assessing the error of these inferred rates and of the error in extrapolated estimates of sediment flux made for streams in the RF1 network. An important implication of the holistic modeling approach adopted in this analysis is that estimates of sediment production and loss ar e based on, and therefore consistent with, measurements of in-stream flux. Other ancillary information, such as direct measurements of long-term sediment storage and release from rese rvoirs (Steffen, 1996), is incorporated into the analysis by specifying additional equations expl aining these ancillary variables. The imposition of cross-equation constraints affords this info rmation a statistically consistent weight in explaining in-stream sediment flux. Thus, the me thodology described here represents a general framework for synthesizing a wide spectrum of available information relevant to the understanding of sediment fate and transport.

Conterminous United States↗

Hydrogeologic Framework in Three Drainage Basins in the New Jersey Pinelands, 2004-06

The U.S. Geological Survey, in cooperation with the New Jersey Pinelands Commission, began a multi-phase hydrologic investigation in 2004 to characterize the hydrologic system supporting the aquatic and wetland communities of the New Jersey Pinelands area (Pinelands). The Pinelands is an ecologically diverse area in the southern New Jersey Coastal Plain underlain by the Kirkwood-Cohansey aquifer system. The demand for ground water from this aquifer system is increasing as local development increases. To assess the effects of ground-water withdrawals on Pinelands stream and wetland water levels, three drainage basins were selected for detailed hydrologic assessments, including the Albertson Brook, McDonalds Branch and the Morses Mill Stream basins. Study areas were defined surrounding the three drainage basins to provide sub-regional hydrogeologic data for the ground-water flow modeling phase of this study. In the first phase of the hydrologic assessments, a database of hydrogeologic information and a hydrogeologic framework model for each of the three study areas were produced. These framework models, which illustrate typical hydrogeologic variations among different geographic subregions of the Pinelands, are the structural foundation for predictive ground-water flow models to be used in assessing the hydrologic effects of increased ground-water withdrawals. During 2004-05, a hydrogeologic database was compiled using existing and new geophysical and lithologic data including suites of geophysical logs collected at 7 locations during the drilling of 21 wells and one deep boring within the three study areas. In addition, 27 miles of ground-penetrating radar (GPR) surface geophysical data were collected and analyzed to determine the depth and extent of shallow clays in the general vicinity of the streams. On the basis of these data, the Kirkwood-Cohansey aquifer system was divided into 7 layers to construct a hydrogeologic framework model for each study area. These layers are defined by their predominant sediment textures as aquifers and leaky confining layers. The confining layer at the base of the Kirkwood-Cohansey aquifer system, depending on location, is defined as one of two distinct clays of the Kirkwood Formation. The framework models are described using hydrogeologic sections, maps of structure tops of layers, and thickness maps showing variations of sediment textures of the various model layers. The three framework models are similar in structure but unique to their respective study areas. The hydraulic conductivity of the Kirkwood-Cohansey aquifer system in the vicinity of the three study areas was determined from analysis of 16 slug tests and 136 well-performance tests. The mean values for hydraulic conductivity in the three study areas ranged from about 84 feet per day to 130 feet per day. With the exception of the basal confining layers, the variable and discontinuous nature of clay layers within the Kirkwood-Cohansey aquifer system was confirmed by the geophysical and lithologic records. Leaky confining layers and discontinuous clays are generally more common in the upper part of the aquifer system. Although the Kirkwood-Cohansey aquifer system generally has been considered a water-table aquifer in most areas, localized clays in the aquifer layers and the effectiveness of the leaky confining layers may act to impede the flow of ground water in varying amounts depending on the degree of confinement and the location, duration, and magnitude of the hydraulic stresses applied. Considerable variability exists in the different sediment textures. The extent to which this hydrogeologic variability can be characterized is constrained by the extent of the available data. Thus, the hydraulic properties of the modeled layers were estimated on the basis of available horizontal hydraulic conductivity data and the range of sediment textures estimated from geophysical and lithologic data.

Scientific Investigations Report↗

Congruent population structure inferred from dispersal behaviour and intensive genetic surveys of the threatened Florida scrub-jay ( Aphelocoma cœrulescens )

The delimitation of populations, defined as groups of individuals linked by gene flow, is possible by the analysis of genetic markers and also by spatial models based on dispersal probabilities across a landscape. We combined these two complimentary methods to define the spatial pattern of genetic structure among remaining populations of the threatened Florida scrub-jay, a species for which dispersal ability is unusually well-characterized. The range-wide population was intensively censused in the 1990s, and a metapopulation model defined population boundaries based on predicted dispersal-mediated demographic connectivity. We subjected genotypes from more than 1000 individual jays screened at 20 microsatellite loci to two Bayesian clustering methods. We describe a consensus method for identifying common features across many replicated clustering runs. Ten genetically differentiated groups exist across the present-day range of the Florida scrub-jay. These groups are largely consistent with the dispersal-defined metapopulations, which assume very limited dispersal ability. Some genetic groups comprise more than one metapopulation, likely because these genetically similar metapopulations were sundered only recently by habitat alteration. The combined reconstructions of population structure based on genetics and dispersal-mediated demographic connectivity provide a robust depiction of the current genetic and demographic organization of this species, reflecting past and present levels of dispersal among occupied habitat patches. The differentiation of populations into 10 genetic groups adds urgency to management efforts aimed at preserving what remains of genetic variation in this dwindling species, by maintaining viable populations of all genetically differentiated and geographically isolated populations.

Florida↗

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗

Monitoring population-level foraging distribution of a marine migratory species from land: Strengths and weaknesses of the isotopic approach on the Northwest Atlantic loggerhead turtle aggregation

Assessing the linkage between breeding and non-breeding areas has important implications for understanding the fundamental biology of and conserving animal species. This is a challenging task for marine species, and in sea turtles a combination of stable isotope analysis (SIA) and satellite telemetry has been increasingly used. The Northwest Atlantic (NWA) loggerhead ( Caretta caretta ) Regional Management Unit, one of the largest sea turtle populations in the world, provides an excellent opportunity to investigate key biological patterns as well as methodological aspects related to the use of stable isotopes to infer spatial distribution of turtles in foraging areas. We provide the first comprehensive assessment of the annual distribution of NWA adult female loggerheads among foraging areas and investigate the efficacy of various analytical approaches as well as the effect of sample size in these types of studies. A total of 5168 individual females were sampled from seven Management Units (MUs) between 2013-2018. We provide the first estimate of the proportion of females originating from each MU that uses each foraging area and show how this proportion varies over time. We also estimate the relative importance (in terms of number of turtles) of each foraging area to the overall loggerhead breeding aggregation nesting in Florida and in the NWA for each year of the study. The foraging area used by reproductively active females differs considerably across MUs. One of these, the Subtropical NWA, is by far the most important foraging area in terms of both number of individuals and genetic diversity, and therefore this region may be considered as a conservation priority. Through simulations, we show that limited sizes of sample groups (unknowns; training; priors) may result in false geographic differentiation and consequently mislead interpretations. We provide thresholds and methodological recommendations for future studies. This study establishes a fundamental baseline for monitoring the annual contribution of foraging area to a terrestrial-based breeding aggregation of a marine animal in a cost-effective way. This type of monitoring allows for early detection of changes in foraging distributions—a possible effect of climate change on marine ecosystems or of area-specific anthropogenic threats.

Frontiers in Marine Science↗

Physical characteristics and fish assemblage composition at site and mesohabitat scales over a range of streamflows in the Middle Rio Grande, New Mexico, winter 2011-12, summer 2012

In winter 2011&ndash;12 and summer 2012, the U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers, Albuquerque District and the U.S. Fish and Wildlife Service New Mexico Fish and Wildlife Conservation Office in Albuquerque, New Mexico, evaluated the physical characteristics and fish assemblage composition of available mesohabitats over a range of streamflows at 15 sites on the Middle Rio Grande in New Mexico. The fish assemblage of the Middle Rio Grande includes several minnow species adapted to hydrologically variable but seasonably predictable rivers, including the Hybognathus amarus (Rio Grande silvery minnow), a federally listed endangered species. Gaining a better understanding of habitat usage by the Rio Grande silvery minnow was the impetus for studying physical characteristics and fish assemblages in the Middle Rio Grande during different streamflow conditions. Data were collected at all 15 sites during winter 2011&ndash;12 (moderate streamflow), and a subset was collected at the 13 most downstream sites in summer 2012 (low streamflow). Sites were grouped into four river reaches separated by diversion dams listed in downstream order (names of the diversion dams are followed by short names of the sites nearest each dam in parentheses, listed in downstream order): (1) Cochiti (Pe&ntilde;a Blanca), (2) Angostura (Bernalillo, La Orilla, Barelas, Los Padillas), (3) Isleta (Los Lunas I, Los Lunas II, Abeytas, La Joya, Rio Salado), and (4) San Acacia (Lemitar, Arroyo del Tajo, San Pedro, Bosque del Apache I, and Bosque del Apache II). Stream habitat was mapped in the field by using a geographic information system in conjunction with a Global Positioning System. Fish assemblage composition was determined during both streamflow regimes, and fish were collected by seining in each mesohabitat where physical characteristic data (depth, velocity, dominant substrate type and size, and percent embeddedness) and water-quality properties (temperature, dissolved oxygen, specific conductance, and pH; during summer 2012 only) were measured. Nineteen species of fish were collected among the 15 sites and four reaches over both sampling periods; 10 of these 19 species are introduced. Fish-species richness (total number of fish species collected at each site during each sampling event) among sites that were sampled during both sampling periods ranged from 6 at Rio Salado to 12 at La Orilla. Fish were most abundant at the Lemitar site (1,786 individuals) and least abundant at the San Pedro site (275 individuals). The native Cyprinella lutrensis (red shiner) was the most abundant species collected among all of the sites, accounting for about 42 percent of fish collected. Fish-species richness and catch per unit effort (CPUE) were higher (or equivalent) at all sites during summer 2012 compared to winter 2011&ndash;12. The relations between fish assemblage composition (that is, total abundance, which refers to the number of individuals of each species that were collected) and selected environmental variables (physical characteristic data collected at the mesohabitat scale [depth, velocity, and substrate particle size], and mesohabitat types) were explored by using canonical correspondence analysis. Environmental variables explained 8 percent (p=0.48) of the variability in the Middle Rio Grande fish assemblage during winter 2011&ndash;12, and Rio Grande silvery minnow were weakly associated with sand substrates, relatively moderate velocities (qualitative descriptors are derived from synthetic gradients extracted from CCAs), and relatively shallow depths. Environmental variables explained 14 percent (p < 0.01) of the variability in the Middle Rio Grande fish assemblage during summer 2012, when Rio Grande silvery minnow were associated with run mesohabitats, relatively high velocities, sand substrates, and relatively moderate depths. The mean fish-species richness was greater in summer 2012 than in winter 2011&ndash;12 for each mesohabitat type, and the overall fish-species richness across all mesohabitat types was 0.62 during winter 2011&ndash;12, compared to 1.49 during summer 2012. The highest mean CPUE during winter 2011&ndash;12 was in isolated pools (54.3 fish per 100 square meters [m 2 ]), whereas the lowest was in flats (18.9 fish per 100 m 2 ). Ranges in CPUE were higher in summer 2012 relative to winter 2011&ndash;12 in each mesohabitat type sampled. As in winter 2011&ndash;12, the highest mean CPUE during summer 2012 was in isolated pools (233 fish per 100 m 2 ), whereas the lowest was in flats (29.6 fish per 100 m 2 ). Overall mean CPUE per mesohabitat across all mesohabitat types was 29.1 fish per 100 m 2 during winter 2011&ndash;12 compared to 85.3 fish per 100 m 2 during summer 2012. Four species of minnows (red shiner, Rio Grande silvery minnow, Pimephales promelas [fathead minnow], and Platygobio gracilis [flathead chub]) were selected to compare preferred mesohabitat characteristics because all are small-bodied minnows and because more than 200 individuals of each of these species were collected. Red shiner were collected across the largest range of depths in both winter 2011&ndash;12 (0.02&ndash;4.31 feet [ft]) and summer 2012 (0.05&ndash;3.4 ft), as well as the largest range of velocities (0.0&ndash;4.31 feet per second [ft/s]) during winter 2011&ndash;12 among the four minnow species of interest. Rio Grande silvery minnow occurred in the narrowest range of depths (0.30&ndash;2.1 ft) during summer 2012, as well as the narrowest range of velocities in both winter 2011&ndash;12 (0.0&ndash;3.18 ft/s) and summer 2012 (0.02&ndash;1.51 ft/s). Water-quality properties were only collected during summer 2012, when low-streamflow conditions existed and water-quality properties were thought to be potentially most limiting to aquatic life. Area-weighted mean water temperatures tended to be higher at the sites that were sampled in August 2012 (25.57 degrees Celsius [&deg;C]) compared to June 2012 (24.61 &deg;C). The highest area-weighted mean water temperature at a given site (29.03 &deg;C) was measured at the Lemitar site on August 7, 2012, coincident with the lowest measured discharge (4.13 cubic feet per second [ft 3 /s]). Area-weighted mean dissolved oxygen concentrations tended to be lower in August (7.46 milligrams per liter [mg/L]) compared to June (8.33 mg/L). The highest area-weighted mean dissolved oxygen concentration (9.13 mg/L) was measured at the Lemitar site on August 7, 2012, and the lowest area-weighted mean dissolved oxygen concentration (6.23 mg/L) was measured at the Los Padillas site on August 10, 2012. Area-weighted specific conductance in the sites upstream from La Joya did not exceed 400 microsiemens per centimeter (&mu;S/cm) at 25 &deg;C, whereas the area-weighted mean specific conductance at La Joya (837 &mu;s/cm at 25 &deg;C), Rio Salado (857 &mu;s/cm at 25 &deg;C), and Lemitar (1,300 &mu;s/cm at 25 &deg;C) were all well above the average of the area-weighted means for the 10 remaining sites (433 &mu;s/cm at 25 &deg;C). Lower area-weighted mean pH values were measured at the 3 sites in and near Albuquerque (La Orilla, Barelas, and Los Padillas&mdash;7.98, 8.08, and 7.81, respectively) compared to any of the 10 remaining sites, which had an overall mean pH of 8.44.

New Mexico↗

Selected data on characteristics of glacial-deposit and carbonate-rock aquifers, Midwestern Basin and Arches region

In 1988, the Geological Survey (USGS) began study to examine the hydrogeologic framework, ground-water-flow systems, water chemistry, and withdrawal response of aquifers in glacial deposits and carbonate rock in the Midwestern Basins and Arches Region in western Ohio and eastern Indiana. As part of this study, data from pumped-well tests and instantaneous-rechange tests (slug tests) of wells completed in the glacial-deposit and carbonate-rock aquifers were compiled from reports and information on file with State agencies, environmental consulting firms, drilling firms, municipalities, universities, and the USGS. The data, from 73 counties in Ohio and Indiana, were entered into a computerized data base in a spreadsheet format and subsequently into a geographic information system (GIS). Aquifer-characteristics data from this compilation include the results of 105 pumped-well tests and 39 slug tests in wells completed in glacial deposits, 174 pumped-well tests in wells completed in the carbonate-rock aquifer, and 4 slug tests in wells completed in limestones and shales of Ordovician age. Transmissivities from the pumped-well tests in wells completed in glacial till and glacial-deposit aquifers (sands and gravels) range from 1.54 to 69,700 feet squared per day. Storage coefficients or specific yields range from 0.00002 to 0.38 at these wells. Horizontal-hydraulic conductivities from the slug tests in wells completed in glacial-deposit aquifers range from 0.33 to 1,000 feet per day. Transmissivities from the pumped-well tests in wells completed in the carbonate-rock aquifer range from 70 to 52,000 feet squared per day. Storage coefficient or specific yields at these wells range from 0.00001 to 0.05. Horizontal hydraulic conductivities from the slug tests in wells completed in limestones and shales of Ordovician age range from 0.0016 to 12 feet per day. These data are summarized in tables and figures within this report. The collection and compilation of selected aquifer-characteristic data for the glacial-deposit and carbonate-rock aquifers within the Midwestern Basin and Arches Region of Shaver (1985) are an essential part of the Midwestern Basins and Arches Regional Aquifer-Systems Analysis (Midwestern Basins and Arches RASA) project of the U.S. Geological Survey (USGS). Specifically, the data are needed to help describe ground-water flow in the regional aquifer system, which isone of the objectives of the Midwestern Basins and Arches RASA project (Bugliosi, 1990). To meet this objective, the Midwestern Basins and Arches RASA began subprojects in the Ohio and Indiana offices of the USGS to collect and compile available aquifer-characteristics data from aquifer tests of the glacial-deposit and carbonate-rock aquifers. The data were not reanalyzed to verify accuracy because of time constraints and insufficient data in many cases.

Illinois, Indiana, Kentucky, Ohio, West Virginia↗

Estimation of low-flow statistics at ungaged sites on streams in the Lower Hudson River Basin, New York, from data in geographic information systems

A variety of watershed properties available in 2015 from geographic information systems were tested in regression equations to estimate two commonly used statistical indices of the low flow of streams, namely the lowest flows averaged over 7 consecutive days that have a 1 in 10 and a 1 in 2 chance of not being exceeded in any given year (7-day, 10-year and 7-day, 2-year low flows). The equations were based on streamflow measurements in 51 watersheds in the Lower Hudson River Basin of New York during the years 1958–1978, when the number of streamflow measurement sites on unregulated streams was substantially greater than in subsequent years. These low-flow indices are chiefly a function of the area of surficial sand and gravel in the watershed; more precisely, 7-day, 10-year and 7-day, 2-year low flows both increase in proportion to the area of sand and gravel deposited by glacial meltwater, whereas 7-day, 2-year low flows also increase in proportion to the area of postglacial alluvium. Both low-flow statistics are also functions of mean annual runoff (a measure of net water input to the watershed from precipitation) and area of swamps and poorly drained soils in or adjacent to surficial sand and gravel (where groundwater recharge is unlikely and riparian water loss to evapotranspiration is substantial). Small but significant refinements in estimation accuracy resulted from the inclusion of two indices of stream geometry, channel slope and length, in the regression equations. Most of the regression analysis was undertaken with the ordinary least squares method, but four equations were replicated by using weighted least squares to provide a more realistic appraisal of the precision of low-flow estimates. The most accurate estimation equations tested in this study explain nearly 84 and 87 percent of the variation in 7-day, 10-year and 7-day, 2-year low flows, respectively, with standard errors of 0.032 and 0.050 cubic feet per second per square mile. The equations use natural values of streamflow and watershed properties; logarithmic transformations yielded less accurate equations inconsistent with some conceptualized relationships.

New York↗

Bundle adjustment using space based triangulation method for improving the Landsat global ground reference

There is an ever-increasing interest and need for accurate geo-registration of remotely sensed data products to a common global geometric reference. Although the geo-registration has improved significantly in the last decade, the lack of an accurate global ground reference dataset poses serious issues for data providers seeking to make geometrically stackable analysis ready data. The existing Global Land Survey 2000 (GLS2000) dataset derived from Landsat 7 images provide global coverage and can be used as a reference dataset, but its accuracy is much lower than what can be attained using the agile and precise pointing capability of the new spacecrafts. The improved position and pointing knowledge of the new spacecrafts such as Landsat 8 can be used to improve the accuracy of the existing global ground control points using a space based triangulation method. This paper discusses the theoretical basis, formulation, and application of the space based triangulation method at a continental scale to improve the accuracy of the GLS-derived ground control points.Our triangulation method involves adjusting the spacecraft position, velocity, attitude, attitude rate, and ground control point locations, iteratively, by linearizing the non-linear viewing geometry, such that the residual errors in the measured image points are minimized. The complexity of the numerical inversion and processing is dealt with in our approach by processing and eliminating the ground points one at a time. This helps to reduce the size of the normal matrix significantly, thereby making the triangulation of a continent-wide scale block feasible and efficient. One of the unique characteristics of our method is the use of a correlation model linking the attitude corrections between images of the same pass, which promotes consistency in the attitude corrections. We evaluated the performance of our triangulation method over the Australian continent using the Australian Geographic Reference Image (AGRI) dataset as a reference. Both a free adjustment, using only the pointing information of the Landsat 8 spacecraft, and a constrained adjustment, using the AGRI as external control were performed and the results compared. The Australian block’s horizontal accuracy improved from 15.4 m to 3.6 m with the use of AGRI controls, and from 15.4 m to 8.8 m without the use of AGRI controls.

Remote Sensing↗

Digital data sets of depth-duration frequency of precipitation for Oklahoma

These geospatial data sets were produced as part of a regional precipitation frequency analysis for Oklahoma. The data sets consist of surface grids of precipitation depths for seven frequencies (expressed as recurrence intervals of 2-, 5-, 10-, 25-, 50-, 100-, and 500-years) and 12 durations (15-, 30-, and 60-minutes; 1-, 2-, 3-, 6-, 12-, and 24-hours; and 1-, 3-, and 7-days). Eighty-four depth-duration-frequency surfaces were produced from precipitation-station data. Precipitation-station data from which the surfaces were interpolated and contour lines derived from each surface also are included. Contour intervals vary from 0.05 to 0.5 inch. Data were used from precipitation gage stations with at least 10 years of record within Oklahoma and a zone extending about 50 kilometers into bordering states. Three different rain gage networks provided the data (15-minute, 1-hour, and 1-day). Precipitation annual maxima (depths) were determined from the station data for each duration for 110 15-minute, 141 hourly, and 413 daily stations. Statistical methods were used to estimate precipitation depths for each duration-frequency at each station. These station depth-duration-frequency estimates were interpolated to produce continuous grids with grid-cell spacing of 2,000 meters. Contour lines derived from these surfaces (grids) were used to produce the maps in the 'Depth-Duration Frequency of Precipitation for Oklahoma,' by R.L. Tortorelli, Alan Rea, and W.H. Asquith, U.S. Geological Survey Water-Resources Investigations Report 99-4232. The geospatial data sets are presented in digital form for use with geographic information systems. These geospatial data sets may be used to determine an interpolated value of depth-duration-frequency of precipitation for any point in Oklahoma.

Oklahoma↗

Age, distribution, and stratigraphic relationship of rock units in the San Joaquin Basin Province, California: Chapter 5 in Petroleum systems and geologic assessment of oil and gas in the San Joaquin Basin Province, California

The San Joaquin Basin is a major petroleum province that forms the southern half of California’s Great Valley, a 700-km-long, asymmetrical basin that originated between a subduction zone to the west and the Sierra Nevada to the east. Sedimentary fill and tectonic structures of the San Joaquin Basin record the Mesozoic through Cenozoic geologic history of North America’s western margin. More than 25,000 feet (>7,500 meters) of sedimentary rocks overlie the basement surface and provide a nearly continuous record of sedimentation over the past ~100 m.y. Further, depositional geometries and fault structures document the tectonic evolution of the region from forearc setting to strike-slip basin to transpressional margin. Sedimentary architecture in the San Joaquin Basin is complicated because of these tectonic regimes and because of lateral changes in depositional environment and temporal changes in relative sea level. Few formations are widespread across the basin. Consequently, a careful analysis of sedimentary facies is required to unravel the basin’s depositional history on a regional scale. At least three high-quality organic source rocks formed in the San Joaquin Basin during periods of sea level transgression and anoxia. Generated on the basin’s west side, hydrocarbons migrated into nearly every facies type in the basin, from shelf and submarine fan sands to diatomite and shale to nonmarine coarse-grained rocks to schist. In 2003, the U.S. Geological Survey (USGS) completed a geologic assessment of undiscovered oil and gas resources and future additions to reserves in the San Joaquin Valley of California (USGS San Joaquin Basin Province Assessment Team, this volume, chapter 1). Several research aims supported this assessment: identifying and mapping the petroleum systems, modeling the generation, migration, and accumulation of hydrocarbons, and defining the volumes of rock to be analyzed for additional resources. To better understand the three dimensional relationships between hydrocarbon source and reservoir rocks, we compiled a database consisting of more than 13,000 well picks and of one-mile resolution seismic grids. Both the well picks and the seismic grids characterize the depths to the top of key stratigraphic units. This database formed the basis of subsequent numerical modeling efforts, including the construction of a three- dimensional geologic model (Hosford Scheirer, this volume, chapter 7) and simulation of the petroleum systems in space and time (Peters, Magoon, Lampe, and others, this volume, chapter 12). To accomplish this modeling, we synthesized the age, geographic distribution, lithology, and petroleum characteristics of hydrocarbon source and reservoir rocks in the basin. The results of that synthesis are presented in this paper in the form of new stratigraphic correlation columns for the northern, central, and southern San Joaquin Valley (fig. 5.1; note that all figures are at the back of this report, following the References Cited). The stratigraphic relationships and ages published here draw heavily on published and unpublished studies of the San Joaquin Basin. The stratigraphy presented in each of the columns necessarily idealizes the subsurface geology over a relatively large area, instead of representing the specific geology at an individual well, oil and gas field, or outcrop. In this paper we present the background rationale for defining the geographic divisions of the basin (inset map, fig. 5.1), the paleontological time scales used for assigning absolute ages to rock units (figs. 5.2 and 5.3), and the supporting maps illustrating the geographic distribution of each rock type included in the stratigraphic column (figs. 5.4 through 5.64).

California↗

Susceptibility index to surface contamination for the Little Cross Creek watershed, Cumberland County, North Carolina

An index of surface-water contamination potential was constructed for the Little Cross Creek Basin, a 9.7-square-mile, water-supply watershed in Cumberland County, North Carolina. The index was developed because previous water-quality investigations raised concerns regarding inputs of bacteria, suspended sediment, and phosphorus from nonpoint sources in the watershed. A geographic information system was used to build map overlays and to categorize and rate three factors that affect the transport of water and contaminants-land-surface slope, distance to water, and land use/land cover. Each factor was weighted to reflect its potential contribution to surface-water contamination; the factors then were combined to estimate susceptibility values for the entire watershed. The numerical susceptibility values were categorized to indicate lowest to highest potential for surface-water contamination, and a map was produced showing the spatial distribution of these categories within the watershed. The susceptibility index for about 17 percent of the Little Cross Creek watershed is rated in the high or highest category. These areas have high slopes, short distances to the nearest surface water, impervious land cover, and land uses that generate contaminants. About 38 percent of the watershed area is rated as having low or lowest susceptibility to contamination. These areas contain flat terrain, greater distances to water, land cover that promotes infiltration, and land uses that pose little risk for generating contaminants. Approximately 43 percent of the watershed is in the moderate category of susceptibility. Open water, which is not rated, accounts for the remaining area. The susceptibility index provides water-resource managers with a tool that can aid in prioritizing areas within the Little Cross Creek Basin for monitoring, protection, and remediation. Previous suspended sediment, total phosphorus, and fecal coliform data collected in the Little Cross Creek watershed support the results of the susceptibility analysis. Although this susceptibility index is specific to the Little Cross Creek Basin, the methods used to develop the index are transferable to other watersheds.

North Carolina↗