USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth System Science Data”

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.

519 records · Page 29Linked to original sources
USGS2019

Physicochemical models of effusive rhyolitic eruptions constrained with InSAR and DEM data: A case study of the 2011-2012 Cordon Caulle eruption

The 9 month long 2011-2012 eruption of Cordon Caulle (Southern Andes, Chile) is the best instrumentally recorded rhyolitic eruption to date and the first time that the effusion of a rhyolitic flow has been observed in detail. We use Interferometric Synthetic Aperture Radar (InSAR), with time-lapse DEMs and numerical models to study the dynamics of coupled magma reservoir deflation and lava effusion. InSAR recorded ~2.2-2.5 m of deflation after the first three days of the eruption, which can be modeled using a spheroidal magma reservoir at a depth of ∼5 km, ∼20 km long, and with a pressure drop of ~20-30 MPa. The source is elongated in the NW-SE direction and its large dimensions imply a large plumbing system spanning neighboring volcanoes and active throughout the eruption, with a slight change halfway through the effusive phase. TanDEM-X DEMs record the extrusion of both the rhyolitic lava flow and the intrusion of a shallow laccolith around the eruptive vent, with a total volume of ~1.2 km3 DRE. The laccolith was emplaced during the first month of the eruption, during both the eruption explosive and effusive stages. Both the reservoir pressure drop and the extruded volume time series follow quasi-exponential trends, and can be explained by a model that couples the reservoir pres- sure decrease, time- and pressured ependent variations in the magma properties inside of the reservoir, and conduit flow. This model predicts both the temporal evolution and amplitude of both time series, and a magma compressibility of ∼10^−10 Pa−1, half the compressibility of the magma of the sub-Plinian explosive phase. Further, we estimate that the reservoir contained 1-3 wt.% dissolved H2O at the onset of lava effusion, with no exsolved CO2 and H2O in the reservoir throughout the eruption. This is in accord with a magma that was significantly degassed after the explosive phase. These remaining volatiles might have been responsible for magma fragmentation, consistent with the hybrid explosive and effusive style observed during the waning of the eruption.

Cordón Caulle↗
USGS2012

Archive of post-Hurricane Charley coastal oblique aerial photographs collected during U.S. Geological Survey field activity 04CCH01 from Marco Island to Fort DeSoto, Florida, August 15, 2004

On August 15, 2004, the U.S. Geological Survey (USGS) conducted an oblique aerial photographic survey off the southwest coast of Florida, from Marco Island to Fort DeSoto, aboard a Navajo Chieftain airplane, tail number N2KK, at an altitude of 500 ft and approximately 1000 ft offshore. These photographs were used to document coastal changes such as beach erosion and overwash caused by Hurricane Charley. They will also be used as baseline data for future coastal change. The oblique photography also served as qualitative ground truthing for the Experimental Advanced Airborne Research Lidar (EAARL) coastal topography and bathymetry data collected on August 16, 2004 (Bonisteel and others, 2009). This report serves as an archive of photographs collected during the August 15, 2004, post-Hurricane Charley coastal oblique aerial survey along with associated flight path maps, KML files, navigation files, digital Field Activity Collection System (FACS) logs, and Federal Geographic Data Committee (FGDC) metadata. Refer to the Acronyms page for expansions of all acronyms and abbreviations used in this report. The USGS St. Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 04CCH01 tells us the data were collected in 2004 for the Coastal Change Hazards (CCH) study and the data were collected during the first field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the ID number. Two separate records of flight navigation were collected during the survey. The first was a continuous ASCII text file from the PLGR that recorded only latitudes, longitudes, and altitudes every 30 sec for the entire flight. No time values were recorded by the PLGR. The second navigation record was recorded by a Trimble Centurion GPS and converted to subtitles on the video, using a Compix Titler unit. The video was shot continuously during the survey. The video subtitles recorded day, month, year, latitude, longitude, and time in hours, minutes, and seconds. In order to produce a digital record of the navigation values that included latitude, longitude, and time, each was manually extracted from the video every 5 min, and these values were matched to the latitude and longitude in the PLGR file. Next, the time was interpolated between these 5-min fixes using Excel to produce time values for each navigation fix recorded in the PLGR file. The location of each photograph taken was determined in the following manner. A Nikon MF-14 data back marks the time each photograph was acquired on the lower right corner of the image in day, hour, and minute format (in UTC). These values were entered from the photographs into an Excel spreadsheet. It is assumed for the purposes of locating the images that the photographs were taken at a constant rate during any given minute of flight. To assign the time value in seconds to each photograph, the number of photographs taken during each minute was evenly distributed across that minute. For example, if 15 photographs were taken during minute 19:00:00, we assume that a picture was taken every 4 sec. The photographs were assigned the time values 19:00:00, 19:00:04, 19:00:08, and so on. The video time navigation file was then merged with the new interpolated photograph file based on time to produce the point of which each photograph was collected. As a result, the positions assigned to each photograph are an estimate of the aircraft position, not the location of the landmark photographed. The photographs provided here are JPEG scanned images of the analog slides. The metadata values for photo creation time, GPS latitude, GPS longitude, GPS position, keywords, credit, artist, caption, copyright, and contact were added to each photograph's EXIF header using EXIFtools (see the Software page). Photographs can be opened directly with any JPEG-compatible image viewer by clicking on a thumbnail on the contact sheet, or, when viewing the Google Earth KML file, by clicking on the marker and then clicking on either the thumbnail or the link below the thumbnail. The KML files were created using the photographic navigation files and the CreateKML JavaScript (see the Software page). To view the survey maps and navigation files, and for more information about these items, see the Navigation page. Figure 1 displays the acquisition geometry. The tables provide detailed information about the assigned location, name, data, and time the photograph was taken along with links to the photo and corresponding 5-min contact sheet. Refer to table 1 and table 2 for details of the northern and southern county photographs, respectively.

Florida↗
USGS2021

Recency of faulting and subsurface architecture of the San Diego Bay pull-apart basin, California, USA

In southern California, plate boundary motion between the North American and Pacific plates is distributed across several sub-parallel fault systems. The offshore faults of the California Continental Borderland (CCB) are thought to accommodate ~10-15% of the total plate boundary motion, but the exact distribution of slip and the mechanics of slip partitioning remain uncertain. The Newport-Inglewood-Rose Canyon fault is the easternmost fault within the CCB whose southern segment splays out into a complex network of faults beneath San Diego Bay. A pull-apart basin model between the Rose Canyon and the offshore Descanso fault has been used to explain prominent fault orientations and subsidence beneath San Diego Bay; however this model does not account for faults in the southern portion of the bay or faulting east of the bay. To investigate the characteristics of faulting and stratigraphic architecture beneath San Diego Bay, we combined a suite of reprocessed legacy airgun multi-channel seismic profiles and high-resolution Chirp data, with age and lithology controls from geotechnical boreholes and shallow sub-surface vibracores. This combined dataset is used to create gridded horizon surfaces, fault maps, and perform a kinematic fault analysis. The structure beneath San Diego Bay is dominated by down-to-the-east motion on normal faults that can be separated into two distinct groups. The strikes of these two fault groups can be explained with a double pull-apart basin model for San Diego Bay. In our conceptual model, the western portion of San Diego Bay is controlled by a right-step between the Rose Canyon and Descanso faults, which matches both observations and predictions from laboratory models. The eastern portion of San Diego Bay appears to be controlled by an inferred step-over between the Rose Canyon and San Miguel-Vallecitos faults and displays distinct fault strike orientations, which kinematic analysis indicates should have a significant component of strike-slip partitioning that is not detectable in the seismic data. The potential of a Rose Canyon-San Miguel-Vallecitos fault connection would effectively cut the stepover distance in half and have important implications for the seismic hazard of the San Diego-Tijuana metropolitan area (population ~3 million people).

California↗
USGS2023

Climate change hotspots and implications for the global subsea telecommunications network

A global network of subsea telecommunications cables underpins our daily lives, enabling >95% of global digital data transfer, $trillions/day in financial trading, and providing critical communications links, particularly to remote, low-income countries. Despite their importance, subsea cables and their landing stations are vulnerable to damage by natural hazards, including storm surges, waves, cyclones, earthquakes, floods, volcanic eruptions , submarine landslides and ice scour. However, the likelihood or recurrence interval of many of these types of events will likely change under future projected climate change scenarios, compounded by sea-level rise; potentially increasing hazard severity, creating previously unanticipated hazards, or hazards may shift to new locations during the 20–30-year operational life of cable systems. To date, no study has assessed the wide-reaching impacts of future climate change on subsea cables and landing stations on a global scale. Here, for the first time we synthesize the current evidence base, based on published peer-reviewed datasets, to fill this crucial knowledge gap, specifically to assess how and where future climate change is likely to impact subsea cables and their shore-based infrastructure. We find that ocean conditions are highly likely to change on a global basis as a result of climate change, but the feedbacks and links between climate change, natural processes and human activities are often complicated, resulting in a high degree of geographic variability. We identify climate change ‘hotspots’ (regions and locations likely to experience the greatest impacts) but find that not all areas will be affected in the same manner, nor synchronously by the same processes. We conclude that cable routes should carefully consider locally-variable drivers of hazard frequency and magnitude. Consideration should be given both to instantaneous events (e.g. landslides, tropical cyclones) as well as longer-term, sustained impacts (e.g. seabed currents that circulate even in deep water). Multiple factors can combine to increase the risk posed to subsea cables, hence a holistic approach is essential to assess the compounded effects of both natural processes and human activities in the future.

Earth Science Reviews↗
USGS2024

Anomalously high relief on Denali, Alaska, caused by tectonic, lithologic, and climatic drivers

We assess the growth of anomalously high relief on Denali, located in the Alaska Range, Alaska, and the tallest mountain in North America (6190 masl). Denali is 3000 m taller than most surrounding peaks. It lies inside a 19° restraining bend in the active Denali fault system that is moving at about 7 mm/yr, and thus there is a tectonic and structural driver for ongoing rock uplift. High relief around Denali is also due, in part, to its granitic rock type and low fracture density relative to adjacent metasediments. Here we show that unique climatic conditions at high elevations also contribute to the growth of relief. We examine 10Be concentrations in 1) three new gravel samples between 3500 and 5200 m elevation from sites unaffected by recent glacial erosion, 2) previously published samples from a sidewall of the Kahiltna Glacier from 2400 to 2800 masl, 3) previously published samples from medial moraines along the length of the Kahiltna Glacier from 1400 to ~500 masl, and 4) previously published alluvial samples along the Kahiltna River at an elevation of ~200 masl. These samples constitute a transect extending more than 5000 vertical meters, and the data establish that erosion rates decrease with elevation and contribute to the growth of relief. Erosion rates for the three new high-elevation samples are calculated to 4.6±0.6 mm/ka at 5.2 km, 28.6±3.7 mm/ka at 4.0 km, and 38±5 mm/ka at 3.5 km. Erosion rates at intermediate elevations, on the sidewall of the Kahiltna Glacier, range between 160-327 mm/ka. Along the medial moraines inferred erosion rates range between 140-537 mm/ka, and basin-wide erosion rates calculated from sediments in the river below the glacier range between 450-896 mm/ka. These differences in erosion rates can create relief of 3 km within 1-10 Ma, well within the estimated period of increase in rock uplift and exhumation on Denali over the last ~6 Ma. Meteorological data from 2130 to 5550 masl at 5 sites show it rarely gets above freezing above 4 km elevation, indicating that frost weathering currently plays a diminished role in erosion at high elevations. The immediate implication of this temperature and erosional correlation is an increase in relief. This is the first study to directly measure a significant decrease in erosion rates at high elevations in the relative absence of frost weathering. The results highlight the combined influence of rock type, glacial erosion, and permanent sub-zero temperatures on erosion rates. In combination with active faulting, the data explain the resultant increase in relief along the southern side of the Alaska Range over the past 100 ka.

Alaska↗
USGS2014

U.S. Geological Survey Karst Interest Group Proceedings, Carlsbad, New Mexico, April 29-May 2, 2014

Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex, and involve biological, chemical, and physical changes. These factors, combined with the diverse climatic regimes under which karst development in these rocks has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs has resulted in the most exceptional karst areas of the United States being designated as national or state parks; commercial caverns and known privately owned caves number in the tens of thousands. Both public and private properties provide access for scientists to study the flow of groundwater in situ. Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and as unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas can impact the ecosystem and water quality of these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Ground-Water Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as other federal and state agencies, and with researchers from academia and institutes. The KIG also encourages younger scientists by participation of students in the poster and oral sessions. To accomplish its mission, the KIG has organized a series of workshops that are held near nationally important karst areas. To date (2014) six KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst features. Proceedings of the workshops are published by the USGS and are available online at http://water.usgs.gov/ogw/karst/kig . The first KIG workshop was held in St. Petersburg, Florida, February 13–16, 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held August 20–22, 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held September 12–15, 2005, in Rapid City, South Dakota, nearby to karst features in evaporites and limestones of the Madison Group in the Black Hills of South Dakota, including Wind Cave National Park and Jewel Cave National Monument. The workshop also included a featured presentation by Thomas Casadevall, Central Region Director, USGS, on the status of earth science at the USGS and evening trips to Jewel Cave led by Mike Wiles, National Park Service (NPS) and Wind Cave led by Rod Horrocks, NPS. The fourth KIG workshop was held May 27–29, 2008, and hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave with Rickard Toomey and Rick Olsen, NPS. The fifth workshop was held April 26–29, 2011, and was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River of the southern Ozarks and a keynote presentation on paleokarst in the United States by Art and Peggy Palmer. This sixth and current 2014 KIG workshop is hosted by the National Cave and Karst Research Institute (NCKRI) in Carlsbad, New Mexico, with Director of NCKRI, George Veni, serving as co-chair of the workshop with Eve Kuniansky, USGS. The session planning committee for this sixth workshop includes Van Brahana, USGS retired and University of Arkansas Professor Emeritus; Tom Byl, USGS and Tennessee State University; Zelda Bailey, former Director of NCKRI and retired Director, National Institute of Standards and Technology, Boulder Laboratory, Colorado; Patrick Tucci, USGS retired; and Mike Bradley, Allan Clark, Geoff Delin, Daniel Doctor, James Kaufmann, Eve Kuniansky, Randy Orndorff, Larry Spangler, and Dave Weary of the USGS. The karst hydrology field trip on Thursday will be led by Lewis Land (NCKRI karst hydrologist) and the optional Friday field trip on the geology of Carlsbad Caverns National Park will be led by George Veni. The keynote speaker is Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at NCKRI, who will address the future of karst research. Additionally, there is a featured presentation “Irish karst and its management,” by Caoimhe Hickey, The Geological Survey of Ireland, preceding a panel discussion on “Collaboration During Times of Limited Resources.” The extended abstracts of USGS authors were peer reviewed and approved for publication by the U.S. Geological Survey. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to our editorial standards or stratigraphic nomenclature and is not research conducted or data collected by the USGS. However, all articles had at a minimum of two peer reviews, and all articles were edited for consistency of appearance in the published Proceedings. The use of trade, firm or product names in any article is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS, Office of Groundwater, provides technical support for the Karst Interest Group website and public availability of the Proceedings from these workshops, and the USGS Groundwater Resources Program funds the publication costs. Finally, the cover illustration is the work of Ann Tihansky, USGS, used since the first KIG workshop in 2000.

Scientific Investigations Report↗
USGS2015

Spatial and seasonal responses of precipitation in the Ganges and Brahmaputra river basins to ENSO and Indian Ocean dipole modes: Implications for flooding and drought

We evaluated the spatial and seasonal responses of precipitation in the Ganges and Brahmaputra basins as modulated by the El Niño Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) modes using Global Precipitation Climatology Centre (GPCC) full data reanalysis of monthly global land-surface precipitation data from 1901 to 2010 with a spatial resolution of 0.5° × 0.5°. The GPCC monthly total precipitation climatology targeting the period 1951–2000 was used to compute gridded monthly anomalies for the entire time period. The gridded monthly anomalies were averaged for the years influenced by combinations of climate modes. Occurrences of El Niño alone significantly reduce (88% of the long-term average (LTA)) precipitation during the monsoon months in the western and southeastern Ganges Basin. In contrast, occurrences of La Niña and co-occurrences of La Niña and negative IOD events significantly enhance (110 and 109% of LTA in the Ganges and Brahmaputra Basin, respectively) precipitation across both basins. When El Niño co-occurs with positive IOD events, the impacts of El Niño on the basins' precipitation diminishes. When there is no active ENSO or IOD events (occurring in 41 out of 110 years), precipitation remains below average (95% of LTA) in the agriculturally intensive areas of Haryana, Uttar Pradesh, Rajasthan, Madhya Pradesh, and Western Nepal in the Ganges Basin, whereas precipitation remains average to above average (104% of LTA) across the Brahmaputra Basin. This pattern implies that a regular water deficit is likely, especially in the Ganges Basin, with implications for the agriculture sector due to its reliance on consistent rainfall for successful production. Historically, major droughts occurred during El Niño and co-occurrences of El Niño and positive IOD events, while major flooding occurred during La Niña and co-occurrences of La Niña and negative IOD events in the basins. This observational analysis will facilitate well-informed decision making in minimizing natural hazard risks and climate impacts on agriculture, and supports development of strategies ensuring optimized use of water resources in best management practice under a changing climate.

Brahmaputra Basin, Ganges Basin↗
USGS2021

Determination of vadose zone and saturated zone nitrate lag times using long-term groundwater monitoring data and statistical machine learning

In this study, we explored the use of statistical machine learning and long-term groundwater nitrate monitoring data to estimate vadose zone and saturated zone lag times in an irrigated alluvial agricultural setting. Unlike most previous statistical machine learning studies that sought to predict groundwater nitrate concentrations within aquifers, the focus of this study was to leverage available groundwater nitrate concentrations and other environmental variables to determine mean regional vertical velocities (transport rates) of water and solutes in the vadose zone and saturated zone (3.50 and 3.75 m yr −1 , respectively). The statistical machine learning results are consistent with two primary recharge processes in this western Nebraska aquifer, namely ( 1 ) diffuse recharge from irrigation and precipitation across the landscape and ( 2 ) focused recharge from leaking irrigation conveyance canals. The vadose zone mean velocity yielded a mean recharge rate (0.46 m yr −1 ) consistent with previous estimates from groundwater age dating in shallow wells (0.38 m yr −1 ). The saturated zone mean velocity yielded a recharge rate (1.31 m yr −1 ) that was more consistent with focused recharge from leaky irrigation canals, as indicated by previous results of groundwater age dating in intermediate-depth wells (1.22 m yr −1 ). Collectively, the statistical machine learning model results are consistent with previous observations of relatively high water fluxes and short transit times for water and nitrate in the primarily oxic aquifer. Partial dependence plots from the model indicate a sharp threshold in which high groundwater nitrate concentrations are mostly associated with total travel times of 7 years or less, possibly reflecting some combination of recent management practices and a tendency for nitrate concentrations to be higher in diffuse infiltration recharge than in canal leakage water. Limitations to the machine learning approach include the non-uniqueness of different transport rate combinations when comparing model performance and highlight the need to corroborate statistical model results with a robust conceptual model and complementary information such as groundwater age.

Nebraska↗
USGS2022

Bedrock depth influences spatial patterns of summer baseflow, temperature and flow disconnection for mountainous headwater streams

In mountain headwater streams, the quality and resilience of summer cold-water habitat is generally regulated by stream discharge, longitudinal stream channel connectivity and groundwater exchange. These critical hydrologic processes are thought to be influenced by the stream corridor bedrock contact depth (sediment thickness), a parameter often inferred from sparse hillslope borehole information, piezometer refusal and remotely sensed data. To investigate how local bedrock depth might control summer stream temperature and channel disconnection (dewatering) patterns, we measured stream corridor bedrock depth by collecting and interpreting 191 passive seismic datasets along eight headwater streams in Shenandoah National Park (Virginia, USA). In addition, we used multi-year stream temperature and streamflow records to calculate several baseflow-related metrics along and among the study streams. Finally, comprehensive visual surveys of stream channel dewatering were conducted in 2016, 2019 and 2021 during summer low flow conditions (124 total km of stream length). We found that measured bedrock depths along the study streams were not well-characterized by soils maps or an existing global-scale geologic dataset where the latter overpredicted measured depths by 12.2 m (mean) or approximately four times the average bedrock depth of 2.9 m. Half of the eight study stream corridors had an average bedrock depth of less than 2 m. Of the eight study streams, Staunton River had the deepest average bedrock depth (3.4 m), the coldest summer temperature profiles and substantially higher summer baseflow indices compared to the other study steams. Staunton River also exhibited paired air and water annual temperature signals suggesting deeper groundwater influence, and the stream channel did not dewater in lower sections during any baseflow survey. In contrast, Paine Run and Piney River did show pronounced, patchy channel dewatering, with Paine Run having dozens of discrete dry channel sections ranging from 1 to greater than 300 m in length. Stream dewatering patterns were apparently influenced by a combination of discrete deep bedrock (20+ m) features and more subtle sediment thickness variation (1–4 m) depending on local stream valley hydrogeology. In combination, these unique datasets show the first large-scale empirical support for existing conceptual models of headwater stream disconnection based on spatially variable underflow capacity and shallow groundwater supply.

Virginia↗
USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗
USGS2024

U.S. Geological Survey Karst Interest Group proceedings, Nashville, Tennessee, October 22-24, 2024

Karst hydrogeologic systems represent challenging and unique conditions to scientists studying groundwater flow and contaminant transport. Karst terrains are characterized by distinct and beautiful landscapes, caverns, and springs, and many of the exceptional karst areas are designated as national or state parks. The range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than any other aquifer type. The U.S. Geological Survey (USGS) Karst Interest Group (KIG), formed in 2000, is a loosely knit, grassroots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst aquifers. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. To accomplish its mission, the KIG has organized a series of workshops. To date (2024), nine KIG workshops, including the workshop documented in this report, have been held. The abstracts and extended abstracts provide a snapshot in time of past and current karst related studies. The USGS Water Availability and Use Science Program funded the workshop and proceedings. The planning committee for the ninth workshop includes Thomas D. Byl (USGS and Tennessee State University), Allan K. Clark (USGS), Laura M. DeMott (USGS), Eve L. Kuniansky (USGS, Emeritus), Benjamin V. Miller (USGS), and Lawrence E. Spangler (USGS, Emeritus). The workshop proceedings are edited by Eve L. Kuniansky and Lawrence E. Spangler. The field trip guide was produced by Benjamin V. Miller and Brian Ham (Tennessee Department of Environment and Conservation) and included in the proceedings from the KIG’s 2021 virtual workshop to be used on the optional field trip held on Thursday, October 24, 2024.

Open-File Report↗
USGS2026

Changing drivers of regional large magnitude avalanche frequency throughout Colorado, USA

Large magnitude snow avalanches (destructive size ≥ D3) impact settlements, transportation corridors, and public safety worldwide. In Colorado, United States, avalanches have killed more people than any other natural hazard since 1950. In March 2019, a large magnitude avalanche cycle occurred throughout the entire mountainous portion of Colorado resulting in more than 1000 reported avalanches during a two-week period. Nearly 200 of these avalanches were size D4 or larger with at least three D5 avalanches. However, placing this 2019 large magnitude avalanche cycle in historic context requires data prior to the instrumental record. Here, we paired tree disturbance data from dendrochronology (1698 to 2020) with meteorological data from the modeled and instrumental record (1901 to 2020) to understand the frequency and climate drivers of large magnitude snow avalanche cycles. The extensive number of downed trees from the 2019 avalanche cycle allowed us to collect 1,188 cross-sections and cores from 1023 individual trees within 24 avalanche paths across the state. From these samples we identified 4135 avalanche-related growth disturbances. We employed a strategic nested sampling design to spatially aggregate avalanche frequency from individual avalanche paths, to counties, to three major sub-regions (i.e., north, central, and south), and across the entire region (i.e., state of Colorado). Over a period spanning more than three centuries (1698 to 2020), we identified 76 avalanche years within 24 individual avalanche paths. Large magnitude avalanche event frequency varied across paths and sub-regions with several notable region-wide avalanche cycles. Both tree-ring and historical written records highlighted 1899 as a year with widespread and large magnitude avalanche activity similar to the March 2019 avalanche cycle. Since the early-20th century (1900 to 2020) regional avalanche probability declined significantly in parallel with decreasing snowpack throughout Colorado. Similarly, dominant avalanche regimes shifted from large magnitude regional cycles driven by above average snowfall years over most of the record, to regional avalanche cycles occurring more commonly in average to low snow years since 1988. In recent decades, a lack of December precipitation and above average March precipitation characterized years with regional large magnitude avalanche activity. Even with declining snow water equivalent, truly extreme regional large magnitude avalanche cycles remain possible – as demonstrated by the 2019 cycle. This underscores that rare but high-impact events are not eliminated by long-term trends. Understanding the changing snow and weather drivers and subsequent behavior of large magnitude avalanche cycles across multiple spatial scales may improve avalanche forecasting and the products and mitigations strategies developed by structural engineers to mitigate avalanche danger. This can decrease the avalanche risk to the public and improve infrastructure design in avalanche terrain.

Colorado↗
USGS2005

Stratigraphy and sedimentology of a dry to wet eolian depositional system, Burns formation, Meridiani Planum, Mars

Outcrop exposures of sedimentary rocks at the Opportunity landing site (Meridiani Planum) form a set of genetically related strata defined here informally as the Burns formation. This formation can be subdivided into lower, middle, and upper units which, respectively, represent eolian dune, eolian sand sheet, and mixed eolian sand sheet and interdune facies associations. Collectively, these three units are at least 7 m thick and define a “wetting-upward” succession which records a progressive increase in the influence of groundwater and, ultimately, surface water in controlling primary depositional processes. The Burns lower unit is interpreted as a dry dune field (though grain composition indicates an evaporitic source), whose preserved record of large-scale cross-bedded sandstones indicates either superimposed bedforms of variable size or reactivation of lee-side slip faces by episodic (possibly seasonal) changes in wind direction. The boundary between the lower and middle units is a significant eolian deflation surface. This surface is interpreted to record eolian erosion down to the capillary fringe of the water table, where increased resistance to wind-induced erosion was promoted by increased sediment cohesiveness in the capillary fringe. The overlying Burns middle unit is characterized by fine-scale planar-laminated to low-angle-stratified sandstones. These sandstones accumulated during lateral migration of eolian impact ripples over the flat to gently undulating sand sheet surface. In terrestrial settings, sand sheets may form an intermediate environment between dune fields and interdune or playa surfaces. The contact between the middle and upper units of the Burns formation is interpreted as a diagenetic front, where recrystallization in the phreatic or capillary zones may have occurred. The upper unit of the Burns formation contains a mixture of sand sheet facies and interdune facies. Interdune facies include wavy bedding, irregular lamination with convolute bedding and possible small tepee or salt-ridge structures, and cm-scale festoon cross-lamination indicative of shallow subaqueous flows marked by current velocities of a few tens of cm/s. Most likely, these currents were gravity-driven, possibly unchannelized flows resulting from the flooding of interdune/playa surfaces. However, evidence for lacustrine sedimentation, including mudstones or in situ bottom-growth evaporites, has not been observed so far at Eagle and Endurance craters. Mineralogical and elemental data indicate that the eolian sandstones of the lower and middle units, as well as the subaqueous and eolian deposits of the Burns upper unit, were derived from an evaporitic source. This indirectly points to a temporally equivalent playa where lacustrine evaporites or ground-water-generated efflorescent crusts were deflated to provide a source of sand-sized particles that were entrained to form eolian dunes and sand sheets. This process is responsible for the development of sulfate eolianites at White Sands, New Mexico, and could have provided a prolific flux of sulfate sediment at Meridiani. Though evidence for surface water in the Burns formation is mostly limited to the upper unit, the associated sulfate eolianites provide strong evidence for the critical role of groundwater in controlling sediment production and stratigraphic architecture throughout the formation.

Earth and Planetary Science Letters↗
USGS2006

Submarine landslides in the Santa Barbara Channel as potential tsunami sources

Recent investigations using the Monterey Bay Aquarium Research Institutes (MBARI) Remotely Operated Vehicles (ROVs) "Ventana" and "Tiburon" and interpretation of MBARI's EM 300 30 kHz multibeam bathymetric data show that the northern flank of the Santa Barbara Basin has experienced massive slope failures. Of particular concern is the large (130 km2) Goleta landslide complex located off Coal Oil Point near the town of Goleta, that measures 14.6-km long extending from a depth of 90 m to nearly 574 m deep and is 10.5 km wide. We estimate that approximately 1.75 km3 has been displaced by this slide during the Holocene. This feature is a complex compound submarine landslide that contains both surfical slump blocks and mud flows in three distinct segments. Each segment is composed of a distinct head scarp, down-dropped head block and a slide debris lobe. The debris lobes exhibit hummocky topography in the central areas that appear to result from compression during down slope movement. The toes of the western and eastern lobes are well defined in the multibeam image, whereas the toe of the central lobe is less distinct. Continuous seismic reflection profiles show that many buried slide debris lobes exist and comparison of the deformed reflectors with ODP Drill Site 149, Hole 893 suggest that at least 200 000 years of failure have occurred in the area (Fisher et al., 2005a). Based on our interpretation of the multibeam bathymetry and seismic reflection profiles we modeled the potential tsunami that may have been produced from one of the three surfical lobes of the Goleta slide. This model shows that a 10 m high wave could have run ashore along the cliffs of the Goleta shoreline. Several other smaller (2 km2 and 4 km2) slides are located on the northern flank of the Santa Barbara Basin, both to the west and east of Goleta slide and on the Concepcion fan along the western flank of the basin. One slide, named the Gaviota slide, is 3.8 km2, 2.6 km long and 1.7 km wide. A distinct narrow scar extends from near the eastern head wall of this slide for over 2 km eastward toward the Goleta slide and may represent either an incipient failure or a remnant of a previous failure. Push cores collected within the main head scar of this slide consisted of hydrogen sulfide bearing mud, possibly suggesting active fluid seepage and a vibra-core penetrated ???50 cm of recent sediment overlying colluvium or landslide debris confirming the age of ???300 years as proposed by Lee et al. (2004). However, no seeps or indications of recent movement were observed during our ROV investigation within this narrow head scar indicating that seafloor in the scar is draped with mud. ?? 2006 Author(s). This work is licensed under a Creative Commons License.

Natural Hazards and Earth System Sciences↗
USGS2025

Calibration of the Stream Salmonid Simulator (S3) model to estimate annual survival, movement, and food consumption by juvenile Chinook salmon (Oncorhynchus tshawytscha) in the restoration reach of the Trinity River, California, 2006–18

Executive Summary The Trinity River is managed in two sections: (1) from the upper 64-kilometer “restoration reach” downstream from Lewiston Dam to the confluence with the North Fork Trinity River, and (2) the 120-kilometer lower Trinity River downstream from the restoration reach. The Stream Salmonid Simulator (S3) has been previously applied to these reaches and the Klamath River. To estimate fish growth, past S3 calibration efforts in the Trinity and Klamath Rivers used maximum likelihood methods that considered only the abundance of juvenile Chinook salmon ( Oncorhynchus tshawytscha ) passing a fish trap to estimate survival and movement parameters, but not fish consumption. Previous calibrations did not estimate the average proportion of maximum consumption ( C y ) when estimating survival ( S y ) and movement ( M 0 y ) parameters across years ( y ) of data, but because no other information was available in the literature a fixed value of C y =0.66 was assumed. Therefore, the goal of this report is to present an alternative approach that calibrates the S3 model to multivariate data (that is, abundance and size), enabling the estimation of the average proportion of maximum consumption, in conjunction with survival and movement parameters for a particular migration year. We fit the S3 model to individual years of weekly trap abundance estimates and mean fish sizes (fork length) at the Pear Tree Gulch (hereafter referred to as Pear Tree) fish trap representing the restoration reach. We used the Earth Mover’s Distance (EMD) as the objective value to be minimized in parameter optimization. This approach estimated survival, movement, and consumption parameters for each migration year. Because we had information on the abundance of natural and hatchery produced juvenile salmon at the fish traps, we estimated survival and movement for natural and hatchery fish. S3 is a deterministic life-stage-structured population model that tracks daily growth, movement, and survival of juvenile Chinook Salmon. A key theme of the model is that river discharge affects habitat availability and capacity, which in turn drives density-dependent population dynamics. To explicitly link population dynamics to habitat quality and quantity, the river environment is constructed as a one-dimensional series of linked habitat units, each of which has an associated daily timeseries of discharge, water temperature, and useable habitat area or carrying capacity. In turn, the physical characteristics of each habitat unit and the number of fish occupying each unit drive survival and growth within each habitat unit and movement of fish among habitat units. The physical template of the restoration reach of the Trinity River was classified into 356 meso-habitat units comprised of runs, riffles, and pools. For each habitat unit, we developed a timeseries of daily discharge, water temperature, amount of available spawning habitat, and fry and parr carrying capacity. Capacity time series were constructed using state-of-the-art models of spatially explicit hydrodynamics and quantitative fish habitat relationships developed for the Trinity River. These variables were then used to drive population dynamics such as egg maturation and survival, and in turn, juvenile movement, growth, and survival. We estimated movement, survival, and consumption parameters by calibrating the model to 12 years of weekly juvenile abundance estimates and fish sizes at the Pear Tree fish trap near the downstream end of the restoration reach. We estimated parameters for 12 years that included a wide range of female spawner abundances (1,414–11,494) and water year types (critically dry–extremely wet). We contrast the estimated parameters to the corresponding number of female spawners and the total annual volume of water discharged for the Trinity River (Trinity River Restoration Program; https://www.trrp.net/restoration/flows/summary/ ). The calibration consisted of replicating historical conditions as closely as possible (for example, discharge; temperature; spawner abundance, spawning location and timing, and hatchery releases), and then running the model to predict weekly abundance passing the trap location from each brood year of adults and subsequent migration year of their juvenile progeny. Because density-dependent movement was favored in past evaluations, we estimated S3 parameters based on density-independent survival and density-dependent movement. Likewise, each year’s estimated survival parameter for natural ( S N y ) and hatchery ( S H y ) fish may be interpreted as the mean daily survival probability from emergence or hatchery release to the Pear Tree fish trap. Under density dependence, the estimated movement parameter for natural ( M 0N y ) and hatchery ( M 0H y ) fish represents the intercept of the Beverton-Holt model; the probability of remaining in a habitat at near-zero abundance. We estimated C y by using EMD and incorporating abundance and fish size into model calibration. Average daily proportions of maximum consumption, , across the years were generally high (=0.640; standard deviation (SD) SD=0.176), suggesting that fish were feeding at about two-thirds of expected maximum consumption rates. This average proportion of maximum consumption,is very similar to what has been assumed (=0.66) in previous Trinity and Klamath River S3 calibration and simulation efforts. In 2017, we estimated the lowest C y , suggesting lower average consumption for juvenile salmon in high-discharge water years. When this high discharge year was excluded, there was no apparent trend in C y with annual water volume. Estimates of survival showed little trend over the range in spawner abundances, but a trend towards higher natural and hatchery fish survival with higher annual volumes of water was apparent. Over the 12 years, the average survival of hatchery fish was =0.888 (SD=0.079) and the average survival natural fish was=0.969 (SD=0.01). With respect to fish movement, we estimated higher M 0N y and M 0H y with higher annual volumes of water in the Trinity River. Higher M N0 y or M H0 y suggest greater probability of remaining in a habitat at low fish densities, with potential for density-dependent processes in movement to occur. The highest M 0N y = 0.676 was estimated during brood year 2012, and the overall average for natural fish was =0.276 (SD=0.188) and for hatchery fish was=0.467 (SD=0.235). Under the Beverton-Holt model, as M 0N y or M 0H y approach zero, there is less capacity for change in fish movement as fish density increases. The S3 model was initialized with only the spatiotemporal distribution of spawners, so it performed well at capturing the essential outmigration features that are ultimately governed by rates of growth, movement, and mortality. We used a new optimization method that could accommodate multivariate data on abundance and fish size collected at the Pear Tree fish trap, enabling the calibration of S3 to estimate five parameters for 12 separate years of data. Incorporating weekly fish size data for each year in our parameter optimization process made the estimation of C y possible and represents a step forward in the fitting of the S3 model to fish trap data for the purposes of parameter calibration and the estimation of growth parameters with respect to annual conditions. We identified lack of fit and adding important effects into the S3 model may improve the S3 estimation and simulation of water scenarios. The Trinity River Restoration Program (TRRP) Science Advisory Board recommended that the TRRP focus on developing core elements of a decision support system (DSS; Buffington and others, 2014). Toward that end, the habitat and S3 models described in this report are both core elements of the DSS. The structure of S3 makes it a particularly useful fish production model for the DSS because population dynamics are sensitive to (1) water temperature, (2) daily discharge management, and (3) habitat quality and quantity. Each of these variables are key management parameters under consideration in the TRRP. As such, the S3 model may provide valuable insights into the potentially variable effects of different management decisions on the Trinity River.

California↗