USGS ScienceSearch

SEARCH · USGS Science

Results for “Article”

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.

979 records · Page 27Linked to original sources

Polystyrene microplastics alter the accumulation and elimination dynamics of silver nanoparticles in Daphnia magna

Rationale Microplastics (MPs) can interact with engineered nanomaterials and alter their environmental fate and bioavailability. However, their influence on the bioaccumulation dynamics of silver nanoparticles (AgNPs) in aquatic filter feeders under environmentally relevant conditions remains poorly understood. We hypothesized that polystyrene microplastics (PS-MPs) alter the uptake, elimination and overall bioaccumulation dynamics of AgNPs in Daphnia magna . Methodology Adult D. magna were exposed to isotopically labeled citrate-coated 109 AgNPs in the presence and absence of PS-MPs. Waterborne uptake, dietary uptake, assimilation efficiency, food ingestion and elimination were quantified experimentally and incorporated into a biodynamic model to predict steady-state silver (Ag) concentrations under environmentally relevant exposure scenarios. Results PS-MPs increased waterborne Ag accumulation by approximately eightfold compared with AgNPs exposure alone. However, Ag elimination was substantially faster, with 98% of accumulated Ag eliminated after five days of depuration compared with 66% without PS-MPs. PS-MPs also reduced food ingestion rates (IRs), while assimilation efficiency remained largely unchanged. Biodynamic modeling predicted that steady-state Ag accumulation was approximately fivefold greater in the presence of PS-MPs, with waterborne exposure becoming the dominant accumulation pathway. Discussion These findings suggest that PS-MPs alter bioaccumulation dynamics of AgNPs in filter-feeding organisms by enhancing organism-associated Ag during waterborne exposure while accelerating Ag elimination, likely through particle-associated transport and gut egestion. Collectively, these results underscore the complexity of Ag bioaccumulation processes when MPs and AgNPs interact and highlight the importance of studies conducted under environmentally relevant conditions.

Environmental Chemistry

Where will the cat cross the road? Comparing camera and GPS-based models for identifying wildlife corridors

Designing effective wildlife corridors is a critical conservation challenge in fragmented landscapes. GPS-based step selection functions strongly predict dispersal corridors and connectivity, but GPS collaring can be expensive and invasive. Camera-based occupancy models are widely used for connectivity analyses but may involve trade-offs in data resolution. Despite widespread use of both approaches, few studies have directly compared them using concurrent datasets. We developed a stacked single-species, single-season occupancy model and a Circuitscape connectivity surface for mountain lions (Puma concolor) on Washington’s Olympic Peninsula, USA, and compared them with a connectivity surface from an existing integrated step selection function. Both models predicted mountain lion GPS locations well, with binned Spearman rank correlations of 1 for Circuitscape and 0.96 for the step selection function, though step selection better identified habitat use by dispersers. Connectivity predictions were moderately correlated across the landscape ( r = 0.26), but agreement was strongest in human-dominated areas most critical for corridor planning. We conclude that GPS-based approaches are advantageous when data collection is feasible and the focus is on dispersal or fine-scale movement. However, camera-based approaches may be preferable for multi-species monitoring, large spatial and temporal scales, noninvasive sampling, when resources are limited, or when fine-scale or dispersal-specific inference is not required.

Washington

Light-dependent activity of deepwater sculpin (Myoxocephalus thompsonii) across substrates with and without predation risk

Light availability strongly influences predator–prey interactions in deepwater ecosystems, where visual constraints shape both foraging success and prey behavior. The behavioral response of deepwater sculpin ( Myoxocephalus thompsonii ) to siscowet lake trout ( Salvelinus namaycush siscowet ) was studied under ecologically relevant light intensities spanning several orders of magnitude typical of daytime downwelling light in the 20–100 m depth range of Lake Superior. Trials were conducted over varying substrates (gravel, sand, and black fabric). Deepwater sculpin showed a significant preference for gravel over sand and black fabric. In the absence of siscowet, sculpin movement frequency increased as light intensity decreased. Sculpin reaction distance to siscowet was influenced by both light and substrate. Reaction distance was shortest at low light intensities, peaked at intermediate light intensities (3.05 × 10⁹ to < 6.0 × 10⁹ photons m⁻2 s⁻1), and declined again at the highest light intensities tested. In the presence of siscowet, sculpin activity was suppressed at the upper end of the light range (≥ 6.0 × 10⁹ photons m⁻2 s⁻1). The greatest increase in movement occurred between 6.0 × 10⁹ and 3.05 × 10⁹ photons m⁻2 s⁻1, which corresponded to the range where siscowet prey capture declined, suggesting sculpin exploit this low-light window to move with reduced risk. Reduced activity in the presence of predators is common among cryptic species, and our findings suggest that sculpin restrict movement at higher light levels to avoid detection by siscowet.

Wisconsin

Relating glassy rind thicknesses to ambient air temperatures at the Lost Jim flow field in the Imuruk Lake volcanic field, Alaska

The Lost Jim flow field, in the Imuruk Lake volcanic field, Alaska, extends west ~34 km from a single vent, crossing subarctic tundra and currently touches several lakes and streams. The weighted mean of five 36Cl cosmogenic exposure ages from the Lost Jim pāhoehoe flow is 7.73 ± 0.37 ka, indicating this eruption occurred substantially after the eruption of the underlying Camille flow, which was emplaced at 39.7 ± 1.3 ka. Paleoclimate records indicate the period when the Lost Jim flow field was emplaced was after deglaciation, and the climate was similar to today. We propose that the emplacement of lava in these cold subarctic conditions can lead to faster cooling of the lava surface compared to lava emplaced in warmer locations such as mid- latitude cold deserts. Glass abundance in the outermost rinds at the Lost Jim flow field was on average 74 % with 6.4 mm thick rims, compared to 60 % with 2.9 mm rims for cold mid-latitude desert samples. We interpret increased glass content as a proxy for rapid cooling likely occurring partly during winter. Glassiness values varied less across vent, margin, and mid-flow locations when compared to the mid-latitude flows suggesting the Lost Jim flow field was broadly impacted by the subarctic climate as opposed to responding to local microclimates. Our results indicate that lava glassiness may be a useful environmental indicator of cooler (in this case subarctic) conditions.

Alaska

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska

Machine learning provides reconnaissance-type estimates of carbon dioxide storage resources in oil and gas reservoirs

Oil and gas reservoirs represent suitable containers to sequester carbon dioxide (CO 2 ) in a supercritical state because they are accessible, reservoir properties are known, and they previously contained stored buoyant fluids. However, planners must quantify the relative magnitude of the CO 2 storage resource in these reservoirs to formulate a comprehensive strategy for CO 2 mitigation. Even reconnaissance-type estimates of CO 2 storage resources of known oil and gas reservoirs may require complicated calculations involving 1) estimates of recoverable oil and gas, 2) reservoir properties (depth, temperature, pressure, etc.), and 3) the physical qualities of the retained fluids. We demonstrate the application of machine learning (ML) algorithms to bypass these computations to yield more rapid estimates of CO 2 storage resources in reservoirs capable of hosting CO 2 in a supercritical state. ML algorithms are computationally efficient because they do not impose the strong assumptions on the data-generating process that standard statistical or engineering procedures require. Further, ML algorithms can capture highly complex, particularly nonlinear, relationships among predictor variables. We demonstrate the application of four different ML algorithms using data from onshore and offshore oil and gas reservoirs in Europe, and show they perform well when predictions are compared to engineering estimates. The proposed methods and models provide an effective and novel way to more rapidly and directly determine the subsurface CO 2 storage capacity of oil and gas reservoirs around the world, information that operators, researchers, and policymakers alike require to meet energy transition and decarbonization goals.

Frontiers in Enviornmental Science

Case study of deep learning image segmentation for the purposes of rapid 2D petrographic analysis in volcanic rocks

Automation using deep learning methods is a useful alternative to manual methods of petrographic segmentation, but often requires user familiarity with coding and/or algorithms. We examine the Dragonfly TM program's deep learning tools for application by users with a variety of skill levels as a method for petrographic image segmentation. An image processing methodology, bimodal image stacking, was created for low-input-data, high-efficacy training of models which can then be applied to varied samples. Using backscatter electron images we show that the resulting model segmentations agree with manual segmentation total and modal crystallinity values within 5%, and calculated plagioclase crystal size distribution (CSD) values within 2σ, despite limitations in discriminating mafic phases. Model creation and training takes <24 hours, 1–3 hours of which are supervised, and the resultant model can then be applied to new uncharacterized samples in <15 minutes per image. This allows for non-experts to create and utilize deep learning models to segment images of variable brightness and texture, at low user-time cost and resulting in size and shape data which are within uncertainty of manual segmentation. While some limitations are noted (for example, sieve-textured phases may need manual correction, and different minerals with similar BSE intensity may not be resolved as separate phases), this methodology can be utilized for general application of models to wide ranges of volcanic crystalline and bubble textures, and to create a library of models for rapid petrological analysis during volcanic eruptions.

Alaska

Multi-site evaluation of a postfire debris-flow runout forecast method

Postfire debris flows pose a hazard to human life, property, and infrastructure when they travel from steep source areas to urbanized alluvial fans or other developed areas. Existing methods for rapid (<1 week) postfire debris-flow hazard assessment document the increase in the likelihood and size of debris flows as the magnitude of high-intensity rain necessary to initiate debris flows increases but do not indicate the extent of downstream debris-flow runout. Although many models for the simulation of debris-flow motion are available, there is no established approach for using these models to delineate locations susceptible to postfire debris-flow runout that (a) is feasible to use at the spatial scale of an entire fire; (b) is appropriate for runout onto unconfined areas; (c) reproduces observed relations between runout and rainfall intensity; and (d) characterizes inherent uncertainty in runout, even without spatiotemporally variable rainfall. We propose and evaluate a method for generating postfire debris-flow runout hazard maps that has all the above qualities. Selection of case studies prioritized events triggered by a range of rainfall intensities, locations within and outside of southern California, and observed runout onto unconfined topography. Qualitative and quantitative assessment of performance for four events indicate that simulation results broadly match observations albeit with some discrepancies at a scale larger than structure or land parcel level (approximately 20-m by 20-m). The method may be used to identify potentially hazardous areas immediately following a fire and to provide approximate runout forecasts when a storm is imminent.

Arizona, California

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Proceedings of the National Academy of Sciences

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth

The Mammoth magnetic anomaly, Pinal County, Arizona

A high-resolution Earth Mapping Resources Initiative airborne geophysical survey was flown in the southwest North American porphyry copper province to improve bedrock geologic maps and to identify areas that have unrecognized critical mineral resource potential. During the review of the aeromagnetic data, a distinctly monopolar-shaped, negative magnetic anomaly was observed at a flight elevation of 200 m above the ground with a maximum amplitude of –9500 nT. We have named this the Mammoth magnetic anomaly (MMA) because it is centered 12 km northeast of the town of Mammoth, Arizona, USA. The total field anomaly (TFA) contour of –500 nT enclosing the MMA defines an elongate shape measuring 2.5 km long by 1 km wide that trends northwest–southeast. Given the striking nature of this negative, monopolar-shaped magnetic anomaly, we conducted a ground campaign in May 2025 to determine its authenticity and potential relationship to critical mineral endowment. The MMA was confirmed on the ground with a TFA approaching –46,000 nT. Total magnetic intensity (TMI) observations routinely fell below the 18,000 nT operating floor of an industry-standard cesium-vapor total field magnetometer, and extremely low TMI measurements were corroborated along coincident traverse lines using two high dynamic range, but lower sensitivity, smartphone vector magnetometers. The lowest TMI values recorded by both smartphone magnetometers were 1000 nT and confirmed with multiple adjacent and crossing lines. Field observations suggest that this magnetic feature is caused by strong remanent magnetization within fine-grained magnetite hosted within locally altered Pinal Schist.

Arizona

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota

Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access

The growing literature on deaf and hard of hearing (DHH+) populations and disasters demonstrates that emergency communication (including alerts) is not reaching global DHH + individuals with dangerous impacts for morbidity and mortality. This is the first research study in the U.S. to qualitatively explore the experiences of DHH + persons with earthquake early warning (EEW) through group-based dialogue sessions. The study investigates eight DHH + university students'past earthquake experiences, access to EEW alerts, and perceptions of ShakeAlert Ⓡ , an EEW system for detecting earthquakes and alerting residents of California, Oregon, and Washington. Findings highlight key gaps in disaster alert usability within four thematic areas: lack of messaging in participants' language(s), unclear alert messaging, deficient message delivery mechanisms for deafblind persons, and insufficient access to earthquake information and training that leads to dependence on informal information networks. Weaknesses identified in these four themes reduce DHH + trust in EEW systems and compromise the capacity of alert recipients to take swift protective action or to mentally prepare before shaking starts. The study also underscores structural factors such as insufficient linguistic representation in disaster language planning and technology design, which ignores the linguistic and sensory access needs of DHH + individuals. Building on disaster language planning frameworks, we recommend involving DHH + populations to co-develop EEW alerts. By centering DHH + perspectives, this research contributes to ongoing efforts to ensure that EEW systems reach everyone.

California

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances