USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth and Space Science”

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

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

At least 361 records · Page 20Linked to original sources
USGS2000

Preliminary report on the 16 October 1999 M 7.1 Hector mine, California, earthquake

The M w 7.1 Hector Mine, California, earthquake occurred at 9:46 GMT on 16 October 1999. The event caused minimal damage because it was located in a remote, sparsely populated part of the Mojave Desert, approximately 47 miles east-southeast of Barstow, with epicentral coordinates 34.59°N 116.27°W and a hypocentral depth of 5 ± 3 km. Twelve foreshocks, M 1.9-3.8, preceded the mainshock during the previous twelve hours. All of these events were located close to the hypocenter of the mainshock. The Hector Mine earthquake occurred within the Eastern California Shear Zone (ECSZ). By virtue of its remote location, the societal impact of the Hector Mine earthquake was, fortunately, minimal in spite of the event's appreciable size. The ECSZ is characterized by high seismicity, a high tectonic strain rate, and a broad, distributed zone of north-northwest-trending faults (ECSZ; Figure 1 ; Dokka and Travis, 1990 ; Sauber et al., 1986 ; Sauber et al., 1994 ; Sieh et al., 1993 ). Data regarding the slip rates of faults within the ECSZ suggest that on the order of 15% of the Pacific-North American plate motion occurs along this zone ( Sauber et al., 1986 ; Wesnousky, 1986 ). Most of the faults in the ECSZ have low slip rates and long repeat times for major earthquakes, on the order of several thousands to tens of thousands of years. The occurrence of the Hector Mine earthquake within seven years and only about 30 km east of the 1992 M w 7.3 Landers earthquake suggests that the closely spaced surface faults in the ECSZ are mechanically related. The Hector Mine event involved rupture on two previously mapped fault zones—the Bullion Fault and an unnamed, more northerly-trending fault that is informally referred to in this paper as the Lavic Lake Fault (Dibblee, 1966 , 1967a , b ). Traces of the Bullion Fault exhibit evidence of Holocene displacement and were zoned as active in 1988 under California's Mquist-Priolo Earthquake Fault Zoning Act ( Hart and Bryant, 1997 ). The pattern of rupture along more than one named fault was also observed from the 1992 Landers earthquake ( Hauksson et al., 1993 ; Sieh et al., 1994). Much of the fault zone that produced the Hector Mine earthquake had been buried by relatively young stream deposits, and the fault scarps in bedrock have a subdued morphology. It appears that these faults have not experienced significant offset for perhaps 10,000 years or more ( Hart, 1987 ). Planned future investigations will refine the age of the last event on these faults. The portion of the Lavic Lake Fault that ruptured between the northern end of the Bullion Mountains and Lavic Lake had not previously been mapped. However, our field investigations have identified ancient, subdued fault scarps along portions of the 1999 rupture zone in this area. It thus appears that the entire segment of the Lavic Lake Fault that was involved in the 1999 event had ruptured in the past. As is typical for most faults within the Eastern California Shear Zone, the rate of movement along the Lavic Lake Fault may be quite slow (<1 mm/yr) and should produce earthquakes only infrequently. This event is a reminder that faults that have ruptured in late Quaternary time, but that lack evidence of Holocene displacement, can still produce earthquakes in this low-slip-rate tectonic setting. Additionally, the Hector Mine earthquake is noteworthy for a couple of other reasons. First, it clearly produced triggered seismicity over much of southern California, from the rupture zone toward the south-southwest in particular. Second, as we will discuss, the event may provide new data and insight into recently developed paradigms concerning earthquake interactions and the role of static stress changes. Questions such as these will, of course, be the subject of extensive detailed analyses in years to come. Fortunately, the Hector Mine sequence will provide one of the best data sets obtained to date for a significant earthquake in the United States. Because it occurred when major upgrades to both the regional seismic network (TriNet) and the regional geodetic network (SCIGN) were well underway, the Earth science community will have abundant high-quality data with which to explore the important and interesting questions that have been raised. In this paper, we present and discuss the basic data and preliminary results from the Hector Mine earthquake.

California↗
USGS2013

A new model for the growth of basaltic shields based on deformation of Fernandina volcano, Galápagos Islands

Space-geodetic measurements of surface deformation produced by the most recent eruptions at Fernandina – the most frequently erupting volcano in the Galápagos Archipelago – reveal that all have initiated with the intrusion of subhorizontal sills from a shallow magma reservoir. This includes eruptions from fissures that are oriented both radially and circumferentially with respect to the summit caldera. A Synthetic Aperture Radar (SAR) image acquired 1–2 h before the start of a radial fissure eruption in 2009 captures one of these sills in the midst of its propagation toward the surface. Galápagos eruptive fissures of all orientations have previously been presumed to be fed by vertical dikes, and this assumption has guided models of the origin of the eruptive fissure geometry and overall development of the volcanoes. Our findings allow us to reinterpret the internal structure and evolution of Galápagos volcanoes and of similar basaltic shields. Furthermore, we note that stress changes generated by the emplacement of subhorizontal sills feeding one type of eruption may control the geometry of subsequent eruptive fissures. Specifically, circumferential fissures tend to open within areas uplifted by sill intrusions that initiated previous radial fissure eruptions. This mechanism provides a possible explanation for the pattern of eruptive fissures that characterizes all the western Galápagos volcanoes, as well as the alternation between radial and circumferential fissure eruptions at Fernandina. The same model suggests that the next eruption of Fernandina will be from a circumferential fissure in the area uplifted by the 2009 sill intrusion, just southwest of the caldera rim.

Galápagos Islands↗
USGS2023

Geochemical evidence for evolving Proterozoic crustal thickness and orogenic styles in southwestern Laurentia

It has long been challenging for researchers to track the crustal thickness and mode(s) of crustal modification in ancient convergent margins, limiting evaluation of the tectonic styles and processes that modify continental crust during orogenesis. We present trace element igneous geochemical crustal thickness proxies that quantitatively track the crustal thickness evolution of the long-lived Proterozoic active margin in the southwestern U.S.A. We integrate these results with geobarometric data to constrain the mode of crustal modification. The data indicate a complex record of crustal thickness change in space and time and evolving orogenic styles. Geochemical proxies at 1.84–1.72 Ga are consistent with 20–40 km thick magmatic arcs that were locally thickened to ∼50 km during ∼1.75 Ga tectonism. During the Yavapai orogeny, 1.72–1.69 Ga, a ∼200-km-wide belt of 50-60 km thick crust extended from southern California to northern Colorado and was rapidly thinned and exhumed by ∼1.68 Ga. Crustal thickening and thinning during the Yavapai orogeny largely occurred by shortening and exhumation, respectively, in the upper 25 km of the crust. Subsequent 1.68–1.60 Ga tectonism involved crustal growth, local crustal thickening, and low-P, high-T metamorphism, consistent with extensional accretionary orogenesis. The 1.47–1.37 Ga Picuris orogeny was associated with 50–60 km thick crust across much of the Southwest and involved crustal shortening with ∼10 km of magmatic underplating. Advective heat from the emplacement of ferroan granites in the mid-crust likely contributed to elevated geothermal gradients and rheologically weakened the crust. Our results suggest evolving orogenic styles in the Southwest from 1.75–1.69 Ga short-lived crustal thickening associated with terrane accretion to 1.69–1.60 Ga largely extensional accretionary orogenesis, and regional, long-lived crustal thickening at 1.47–1.37 Ga involving extensive basaltic underplating. Contrasting with some recent hypotheses, our data document a complex middle Proterozoic record for the Southwest that was not orogenically quiescent or tectonically stagnant but involved complex mountain building styles.

Earth and Planetary Science Letters↗
USGS2023

Using simulated GEDI waveforms to evaluate the effects of beam sensitivity and terrain slope on GEDI L2A relative height metrics over the Brazilian Amazon Forest

The vertical structure of forests provides important parameters for estimating aboveground biomass (AGB) and it can be measured by LiDAR sensors. The Global Ecosystem Dynamics Investigation (GEDI) full-waveform LiDAR sensor collects data systematically over the Earth’s surface from the International Space Station. Since GEDI became operational, it has collected billions of ~25 m diameter footprints. This massive dataset has been used to create higher level gridded and non-grided products. However, GEDI’s ~25 m footprints can be subject to errors associated with effects of geolocation, terrain slope, and beam sensitivity, among others, which are likely transferred to the downstream products. This study aims to (1) evaluate the effect of beam sensitivity and terrain slope on the accuracy of relative heights (RH) of GEDI product L2A version 2 through comparison with discrete-return airborne LiDAR data collected over transects in the Brazilian Amazon Forest biome, (2) assess GEDI’s geolocation uncertainty and investigate its combined effects with beam sensitivity and terrain slope, and (3) re-evaluate beam sensitivity and terrain slope effects on the GEDI L2A RHs using footprints that were geolocation-adjusted through a simple novel approach. The analysis separates GEDI footprints by acquisition time, i.e., daytime, nighttime, and combined (all-data). The discrete-return airborne LiDAR point clouds are used to derive terrain slope within the GEDI footprints, and to simulate GEDI waveforms and derive RHs comparable to the GEDI L2A product. Results indicate that terrain slope only causes significant effects on GEDI data collected during daytime because solar radiation affects waveform signal-to-noise ratio. Beam sensitivity causes significant effects on nighttime and all-data GEDI L2A RHs. If geolocation uncertainty is considered, the effects of beam sensitivity and terrain slope have just minor changes. Geolocation-adjusted data continue showing significant effects on nighttime and all-data RH differences caused by beam sensitivity but produce unbiased results. This study improves the understanding of how beam sensitivity, terrain slope and their combined effect with geolocation uncertainty may affect the GEDI L2A RH collected over the Brazilian Amazon forest during daytime and nighttime.

Brazilian Amazon Forest↗
USGS2023

Summary of the history and research of the U.S. Geological Survey gas hydrate properties laboratory in Menlo Park, California, active from 1993 to 2022

The U.S. Geological Survey (USGS) Clathrate Hydrate Properties Project, active from 1993 to 2022 in Menlo Park, California, stemmed from an earlier project on the properties of planetary ices supported by the National Aeronautics and Space Administration’s (NASA’s) Planetary Geology and Geophysics Program. We took a material science approach in both projects, emphasizing chemical purity of samples, having controlled grain size and grain texture, and having verified crystal structures and phase relations. A foundational contribution from our USGS Gas Hydrate Properties Laboratory (GHPL) was in demonstrating the ability to reproducibly create such pure clathrate hydrate samples for study. Clathrate sample synthesis was achieved by heating sieved and weighed pure granular water ice in the presence of cold clathrate-forming gas or liquid. During heating, the ice melts at the grain scale and reacts with the gas to form clathrate. The resulting material has the desired uniformity and purity, with known intergranular porosity; our subsequent measurements showed that these clathrates exhibited the established clathrate structures and phase relations. This novel synthesis method was successful in creating clathrates of pure methane, ethane, propane, carbon dioxide, and multi-component gases. By mixing sand or silt with granular ice, we were also able to make clathrate-sediment aggregates with controlled grain textures. This simple method, adopted by many others in the community, permitted us to measure the physical and chemical properties of well-characterized and well-crystallized clathrates and clathrate/sediment aggregates. At about the same time, we adapted conventional scanning electron microscopy to cryogenic conditions for analysis of grain-scale characteristics of clathrates made in the GHPL as well as those collected from nature by drill core. The uniformity and reproducibility of our samples also allowed us to investigate how clathrates respond to environmental changes in chemistry, temperature, and pressure: we measured chemical exchange rates with dissolved gas species—such as noble gases and chlorofluorocarbons—as well as rates of clathrate dissolution and decomposition. These advances include the first accurate mapping of the conditions that promote the remarkable process of “anomalous preservation” at room pressure, a metastability that offers potential application for low-cost and safe transportation of natural gas from gas fields far from pipelines. Another advancement stemming from the GHPL was the compaction of as-synthesized porous clathrates to nearly full density by applying external pressure using three different techniques. Compaction allows for high-accuracy measurements of many fundamental physical and chemical properties of these materials, such as elastic wavespeeds and moduli, complete thermal properties, decomposition rates, thermal expansion, and clathrate equations of state. These properties and others, in turn, have helped USGS scientists to interpret geophysical well logs and active geophysical surveys, as well as model the rates of gas production from hydrate deposits in nature. Studying this class of icy minerals that occur in abundance on Earth and in the outer solar system has been a fascinating laboratory journey. Here, we summarize the history and major findings of the USGS GHPL in Menlo Park, including both in-house research as well as findings from the synergistic collaborations with other agencies and institutes that were key to the success of our laboratory. The Menlo Park GHPL was more formally incorporated within the USGS Gas Hydrates Project, a collaboration among multiple USGS Science Centers, in the early 2000s under the leadership of Deborah Hutchinson, and now under the leadership of Carolyn Ruppel and Timothy Collett.

Open-File Report↗
USGS2023

ENSO and NAO linkages to interannual salinity variability in north central Gulf of Mexico estuaries through teleconnections with precipitation

Though the importance of Earth's internal climate modes such as the El Niño-Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO) to regional-scale climate variability is well recognized, the degree to which these oscillations are reflected by spatio-temporal salinity variability over interannual timescales in estuaries is less understood. Here an 11-year continuous salinity monitoring dataset spanning 223 stations across Louisiana's coastal wetlands along the northern Gulf of Mexico is examined with empirical orthogonal function (EOF) analysis to identify dominant modes of interannual variability in the salinity field. The first EOF mode accounts for 72% of the variance in the salinity field and captures a domain-wide pattern where salinities vary in-phase through space in response to local precipitation anomalies occurring in the vicinity of the study area. This local precipitation anomaly is positively correlated with ENSO (Nino3.4 index), consistent with the El Niño – wet (La Niña – dry) precipitation teleconnection that is prevalent throughout the northern Gulf of Mexico coast. The second EOF mode, which accounts for 13% of the variance in the salinity field, is expressed primarily in the marshes across the lower reaches of the Mississippi River deltaic plain (MRDP). EOF2 is anticorrelated with annual Mississippi River discharge anomaly such that salinities in the lower MRDP decrease as discharge increases, pointing to enhanced advection of fresh river plume waters over the shelf into the estuary via estuary-ocean exchange during years of anomalously high river discharge. Mississippi River discharge anomaly is positively correlated with the NAO at a one-year time lag, through a teleconnection with precipitation throughout much of the central region of the Mississippi River drainage basin. Together, these findings indicate that most of the interannual salinity variability across Louisiana's coastal wetlands can be linked to climate variability through teleconnections with precipitation. Incorporating these dynamics into restoration planning, monitoring, and adaptive management efforts may help constrain background environmental variation and better isolate restoration effects.

Louisiana↗
USGS1996

Channels and valleys on Mars: Cold climate features formed as a result of a thickening cryosphere

Large flood channels, valley networks, and a variety of features attributed to the action of ground ice indicate that Mars emerged from heavy bombardment around 3.8Gyr ago, with an inventory of water at the surface equivalent to at least a few hundred meters spread over the whole planet, as compared with 3 km for the Earth. The surface water resided primarily in a porous, kilometers thick, megaregolith created by the high impact rates. At the end of heavy bombardment a rapid decline in erosion rates by a factor of 1000 suggests a major change in the global climate. It is proposed that at this time the climate became similar to today's and that this climate has been maintained throughout the rest of Mars' history. The various drainage features represent an adjustment of the distribution of water to the surface relief inherited from the period of heavy bombardment and to a thickening of the cryosphere as the heat flow declined. The valley networks formed mostly at the end of heavy bombardment when erosion rates were high and climatic conditions permitted an active water cycle. They continued to form after heavy bombardment when the cryosphere started to form by a combination of episodic flooding and mass-wasting aided by the presence of liquid water at shallow depths. As the cryosphere thickened with declining heat flow, water could no longer easily access the surface and the rate of valley formation declined. Hydrostatic pressures built below the cryosphere. Eruptions of groundwater became more catastrophic and massive floods resulted, mainly in upper Hesperian time. Flood sources were preferentially located in low-lying, low-latitude areas where the cryosphere was thin, or near volcanoes where a thinner than typical cryosphere is also expected. Floods caused a drawdown in the global water table so that few formed in the second half of Mars' history. The floodwaters pooled in low-lying areas, mostly in the northern plains. Some of the water may still be present as thick ice deposits, some has been lost to space, particularly during periods of high obliquity.

Planetary and Space Science↗
USGS2013

Current research issues related to post-wildfire runoff and erosion processes

Research into post-wildfire effects began in the United States more than 70 years ago and only later extended to other parts of the world. Post-wildfire responses are typically transient, episodic, variable in space and time, dependent on thresholds, and involve multiple processes measured by different methods. These characteristics tend to hinder research progress, but the large empirical knowledge base amassed in different regions of the world suggests that it should now be possible to synthesize the data and make a substantial improvement in the understanding of post-wildfire runoff and erosion response. Thus, it is important to identify and prioritize the research issues related to post-wildfire runoff and erosion. Priority research issues are the need to: (1) organize and synthesize similarities and differences in post-wildfire responses between different fire-prone regions of the world in order to determine common patterns and generalities that can explain cause and effect relations; (2) identify and quantify functional relations between metrics of fire effects and soil hydraulic properties that will better represent the dynamic and transient conditions after a wildfire; (3) determine the interaction between burned landscapes and temporally and spatially variable meso-scale precipitation, which is often the primary driver of post-wildfire runoff and erosion responses; (4) determine functional relations between precipitation, basin morphology, runoff connectivity, contributing area, surface roughness, depression storage, and soil characteristics required to predict the timing, magnitudes, and duration of floods and debris flows from ungaged burned basins; and (5) develop standard measurement methods that will ensure the collection of uniform and comparable runoff and erosion data. Resolution of these issues will help to improve conceptual and computer models of post-wildfire runoff and erosion processes.

Earth-Science Reviews↗
USGS1977

Voyager imaging experiment

The overall objective of this experiment is exploratory reconnaissance of Jupiter, Saturn, their satellites, and Saturn's rings. Such reconnaissance, at resolutions and phase angles unobtainable from Earth, can be expected to provide much new data relevant to the atmospheric and/or surface properties of these bodies. The experiment also has the following specific objectives: Observe and characterize the global circulation of the atmospheres of Jupiter and Saturn; Determine the horizontal and vertical structure of the visible clouds and establish their relationship to the belted appearance and dynamical properties of the planetary atmospheres; Determine the vertical structure of high, optically-thin, scattering layers on Jupiter and Saturn; Determine the nature of anomalous features such as the Great Red Spot, South Equatorial Belt disturbances, etc.; Characterize the nature of the colored material in the clouds of Jupiter and Saturn, and identify the nature and sources of chromophores on Io and Titan; Perform comparative geologic studies of many satellites at less than 15-km resolution; Map and characterize the geologic structure of several satellites at high resolution (???1 km); Investigate the existence and nature of atmospheres on the satellites; Determine the mass, size, and shape of many of the satellites by direct measurement; Determine the direction of the spin axes and periods of rotation of several satellites, and establish coordinate systems for the larger satellites; Map the radial distribution of material in Saturn's rings at high resolution; Determine the optical scattering properties of the primaries, rings, and satellites at several wavelengths and phase angles; Search for novel physical phenomena, e.g., phenomena associated with the Io flux tube, meteors, aurorae, lightning, or satellite shadows. ?? 1977 D. Reidel Publishing Company.

Space Science Reviews↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2022

Science mission requirements for a globally ranging, riserless drilling vessel for U.S. Scientific Ocean Drilling

Through the collection and analysis of shallow and deep subseafloor sediments, rocks, fluids, and life, scientific ocean drilling has enriched our understanding of the complex Earth system. Among other achievements, scientific ocean drilling has documented the history of Earth’s climate, the waxing and waning of polar ice sheets, the past changes in ocean and atmospheric circulation, the existence and function of microbial life in the subseafloor, the compositional variations in Earth’s crust and underlying mantle, and the physical and chemical processes acting at subduction zones, including those associated with tsunamigenic earthquakes. Over the decades, more than 12,000 articles that depend on analyses of scientific ocean drilling samples and geophysical data have been published, many detailing breakthrough contributions to global knowledge about the Earth system. Approximately 45% of these publications were led by U.S.-affiliated authors (International Ocean Discovery Program Publication Services, 2021). Since the mid-1980s, the workhorse of this multidisciplinary, international research effort has been the riserless D/V JOIDES Resolution, operated by Texas A&M University with funding from the U.S. National Science Foundation (NSF). D/V JOIDES Resolution has conducted the vast majority of scientific ocean drilling expeditions and collected most of the scientific cores over that period, including 82% of the expeditions and 93% of the cores in the last decade alone, despite being one of three platforms that is operated within the International Ocean Discovery Program. However, D/V JOIDES Resolution is approaching the end of its useful life. With a strong commitment to continue scientific ocean drilling beyond the end of the current phase, the community developed a document outlining the research frontiers that should be pursued. Exploring Earth by Scientific Ocean Drilling: 2050 Science Framework (Koppers and Coggon, 2020) describes seven scientific strategic objectives that focus on understanding interconnections within the Earth system and five flagship initiatives that integrate these objectives into long-term research efforts that address issues facing society. Additional elements in the 2050 Science Framework, including STEM education, workforce development, technology development, and innovative applications of data analytics, will advance the goals of scientific ocean drilling. Addressing the 2050 Science Framework also requires building partnerships with allied U.S. and international science programs and strengthening existing ones. To implement a significant portion of the 2050 Science Framework, the U.S. scientific community seeks to lease or acquire a newly built, globally ranging, state-of-the art, riserless drilling vessel. The many and varied technical and human resources requirements for successful accomplishment of scientific and educational goals summarized in this document and described in detail in the 2050 Science Framework require broad community input and careful consideration. Following receipt of NSF’s formal Request for Assistance to the United States Science Support Program (USSSP), the U.S. scientific ocean drilling community conducted a one-year exercise to identify its national scientific needs and priorities in order to determine the Science Mission Requirements (SMRs) presented here. This community effort included: (1) a U.S. community-wide survey to identify the specific operational and technical capabilities critical to addressing science in the 2050 Science Framework; (2) a series of online workshops focusing on critical capabilities identified by the survey; and (3) a large in-person workshop to synthesize the results of the survey and the virtual workshops (Appendix 1). The approach was designed to reach as many participants as possible. Overall, 278 survey responses were received from U.S. community members, representing 104 unique institutions from 39 states and the District of Columbia, and 137 unique individuals participated in the workshops (Appendix 2). The results of this effort comprise two classes of SMRs: Foundational Science Mission Requirements and Primary Science Mission Requirements. Foundational SMRs define minimum criteria for a new riserless drilling vessel that can address significant portions of the 2050 Science Framework. Primary SMRs build upon the Foundational SMRs and will create more robust science opportunities and data collection capabilities, will increase progress in addressing the 2050 Science Framework objectives, and will provide more real-time ship-to-shore interaction to improve science productivity, engagement, and outreach. Modern safety and environmental standards, including meeting standards to access protected waters such as exclusive economic zones, extended continental shelves, or high latitudes, while being cognizant of the vessel’s environmental footprint. Safe and efficient operations in global locations and in water depths from 70 m to 6000+ m, with total drill string length of at least 7000 m. High-quality core and data collection from a range of key subseafloor environments. Advanced heave compensation, dynamic positioning, and drill pipe stability. Modern mud and cement/casing systems. Critical onboard measurements for safety, operational decision-making, documentation of ephemeral properties, mission-specific science, and long-term science goals that extend beyond a single expedition. Designated and appropriate space for sample and data preservation. Highly skilled onboard personnel, including technical staff for curation and core handling; support for safety, time-sensitive, and critical shipboard measurements; computer support; equipment and instrument repair; application support; and data assurance. Primary Science Mission Requirements include: Flexible shipboard space for laboratories and on deck to ensure safe, successful implementation of diverse science objectives and operations. Minimizing contamination of recovered samples. Over-the-side capabilities for science-supporting technology (e.g., remotely operated vehicles, water- column sampling, sediment-water interface sampling). Downhole logging and measurements. Expanded borehole observatory capabilities. Reliable and consistent ship-to-shore communications. NSF’s investment in a new globally ranging, riserless drilling vessel will have a powerful economic multiplier effect, including the infusion of additional science support funds in the United States for training and research, the development of new technologies and tools, and the associated scientific and technical workforce development. The skills and knowledge gained through scientific ocean drilling are translatable to careers in fields such as sustainable energy development (e.g., geothermal and offshore wind), carbon sequestration, data management and cyberinfrastructure, biotechnology, communications, science education, policy, hazard mitigation, and environmental management. The United States is a leader in a well-established and internationally collaborative scientific ocean drilling community. A modern, globally ranging, riserless drilling vessel will allow the United States to expand its leadership position, address broad scientific questions that current capabilities preclude, and cultivate equitable international, multidisciplinary collaborations that will ensure scientific ocean drilling’s future success.

Report↗
USGS2021

The making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

Frontiers in Earth Science↗
USGS2019

Size distributions of Arctic waterbodies reveal consistent relations in their statistical moments in space and time

Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km 2 (100 m 2 ) to 1 km 2 . We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size ( R 2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.

Alaska↗
USGS2016

Large along-strike variations in the onset of Subandean exhumation: Implications for Central Andean orogenic growth

Plate tectonics drives mountain building in general, but the space-time pattern and style of deformation is influenced by how climate, geodynamics, and basement structure modify the orogenic wedge. Growth of the Subandean thrust belt, which lies at the boundary between the arid, high-elevation Central Andean Plateau and its humid, low-elevation eastern foreland, figures prominently into debates of orogenic wedge evolution. We integrate new apatite and zircon (U-Th)/He thermochronometer data with previously published apatite fission-track data from samples collected along four Subandean structural cross-sections in Bolivia between 15° and 20°S. We interpret cooling ages vs. structural depth to indicate the onset of Subandean exhumation and signify the forward propagation of deformation. We find that Subandean growth is diachronous south (11 ± 3 Ma) vs. north (6 ± 2 Ma) of the Bolivian orocline and that Subandean exhumation magnitudes vary by more than a factor of two. Similar north-south contrasts are present in foreland deposition, hinterland erosion, and paleoclimate; these observations both corroborate diachronous orogenic growth and illuminate potential propagation mechanisms. Of particular interest is an abrupt shift to cooler, more arid conditions in the Altiplano hinterland that is diachronous in southern Bolivia (16-13 Ma) vs. northern Bolivia (10-7 Ma) and precedes the timing of Subandean propagation in each region. Others have interpreted the paleoclimate shift to reflect either rapid surface uplift due to lithosphere removal or an abrupt change in climate dynamics once orographic threshold elevations were exceeded. These mechanisms are not mutually exclusive and both would drive forward propagation of the orogenic wedge by augmenting the hinterland backstop, either through surface uplift or spatially variable erosion. In summary, we suggest that diachronous Subandean exhumation was driven by piecemeal hinterland uplift, orography, and the outward propagation of deformation.

Andes↗
USGS2008

EAARL coastal topography - Northern Gulf of Mexico

These remotely sensed, geographically referenced elevation measurements of Lidar-derived coastal topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. One objective of this research is to create techniques to survey areas for the purposes of geomorphic change studies following major storm events. The USGS Coastal and Marine Geology Program's National Assessment of Coastal Change Hazards project is a multi-year undertaking to identify and quantify the vulnerability of U.S. shorelines to coastal change hazards such as effects of severe storms, sea-level rise, and shoreline erosion and retreat. Airborne Lidar surveys conducted during periods of calm weather are compared to surveys collected following extreme storms in order to quantify the resulting coastal change. Other applications of high-resolution topography include habitat mapping, ecological monitoring, volumetric change detection, and event assessment. The purpose of this project is to provide highly detailed and accurate datasets of the northern Gulf of Mexico coastal areas, acquired on September 19, 2004, immediately following Hurricane Ivan. The datasets are made available for use as a management tool to research scientists and natural resource managers. An innovative airborne Lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Airborne Advanced Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532 nanometer) Lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive Lidar, a down-looking RGB (red-green-blue) digital camera, a high-resolution multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers and an integrated miniature digital inertial measurement unit which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a Lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system on September 19, 2004. The survey resulted in the acquisition of 3.2 gigabytes of data. The data were processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of Lidar data in an interactive or batch mode. Modules for pre-survey flight line definition, flight path plotting, Lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is routinely used to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of 'last return' elevations.

Alabama, Florida, Mississippi↗
USGS2018

Wetlands inform how climate extremes influence surface water expansion and contraction

Effective monitoring and prediction of flood and drought events requires an improved understanding of how and why surface water expansion and contraction in response to climate varies across space. This paper sought to (1) quantify how interannual patterns of surface water expansion and contraction vary spatially across the Prairie Pothole Region (PPR) and adjacent Northern Prairie (NP) in the United States, and (2) explore how landscape characteristics influence the relationship between climate inputs and surface water dynamics. Due to differences in glacial history, the PPR and NP show distinct patterns in regards to drainage development and wetland density, together providing a diversity of conditions to examine surface water dynamics. We used Landsat imagery to characterize variability in surface water extent across 11 Landsat path/rows representing the PPR and NP (images spanned 1985–2015). The PPR not only experienced a 2.6-fold greater surface water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change in surface water extent between drought and deluge conditions. The relationship between surface water extent and accumulated water availability (precipitation minus potential evapotranspiration) was quantified per watershed and statistically related to variables representing hydrology-related landscape characteristics (e.g., infiltration capacity, surface storage capacity, stream density). To investigate the influence stream connectivity has on the rate at which surface water leaves a given location, we modeled stream-connected and stream-disconnected surface water separately. Stream-connected surface water showed a greater expansion with wetter climatic conditions in landscapes with greater total wetland area, but lower total wetland density. Disconnected surface water showed a greater expansion with wetter climatic conditions in landscapes with higher wetland density, lower infiltration and less anthropogenic drainage. From these findings, we can expect that shifts in precipitation and evaporative demand will have uneven effects on surface water quantity. Accurate predictions regarding the effect of climate change on surface water quantity will require consideration of hydrology-related landscape characteristics including wetland storage and arrangement.

Prairie Pothole Regino↗
USGS2015

Coral 13C/12C records of vertical seafloor displacement during megathrust earthquakes west of Sumatra

The recent surge of megathrust earthquakes and tsunami disasters has highlighted the need for a comprehensive understanding of earthquake cycles along convergent plate boundaries. Space geodesy has been used to document recent crustal deformation patterns with unprecedented precision, however the production of long paleogeodetic records of vertical seafloor motion is still a major challenge. Here we show that carbon isotope ratios ( ) in the skeletons of massive Porites corals from west Sumatra record abrupt changes in light exposure resulting from coseismic seafloor displacements. Validation of the method is based on the coral response to uplift (and subsidence) produced by the March 2005 M w 8.6 Nias&ndash;Simeulue earthquake, and uplift further south around Sipora Island during a M&sim;8.4 megathrust earthquake in February 1797. At Nias, the average step-change in coral was 0.6&plusmn;0.1&permil;/m for coseismic displacements of +1.8 m and &minus;0.4 m in 2005. At Sipora, a distinct change in Porites microatoll growth morphology marks coseismic uplift of 0.7 m in 1797. In this shallow water setting, with a steep light attenuation gradient, the step-change in microatoll is 2.3&permil;/m , nearly four times greater than for the Nias Porites . Considering the natural variability in coral skeletal , we show that the lower detection limit of the method is around 0.2 m of vertical seafloor motion. Analysis of vertical displacement for well-documented earthquakes suggests this sensitivity equates to shallow events exceeding M w &sim;7.2 in central megathrust and back-arc thrust fault settings. Our findings indicate that the coral paleogeodesy technique could be applied to convergent tectonic margins throughout the tropical western Pacific and eastern Indian oceans, which host prolific coral reefs, and some of the world's greatest earthquake catastrophes. While our focus here is the link between coral , light exposure and coseismic crustal deformation, the same principles could be used to characterize interseismic strain during earthquake cycles over the last several millennia.

Earth and Planetary Science Letters↗
USGS2024

Automatic identification and quantification of volcanic hotspots in Alaska using HotLINK: The hotspot learning and identification network

An increase in volcanic thermal emissions can indicate subsurface and surface processes that precede, or coincide with, volcanic eruptions. Space-borne infrared sensors can detect hotspots—defined here as localized volcanic thermal emissions—in near-real-time. However, automatic hotspot detection systems are needed to efficiently analyze the large quantities of data produced. While hotspots have been automatically detected for over 20 years with simple thresholding algorithms, new computer vision technologies, such as convolutional neural networks (CNNs), can enable improved detection capabilities. Here we introduce HotLINK: the Hotspot Learning and Identification Network, a CNN trained to detect hotspots with a dataset of −3,800 satellite-based, Visible Infrared Imaging Radiometer Suite (VIIRS) images from Mount Veniaminof and Mount Cleveland volcanoes, Alaska. We find that our model achieves an accuracy of 96% (F1-score 0.92) when evaluated on −1,700 unseen images from the same volcanoes, and 95% (F1-score 0.67) when evaluated on −3,000 images from six additional Alaska volcanoes (Augustine Volcano, Bogoslof Island, Okmok Caldera, Pavlof Volcano, Redoubt Volcano, Shishaldin Volcano). In comparison with an existing threshold-based hotspot detection algorithm, MIROVA (Coppola et al., Geological Society, London, Special Publications, 2016, 426, 181–205), our model detects 22% more hotspots and produces 12% fewer false positives. Additional testing on −700 labeled Moderate Resolution Imaging Spectroradiometer (MODIS) images from Mount Veniaminof demonstrates that our model is applicable to this sensor’s data as well, achieving an accuracy of 98% (F1-score 0.95). We apply HotLINK to 10 years of VIIRS data and 22 years of MODIS data for the eight aforementioned Alaska volcanoes and calculate the radiative power of detected hotspots. From these time series we find that HotLINK accurately characterizes background and eruptive periods, similar to MIROVA, but also detects more subtle warming signals, potentially related to volcanic unrest. We identify three advantages to our model over its predecessors: 1) the ability to detect more subtle volcanic hotspots and produce fewer false positives, especially in daytime images; 2) probabilistic predictions provide a measure of detection confidence; and 3) its transferability, i.e., the successful application to multiple sensors and multiple volcanoes without the need for threshold tuning, suggesting the potential for global application.

Alaska↗