SEARCH · USGS Science

Results for “Final Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76Linked to original sources

USGS2006

Recent volcanic history of Irazu volcano, Costa Rica: alternation and mixing of two magma batches, and pervasive mixing

40 Ar/ 39 Ar dates, field observations, and geochemical data are reported for Irazú volcano, Costa Rica. Volcanism dates back to at least 854 ka, but has been episodic with lava shield construction peaks at ca. 570 ka and 136–0 ka. The recent volcanic record on Irazú volcano comprises lava flows and a variety of Strombolian and phreatomagmatic deposits, with a long-term trend toward more hydrovolcanic deposits. Banded scorias and hybridized rocks reflect ubiquitous magma mixing and commingling. Two distinct magma batches have been identified. One magma type or batch, Haya, includes basalt with higher high field strength (HFS) and rare-earth element contents, suggesting a lower degree melt of a subduction modified mantle source. The second batch, Sapper, has greater enrichment of large ion lithophile elements (LILE) relative to HFS elements and rare-earth elements, suggesting a higher subduction signature. The recent volcanic history at Irazú records two and one half sequences of the following pattern: eruptions of the Haya batch; eruptions of the Sapper batch; and finally, an unusually clear unconformity, indicating a pause in eruptions. In the last two sequences, strongly hybridized magma erupted after the eruption of the Haya batch. The continuing presence of two distinct magma batches requires two active magma chambers. The common occurrence of hybrids is evidence for a small, nearer to the surface chamber for mixing the two batches. Estimated pre-eruptive temperatures based on two-pyroxene geothermometry range from &sim;1000–1176 °C in basalts to 922 °C in hornblende andesites. Crystallization occurred mainly between 4.6 and 3 kb as measured by different geobarometers. Hybridized rocks show intermediate pressures and temperatures. High silica magma occurs in very small volumes as banded scorias but not as lava flows. Although eruptions at Irazú are not often very explosive, the pervasiveness of magma mixing presents the danger of larger, more explosive hybrid eruptions.

Irazu volcano↗

USGS2004

Geologic map of the Reyes Peak quadrangle, Ventura County, California

New 1:24,000-scale geologic mapping in the Cuyama 30' x 60' quadrangle, in support of the USGS Southern California Areal Mapping Project (SCAMP), is contributing to a more complete understanding of the stratigraphy, structure, and tectonic evolution of the complex junction area between the NW-trending Coast Ranges and EW-trending western Transverse Ranges. The 1:24,000-scale geologic map of the Reyes Peak quadrangle, located in the eastern part of the Cuyama map area, is the final of six contiguous 7 ?' quadrangle geologic maps compiled for a more detailed portrayal and reevaluation of geologic structures and rock units shown on previous maps of the region (Carman, 1964; Dibblee, 1972; Vedder and others, 1973). SCAMP digital geologic maps of the five other contiguous quadrangles have recently been published (Minor, 1999; Kellogg, 1999, 2003; Stone and Cossette, 2000; Kellogg and Miggins, 2002). This digital compilation presents a new geologic map database for the Reyes Peak 7?' quadrangle, which is located in southern California about 75 km northwest of Los Angeles. The map database is at 1:24,000-scale resolution.

Open-File Report↗

USGS2010

June and August median streamflows estimated for ungaged streams in southern Maine

Methods for estimating June and August median streamflows were developed for ungaged, unregulated streams in southern Maine. The methods apply to streams with drainage areas ranging in size from 0.4 to 74 square miles, with percentage of basin underlain by a sand and gravel aquifer ranging from 0 to 84 percent, and with distance from the centroid of the basin to a Gulf of Maine line paralleling the coast ranging from 14 to 94 miles. Equations were developed with data from 4 long-term continuous-record streamgage stations and 27 partial-record streamgage stations. Estimates of median streamflows at the continuous-record and partial-record stations are presented. A mathematical technique for estimating standard low-flow statistics, such as June and August median streamflows, at partial-record streamgage stations was applied by relating base-flow measurements at these stations to concurrent daily streamflows at nearby long-term (at least 10 years of record) continuous-record streamgage stations (index stations). Weighted least-squares regression analysis (WLS) was used to relate estimates of June and August median streamflows at streamgage stations to basin characteristics at these same stations to develop equations that can be used to estimate June and August median streamflows on ungaged streams. WLS accounts for different periods of record at the gaging stations. Three basin characteristics-drainage area, percentage of basin underlain by a sand and gravel aquifer, and distance from the centroid of the basin to a Gulf of Maine line paralleling the coast-are used in the final regression equation to estimate June and August median streamflows for ungaged streams. The three-variable equation to estimate June median streamflow has an average standard error of prediction from -35 to 54 percent. The three-variable equation to estimate August median streamflow has an average standard error of prediction from -45 to 83 percent. Simpler one-variable equations that use only drainage area to estimate June and August median streamflows were developed for use when less accuracy is acceptable. These equations have average standard errors of prediction from -46 to 87 percent and from -57 to 133 percent, respectively.

Scientific Investigations Report↗

USGS2017

Multiscale guidance and tools for implementing a landscape approach to resource management in the Bureau of Land Management

The Bureau of Land Management (BLM) is implementing a landscape approach to resource management (hereafter, landscape approach) to more effectively work with partners and understand the effects of management decisions. A landscape approach is a set of concepts and principles used to guide resource management when multiple stakeholders are involved and goals include diverse and sustainable social, environmental, and economic outcomes. Core principles of a landscape approach include seeking meaningful participation of diverse stakeholders, considering diverse resource values in multifunctional landscapes, acknowledging the tradeoffs needed to meet diverse objectives in the context of sustainable resource management, and addressing the complexity of social and ecological processes by embracing interdisciplinarity and considering multiple and broad spatial and temporal perspectives. In chapter 1, we outline the overall goal of this report: to provide a conceptual foundation and framework for implementing a landscape approach to resource management in the BLM, focusing on the role of multiscale natural resource monitoring and assessment information. In chapter 2, we describe a landscape approach to resource management. BLM actions taken to implement a landscape approach include a major effort to compile broad-scale data on natural resource status and condition across much of the west. These broadscale data now provide a regional context for interpreting monitoring data collected at individual sites and informing decisions made for local projects. We also illustrate the utility of using multiscale data to understand potential effects of different resource management decisions, define relevant terms in landscape ecology, and identify spatial scales at which planning and management decisions may be evaluated. In chapter 3, we describe how the BLM Rapid Ecoregional Assessment program and Assessment, Inventory and Monitoring program may be integrated to provide the multiscale monitoring data needed to inform a landscape approach. We propose six core, broad-scale indicators of natural resource status and condition: the amount, spatial distribution, patch size and connectivity of ecosystems and wildlife habitats, and the pattern of existing development across the landscape. Additional supplemental broad-scale indicators may include fire return intervals, distributions of invasive species, and vulnerability of ecosystems to a changing climate. Landscape intactness is an additional derived indicator that is calculated from one or more of the core and supplemental broad-scale indicators. We then outline a process for assessing broad-scale indicators that is consistent with the overall BLM Assessment, Inventory, and Monitoring process, facilitating development of a multiscale natural resource monitoring program. Finally, we describe how broad-scale indicators of natural resource status and condition may guide field monitoring implemented through the BLM Assessment, Inventory and Monitoring program and help address complex management questions. In chapter 4, we consider the specific question of assessing the ecological integrity of rangelands across the western United States. We first define ecological integrity and its relation to land health. We then suggest that a combination of six local-scale indicators collected through field sampling at individual sites and five complementary broad-scale indicators together provide information on the composition, structure, and function of rangelands. The terrestrial monitoring indicators collected at the level of individual field sites are the amount of bare ground, vegetation composition (including invasive plants and plants of management concern), vegetation height, and the proportion of the soil surface in large intercanopy gaps. The broad-scale indicators are vegetation amount, distribution, patch size, connectivity, and productivity, along with the pattern of terrestrial development. Our suggested approach to quantifying ecological integrity focuses specifically on informing management of public lands for multiple resource uses, and illustrates how existing data collected through BLM monitoring and assessment programs may be used together to provide multiscale information on land condition across broad extents. In chapter 5, we develop a method for quantifying landscape intactness and apply this method to the western United States. Our multiscale index of landscape intactness is designed to be defensible, decomposable, and easy to understand. The foundation of the multiscale index of landscape intactness is the surface disturbance footprint of anthropogenic development, including energy and urban development, roads and railroads, cultivated croplands, surface mines and quarries, and energy transmission lines and pipelines. The index represents a gradient of anthropogenic influence as represented by development summarized at two spatial scales of analysis: 2.5 and 20 kilometers. We provide several example applications of the index, illustrating how these data may inform natural resource decisions at the spatial extent of BLM field and district offices, states, ecoregions, and the western United States. We find that 19.2 percent of lands managed by the BLM across the 17 western states of the conterminous United States had the highest landscape intactness. The largest intact areas occur on public lands at high elevations or in the Great Basin. We believe the frameworks, processes, and analyses provided in this report will improve the ability of the BLM to identify and evaluate potential direct and indirect effects of management actions (such as habitat restoration and renewable energy development), and assist the BLM in further implementing a landscape approach to resource management.

Open-File Report↗

USGS2023

Seismic source characterization in central and eastern United States

This report documents earthquake geology studies conducted in the central and eastern United States (CEUS) between 2018 and 2023. The overarching goal of the earthquake geology studies was to identify and characterize fault sources in the New Madrid seismic zone, the eastern Tennessee seismic zone, and the Charleston seismic zone. These studies primarily relied on the interpretation of high-resolution topography, morphotectonic analysis of 10-meter (m) digital elevation models, field reconnaissance, and integration with new and existing subsurface datasets. Results from each study are placed into a seismic hazard framework. The first two chapters focus on the New Madrid seismic zone, which is the most seismically active region in the CEUS, to better characterize the southern extent of the Reelfoot fault in the Obion River valley, and identify and characterize proposed faults along the margins of Crowleys Ridge in the center of the Mississippi River embayment. Crowleys Ridge is a ~320-kilometer (km) long landform that has been proposed to be formed by intermittent faulting or erosion from the Mississippi River. Legacy seismic reflection data documented Eocene, Mesozoic, and older strata offset by faults on the margins of the ridge (Van Arsdale et al., 1999). A morphotectonic analysis of the topography along Crowleys Ridge indicated southward increases in catchment hypsometric integral and slope. Neotectonic mapping on high-resolution lidar data revealed a series of subparallel linear fault and fold scarps along the margins of Crowleys Ridge. These scarps offset previously mapped fluvial and alluvial surfaces that were <56 thousand years ago (ka) (Rittenour et al., 2007). Integration of new seismic reflection and airborne electromagnetic data, legacy seismic data, with the landscape morphotectonic analysis and neotectonic mapping supports the interpretation that Crowleys Ridge is bound by faults that have been active in the late Quaternary but do not experience modern seismicity. By placing these newly characterized faults into the regional framework, the Crowleys Ridge faults can be interpreted as thrust faults associated with stepovers in the overall dextral fault system, similar to the Reelfoot fault. The results of these studies indicate that seismicity has migrated eastward during the late Quaternary, with the Reelfoot fault accommodating much of the modern-day strain. The Reelfoot fault ruptured in the 1811-1812 New Madrid earthquake sequence (Fuller, 1912) and has a record of past large earthquakes (Tuttle et al., 2002; 2019). Whereas the northern extent of the fault intersects the Mississippi River and has well-expressed surface deformation, questions remained regarding the southern limit of seismicity, deformation, and overall fault length. Analysis of <25 ka river terraces along the Obion River valley in western Tennessee revealed subtle fold scarps, with increased magnitudes of deformation on progressively older terraces, indicative of a longer record of fault movement than has been documented by paleoseismic and paleoliquefaction studies alone (e.g., Kelson et al., 1996; Gold et al., 2019; Tuttle et al., 2019). The third chapter integrates surface and subsurface data to assess the landscape record of surface deformation in the eastern Tennessee seismic zone (ETSZ), which is the second most seismically active region in the CEUS. Seismicity in the ETSZ is deep (5-26 km) and focused in the Proterozoic rock below the Paleozoic detachment. Although the region has experienced moderate (≤ Mw 4.8) seismicity in the instrumental record, the lack of a large historical earthquake and limited paleoseismic and paleoliquefaction evidence of prehistoric large ground shaking events have made seismic hazard characterization of the region challenging. Two models have emerged to predict the style and orientation of active surface deformation in the ETSZ. One suggests that surface rupturing earthquakes would exploit pre-existing faults along the regional structural grain (Cox et al., 2022). The second model advocates for primarily strike-slip motion on east-west or north-south faults, following analyses of recent seismicity (Chapman et al., 1997; Dunn and Chapman, 2006; Daniels and Peng, 2022) and the modern stress field (Levandowski et al., 2018). In our analysis, neotectonic mapping of high-resolution lidar data reveals a concentration of lineaments with east-west orientations, with some corresponding to previously mapped east-west faults that crosscut the northeast-southwest regional structural grain from the Paleozoic orogeny. A morphotectonic analysis of catchments and river segments within the same lithology indicates subtle differences that may be indicative of a longer-term tectonic uplift signal. These changes in morphotectonic metrics spatially correspond to newly mapped lineaments and previously mapped east-west trending faults. Within a regional framework, we suggest that diffuse surface deformation associated with deep seismicity is accommodated on a network on east-west faults. However, further work is needed to better understand the potential late Pleistocene fault activity of these lineaments, and we propose that the ETSZ is still best characterized as an area source in seismic hazard models. The final chapter presents preliminary work to create a geographic information system (GIS) database of recent studies analyzing surface and subsurface datasets and interpretations in the Charleston seismic zone, South Carolina. This GIS database will serve as a foundation for future work to analyze new quality level 1 (QL1; <0.5 m resolution) lidar data over the 1886 magnitude (M) 7 Charleston epicentral region to identify possible fault source sources responsible for the 1886 or other surface-rupturing events. This chapter briefly summarizes three new subsurface datasets, including seismic reflection (Pratt et al., 2022; Liberty, 2022) and seismicity data (Chapman et al., 2016), and one new surface dataset (Marple and Hurd, 2020), and their interpretations of potentially active lineaments and faults in the region.

central and eastern United States↗

USGS2002

Crustal structure of the coastal and marine San Francisco Bay region, California

As of the time of this writing, the San Francisco Bay region is home to about 6.8 million people, ranking fifth among population centers in the United States. Most of these people live on the coastal lands along San Francisco Bay, the Sacramento River delta, and the Pacific coast. The region straddles the tectonic boundary between the Pacific and North American Plates and is crossed by several strands of the San Andreas Fault system. These faults, which are stressed by about 4 cm of relative plate motion each year, pose an obvious seismic hazard. We have many ways to study earthquake faults. Where faults break the land surface, we may learn valuable information needed for hazard assessment, such as cumulative offset, slip rate, and earthquake history. However, many of the major faults in the region are partly submerged beneath San Francisco and Monterey Bays. Although this situation poses problems in gathering observational data for hazard assessment, bay-region waterways provide an opportunity to study faultzone structure by using marine subsurface-imaging techniques, which are easier and cheaper than equivalent studies on land. In 1993, the U.S. Geological Survey (USGS) launched a 5-year project aimed at unearthing the basic science of the submerged San Andreas strike-slip fault system in the San Francisco Bay region with its many interacting strands. Primary project goals were structural, such as to discover how the San Andreas and Hayward Faults are connected or related at depth, to learn how the complex of faults in the San Andreas stepover zone on the Golden Gate platform functions, and to locate previously unknown faults. This volume thus contains mostly structural information about the San Francisco Bay region, much of it gathered through exploratory geophysical experiments. The volume is organized “top down,” from studies in the shallowest crust to the base of the crust. The first three chapters are linked through their use of novel geophysical techniques to study earthquake effects, coseismic slip, and shallow stratigraphy. Kayen and others examine crustal structure at very high resolution and demonstrate the use of ground-penetrating-radar tomography to measure the liquefaction potential of coastal sedimentary deposits. McGann and others use microfossils from drill cores along the San Francisco-Oakland Bay Bridge to determine a more detailed late Pleistocene stratigraphy of San Francisco Bay than was previously available. Geist and Zoback use the historical record of a small local tsunami generated by the great 1906 San Francisco earthquake to model the rupture process of that earthquake. The last four chapters are dedicated to studies of fault related structure of the seismogenic crust in and around the San Andreas Fault system in the San Francisco Bay region. Jachens and others compile an aeromagnetic anomaly map from new high-resolution flights across the bay region. Some of these anomalies mark the positions of offshore faults, and others are offset by faults, providing constraints on cumulative slip. Hart and others concisely summarize the marine seismic data recorded in and around San Francisco Bay, map the coverage, and provide archival information for those interested in acquiring data. The last two chapters present the results of the seismic data that have been analyzed. Bruns and others present their analysis of high-quality intermediate-resolution (~5-km penetration) seismic-reflection data gathered over the complex San Andreas-San Gregorio Fault junction. This junction, which is thought to be where the 1906 San Francisco earthquake originated (see Geist and Zoback, this volume), contains an apparent extensional right stepover in the San Andreas Fault. Finally, Parsons and others review and summarize the results of deep-crustal seismic-reflection experiments and local-earthquake tomographic studies, including previously unpublished data, and provide additional support and discussion for already-published studies. In summary, these studies were carried out in an environment where background information on faults in the San Francisco Bay region was sought. Much of the structural information presented here comes from experiments of a style unlikely to be conducted by the USGS in the near future. Together, the chapters in this volume provide a structural framework for a major part of a complex strike-slip fault system.

California↗

USGS2010

Determination of baseline periods of record for selected streamflow-gaging stations in and near Oklahoma for use in modeling applications

Use of historical streamflow data from a least-altered period of record can be used in calibration of various modeling applications that are used to characterize least-altered flow and predict the effects of proposed streamflow alteration. This information can be used to enhance water-resources planning. A baseline period of record was determined for selected streamflow-gaging stations that can be used as a calibration dataset for modeling applications. The baseline period of record was defined as a period that is least-altered by anthropogenic activity and has sufficient streamflow record length to represent extreme climate variability. Streamflow data from 171 stations in and near Oklahoma with a minimum of 10 complete water years of daily streamflow record through water year 2007 and drainage areas that were less than 2,500 square miles were considered for use in the baseline period analysis. The first step to determine the least-altered period of record was to evaluate station information by using previous publications, historical station record notes, and information gathered from oral and written communication with hydrographers familiar with selected stations. The second step was to indentify stations that had substantial effects from upstream regulation by evaluating the location and extent of dams in the drainage basin. The third step was (a) the analysis of annual hydrographs and included visual hydrograph analysis for selected stations with 20 or more years of streamflow record, (b) analysis of covariance of double-mass curves, and (c) Kendall's tau trend analysis to detect statistically significant trends in base flow, runoff, total flow, and base-flow index related to anthropogenic activity for selected stations with 15 or more years of streamflow record. A preliminary least-altered period of record for each stream was identified by removing the period of streamflow record when streams were substantially affected by anthropogenic activity. After streamflow record was removed from designation as a least-altered period, stations that did not have at least 10 years of remaining continuous streamflow record were considered to have an insufficient baseline period for modeling applications. An optimum minimum period of record was determined for each of the least-altered periods for each station to ensure a sufficient streamflow record length to provide a representative sample of annual climate variability. An optimum minimum period of 10 years or more was evaluated by analyzing the variability of annual precipitation for selected 5-, 10-, 15-, 25-, and 35-year periods for each of 20 climate divisions that contained stations used in the baseline period analysis. The distribution of annual precipitation was compared for each consecutive overlapping 5-year period to the period 1925-2007 by using a Wilcoxon rank-sum test. The least-altered period of record for stations was also compared to the period 1925-2007 by using a Wilcoxon rank-sum test. The results of this analysis were used to determine how many years of annual precipitation data were needed for the selected period to be statistically similar to the distribution of annual precipitation data for a long-term period, 1925-2007. Minimum optimum periods ranged from 10 to 35 years and varied by climate division. A final baseline period was determined for 111 stations that had a baseline period of at least 10 years of continuous streamflow record after the record-elimination process. A suitable baseline period of record for use in modeling applications could not be identified for 58 of the initial 171 stations because of substantial anthropogenic alteration of the stream or drainage basin and for 2 stations because the least-altered period of record was not representative of annual climate variability. The baseline period for each station was rated ?excellent?, ?good?, ?fair?, ?poor?, or ?no baseline period.? This rating was based on a qualitative evaluation of t

Scientific Investigations Report↗

USGS2021

Exploration of the 2016 Yellowstone River fish kill and proliferative kidney disease in wild fish populations

Proliferative kidney disease (PKD) is an emerging disease that recently resulted in a large mortality event of salmonids in the Yellowstone River (Montana, USA). Total PKD fish mortalities in the Yellowstone River were estimated in the tens of thousands, which resulted in a multi‐week river closure and an estimated economic loss of US$500,000. This event shocked scientists, managers, and the public, as this was the first occurrence of the disease in the Yellowstone River, the only reported occurrence of the disease in Montana in the past 25 yr, and arguably the largest wild PKD fish kill in the world. To understand why the Yellowstone River fish kill occurred, we used molecular and historical data to evaluate evidence for several hypotheses: Was the causative parasite Tetracapsuloides bryosalmonae a novel invader, was the fish kill associated with a unique parasite strain, and/or was the outbreak caused by unprecedented environmental conditions? We found that T. bryosalmonae is widely distributed in Montana and have documented occurrence of this parasite in archived fish collected in the Yellowstone River prior to the fish kill. T. bryosalmonae had minimal phylogeographic population structure, as the DNA of parasites sampled from the Yellowstone River and distant water bodies were very similar. These results suggest that T. bryosalmonae could be endemic in Montana. Due to data limitations, we could not reject the hypothesis that the fish kill was caused by a novel and more virulent genetic strain of the parasite. Finally, we found that single‐year environmental conditions are insufficient to explain the cause of the 2016 Yellowstone River PKD outbreak. Other regional rivers where we documented T. bryosalmonae had similar or even more extreme conditions than the Yellowstone River and similar or more extreme conditions have occurred in the Yellowstone River in the recent past, yet mass PKD mortalities have not been documented in either instance. We conclude by placing these results and unresolved hypotheses into the broader context of international research on T. bryosalmonae and PKD, which strongly suggests that a better understanding of bryozoans, the primary host of T. bryosalmonae , is required for better ecosystem understanding.

Montana↗

USGS2016

Solving for source parameters using nested array data: A case study from the Canterbury, New Zealand earthquake sequence

The seismic spectrum can be constructed by assuming a Brune spectral model and estimating the parameters of seismic moment ( M 0 ), corner frequency ( f c ), and high-frequency site attenuation ( κ ). Using seismic data collected during the 2010–2011 Canterbury, New Zealand, earthquake sequence, we apply the non-linear least-squares Gauss–Newton method, a deterministic downhill optimization technique, to simultaneously determine the M 0 , f c , and κ for each event-station pair. We fit the Brune spectral acceleration model to Fourier-transformed S-wave records following application of path and site corrections to the data. For each event, we solve for a single M 0 and f c , while any remaining residual kappa, κ r "> κ r κr , is allowed to differ per station record to reflect varying high-frequency falloff due to path and site attenuation. We use a parametric forward modeling method, calculating initial M 0 and f c values from the local GNS New Zealand catalog M w, GNS magnitudes and measuring an initial κ r "> κ r κr using an automated high-frequency linear regression method. Final solutions for M 0 , f c , and κ r "> κ r κr are iteratively computed through minimization of the residual function, and the Brune model stress drop is then calculated from the final, best-fit f c . We perform the spectral fitting routine on nested array seismic data that include the permanent GeoNet accelerometer network as well as a dense network of nearly 200 Quake Catcher Network (QCN) MEMs accelerometers, analyzing over 180 aftershocks M w,GNS ≥ 3.5 that occurred from 9 September 2010 to 31 July 2011. QCN stations were hosted by public volunteers and served to fill spatial gaps between existing GeoNet stations. Moment magnitudes determined using the spectral fitting procedure ( M w,SF ) range from 3.5 to 5.7 and agree well with M w,GNS , with a median difference of 0.09 and 0.17 for GeoNet and QCN records, respectively, and 0.11 when data from both networks are combined. The majority of events are calculated to have stress drops between 1.7 and 13 MPa (20th and 80th percentile, correspondingly) for the combined networks. The overall median stress drop for the combined networks is 3.2 MPa, which is similar to median stress drops previously reported for the Canterbury sequence. We do not observe a correlation between stress drop and depth for this region, nor a relationship between stress drop and magnitude over the catalog considered. Lateral spatial patterns in stress drop, such as a cluster of aftershocks near the eastern extent of the Greendale fault with higher stress drops and lower stress drops for aftershocks of the 2011 M w,GNS 6.2 Christchurch mainshock, are found to be in agreement with previous reports. As stress drop is arguably a method-dependent calculation and subject to high spatial variability, our results using the parametric Gauss–Newton algorithm strengthen conclusions that the Canterbury sequence has stress drops that are more similar to those found in intraplate regions, with overall higher stress drops that are typically observed in tectonically active areas.

Christchurch↗

USGS2023

Latency and geofence testing of wireless emergency alerts intended for the ShakeAlert® earthquake early warning system for the West Coast of the United States of America

ShakeAlert, the earthquake early warning (EEW) system for the West Coast of the United States, attempts to provides crucial warnings before strong shaking occurs. However, because the alerts are triggered only when an earthquake is already in progress, and the alert latencies and delivery times are platform dependent, the time between these warnings and the arrival of shaking is variable. The ShakeAlert system uses, among other public alerting platforms like a mobile phone operating system, smartphone apps, and the Federal Emergency Management Agency Integrated Public Alert & Warning System (IPAWS). IPAWS sends Wireless Emergency Alerts (WEAs) informing people via their smartphones and other mobile devices about various events, such as natural hazards, child abductions, or public health information about COVID-19. However, little is known about the IPAWS delivery latencies. Given that people may have only a few seconds of notice after they receive an alert to take a protective action before they feel earthquake shaking, quantifying latencies is critical to understanding whether the IPAWS system is useful for EEW. In this study, we developed new methods to test the IPAWS distribution system's performance, both with devices in a controlled environment and as well as with a 2019 community-based feedback form, in Oakland and San Diego County, California, respectively. The controlled environment test used mobile phones (including smart and non-smart phones) and associated devices to determine alert receipt times; the community research form had participants self-report their receipt times. By triangulating the data between the controlled test environment and the community research, we determined the latency statistics as well as whether the geofence (the geographic area where the alert was intended to be sent) held broadly. We found that the latencies were similar between the two tests despite the large differences in population sizes. WEA messages were received within a median time frame of 6–12 s, and the geofence held with only a few exceptions. We use this latency to assess how the system would have performed in two large earthquakes, the 1989 M6.9 Loma Prieta and 2019 M7.1 Ridgecrest earthquakes, which both occurred near our WEA test locations. Our analysis revealed that had IPAWS been available during those earthquakes, particularly Loma Prieta, it would have provided crucial seconds of notice that damaging shaking was imminent in some locations relatively far from the epicenter. Further, we find affordable non-smart phones can receive WEAs as fast as smartphones. Finally, our new method can be used for latency and geospatial testing going forward for IPAWS and other similar alerting systems.

California↗

USGS2009

Terrestrial Ecosystems - Land Surface Forms of the Conterminous United States

As part of an effort to map terrestrial ecosystems, the U.S. Geological Survey has generated land surface form classes to be used in creating maps depicting standardized, terrestrial ecosystem models for the conterminous United States, using an ecosystems classification developed by NatureServe . A biophysical stratification approach, developed for South America and now being implemented globally, was used to model the ecosystem distributions. Since land surface forms strongly influence the differentiation and distribution of terrestrial ecosystems, they are one of the key input layers in this biophysical stratification. After extensive investigation into various land surface form mapping methodologies, the decision was made to use the methodology developed by the Missouri Resource Assessment Partnership (MoRAP). MoRAP made modifications to Hammond's land surface form classification, which allowed the use of 30-meter source data and a 1-km2 window for analyzing the data cell and its surrounding cells (neighborhood analysis). While Hammond's methodology was based on three topographic variables, slope, local relief, and profile type, MoRAP's methodology uses only slope and local relief. Using the MoRAP method, slope is classified as gently sloping when more than 50 percent of the area in a 1-km2 neighborhood has slope less than 8 percent, otherwise the area is considered moderately sloping. Local relief, which is the difference between the maximum and minimum elevation in a neighborhood, is classified into five groups: 0-15 m, 16-30 m, 31-90 m, 91-150 m, and >150 m. The land surface form classes are derived by combining slope and local relief to create eight landform classes: flat plains (gently sloping and local relief =< 15 m), smooth plains (gently sloping and 15 m < local relief =< 30 m), irregular plains (gently sloping and 30 m < local relief =< 90 m), escarpments (gently sloping and local relief > 90 m), low hills (not gently sloping and local relief =< 30 m), hills (not gently sloping and 30 m < local relief =< 90 m), breaks/foothills (not gently sloping and 90 m < local relief =< 150 m), and low mountains (not gently sloping and local relief > 150 m). However, in the USGS application of the MoRAP methodology, an additional local relief group was used (> 400 m) to capture additional local topographic variation. As a result, low mountains were redefined as not gently sloping and 151 m < local relief < 400 m, and a new land surface form class, high mountains/deep canyons, was identified as not gently sloping and local relief > 400 m. The final application of the MoRAP methodology was implemented using the USGS 30-meter National Elevation Dataset and an existing USGS slope dataset that had been derived by calculating the slope from the NED in Universal Transverse Mercator (UTM) coordinates in each UTM zone, and then combining all of the zones into a national dataset. This map shows a smoothed image of the nine land surface form classes based on MoRAP's methodology. Additional information about this map and any data developed for the ecosystems modeling of the conterminous United States is available online at http://rmgsc.cr.usgs.gov/ecosystems/.

Scientific Investigations Map↗

USGS2023

Estimating flood magnitude and frequency for unregulated streams in Wisconsin

Flood frequency characteristics and estimated flood discharges for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities were computed at 299 streamgaged locations in Wisconsin. The State was divided into four flood frequency regions using a cluster analysis to produce regions which are homogeneous with respect to physical basin characteristics. Regression equations relating flood discharges to basin characteristics within each region were developed and can be used to estimate flood discharges at ungaged locations in Wisconsin. Basin characteristics included in the final regression equations include drainage area, saturated hydraulic conductivity, percent forest, percent herbaceous upland, percent open water, and the maximum 24-hour precipitation with a 10-year recurrence interval. The standard error of prediction for regression equations ranges between 40 and 71 percent, and the pseudo coefficient of determination ranges between 0.8 and 0.95. Nonmonotonic trends in the annual peak flow time series in the southwest part of the State are producing bias in some flood discharge estimates at streamgages with shorter (less than 20 years) periods of record. This bias increases the uncertainty in regression equations in this flood frequency region.

Wisconsin↗

USGS1986

Three-dimensional steady-state simulation of flow in the sand-and-gravel aquifer, southern Escambia County, Florida

The sand-and-gravel aquifer is the only freshwater aquifer in southern Escambia County, Florida and is the source of public water supply for the area, including the City of Pensacola. The aquifer was simulated by a two-layer, digital model to provide hydrologic information for water resource planning. The lower layer represents the main-producing zone; the upper layer represents all of the aquifer above the main-producing zone including an unconfined zone and discontinuous perched, confined , and confining zones. The model was designed for steady-state simulation and predicts the response of the aquifer (changes in water levels) to groundwater pumping where steady-state conditions have been reached. Input to the model includes matrices representing constant-head nodes, starting head, transmissivity of layer 1, leakance between layers 1 and 2, lateral hydraulic conductivity of layer 2, and altitude of the base layer 2. The sources of water to the model are from recharge by infiltrated precipitation (estimated from base runoff), inflow across boundaries, and induced recharge from river leakance in periods of prolonged groundwater pumping. Model output includes final head and drawdown for each layer and total values for discharge and recharge in the model area. The model was calibrated for 1972 pumping and tested by simulating pumpages during 1939-40, 1958, and 1977. Sensitivity analyses showed water levels in both layers were most sensitive to changes in the recharge matrix and least sensitive to river leakage. Suggestions for further development of the model include subdivision and expansion of the grid, assignment of storage coefficients for transient simulations, more intensive study of the stream-aquifer relations, and consideration of the effects of infiltration basins on recharge. (Author 's abstract)

Water-Resources Investigations Report↗

USGS1985

Selected papers in the hydrologic sciences, 1985; May 1985

The University of Minnesota, the Minnesota Geological Survey, and the U.S. Geological Survey are studying the feasibility of storing water at a temperature of 150 degrees Celsius in the Franconia-Ironton-Galesville aquifer. The Aquifer Thermal-Energy Storage project has a doublet-well design with a well spacing of approximately 250 meters. One well will be used for cool-water supply, and, the other, for hot-water injection. The U.S. Geological Survey is constructing a model of ground-water flow and thermal-energy transport to aid in determining the efficiency of the Aquifer Thermal Energy Storage system. A preliminary model of radial flow and thermal-energy transport was constructed, based on hydraulic and thermal properties of the Franconia-Ironton-Galesville aquifer determined in previous studies. The model was used to investigate the sensitivity of model results to various hydraulic and thermal properties and to study the potential for buoyancy flow within the aquifer and the effect of various cyclic injection-withdrawal schemes on the relative thermal efficiency of the aquifer. Sensitivity analysis was performed assuming 8 days of injection of 150-degree-Celsius water at 18.9 liters per second, 8 days of storage, and 8 days of withdrawal of hot water at 18.9 liters per second. The analysis indicates that, for practical ranges of hydraulic and thermal properties, rock-heat capacity is the least important property and thermal dispersivity is the most important property used to compute temperature and aquifer thermal efficiency. The amount of buoyancy flow was examined for several values of hydraulic conductivity and ratios of horizontal to vertical hydraulic conductivities. For the assumed base values of hydraulic and thermal properties, buoyancy flow was negligible. The greatest simulated buoyancy flow resulted from simulations in which horizontal hydraulic conductivity was increased to 10 times the base value, and the vertical hydraulic conductivity was set equal to the horizontal hydraulic conductivity. The effects of various injection-withdrawal rates and durations on computed values of aquifer relative thermal efficiency and final well-bore temperature were studied for five 1-year hypothetical test cycles of injection and withdrawal. The least efficient scheme was 8 months injection of 150-degree-Celsius water and 4 months of withdrawal of hot water at 18.9 liters per second. The most efficient scheme was obtained with 6 months of injection of 150-degree-Celsius water at 18.9 liters per second and 6 months of withdrawal of hot water at 37.8 liters per second. The hypothetical simulations indicate that the subsequent calibrated model of the doublet-well system will be a valuable tool in determining the most efficient system operation.

Water Supply Paper↗

USGS1990

Conceptualization and analysis of ground-water flow system in the Coastal Plain of Virginia and adjacent parts of Maryland and North Carolina

This report presents the results of a study of the ground-water flow system in the Coastal Plain of Virginia and adjacent parts of Maryland and North Carolina. The ground-water flow system consists of water-table aquifer and an underlying sequence of confined aquifers and intervening confining units composed of unconsolidated sand and clay .Water levels have declined steadily, and cones of depression have expanded and coalesced around major ground-water withdrawal centers. A digital flow model was developed to enhance knowledge of the behavior of the ground-water flow system in response to its development. Transmissivity and vertical leakance maps were developed for each aquifer and confining unit. The model was calibrated to simulate ground-water flow within the system under both prepumping and pumping conditions. Simulated prepumping potentiometric-surface maps indicate that regional movement of ground water was from the Fall Line toward coastal areas and that local movement of ground water was from interfluves toward major river valleys. Maps of simulated prepumping flow across confining units show that most recharge occurred in narrow bands approximately parallel to the Fall Line and under interfluves and that discharge was toward major river valleys and coastal water. Simulated prepumping rates of recharge into the confined aquifer system from the water-table aquifer varied up to 3 .2 inches per year (in/yr), and rates of discharge out of the confined system varied up to 2.8 in/yr. Ten pumping periods covering 90 years (yr) of withdrawal simulated the history of ground-water development. Simulated potentiometric-surface maps for 1980 show lowered water levels and the development of coalescing cones of depression around the cities of Franklin, Suffolk, and Williamsburg and the town of West Point, all in Virginia. The largest simulated decline in water level, about 210 feet (ft) ,was near Franklin. Water budgets indicate that over the period of simulation (1891-1980) (1) pumpage from the model area increased by about 105 million gallons per day (Mgal/d), (2) lateral boundary outflow increased by about 5 Mgal/d, (3) ground-water flow to streams and coastal water decreased by about 107.5 Mgal/d, (4) lateral boundary inflow increased by about0.7 Mgal/d, and (5) water released from aquifer storage increased by about 1 .6 Mgal/d. The difference between total inflow and total outflow is the numerical truncation error of the digital simulation. Analysis of water budgets for individual confined aquifers shows that the major source of water supplied to wells was vertical leakage induced through confining units by pumping. Simulated rates of recharge into the confined aquifer system at the end of the final pumping period (1980) varied up to 3.8in/yr,and simulated rates of discharge out of the confined system varied up to 2.2in/yr. Results of simulations show an increase of about 110Mgal/d into the confined system from the unconfined system over the period of simulation. This increase in flow into the confined system affected local discharge of ground water to streams and regional discharge to coastal water. Withdrawal of groundwater from the confined aquifers also induced brackish water from Chesapeake Bay into the confined system. Results of sensitivity analyses indicate that simulated water levels are more sensitive to decreases in aquifer transmissivity and confining unit vertical hydraulic conductivity than to increases in these properties. Lowering the storage coefficient of an aquifer had minimal effect on simulated water levels, whereas increasing the storage coefficient had a much more significant effect. The effect of confining unit storage is shown to be insignificant if it is assumed that the water released from confining unit storage is attributable to the compressibility of water only.

Maryland, North Carolina, Virginia↗

USGS2010

Baseline survey for rare plant species and native plant communities within the Kamehameha Schools 'Lupea Safe Harbor Planning Project Area, North Kona District, Island of Hawai'i

Kamehameha Schools, in conjunction with several federal, state, and private organizations, has proposed to conduct conservation management on approximately 5,340 ha (~13,200 acres) of land they own in the vicinity of Kīpukalupea in the North Kona District on the island of Hawai'i. The goal of this program is to restore and enhance the habitat to benefit native plant and animal populations that are currently, or were formerly, found in this site. The initial phase of this project has been focused on various activities including conducting baseline surveys for bird and plant species so Kamehameha Schools could develop a Safe Harbor Agreement (SHA) for the proposed project lands relative to the habitat management and species reintroduction efforts they would like to conduct in the Lupea Project area. This report summarizes methods that were used to collect field data on plant species and communities within the project area, and the results of that initial survey. The information was used to calculate baseline values for all listed threatened or endangered plant species found, or expected to be found, within the project area, and to design a monitoring program to assess changes in plant communities and rare plant species relative to management activities over the duration of the SHA. The Lupea Project area contains excellent examples of several high elevation native plant communities including montane dry forest and woodland, native subalpine shrubland, and native grassland. Between November 2003 and January 2004 we sampled plant communities and species along seven transects established through the project area. A total of 109 plant species were found during this survey, within the transect grid and in nearby areas. Forty-four of these plants are endemic species, 21 are indigenous species, 43 are introduced, and one species is believed to have been introduced to Hawai&bdquo;i by early Polynesian settlers. Only one federally listed Endangered plant, Asplenium peruvianum var. insulare , was found within the survey area. Additionally, we found one immature plant that may be Sicyos macrophyllus , a candidate species for listing. However, we were not able to make a definite determination of this species‟ identity since it did not have fruits or flowers. Finally, we documented four plant species within the survey area that have no official status designation but are considered to be rare and informally recognized as “species of concern” (SOC) as they appear to be declining in distribution and abundance statewide. These included Chamaesyce olowaluana, Eragrostis deflexa, Sisyrinchium acre, and Tetramolopium consanguineum . In addition to conducting field surveys, we performed a query on a spatial database developed by Dr. Jonathan Price of the University of Hawai&bdquo;i at Hilo which models the potential range of all native Hawaiian plant species based on historic observations and a set of environmental parameters. The potential species list for the Lupea Project area includes 47 taxa that we did not find during our surveys, as well as three other listed species that were not modeled by Price, but known from historic records in adjacent habitats. Some of these species are extremely rare or, in some cases have been locally extirpated. However, most of the plants that were predicted but not found during our surveys are expected to be located with additional searching, or may potentially recolonize the area following the elimination of ungulates and initiation of other restoration efforts. Forty-four introduced plant species were found within the survey area, seven of which are considered to be highly invasive. These include the grasses Pennisetum clandestinum and Pennisetum setaceum , vines Delairea odorata and Passiflora tarminiana , herbs Senecio madagascariensis and Verbascum thapsus , and the shrub Rubus niveus . Non-zero baseline values are proposed for the one listed plant species found within the Lupea Project area, one species that is a candidate for listing, and the four other rare species we found that may be considered for listing in the future. Additionally, a zero baseline is proposed for 23 other species that were predicted, but not found within the project area. These include 14 Endangered species, one Threatened species, two candidates for listing, and six species of concern. Subsequent monitoring of the site will be necessary to determine if the populations of these species have increased or decreased relative to their baseline values. It is presumed that the management activities Kamehameha Schools has proposed for this area, particularly removal of the ungulates and weed control, will provide a benefit to the habitat as a whole and allow for natural regeneration and maintenance of the all elements of the plant communities found there.

Hawaii↗

USGS2009

A landscape indicator approach to the identification and articulation of the consequences of land-cover change in the Mid-Atlantic Region, 1973-2001

Landscape indicators, derived from land-use and land-cover data, hydrology, nitrate deposition, and elevation data, were used by Jones and others (2001a) to calculate the ecological consequences of land-cover change. Nitrate loading and physical bird habitat were modeled from 1973 and 1992 land-cover and other spatial data for the Mid-Atlantic region. Utilizing the same methods, this study extends the analysis another decade with the use of the 2001 National Land Cover Dataset. Land-cover statistics and trends are calculated for three time periods: 1973-1992, 1992-2001 and 1973-2001. In addition, high-resolution aerial photographs (1 meter or better ground-sample distance) were acquired and analyzed for thirteen pairs of adjacent USGS 7.5 minute quadrangle maps in areas where distinct positive or negative changes to nitrogen loading and bird habitat were previously calculated. During the entire 30 year period, the data show that there was extensive loss of agriculture and forest area and a major increase in urban land-cover classes. However, the majority of the conversion of other classes to urban occurred during the 1992-2001 period. During the 1973-1992 period, there was only moderate increase in urban area, while there was an inverse relationship between agricultural change and forest change. In general, forest gain and agricultural loss was found in areas of improving landscape indicators, and forest loss and agricultural gain was found to occur in areas of declining indicators related to habitat and nitrogen loadings, which was generally confirmed by the aerial photographic analysis. In terms of the specific model results, bird habitat, which is mainly related to the extent of forest cover, declined overall with forest extent, but was also affected more in the decline of habitat quality. Nitrate loading, which is mainly related to agricultural land cover actually improved from 1992-2001, and in the overall study, mainly due to the conversion of agriculture to forests and urban. The high-resolution imagery analysis was significant in that it confirmed, at a very local level, the specific land-cover changes that were driving the landscape metrics and model results that were calculated from moderate resolution land-cover data and models. These were generally subtle changes in patch size of agriculture, forest, and urban areas, but had substantial effects on bird habitat and nitrogen loadings. This analysis of high-resolution imagery demonstrates and confirms the important ability of moderate-resolution land-cover data to capture significant landscape-level activity that is directly related to specific metrics of ecological significance. It also demonstrates consistent landscape-scale relationships between data derived from high-resolution, moderate-resolution and landscape-model sources. Finally, many of the areas of improvement and decline in bird habitat and nitrogen loadings appear to be potentially regional in nature and likely reflect some local trend in landscape activity. Although the use of ecoregions as sampling units has been criticized in recent years, these results show that basic changes in Level 1 land-cover categories, such as forest and agriculture, may still reflect ecoregional patterns and considerations at some scale of mapping and analysis. This is a potentially important area for future landscape-indicator research. This and other follow-on research opportunities are discussed.

Open-File Report↗

USGS2017

Methods for estimating selected low-flow frequency statistics and mean annual flow for ungaged locations on streams in North Georgia

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources, Environmental Protection Division, developed regional regression equations for estimating selected low-flow frequency and mean annual flow statistics for ungaged streams in north Georgia that are not substantially affected by regulation, diversions, or urbanization. Selected low-flow frequency statistics and basin characteristics for 56 streamgage locations within north Georgia and 75 miles beyond the State’s borders in Alabama, Tennessee, North Carolina, and South Carolina were combined to form the final dataset used in the regional regression analysis. Because some of the streamgages in the study recorded zero flow, the final regression equations were developed using weighted left-censored regression analysis to analyze the flow data in an unbiased manner, with weights based on the number of years of record. The set of equations includes the annual minimum 1- and 7-day average streamflow with the 10-year recurrence interval (referred to as 1Q10 and 7Q10), monthly 7Q10, and mean annual flow. The final regional regression equations are functions of drainage area, mean annual precipitation, and relief ratio for the selected low-flow frequency statistics and drainage area and mean annual precipitation for mean annual flow. The average standard error of estimate was 13.7 percent for the mean annual flow regression equation and ranged from 26.1 to 91.6 percent for the selected low-flow frequency equations. The equations, which are based on data from streams with little to no flow alterations, can be used to provide estimates of the natural flows for selected ungaged stream locations in the area of Georgia north of the Fall Line. The regression equations are not to be used to estimate flows for streams that have been altered by the effects of major dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, or wastewater discharges. The regression equations should be used only for ungaged sites with drainage areas between 1.67 and 576 square miles, mean annual precipitation between 47.6 and 81.6 inches, and relief ratios between 0.146 and 0.607; these are the ranges of the explanatory variables used to develop the equations. An attempt was made to develop regional regression equations for the area of Georgia south of the Fall Line by using the same approach used during this study for north Georgia; however, the equations resulted with high average standard errors of estimates and poorly predicted flows below 0.5 cubic foot per second, which may be attributed to the karst topography common in that area. The final regression equations developed from this study are planned to be incorporated into the U.S. Geological Survey StreamStats program. StreamStats is a Web-based geographic information system that provides users with access to an assortment of analytical tools useful for water-resources planning and management, and for engineering design applications, such as the design of bridges. The StreamStats program provides streamflow statistics and basin characteristics for U.S. Geological Survey streamgage locations and ungaged sites of interest. StreamStats also can compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in Georgia.

Georgia↗