USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geosciences”

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 271 records · Page 15Linked to original sources

Altitude-age relationships of the lunar maria

Altitudes and relative ages of mare surface units were compared to test if a systematic correlation in height of lava eruption surfaces and age might reflect a corresponding increase in depth of the magma chamber with time; in addition the altitudes were studied to shed light on the time and place of warping of mare surfaces. The laser altimeter data from the Apollo missions and relative age data based on crater erosion models and crater counts were used for the study. The data were correlated by using the image data bank of the Lunar Geoscience Consortium. Results of the first part of the study are inconclusive as no systematic increase in height of lavas with time could be shown. The data of the second part of the study support the conclusion that mare surfaces may have warped throughout most of the time represented by sample ages from the Apollo missions, and that the lithosphere may have become strong enough to remain stable after the time of Apollo 12 samples around 3 b.y. ago.

Conference Paper↗

The Appalachian Geo-STEM Camp: Learning about geology through experiential adventure recreation

The inaugural Appalachian Geo-STEM Camp (AGC) was a partnership between West Virginia University (WVU), the U.S. Geological Survey (USGS) and the West Virginia Geological and Economic Survey (WVGES). Designed to engage high school students in geoscience-oriented Science, Technology, Engineering and Mathematics (STEM) activities through adventure-based outdoor recreation, the inaugural AGC took place in June 2018, with its base operations at the WVU Natural Resources Center (NRC), located northeast of Morgantown, West Virginia. The goals of the AGC are to increase the knowledge of the teenaged campers about the geological formations and biodiversity in the region, to acquaint them with geologic mapping technology used by USGS, WVGES, and WVU, and to foster interest in STEM-based careers. Nine students participated, with a cadre from the USGS, WVGES, and WVU teaching lessons in local geology and ecology. Inaugural-year efforts were focused on development and logistics of the camp and what activities best complimented the STEM research. Post-evaluations by the participants were generally favorable. Year-two goals are to fully develop a curriculum, and conduct a thorough pre-camp and post-camp participant survey to quantify learning objectives and guide the sustainability of the effort.

West Virginia↗

Opportunities to enhance seismic demand parameters for future editions of the AS1170.4

Geoscience Australia has recently released its 2018 National Seismic Hazard Assessment (NSHA18). Results from the NSHA18 indicate significantly lower seismic hazard across almost all Australian localities at the 1/500 annual exceedance probability level relative to the factors adopted for the current Australian Standard AS1170.4–2007 (R2018). These new hazard estimates, coupled with larger probability factors (kp) for long return periods, have challenged notions of seismic hazard in Australia in terms of the recurrence of damaging ground motions. As a consequence, the new hazard estimates have raised questions over the appropriateness of the prescribed National Construction Code probability level as used in the AS1170.4 to determine appropriate seismic demands for the design of ordinary-use structures. Therefore, it is suggested that the ground-motion exceedance probability used in the current AS1170.4 be reviewed in light of the recent hazard assessment and the expected performance of modern buildings for rarer ground motions. Whilst adjusting the AS1170.4 exceedance probability level would be a major departure from previous earthquake loading standards, it would bring it into line with other international building codes in similar tectonic environments. Additionally, it would offer opportunities to further modernise how seismic demands are considered in Australian building design. In particular, the authors highlight the following additional opportunities: 1) the use of uniform hazard spectra to replace and simplify the spectral shape factors, which do not deliver uniform hazard across all natural periods; 2) updated site amplification factors to ensure continuity with modern ground-motion models, and; 3) the potential to define design ground motions in terms of uniform collapse risk rather than uniform hazard. Estimation of seismic hazard at any location is an uncertain science. However, as our knowledge improves, our estimates of the hazard will converge more closely to the actual – but unknowable – (time independent) hazard. It is therefore prudent to regularly update the estimates of the seismic demands in our building codes using the best available evidence-based methods and models.

Conference Paper↗

EMD Gas Hydrates Committee annual report

Global research on the potential commercial viability of gas extraction from gas hydrates is continuing, predominantly in Asia and in the United States, where recent efforts have focused on the exploration and characterization of gas hydrate petroleum systems and conducting controlled production tests of gas hydrate deposits hosted in mostly sand-rich reservoir systems. In 2017, gas hydrate production testing was conducted in marine settings in the offshore of Japan and China. An additional test was conducted in the Shenhu region of the South China Sea in 2020 to further evaluate the effectiveness of various well completion technologies to produce gas from gas hydrates. Late in 2018, a research partnership led by the U.S. Department of Energy National Energy Technology Laboratory and the Japan Oil, Gas and Metals National Corporation established a new Alaska North Slope test site for an extended gas hydrate production pilot with the completion of a stratigraphic test well. With the successful completion of the Prudhoe Bay Unit Hydrate-01 stratigraphic test well, planning is now underway for drilling of three additional wells at the Alaska test site, which will include a geoscience data well and two production test wells. In 2019, the Government of India released the results of one of the most comprehensive gas hydrate scientific drilling investigations ever undertaken, the leadership of the Indian National Gas Hydrate Program Expedition 02 reported on the discovery of several significant gas hydrate accumulations that were considered suitable for future gas hydrate production testing. In September of 2019, the U.S. Geological Survey released a new assessment of the technically recoverable gas hydrate resources on the Alaska North Slope. In 2018, the University of Texas at Austin released a comprehensive report describing the operational and initial technical results of the UT-GOM2-1 Hydrate Pressure Coring Expedition, which drilled and partially cored two wells in Green Canyon Block 955 in the deepwater Gulf of Mexico. The results of these recently complete gas hydrate geologic and production studies have been further reviewed in this AAPG-EMD Gas Hydrate Committee report.

Alaska↗

Geologic map of the Butte City 7.5' Quadrangle, Butte County, Idaho

The geologic map of the Butte City 7.5’ quadrangle is based on mapping summarized in the 1:100,000 scale map of the Idaho National Laboratory, U.S. Geological Survey Miscellaneous Investigations Map I-2330, by Kuntz and others, 1994. New surficial geologic mapping was completed by National Association of Geoscience Teachers (NAGT) interns, Evan Martin (2015) and Samuel Helmuth (2017). Previously published maps and associated data, by various authors, including Kuntz, M.A., Betty Skipp, M.A. Lanphere, W.E. Scott, K.L. Pierce, G.B. Dalrymple, L.A. Morgan, D.E. Champion, G.F. Embree, W.R. Page, R.P. Smith, W.R. Hackett, and D.W. Rodgers from 1994 to 2017, were incorporated into this project. The U. S. Geological Survey (USGS) Idaho National Laboratory Project Office, a field office of the Idaho Water Science Center, conducts subsurface investigations to study the Snake River aquifer at the Idaho National Laboratory. Many subsurface basalts erupted from vents exposed at the surface. Tracing the flows from the surface to the subsurface provides data for numerical models of groundwater flow and contaminant transport.

Idaho↗

Toward an effective global green economy: The Critical Minerals Mapping Initiative (CMMI)

Global population growth, economic development and the accelerating pace of technological innovation are driving increased demand for non-fuel mineral commodities that are vital for emerging and low-carbon technologies. Examples of such commodities include cobalt and graphite for rechargeable batteries, tellurium in thin-film solar photovoltaics and rare earth elements (REE) in permanent magnets, electronics and medical technologies. These commodities are known as critical elements and/or minerals (collectively referred to as critical minerals), with the term critical used to define not only their importance for new technologies, but also their demand and vulnerability to supply disruption. The demand for critical minerals is likely to continue to grow, but supply is not assured. Therefore, national strategies in Australia, Canada, United States and elsewhere (e.g. Europe: Wittenberg et al. 2021) are being developed to encourage exploration and production, including resolving the geological processes responsible for their enrichrichment into viable ore deposits. In 2019, Geoscience Australia (GA), the Geological Survey of Canada (GSC) and the United States Geological Survey (USGS) formed the Critical Minerals Mapping Initiative (CMMI) to undertake research to develop a better understanding of critical mineral resources in known deposits, determine the geological controls on critical mineral distribution for deposits currently producing by-products, and identify new sources of supply through mineral prospectivity mapping and resource assessment.

SGA News↗

Data-driven prospectivity modelling of sediment-hosted mineral systems

Mississippi Valley-type (MVT) and clastic-dominated (CD) deposits are important sources for Zn, Pb, Ag, and Cd as well as the critical elements Ga, Ge, In, and Sb. However, mapping the drivers, sources, pathways, and traps of MVT and CD deposits within the much larger and mostly unmineralized sedimentary basins remain some of the least understood aspects of these mineral systems. Herein we address those knowledge gaps by integrating public geoscience datasets from Canada, the United States of America, and Australia using a discrete global grid system to map the continent-scale footprints of MVT and CD deposits.

Conference Paper↗

City-scale geothermal energy everywhere to support renewable resilience – A transcontinental cooperation

Cities have important and varying incentives to transform their energy sector to all-electric with low carbon emissions. However, they often encounter a number of impediments when attempting to implement such a change. For example, while urban areas have the highest energy demand-density, cities often lack the space for installing additional energy generation and/or long-duration energy storage systems. Cities also have existing environmental issues from energy sources (e.g., pollution from dust, waste heat or noise) that make residents sensitive to energy infrastructure development. Utilizing power from conventional sources, such as natural gas, biomass and hydropower, which usually are distanced from the urban areas, also make cities more vulnerable to supply disruptions. One promising de-carbonizing energy option for cities focuses on their heating and cooling needs, which constitutes around 1/3 of U.S. and 1/2 of European energy consumption (including industrial processes like drying, pasteurization, etc.; Jadun and others, 2017; EU Commission 2022). If heating and cooling loads can be met by geothermal direct-use technologies, then the need for new electric sources can be greatly lessened. Despite the proven efficacy of geothermal energy as a city/community-scale heating and cooling resource, it is currently only a niche resource in the heating and cooling sector, though has significant potential for future growth. Historically, emphasis has been placed on geothermal electricity generation potential that requires higher temperature (greater than 90 °C) resources at drillable depths, but potentially viable areas are geographically limited and typically well removed from urban centers. Key drivers for investments were represented by greater political interest in renewable electricity production, higher revenues and less effort in distributing the produced energy via grids. In contrast, low-temperature (less than 90 °C) geothermal resources can be used directly for heating and cooling almost everywhere and are cost-effective in urban/suburban settings. In addition, the increased prominence of renewable electricity sources, such as wind and solar onto city-scale electric grids, has led to new urgency around questions of energy storage. Underground thermal energy storage (UTES), wherein surplus or waste heat is stored underground for later use, could present a long-duration energy storage solution. From October 2022 through September 2024, a transcontinental consortium consisting of geological surveys, geoscience organizations, industry representatives and universities aims to develop an understanding of the global potential for city-scale geothermal, proposing guidelines to aid in promoting the economic utilization of low temperature geothermal resources. Efforts will focus on providing city managers and other decision makers with the information needed to evaluate and implement suitable city/community-scale geothermal technologies. Funded by the U.S. Geological Survey’s John Wesley Powell Center for Analysis and Synthesis, this interdisciplinary consortium will showcase tools, datasets, and scientific recommendations to accelerate the broader understanding and adoption of renewable energy systems that access geothermal resources. The collaborative research activities include standardization of nomenclatures, resource description and characterization strategies globally. The results from these activities will be combined with a preliminary climate-driven, city-based energy needs related analysis to perform energy supply/demand matching analysis. The identification of city-specific applications that would benefit from the geothermal technologies provides the basis to up-scale city-specific determinations to regional and national assessments of resource estimates. The city-scale geothermal energy research initiative will ultimately provide the synergies and management analysis that can address benefits, environmental impacts, regulatory frameworks, sustainability, and suitability in retrofitted buildings or new as well as existing heating networks.

Conference Paper↗

U.S. Geological Survey geoheritage sites of the nation: Broadening the relevance of geologic mapping through the intersection of geology and society

Geoheritage lies at the intersection of geology and society (NASEM 2021). It is a value-laden term applied to geologically significant features and landforms to highlight the scientific, educational, cultural, economic, and aesthetic values that society places on them (Brilha 2018, GSA 2017, NPS and AGI 2015). The U.S. Geological Survey (USGS) National Cooperative Geologic Mapping Program (NCGMP) proposed the development of a geoheritage geospatial product (GGP) to broaden the relevance of geologic mapping to stakeholders within the scientific community and the public (Shelton et al. 2022). To accomplish this, the USGS Geoheritage Sites of the Nation project was initiated to design the GGP and build a geoheritage site inventory to showcase geodiversity and natural heritage throughout the United States and its territories to accompany the seamless 2D and 3D geologic map and model being developed by the NCGMP and its partners (DeVera et al. 2022). The GGP serves as an educational and outreach tool to 1) raise awareness of significant geologic sites, 2) communicate the role of geology in our natural heritage, and 3) increase accessibility and relevance of geologic mapping and geoscience topics more broadly to the public (DeVera et al. 2023b).

Conference Paper↗

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Computers & Geosciences↗

Methane hydrate-bearing seeps as a source of aged dissolved organic carbon to the oceans

Marine sediments contain about 500–10,000 Gt of methane carbon1, 2, 3, primarily in gas hydrate. This reservoir is comparable in size to the amount of organic carbon in land biota, terrestrial soils, the atmosphere and sea water combined1, 4, but it releases relatively little methane to the ocean and atmosphere5. Sedimentary microbes convert most of the dissolved methane to carbon dioxide6, 7. Here we show that a significant additional product associated with microbial methane consumption is methane-derived dissolved organic carbon. We use Δ14C and δ13C measurements and isotopic mass-balance calculations to evaluate the contribution of methane-derived carbon to seawater dissolved organic carbon overlying gas hydrate-bearing seeps in the northeastern Pacific Ocean. We show that carbon derived from fossil methane accounts for up to 28% of the dissolved organic carbon. This methane-derived material is much older, and more depleted in 13C, than background dissolved organic carbon. We suggest that fossil methane-derived carbon may contribute significantly to the estimated 4,000–6,000 year age of dissolved organic carbon in the deep ocean8, and provide reduced organic matter and energy to deep-ocean microbial communities.

Nature Geoscience↗

Broad accommodation of rift-related extension recorded by dyke intrusion in Saudi Arabia

The extensive harrat lava province of Arabia formed during the past 30 million years in response to Red Sea rifting and mantle upwelling. The area was regarded as seismically quiet, but between April and June 2009 a swarm of more than 30,000 earthquakes struck one of the lava fields in the province, Harrat Lunayyir, northwest Saudi Arabia. Concerned that larger damaging earthquakes might occur, the Saudi Arabian government evacuated 40,000 people from the region. Here we use geologic, geodetic and seismic data to show that the earthquake swarm resulted from magmatic dyke intrusion. We document a surface fault rupture that is 8 km long with 91 cm of offset. Surface deformation is best modelled by the shallow intrusion of a north-west trending dyke that is about 10 km long. Seismic waves generated during the earthquakes exhibit overlapping very low- and high-frequency components. We interpret the low frequencies to represent intrusion of magma and the high frequencies to represent fracturing of the crystalline basement rocks. Rather than extension being accommodated entirely by the central Red Sea rift axis, we suggest that the broad deformation observed in Harrat Lunayyir indicates that rift margins can remain as active sites of extension throughout rifting. Our analyses allowed us to forecast the likelihood of a future eruption or large earthquake in the region and informed the decisions made by the Saudi Arabian government to return the evacuees.

Nature Geoscience↗

Crop type classification, trends, and patterns of central California agricultural fields from 2005 to 2020

California produces many key agricultural products in the United States. Current geospatial agricultural datasets are limited in mapping accuracy, spatial context, or observation period. This study uses machine learning and high-resolution imagery to produce a time series of crop maps to assess crop type trends and patterns across central California from 2005 to 2020. National Agriculture Imagery Program and Landsat imagery were used to classify nine crop types that are common in the study region: grain crops, field crops, rice, citrus and subtropical, deciduous fruit and nut, vineyard, berry and vegetable, pasture, and fallow/young perennial crop types. To create labeled data, we sampled 1253 fields and manually identified crop types for each examined year using high-resolution imagery and Landsat normalized difference vegetation index time series. We applied a random forest machine learning algorithm in Google Earth Engine. Results show that the mean overall classification accuracy of the nine-class map was 93.1%, with individual accuracies ranging from 99.3% (rice) to 89.5% (fallow/young perennial). Mann–Kendall trend tests showed significant ( p less than 0.05) declines in field crop and pasture area during the study period, while deciduous fruit and nut, citrus and subtropical, and fallow/young perennial crop types experienced significant increases. At an aggregate level, there was a general shift from annual crop types to perennial crop types. These data provide a 16-year time span of spatially explicit crop type classifications, trends, and patterns in central California that can be used to aid managers and decision makers for resource planning or hazard mitigation.

California↗

Patterns of water-extractable soil organic matter in the US Great Plains: Insights from the Haas Soil Archive

Novel approaches that are fast and sensitive are needed to evaluate soil change and integrate soil ecosystem properties. Carbon (C) and nitrogen (N) extracted from soil with water are associated with plant nutrients and microbial activity but information about change over time in the US Great Plains is sparse. We used cool (20°C) and hot (80°C) water extracts from historic (1947) and contemporary (2018) soil samples collected at Moccasin, MT; Akron, CO; and Big Spring, TX; to examine changes to labile C and N and optical properties after 71 years of dryland cropping. Concentrations of C and N extracted with cool water decreased between 1947 and 2018 in surface (0–15.2 cm) samples from Moccasin, by 52% and 35%, and Big Spring, by 37% and 32%, but remained unchanged at Akron. Conversely, net (hot−cool) extractable C did not change at Moccasin or Big Spring but increased at Akron by 26%. Net extractable N decreased at Moccasin by 22% but did not change elsewhere. Sequential principal component analysis and stepwise discriminant analysis identified three important optical properties. Values of SUVA 254 (where SUVA 254 is the specific ultraviolet absorbance at 254 nm) in extracts did not change at Moccasin between 1947 and 2018 but increased at Akron, indicating increased aromaticity. Conversely, SUVA 254 decreased at Big Spring. Values for Sag 350–400 (where Sag 350–400 is the slope from a nonlinear fit of an exponential function to the absorption spectrum over the wavelength range from 350 to 400 nm), inversely related to extract molecular weight and aromaticity, decreased at Moccasin but not elsewhere. The proportion of recalcitrant to labile compounds, C:T (where C:T is the ratio of fluorescence intensity from Peak C [ex340/em440] to Peak T [ex275/em340]), increased in extracts from all sites but especially at Akron. Together, these methods provided insights into soil change while conserving samples.

Colorado, Kansas, Montana, Nebraska, New Mexico, N↗

Analytical model for screening potential CO2 repositories

Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small.

Computational Geosciences↗

Linear complementarity formulation for 3D frictional sliding problems

Frictional sliding on quasi-statically deforming faults and fractures can be modeled efficiently using a linear complementarity formulation. We review the formulation in two dimensions and expand the formulation to three-dimensional problems including problems of orthotropic friction. This formulation accurately reproduces analytical solutions to static Coulomb friction sliding problems. The formulation accounts for opening displacements that can occur near regions of non-planarity even under large confining pressures. Such problems are difficult to solve owing to the coupling of relative displacements and tractions; thus, many geomechanical problems tend to neglect these effects. Simple test cases highlight the importance of including friction and allowing for opening when solving quasi-static fault mechanics models. These results also underscore the importance of considering the effects of non-planarity in modeling processes associated with crustal faulting.

Computational Geosciences↗

Bayesian analysis of the impact of rainfall data product on simulated slope failure for North Carolina locations

In the past decades, many different approaches have been developed in the literature to quantify the load-carrying capacity and geotechnical stability (or the Factor of Safety, F_s) of variably saturated hillslopes. Much of this work has focused on a deterministic characterization of hillslope stability. Yet, simulated F_s values are subject to considerable uncertainty due to our inability to characterize accurately the soil mantle’s properties (hydraulic, geotechnical and geomorphologic) and spatiotemporal variability of the moisture content of the hillslope interior. This is particularly true at larger spatial scales. Thus, uncertainty-incorporating analyses of physically based models of rain-induced landslides are rare in the literature. Such landslide modeling is typically conducted at the hillslope scale using gauge-based rainfall forcing data with rather poor spatiotemporal coverage. For regional landslide modeling, the specific advantages and/or disadvantages of gauge-only, radar-merged and satellite-based rainfall products are not clearly established. Here, we compare and evaluate the performance of the Transient Rainfall Infiltration and Grid-based Regional Slope-stability analysis (TRIGRS) model for three different rainfall products using 112 observed landslides in the period between 2004 and 2011 from the North Carolina Geological Survey database. Our study includes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis Version 7 (TMPA V7), the North American Land Data Assimilation System Phase 2 (NLDAS-2) analysis, and the reference ‘truth’ Stage IV precipitation. TRIGRS model performance was rather inferior with the use of literature values of the geotechnical parameters and soil hydraulic properties from ROSETTA using soil textural and bulk density data from SSURGO (Soil Survey Geographic database). The performance of TRIGRS improved considerably after Bayesian estimation of the parameters with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm using Stage IV precipitation data. Hereto, we use a likelihood function that combines binary slope failure information from landslide event and ‘null’ periods using multivariate frequency distribution-based metrics such as the False Discovery and False Omission Rates. Our results demonstrate that the Stage IV-inferred TRIGRS parameter distributions generalize well to TMPA and NLDAS-2 precipitation data, particularly at sites with considerably larger TMPA and NLDAS-2 rainfall amounts during landslide events than null periods. TRIGRS model performance is then rather similar for all three rainfall products. At higher elevations, however, the TMPA and NLDAS-2 precipitation volumes are insufficient and their performance with the Stage IV-derived parameter distributions indicate their inability to accurately characterize hillslope stability.

North Carolina↗