USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research Earth Surface”

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 37 records · Page 2Linked to original sources
USGS2020

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗
USGS2019

Combining numerical and statistical models to predict storm-induced dune erosion

Dune erosion is an important aspect to consider when assessing coastal flood risk, as dune elevation loss makes the protected areas more susceptible to flooding. However, most advanced dune erosion numerical models are computationally expensive, which hinders their application in early-warning systems. Based on a combination of probabilistic and process-based numerical modeling, we develop an efficient statistical tool to predict dune erosion during storms. The analysis focuses on Dauphin Island, AL in the northern Gulf of Mexico, where we combine synthetic sea storms with a calibrated and validated XBeach model to develop and test a range of different surrogate models for their ability to predict barrier-island geometric parameters under storm conditions. Surrogate models are developed by combining the oceanographic forcing from 100 optimally sampled sea-storm events covering the entire multivariate parameter space (used as XBeach input) and associated changes in the dune system (XBeach output). We test four surrogate models using a k-fold approach for validation. All models perform well in predicting changes in dune elevation, barrier-island area, and width, but are less accurate in predicting alterations in the cross-shore locations of dune morphological features. Multivariate adaptive regression splines are identified as the best surrogate model based on their fast development and good performance, attaining a modified Mielke index of 0.81 for dune crest height. As demonstrated at Dauphin Island, our approach shows potential to be used in an operational framework to predict dune response (in particular crest elevation change) when water level and wave forecasts are available.

Alabama↗
USGS2019

Widespread initiation, reactivation, and acceleration of landslides in the northern California Coast Ranges due to extreme rainfall

Episodically to continuously active slow-moving landslides are driven by precipitation. Climate change, which is altering both the frequency and magnitude of precipitation world21 wide, is therefore predicted to have a major impact on landslides. Here we examine the behavior of hundreds of slow-moving landslides in northern California in response to large changes in annual precipitation that occurred between 2016 and 2018. We quantify the landslide displacement using repeat-pass radar interferometry and pixel offset tracking techniques on a novel dataset from the airborne NASA/JPL Uninhabited Aerial Vehicle Synthetic Aperture Radar. We found that 312 landslides were moving due to extreme rainfall during 2017, compared to 119 during 2016, which was the final year of a historic multi-year drought. However, with a return to below average rainfall in 2018, only 146 landslides remained in motion. The increased number of landslides during 2017 was primarily accommodated by landslides that were smaller than the landslides that remained active between 2016 and 2018. Furthermore, by examining a subset of 51 landslides, we found that 49 had increased velocities during 2017 when compared to 2016. Our results show that slow-moving landslides are sensitive to large changes in annual precipitation, particularly the smaller and thinner landslides that likely experience larger basal pore water pressure changes. Based on climate model predictions for the next century in California, which include increases in annual precipitation and increases in the frequency of dry-to-wet extremes, we hypothesize that there will be an overall increase in landslide activity.

California↗
USGS2020

Causes of variability in suspended‐sand concentration evaluated using measurements in the Colorado River in Grand Canyon

Rivers commonly exhibit substantial variability in suspended‐sand concentration, even at constant water discharge. Here we derive an approach for evaluating how much of this variability arises from mean bed‐sand grain size. We apply this approach to the Colorado River in Grand Canyon, where discharge‐independent concentration of suspended sand varies by more than a factor of 23 ( N = 1.4 × 10 6 ). Theory predicts that where concentration is controlled by bed‐sand grain size, concentration and grain size in suspension will be inversely correlated (i.e., coarsening of the bed causes suspended sand to become coarser in grain size and lower in concentration). Although the observed correlation is negative, riverbed grain size accounts for only 40% of the variability in concentration. The residuals vary by an order of magnitude; they arise from other processes, such as changes in topography or distribution of sand that cause shear stress to change at constant discharge, changes in the fine tail of bed‐sand grain sizes or changing bedforms. Both bed sand and the other factors influence concentration for durations from less than 1 day to several years. Predictions of concentration based on bed‐sand grain size ( N = 4 × 10 4 ) are less accurate than predictions based on suspended‐sand grain size, probably because suspended sand is a natural integrator of sand‐transporting processes, giving more weight to those areas of the bed that exchange more sand with the flow. Although the causes of variability vary from one river to another, the approach illustrated here is applicable to any river in which concentration varies at constant water discharge.

Arizona↗
USGS2020

Does channel narrowing by floodplain growth necessarily indicate sediment surplus? Lessons from sediment‐transport analyses in the Green and Colorado rivers, Canyonlands, Utah

Analyses of suspended sediment transport provide valuable insight into the role that sediment supply plays in causing geomorphic change. The sediment supply within a river system evolves depending on the discharge, flood frequency and duration, changes in sediment input, and ecohydraulic conditions that modify sediment transport processes. Changes in supply can be evaluated through analyses of coupled changes in suspended sediment concentration and grain size. The concentration of sand in transport in the Green and Colorado Rivers is most strongly controlled by discharge and the bed sand grain size distribution. Since the 1950s, sand loads have decreased in response to declines in peak discharge in the Green River and coarsening of the bed sand in the Colorado River. However, changes in the bed sand grain size distribution are associated with large changes in suspended sand concentration in both rivers; concentration varies by a factor of ~3 in the Green River and a factor of ~8 in the Colorado River, depending on the bed sand grain size distribution. Analyses of hysteresis in suspended sediment measurements show that sediment depletion during annual floods is most strongly controlled by flood duration, with peak discharge being nearly equally important in the Green River. Despite channel narrowing in both rivers, periods of bed sand coarsening and sediment depletion during annual floods indicate that these rivers are not necessarily in sediment surplus. Channel narrowing appears to be strongly controlled by short‐term declines in flood magnitude and the ecohydraulic effects of vegetation and may not be indicative of the long‐term sediment budget.

Utah↗
USGS2020

Atmospheric dust deposition varies by season and elevation in the Colorado Front Range, USA

As atmospheric dust deposition continues to increase across the southwestern United States, it has the potential to alter ecosystem productivity and structure by delivering nutrients, base cations, and pollutants to remote mountain sites. Due to the sparse distribution of dust monitoring sites, open questions remain about the spatial and temporal variability of dust fluxes and composition across mountainous terrain. We present a 1 year (November 2017 to November 2018) record of seasonal dust fluxes and composition from an elevation transect across the Colorado Front Range extending from the urban plains to the remote alpine. At all nine sites, dust was enriched in the essential nutrient phosphorus and the metals copper, zinc, lead, and cadmium, elements that are enriched in dust deposited at sites across the Rocky Mountain West. We observed a seasonal pattern in dust composition, with the highest concentrations of zinc and cadmium during the summer, when back trajectory modeling suggested a greater contribution of dust from local urban and agricultural regions to the east of the collection sites. During the summer, there was also a trend of higher dust fluxes at lower elevations; dust fluxes ranged from 18.9 ± 0.1 g m −2 yr −1 on the plains to 5.9 ± 0.2 g m −2 yr −1 in the alpine. Our results suggest that urban and agricultural land east of the Colorado Front Range is an important source of nutrients and pollutants to all elevations of the mountain range.

Colorado↗
USGS2020

Sea‐level rise will drive divergent sediment transport patterns on fore reefs and reef flats, potentially causing erosion on atoll islands

Atoll reef islands primarily consist of unconsolidated sediment, and their ocean‐facing shorelines are maintained by sediment produced and transported across their reefs. Changes in incident waves can alter cross‐shore sediment exchange and, thus, affect the sediment budget and morphology of atoll reef islands. Here we investigate the influence of sea level rise and projected wave climate change on wave characteristics and cross‐shore sediment transport across an atoll reef at Kwajalein Island, Republic of the Marshall Islands. Using a phase‐resolving model, we quantify the influence on sediment transport of quantities not well captured by wave‐averaged models, namely, wave asymmetry and skewness and flow acceleration. Model results suggest that for current reef geometry, sea level, and wave climate, potential bedload transport is directed onshore, decreases from the fore reef to the beach, and is sensitive to the influence of flow acceleration. We find that a projected 12% decrease in annual wave energy by 2100 CE has negligible influence on reef flat hydrodynamics. However, 0.5–2.0 m of sea level rise increases wave heights, skewness, and shear stress on the reef flat and decreases wave skewness and shear stress on the fore reef. These hydrodynamic changes decrease potential sediment inputs onshore from the fore reef where coral production is greatest but increase potential cross‐reef sediment transport from the outer reef flat to the beach. Assuming sediment production on the fore reef remains constant or decreases due to increasing ocean temperatures and acidification, these processes have the potential to decrease net sediment delivery to atoll islands, causing erosion.

Journal of Geophysical Research – Earth Surface↗
USGS2021

Self-limitation of sand storage in a bedrock-canyon river arising from the interaction of flow and grain size

Bedrock-canyon rivers tend to be supply limited because they are efficient transporters of sediment and not because the upstream supply of sediment is small. A byproduct of this supply limitation is that the finer alluvium stored in these rivers has shorter residence times and smaller volumes than in alluvial rivers. To improve our understanding of disequilibrium sediment transport and its effect on sand storage in bedrock-canyon rivers, we undertook a 20-year study, synthesized herein, of the Colorado River in Grand Canyon. Despite the large loads for which it was renowned, this river exhibited evidence of natural sand-supply limitation and became the perfect natural laboratory for studying sand transport in a bedrock canyon after upstream dam construction exacerbated this supply limitation. During our study, we made and analyzed an unprecedented ∼2.5 million measurements of the suspended and bed sediment. Results indicate that sand storage in this bedrock-canyon river is self-limiting owing to the physical controls of flow and grain size causing negative feedbacks that likely also operate in other bedrock-canyon rivers. Following episodic tributary floods that supply finer sand, sand migrates quickly downstream in the form of a wave in which large systematic changes in bed-sand grain size occur. These grain-size changes cause discharge-independent systematic changes in suspended-sand concentration in excess of a factor of 20. Although the tributary supply of sand increases the amount of sand storage, it also greatly increases the downstream sand transport by causing bed-sand fining, thus limiting the residence time and volume of sand storage.

Arizona↗
USGS2020

Measuring basal force fluctuations of debris flows using seismic recordings and empirical green's functions

We present a novel method for measuring the fluctuating basal normal and shear stresses of debris flows by using along‐channel seismic recordings. Our method couples a simple parameterization of a debris flow as a seismic source with direct measurements of seismic path effects using empirical Green's functions generated with a force hammer. We test this method using two large‐scale (8 and 10 m 3 ) experimental flows at the U.S. Geological Survey debris‐flow flume that were recorded by dozens of three‐component seismic sensors. The seismically derived basal stress fluctuations compare well in amplitude and timing to independent force plate measurements within the valid frequency range (15–50 Hz). We show that although the high‐frequency seismic signals provide band‐limited forcing information, there are systematic relations between the fluctuating stresses and independently measured flow properties, especially mean basal shear stress and flow thickness. However, none of the relationships are simple, and since the flow properties also correlate with one another, we cannot isolate a single factor that relates in a simple way to the fluctuating forces. Nevertheless, our observations, most notably the gradually declining ratio of fluctuating to mean basal stresses during flow passage and the distinctive behavior of the coarse, unsaturated flow front, imply that flow style may be a primary control on the conversion of translational to vibrational kinetic energy. This conversion ultimately controls the radiation of high‐frequency seismic waves. Thus, flow style may provide the key to revealing the nature of the relationship between fluctuating forces and other flow properties.

Journal of Geophysical Research Earth Surface↗
USGS2020

Sediment connectivity: A framework for analyzing coastal sediment transport pathways

Connectivity provides a framework for analyzing coastal sediment transport pathways, building on conceptual advances in graph theory from other scientific disciplines. Connectivity schematizes sediment pathways as a directed graph (i.e., a set of nodes and links). This study presents a novel application of graph theory and connectivity metrics like modularity and centrality to coastal sediment dynamics, exemplified here using Ameland Inlet in the Netherlands. We divide the study site into geomorphic cells (i.e., nodes) and then quantify sediment transport between these cells (i.e., links) using a numerical model. The system of cells and fluxes between them is then schematized in a network described by an adjacency matrix. Network metrics like link density, asymmetry, and modularity quantify system‐wide connectivity. The degree, strength, and centrality of individual nodes identify key locations and pathways throughout the system. For instance, these metrics indicate that under strictly tidal forcing, sand originating near shore predominantly bypasses Ameland Inlet via the inlet channels, whereas sand on the deeper foreshore mainly bypasses the inlet via the outer delta shoals. Connectivity analysis can also inform practical management decisions about where to place sand nourishments, the fate of nourishment sand, or how to monitor locations vulnerable to perturbations. There are still open challenges associated with quantifying connectivity at varying space and time scales and the development of connectivity metrics specific to coastal systems. Nonetheless, connectivity provides a promising technique for predicting the response of our coasts to climate change and the human adaptations it provokes.

Ameland Inlet↗
USGS2020

Linking mesoscale meteorology with extreme landscape response: Effects of narrow cold frontal rainbands (NCFR)

Landscapes evolve in response to prolonged and/or intense precipitation resulting from atmospheric processes at various spatial and temporal scales. Whereas synoptic (large‐scale) features (e.g., atmospheric rivers and hurricanes) govern regional‐scale hydrologic hazards such as widespread flooding, mesoscale features such as thunderstorms or squall lines are more likely to trigger localized geomorphic hazards such as landslides. Thus, to better understand relations between hydrometeorological drivers and landscape response, a knowledge of mesoscale meteorology and its impacts is needed. Here we investigate the extreme geomorphic response associated with one type of mesoscale meteorological feature, the narrow cold frontal rainband (NCFR). Resulting from low‐level convergence and shallow convection along a cold front, NCFRs are narrow bands of high‐intensity rainfall that occur in midlatitude areas of the world. Our study examines an NCFR impacting the Sierra Nevada foothills (California, USA) that initiated over 500 landslides, mobilized ~360,000 metric tons of sediment to the fluvial system (as much as 16 times the local annual sediment yield), and severely damaged local infrastructure and regional water transport facilities. Coupling geomorphological field investigations with meteorological analyses, we demonstrate that precipitation associated with the NCFR was both intense (maximum 15 min intensity of 70 mm/hr) and localized, resulting in a highly concentrated band of shallow landsliding. This meteorological phenomenon likely plays an important role in landscape evolution and hazard initiation. Other types of mesoscale meteorological features also occur globally and offer new avenues for understanding the effects of storms on landscapes.

California↗
USGS2021

Role of future reef growth on morphological response of coral reef islands to sea-level rise

Coral reefs are widely recognised for providing a natural breakwater effect that modulates erosion and flooding hazards on low‐lying sedimentary reef islands. Increased water depth across reef platforms due sea‐level rise (SLR) can compromise this breakwater effect and enhance island exposure to these hazards, but reef accretion in response to SLR may positively contribute to island resilience. Morphodynamic studies suggest that reef islands can adjust to SLR by maintaining freeboard (island crest elevation above still water level) through overwash deposition and island accretion, but the impact of different future reef accretion trajectories on the morphological response of islands remain unknown. Here we show, using a process‐based morphodynamic model, that, although reef growth significantly affects wave transformation processes and island morphology, it does not lead to decreased coastal flooding and island inundation. According to the model, reef islands evolve during SLR by attuning their elevation to the maximum wave runup and islands fronted by a growing reef platform attain lower elevations than those without reef growth, but have similar overwash regimes. The mean overwash discharge Q over across the island crest plays a key role in the ability of islands to keep up with SLR and maintain freeboard, with a Q over value of O (10 l m ‐1 s ‐1 ) separating island construction from destruction. Islands, therefore, can grow vertically to keep up with SLR via flooding and overwash if specific forcing and sediment supply conditions are met, offering hope for uninhabited and sparely populated islands. However, this physical island response will negatively impact infrastructure and assets on developed islands.

Journal of Geophysical Research--Earth Surface↗
USGS2021

Identifying geomorphic process domains in the synthetic landscapes of West Virginia, USA

Human activities such as mining, agriculture, and urbanization have resulted in severe, large‐scale alteration to landform organization and associated geomorphic processes. The mountaintop mining (MTM) region of West Virginia, USA has experienced dramatic topographic alteration, by removing steep slopes and introducing plateau‐like areas at ridgelines and benches on valley fills. The resulting engineered landforms create synthetic landscapes, disconnected from previous geomorphic processes. Invoking the process domain concept, we compare differences in slope‐area relations, cumulative area distributions (CADs), elevation, slope, upslope accumulated area, and a slope*area product before and after mining to adjacent unmined sub‐catchments in five study basins. Differences in the slope‐area relation include a 42% slope reduction in low drainage areas, corresponding to hillslopes, unchanneled valleys, and debris flow dominated channels, which may fall below thresholds required for debris flow processes. The curved slope‐area relation that represents valley incision by debris flows is replaced by slope‐area relations that resemble basins where gullying and the stream power law dominate. Extremely high chemical weathering of unconsolidated valley fills materials may facilitate process domain shifts from debris flows to gullying and fluvial erosion. The characteristic power law scaling break in CADs that represents the headward limit of the channelized network is subdued in post‐mined sites and may reflect headward channelized network extension in mined basins. Slope‐area relations and CADs present a unique topographic signature of MTM activity, potentially providing an analytical approach to assess impacts on underlying geomorphic processes for other synthetic landscapes such as cities or large‐scale agricultural production.

West Virginia↗
USGS2021

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗
USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗
USGS2021

Multi-model comparison of computed debris flow runout for the 9 January 2018 Montecito, California post-wildfire event

Hazard assessment for post-wildfire debris flows, which are common in the steep terrain of the western United States, has focused on the susceptibility of upstream basins to generate debris flows. However, reducing public exposure to this hazard also requires an assessment of hazards in downstream areas that might be inundated during debris flow runout. Debris flow runout models are widely available, but their application to hazard assessment for post-wildfire debris flows has not been extensively tested. Necessary inputs to these models include the total volume of the mobilized flow, flow properties (either inherent material properties or calibration coefficients), and site topography. Estimates of volume are possible in post-event (“back calculation”) studies, yet before an event, volume is an uncertain quantity. We simulated debris flow runout for the well-constrained 9 January 2018 Montecito event using three models (RAMMS, FLO2D, and D-Claw) to determine the relative importance of volume and flow properties. We broke the impacted area into three domains, and for each model-domain combination, we performed a numerical sampling study in which volume and flow properties varied within a wide, but plausible range. We assessed model performance based on inundation patterns and peak flow depths. We found all models could simulate the event with comparable results. Simulation performance was most sensitive to flow volume and less sensitive to flow properties. Our results emphasize the importance of reducing uncertainty in pre-event estimates of flow volume for hazard assessment.

California↗
USGS2022

Multi-stage soil-hydraulic recovery and limited ravel accumulations following the 2017 Nuns and Tubbs wildfires in Northern California

Wildfire can impact soil-hydraulic properties by reducing saturated hydraulic conductivity and sorptivity, making recently burned landscapes prone to debris flows and flash floods. The post-fire hazard window can range from years to decades. In Northern California, where wildfire frequency is steadily increasing, the impact and soil-hydraulic recovery from wildfires is unknown. Following the October 2017 Nuns and Tubbs fires in the Northern Bay Area of California, we established 41 monitoring sites for repeat tension-disc infiltrometer measurements of field-saturated hydraulic conductivity ( K fs ) over 3.5 years. Our site arrays, which encompass grasslands, chaparral, and oak and conifer forests across a range in lithology, show a marked decrease in K fs following the wildfires and a swift partial recovery following the initial post-fire rainy season. Our time series reveals a complex path to soil-hydraulic recovery marked by distinct seasonal stages. Analysis of changing K fs , sorptivity, and infiltration model residuals collectively suggests that these stages are related to transitions between soil-hydraulic processes like structural soil sealing from rainsplash, thermal cracking of bare soil, and vegetation regrowth. While soil infiltration rates were strongly impacted by the 2017 fires, dry ravel estimates are an order of magnitude less for similar slopes than the 2009 Station fire in the San Gabriel mountains of Southern California, suggesting that limited ravel flux may insufficiently load channels for debris flows that initiate from within-channel failure. Our analysis suggests that burned landscapes in the Northern Bay Area of California may experience rapid soil-hydraulic recovery and limited pathways toward post-fire debris flow initiation.

California↗
USGS2022

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface↗