USGS ScienceSearch

SEARCH · USGS Science

Results for “Scientific Reports”

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.

751 records · Page 39Linked to original sources

Black rails in Baja California

Summarizing records from 1905 to 2025, we document the presence and distribution of the California black rail ( Laterallus jamaicensis coturniculus ) in northwestern México, particularly Baja California and the Colorado River Delta. This subspecies has a restricted distribution in the western U.S. and northwestern México. We found few records over the past century: sparse detections in the early 1900s and an increase in reported records after 2000. Over 120 years, we compiled 35 documented records of California black rail occurrence in Baja California with a minimum reported count of 56 rails. Our call-broadcast survey methods reconfirmed previously occupied sites and documented the species at a previously unreported site. Surveys in the Baja California portion of the Colorado River Delta yielded 20 detections across nine survey events between 2000 and 2022. In April 2003, standardized surveys detected 12 California black rails: 7 in Bahía de San Quintín, 3 in Arroyo San Telmo, and 2 in Arroyo El Rosario. Our recent survey work (2022–2025) recorded maximum monthly counts of 31 detections at Bahía de San Quintín and 15 detections at Estero de Punta Banda, the latter representing a previously undocumented location for the subspecies. These surveys also yielded the first photographic and acoustic documentation of California black rail in México. The contribution of spatial and temporal data provided by citizen science platforms increased our understanding of the subspecies’ distribution. Enforceable protection of this subspecies’ habitat is needed to safeguard its persistence.

Baja California

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California

Spatio-temporal modeling for assessing geoenergy resources: A workflow applied to gas in place variation in coal beds

The ability to estimate spatio-temporal changes in hydrocarbon reservoir properties and energy resources within pore volumes is essential for optimizing production, reservoir management, geologic energy storage, and safety in underground mining operations. In coal seams, predicting remaining methane gas-in-place (GIP) is critical for quantifying producible gas and improving mine safety and productivity through effective ventilation planning. Although such changes are commonly evaluated using physics-based numerical simulation models, these approaches often require extensive data, calibration effort, and time. This study presents a spatio-temporal geostatistical modeling approach that bridges the gap between purely spatial models and full numerical simulations. The method is applied to a case study of coal seam degasification in the Mary Lee coal group, Black Warrior Basin, Alabama, USA, to estimate GIP evolution over time within a selected mining district. The analysis uses published data from prior natural gas production history-matching of degasification using vertical wells. Empirical spatial and temporal statistics were calculated for reservoir pressure and water saturation, and spatio-temporal variogram models were fitted to experimental variograms. These models provided the structural basis for spatio-temporal kriging, integrated with spatial estimates of time-invariant parameters (porosity, density, and thickness) to estimate GIP. This approach enabled estimation of GIP changes over time, including periods without data. Boxplots of GIP estimates indicated systematic depletion and decreasing spatial variability, reflecting the impacts of degasification. Comparison with cumulative gas production from empirical well records showed approximately 85% agreement based on a relative similarity metric. Spatio-temporal GIP estimates were also used to estimate methane emissions to longwall ventilation systems and compared with reported emissions from the U.S. EPA Greenhouse Gas Reporting Program, showing similar distributions (≈80%) given data limitations. Overall, this integrated modeling approach provides time-dependent GIP estimates with broader implications for resource assessment applications.

Alabama

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California

Estimating agricultural irrigation water consumption for the High Plains aquifer region with integrated energy- and water-balance evapotranspiration modeling approaches

Estimation of irrigation water use provides essential information for the management and conservation of agricultural water resources. Conventionally, water use data are created based on reports and surveys from water users, whereas manual records may not be complete due to lacking flow meters, measurement gaps, inconsistent methods across regions, and time- and cost-consuming data processing. Alternatively, spatially explicit estimation of irrigation water use can be conducted efficiently using remote sensing evapotranspiration (ET) modeling approaches. In this study, we created a gridded blue water evapotranspiration (BWET) dataset to estimate historical irrigation water consumption (1986 – 2020) in the croplands across the United States High Plains aquifer region. The BWET data were generated by integrating an energy-balance ET model [Operational Simplified Surface Energy Balance model (SSEBop)] and a water-balance ET model [Vegetation ET model (VegET)]. BWET in croplands indicates crop consumptive use of irrigation water extracted from surface water and groundwater resources. The BWET estimates were compared with reported irrigation water use data for all counties within the aquifer region. The results revealed high agreement between growing season (May – September) BWET and annual water withdrawal at county level. Specifically, correlation coefficients of volumetric BWET and water withdrawal were 0.90 and 0.96, respectively, for the entire aquifer region and western Kansas. The timeseries of BWET and water withdrawal showed similar temporal trends and high covariations. The BWET estimates were systematically lower than the water withdrawal measurements, which was primarily attributed to blue water losses in the irrigation system. The irrigation efficiency, calculated as the ratio of BWET to water withdrawal depth, was 0.57 and 0.74 for the entire aquifer region and western Kansas, respectively. This study demonstrates the capability of using satellite-based ET models (e.g., SSEBop and VegET) to efficiently estimate crop water consumption and evaluate irrigation efficiency at landscape, county, and regional scales.

Colorado, Kansas, Nebraska, New Mexico, Oklahoma,

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan

Measurement of in situ-produced cosmogenic nuclides in fine-grained quartz from shale

In situ -produced 10 Be in quartz is widely used to constrain exposure ages and denudation rates, traditionally measured in sand-sized grains. Here we report a new method for isolating fine-grained quartz from shale and demonstrate its reliability for grain sizes down to single microns. Sequential dissolution tests and analyses of grain size separates show that meteoric 10 Be is eliminated and that recoil losses during cosmogenic nuclide production are balanced by implantation from other mineral grains. High-temperature combustion leads to diffusion of meteoric 10 Be into fine-grained quartz and must be used with caution.

Nuclear Instruments and Methods in Physics Researc

ShakeAlert®—Communication, education, outreach and technical engagement strategic vision

Executive Summary In 2006, the U.S. Geological Survey (USGS) began directly supporting ShakeAlert ® research and in 2012 the ShakeAlert demonstration system began testing ( Given and others, 2018 ). The ShakeAlert earthquake early warning (EEW) system is a partnership between the U.S. Geological Survey (USGS) and the three West Coast States (Washington, Oregon, and California) served by the ShakeAlert System, which is part of the larger Advanced National Seismic System (ANSS). With more than 143 million people exposed to potentially damaging shaking in the United States ( Jaiswal and others, 2015 ), earthquakes are a national hazard. Most of our Nation’s earthquake risk is concentrated in the highly populated areas on the active plate tectonic boundaries on the West Coast of the conterminous United States. ShakeAlert is the first public alert system in the United States to provide rapid mass notification of earthquake detection, potentially offering seconds of warning before strong shaking arrives. A few seconds may not seem like much time, but the information in ShakeAlert Messages can be used to trigger automated actions that can prevent injury or death, reduce immediate damage, and speed recovery from earthquakes. The information product issued by the ShakeAlert system is called a ShakeAlert Message and is one of the information products and tools of the ANSS. The ShakeAlert System includes the USGS component, plus the pathways by which ShakeAlert-powered products and (or) services are delivered to end users. Alerts can be delivered to cell phones or be used to trigger automated systems to protect equipment, facilities, and infrastructure, such as slowing or stopping a train. ShakeAlert-powered automated actions can include fire house doors that can be opened to prevent jamming, heavy equipment (for example, trains, elevators, and cranes) that can be automatically stopped or parked in safe positions, and pipeline valves that can be closed to prevent surges and spills. A few seconds of warning also may be sufficient for people to take protective actions, such as drop, cover, and hold on or modified protective actions for a broad range of populations. Advance training may increase the benefit of a speedy response to an alert. Outreach and education about EEW may raise awareness of the overall earthquake threat and how people can best react when they receive an alert or feel shaking. ShakeAlert communication, education, outreach, and technical engagement (CEO&TE) efforts are highly collaborative and essential for the success of the ShakeAlert System. This strategic vision informs how the vast ShakeAlert CEO&TE Community operates and works together. The CEO&TE Community delineates a strategic framework that is intended to set the path for a long-term, sustainable approach to CEO&TE through three focus areas and five priorities. Focus Areas Technical engagement.— The objective of this focus area is to expand and broaden ShakeAlert technical engagement and implementation. External engagement.— This focus area targets engagement in long-term communication, education, and outreach planning, implementation, and evaluation. Internal engagement.— The purpose of this focus area is to build and sustain a robust infrastructure to optimize collaboration, information sharing, and project planning among the community of stakeholders charged with implementing ShakeAlert CEO&TE. Priorities The enumeration of the five priorities listed below does not suggest priority ranking. Public safety, preparedness, and resilience; Technical implementation and engagement; Consistent messaging and communication; Integration with other Federal and State earthquake hazards products; and Educational resources development and dissemination. This strategic vision is a tangible outcome of collaboration among many stakeholders beginning in July 2016. Since then, the work of the ShakeAlert CEO&TE Community has grown into an international effort. The USGS has developed, tested, and implemented a broad spectrum of communication, education, and outreach tools and resources—all of which recognize that seconds matter when it comes to safety and mitigating harm from earthquake hazards. The CEO&TE social science research effort has provided invaluable insights into the ShakeAlert System’s human interface. USGS-licensed technical partners develop, test, and implement real-world applications using ShakeAlert Messages. The success of ShakeAlert CEO&TE efforts is predicated on robust collaboration across numerous agencies, organizations, and groups. As such, this strategic vision outlines a “partnership model” that delineates roles and responsibilities to ensure alignment with focus areas and priorities. The partnership model includes the CEO&TE lead agency (USGS); its principal partners (State agencies and university partners); its implementation partners (for example, technical partners who build systems to deliver ShakeAlert-powered products and (or) services [focus area one]), earthquake education partners who work to increase public preparedness for seismic events (focus area two); and other organizations that work together to enhance the adoption and effectiveness of the ShakeAlert System. These partners collaborate and convene through a variety of working groups and forums, which are also described in this strategic vision and align with focus area three (internal engagement). The CEO&TE Community collaboratively developed its operating principles and a consensus-based, decision-making strategic framework to guide its collective work. Performance metrics are used to continually measure success. Ultimately, the USGS and ShakeAlert CEO&TE Community are advancing the ShakeAlert System that as of the publication of this strategic vision to “provide earthquake early warning for all” serves more than 50 million people.

Circular

U.S. Geological Survey research and assessments supporting carbon dioxide removal

Both carbon capture and storage (CCS) and carbon dioxide removal (CDR) are methods to limit future global temperature rise and ocean acidification. CCS sequesters (stores) carbon dioxide (CO2) captured from industrial sources thereby preventing the CO2 from reaching the atmosphere. The captured CO2 is injected into underground geologic reservoirs or is converted into stable industrial products. CDR enhances natural processes or uses technical approaches to remove excess CO2 from the atmosphere. Carbon dioxide is sequestered in natural sinks such as vegetation and soil, or injected into underground reservoirs or converted into durable products. The U.S. Geological Survey (USGS) conducts research and resource assessments that support biologic and geologic CDR approaches. This report will review some of the ongoing USGS CDR research and assessment efforts.

Conference Paper

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida

Macroseismology

In this chapter I discuss the use of so-called macroseismic data, i.e., reports of damage and other effects of shaking on humans and the built environment, to improve the characterization of earthquakes and the ground motions they produce. Macroseismic data are critical not only to investigate earthquakes that occurred before the start of the instrumental era in seismology, but are also, by virtue of the spatially rich nature of the data, of tremendous potential importance to investigate early instrumental and even recent earthquakes. Fueled by both the advent of on-line systems that now produce unprecedented volumes of macroseismic intensity data and by recognition of the importance of understanding key historical earthquakes, there has been a growing appreciation for the value of macroseismic data. It is, however, critical to understand both the unique challenges and the unique opportunities associated with these data. I review the evolution of thought and practices regarding analysis of macroseismic data, and the use of such data to improve the characterization of historic and early instrumental earthquakes.

Book chapter

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas

3-D mapping of the conterminous U.S. within the USGS National Cooperative Geologic Mapping Program: Progress and future prospects

The U.S. Geological Survey (USGS) National Cooperative Mapping Program (NCGMP) is bringing together subsurface and three-dimensional information at multiple scales for the conterminous United States from data produced throughout the USGS and by federal and state partners. Components of this work include data inventory and catalog development, data integration and database development, and construction of subregional and basin-scale three-dimensional (3-D) geologic models. At small scales, the NCGMP is compiling subsurface data to create a limited number of subsurface layers that map the majority of the lower 48 States. At larger scales, 3-D models of stratigraphic units are intended to be the subsurface analog and extension of the surface geologic maps made within the Program. The Mapping Program is developing methods of storing, visualizing, and distributing 3-D models and subsurface data from multiple sources in easily shareable, queryable, non-proprietary format(s). We report on progress in 3-D mapping of the conterminous U.S. and on future prospects.

conterminous United States

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada

The influence of nodule versus crust morphology on the composition of seamount-hosted ferromanganese minerals

Ferromanganese mineral precipitation in the global oceans is ubiquitous, occurring in the form of both crusts and nodules at a broad range of depths and seafloor terrains. Although ferromanganese crusts and nodules are both composed of ferromanganese minerals, mineralogy and mean element concentrations compiled for regional crust versus nodule occurrences differ. Notably, most published compilations compare nodules from abyssal plain sediments to crusts forming on seamounts, and do not address the question as to whether mineral morphology may affect the composition of crusts and nodules in addition to environmental factors. To address this, we have compiled a dataset of co-located hydrogenetic (seawater sourced) open-ocean ferromanganese crusts and seamount-hosted nodules. The preliminary data presented here support our hypothesis that co-located crusts and nodules exhibit similar compositions, both in terms of their mineralogy and mean element concentrations, as well as overall variance in concentrations for any given element. These data indicate that local processes influence sample-to-sample variation in both crusts and nodules, in contrast to broader ferromanganese variation reported between open-ocean seamounts- hosted crusts and abyssal plain-hosted nodules.

Conference Paper