USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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 739 records · Page 41Linked to original sources

Reoccupation of floodplains by rivers and its relation to the age structure of floodplain vegetation

River channel dynamics over many decades provide a physical control on the age structure of floodplain vegetation as a river occupies and abandons locations. Floodplain reoccupation by a river, in particular, determines the interval of time during which vegetation can establish and mature. A general framework for analyzing floodplain reoccupation and a time series model are developed and applied to five alluvial rivers in the United States. Channel dynamics in these rivers demonstrate time-scale dependence with short-term oscillation in active channel area in response to floods and subsequent vegetation growth and progressive lateral movement that accounts for much of the cumulative area occupied by the rivers over decades. Rivers preferentially reoccupy locations recently abandoned causing a decreasing probability of reoccupation with time since abandonment. For a typical case, a river is 10 times more likely to reoccupy an area it abandoned in the past decade than it is to reoccupy an area it abandoned 30 yrs ago. The decreasing probability of reoccupation over time is consistent with observations of persistent stands of late seral stage floodplain forest. A power function provides a robust approach for estimating the cumulative area occupied by a river and the age structure of riparian forests resulting from a specific historical sequence of streamflow in comparison to either linear or exponential alternatives.

Journal of Geophysical Research↗

Dissolved organic carbon and chromophoric dissolved organic matter properties of rivers in the USA

Dissolved organic carbon (DOC) concentration and chromophoric dissolved organic matter (CDOM) parameters were measured over a range of discharge in 30 U.S. rivers, covering a diverse assortment of fluvial ecosystems in terms of watershed size and landscape drained. Relationships between CDOM absorption at a range of wavelengths (a 254 , a 350 , a 440 ) and DOC in the 30 watersheds were found to correlate strongly and positively for the majority of U.S. rivers. However, four rivers (Colorado, Colombia, Rio Grande and St. Lawrence) exhibited statistically weak relationships between CDOM absorption and DOC. These four rivers are atypical, as they either drain from the Great Lakes or experience significant impoundment of water within their watersheds, and they exhibited values for dissolved organic matter (DOM) parameters indicative of autochthonous or anthropogenic sources or photochemically degraded allochthonous DOM and thus a decoupling between CDOM and DOC. CDOM quality parameters in the 30 rivers were found to be strongly correlated to DOM compositional metrics derived via XAD fractionation, highlighting the potential for examining DOM biochemical quality from CDOM measurements. This study establishes the ability to derive DOC concentration from CDOM absorption for the majority of U.S. rivers, describes characteristics of riverine systems where such an approach is not valid, and emphasizes the possibility of examining DOM composition and thus biogeochemical function via CDOM parameters. Therefore, the usefulness of CDOM measurements, both laboratory-based analyses and in situ instrumentation, for improving spatial and temporal resolution of DOC fluxes and DOM dynamics in future studies is considerable in a range of biogeochemical studies.

Journal of Geophysical Research G: Biogeosciences↗

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research↗

Problem of the Love‐Gannon relation between the asymmetric disturbance field and Dst

Love and Gannon (2009) discovered that statistically, over a fifty year period the difference in the dawn and dusk disturbance‐field H component at low latitudes (hourly averaged) is linearly proportional to Dst. If the difference is designated by δ DD in units of nT/R E , then the Love‐Gannon (L‐G) relation is δ DD = − 0.2 Dst. At any time departures from the relation can be large. Nonetheless, the relation is evident for all values of Dst and persists throughout magnetic storms, both the main phase and the recovery phase. The Love‐Gannon discovery presents a problem to current understanding of the relation between the causes of δ DD and Dst because the dawn dusk asymmetry in the disturbance field is presumably governed by a long‐established magnetosphere‐ionosphere coupling theory which predicts a characteristic time scale (the shielding time) of less than an hour whereas the characteristic time scale for Dst (the ring current decay time) is more like ten hours. Thus, without forcing both time scales toward each other to the limits of their ranges, a linear proportionality between δ DD and Dst cannot be derived from the current understanding of the causes of the asymmetry and the ring current. This conclusion is the paper's main contribution. In addition, we attempt to get around the conflict of time scales by looking at other possibilities for generating δ DD that depend directly on the ring current. The most promising of these is the possibility that the ring current decay mechanism creates a quasi‐permanent, local‐time modification of the ring current compared to what it would be in the absence of the decay mechanism and that this modification causes a field‐aligned current that closes through the ionosphere and generates the asymmetry δ DD . This idea has the virtue of coupling the asymmetry directly to the ring current and of accounting for the persistence of the L‐G proportionality through the recovery phase of magnetic storms.

Journal of Geophysical Research A: Space Physics↗

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema&lsquo;uma&lsquo;u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema&lsquo;uma&lsquo;u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 &mu; m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 &mu; m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗

Clustering of GPS velocities in the Mojave Block, southeastern California

We find subdivisions within the Mojave Block using cluster analysis to identify groupings in the velocities observed at GPS stations there. The clusters are represented on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. The most significant representation as judged by the gap test involves 4 clusters within the Mojave Block. The fault systems bounding the clusters from east to west are 1) the faults defining the eastern boundary of the Northeast Mojave Domain extended southward to connect to the Hector Mine rupture, 2) the Calico-Paradise fault system, 3) the Landers-Blackwater fault system, and 4) the Helendale-Lockhart fault system. This division of the Mojave Block is very similar to that proposed by Meade and Hager. However, no cluster boundary coincides with the Garlock Fault, the northern boundary of the Mojave Block. Rather, the clusters appear to continue without interruption from the Mojave Block north into the southern Walker Lane Belt, similar to the continuity across the Garlock Fault of the shear zone along the Blackwater-Little Lake fault system observed by Peltzer et al. Mapped traces of individual faults in the Mojave Block terminate within the block and do not continue across the Garlock Fault [Dokka and Travis, ].

California↗

Currents, drag, and sediment transport induced by a tsunami

We report observations of water surface elevation, currents, and suspended sediment concentration (SSC) from a 10-m deep site on the inner shelf in northern Monterey Bay during the arrival of the 2010 Chile tsunami. Velocity profiles were measured from 3.5 m above the bed (mab) to the surface at 2 min intervals, and from 0.1 to 0.7 mab at 1 Hz. SSC was determined from the acoustic backscatter of the near-bed profiler. The initial tsunami waves were directed cross shore and had a period of approximately 16 min. Maximum wave height was 1.1 m, and maximum current speed was 0.36 m/s. During the strongest onrush, near-bed velocities were clearly influenced by friction and a logarithmic boundary layer developed, extending more than 0.3 mab. We estimated friction velocity and bed shear stress from the logarithmic profiles. The logarithmic structure indicates that the flow can be characterized as quasi-steady at these times. At other phases of the tsunami waves, the magnitude of the acceleration term was significant in the near-bed momentum equation, indicating unsteady flow. The maximum tsunami-induced bed shear stress (0.4 N/m 2 ) exceeded the critical shear stress for the medium-grained sand on the seafloor. Cross-shore sediment flux was enhanced by the tsunami. Oscillations of water surface elevation and currents continued for several days. The oscillations were dominated by resonant frequencies, the most energetic of which was the fundamental longitudinal frequency of Monterey Bay. The maximum current speed (hourly-timescale) in 18 months of observations occurred four hours after the tsunami arrived.

California↗

Bathymetric controls on sediment transport in the Hudson River estuary: Lateral asymmetry and frontal trapping

Analyses of field observations and numerical model results have identified that sediment transport in the Hudson River estuary is laterally segregated between channel and shoals, features frontal trapping at multiple locations along the estuary, and varies significantly over the spring-neap tidal cycle. Lateral gradients in depth, and therefore baroclinic pressure gradient and stratification, control the lateral distribution of sediment transport. Within the saline estuary, sediment fluxes are strongly landward in the channel and seaward on the shoals. At multiple locations, bottom salinity fronts form at bathymetric transitions in width or depth. Sediment convergences near the fronts create local maxima in suspended-sediment concentration and deposition, providing a general mechanism for creation of secondary estuarine turbidity maxima at bathymetric transitions. The lateral bathymetry also affects the spring-neap cycle of sediment suspension and deposition. In regions with broad, shallow shoals, the shoals are erosional and the channel is depositional during neap tides, with the opposite pattern during spring tides. Narrower, deeper shoals are depositional during neaps and erosional during springs. In each case, the lateral transfer is from regions of higher to lower bed stress, and depends on the elevation of the pycnocline relative to the bed. Collectively, the results indicate that lateral and along-channel gradients in bathymetry and thus stratification, bed stress, and sediment flux lead to an unsteady, heterogeneous distribution of sediment transport and trapping along the estuary rather than trapping solely at a turbidity maximum at the limit of the salinity intrusion.

Journal of Geophysical Research↗

Along-shelf current variability on the Catalan inner-shelf (NW Mediterranean)

We examine the circulation over the inner shelf of the Catalan Sea using observations of currents obtained from three ADCPs within the inner-shelf (24 and 50 m depth) during March-April 2011. The along-shelf current fluctuations during that period are mainly controlled by the local wind stress on short time scales and by remote pressure gradients on synoptic time scales. Different forcing mechanisms are involved in the along-shelf momentum balance. During storm conditions, wind stress, sea level gradients and the non-linear terms dominate the balance. During weak wind conditions, the momentum balance is controlled by the pressure gradient, while during periods of moderate wind in the presence of considerable stratification, the balance is established between the Coriolis and wind stress terms. Vertical variations of velocity are affected by the strong observed density gradient. The increased vertical shear is accompanied by the development of stratified conditions due to local heating when the wind is not able to counteract (and destroy) stratification. The occasional influence of the Besòs river plume is observed in time scales of hours to days in a limited area in front of Barcelona. The area affected by the plume depends on the vertical extend of the fresher layer, the fast river discharge peak, and the relaxation of cross-shore velocities after northeast storm events. This contribution provides a first interpretation of the inner-shelf dynamics in the Catalan Sea.

Journal of Geophysical Research↗

Spatial distribution of damage around faults in the Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah: A mechanical analog for faulting in pyroclastic deposits on Mars

Volcanic ash is thought to comprise a large fraction of the Martian equatorial layered deposits and much new insight into the process of faulting and related fluid flow in these deposits can be gained through the study of analogous terrestrial tuffs. This study identifies a set of fault-related processes that are pertinent to understanding the evolution of fault systems in fine-grained, poorly indurated volcanic ash by investigating exposures of faults in the Miocene-aged Joe Lott Tuff Member of the Mount Belknap Volcanics, Utah. The porosity and granularity of the host rock are found to control the style of localized strain that occurs prior to and contemporaneous with faulting. Deformation bands occur in tuff that was porous and granular at the time of deformation, while fractures formed where the tuff lost its porous and granular nature due to silicic alteration. Non-localized deformation of the host rock is also prominent and occurs through compaction of void space, including crushing of pumice clasts. Significant off-fault damage of the host rock, resembling fault pulverization, is recognized adjacent to one analog fault and may reflect the strain rate dependence of the resulting fault zone architecture. These findings provide important new guidelines for future structural analyses and numerical modeling of faulting and subsurface fluid flow through volcanic ash deposits on Mars.

Utah↗

Soil genesis on the island of Bermuda in the Quaternary: the importance of African dust transport and deposition

The origin of terra rossa, red or reddish-brown, clay-rich soils overlying high-purity carbonate substrates, has intrigued geologists and pedologists for decades. Terra rossa soils can form from accumulation of insoluble residues during dissolution of the host limestones, addition of volcanic ash, or addition of externally derived, long-range-transported (LRT) aeolian particles. We studied soils and paleosols on high-purity, carbonate aeolianites of Quaternary age on Bermuda, where terra rossa origins have been debated for more than a century. Potential soil parent materials on this island include sand-sized fragments of local volcanic bedrock, the LRT, fine-grained (<20 &mu;m) component of distal loess from the lower Mississippi River Valley, and LRT dust from Africa. These parent materials can be characterized geochemically using trace elements that are immobile in the soil-forming environment. Results indicate that local volcanic bedrock on Bermuda has Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that can be distinguished from African dust and lower Mississippi River valley loess. Bermuda soils have Sc-Th-La, Cr-Ta-Nd, and Eu/Eu*, La N /Yb N , Gd N /Yb N that indicate derivation from a combination of LRT dust from Africa and local volcanic bedrock. Our results indicate that soils on islands in a very broad latitudinal belt of the western Atlantic margin have been influenced by African LRT dust inputs over much of the past &ndash;500 ka.

Journal of Geophysical Research↗

Nitrate removal in deep sediments of a nitrogen-rich river network: A test of a conceptual model

Many estimates of nitrogen removal in streams and watersheds do not include or account for nitrate removal in deep sediments, particularly in gaining streams. We developed and tested a conceptual model for nitrate removal in deep sediments in a nitrogen-rich river network. The model predicts that oxic, nitrate-rich groundwater will become depleted in nitrate as groundwater upwelling through sediments encounters a zone that contains buried particulate organic carbon, which promotes redox conditions favorable for nitrate removal. We tested the model at eight sites in upwelling reaches of lotic ecosystems in the Waupaca River Watershed that varied by three orders of magnitude in groundwater nitrate concentration. We measured denitrification potential in sediment core sections to 30 cm and developed vertical nitrate profiles to a depth of about 1 m with peepers and piezometer nests. Denitrification potential was higher, on average, in shallower core sections. However, core sections deeper than 5 cm accounted for 70%, on average, of the depth-integrated denitrification potential. Denitrification potential increased linearly with groundwater nitrate concentration up to 2 mg NO 3 -N/L but the relationship broke down at higher concentrations (> 5 mg NO 3 -N/L), a pattern that suggests nitrate saturation. At most sites groundwater nitrate declined from high concentrations at depth to much lower concentrations prior to discharge into the surface water. The profiles suggested that nitrate removal occurred at sediment depths between 20 and 40 cm. Dissolved oxygen concentrations were much higher in deep sediments than in pore water at 5 cm sediment depth at most locations. The substantial denitrification potential in deep sediments coupled with the declines in nitrate and dissolved oxygen concentrations in upwelling groundwater suggest that our conceptual model for nitrate removal in deep sediments is applicable to this river network. Our results suggest that nitrate removal rates can be high in deep sediments of upwelling stream reaches, which may have implications for efforts to understand and quantify nitrogen transport and removal at larger scales.

Journal of Geophysical Research: Biogeosciences↗

Climatic limits on foliar growth during major droughts in the Southwestern U.S.A.

Pronounced droughts during the 1950s and 2000s in the Southwestern U.S.A. (SW) provide an opportunity to compare mesoscale ecosystem responses to anomalously dry conditions before and during the regional warming that started in the late 1970s. This year-round warming has produced fewer cool season freezes, losses in regional snowpack, an 8-10 day advance in spring onset, and hotter summers, all of which should affect vegetation differently across seasons and elevations. Here, we examine indices that represent climatic limits on foliar growth for both drought periods, and evaluate these indices for areas that experienced tree mortality during the 2000s drought. Relative to the 1950s drought, warmer conditions during the 2000s drought decreased the occurrence of temperatures too low for foliar growth at lower elevations in winter and higher elevations in summer. Higher vapor pressure deficits (VPDs) largely driven by warmer temperatures in the more recent drought were more limiting to foliar growth from spring through summer at lower and middle elevations. At many locations where tree mortality occurred during the 2000s drought, low-temperature constraints on foliar growth were extremely unlimiting, whereas VPD constraints were extremely limiting from early spring through late autumn. Our analysis shows that in physiographically complex regions like the SW, seasonality and elevational gradients are important for understanding vegetative responses to warming. It also suggests that continued warming will increase the degree to which VPD limits foliar growth during future droughts, and expand its reach to higher elevations and other seasons.

Journal of Geophysical Research G: Biogeosciences↗

Quantifying soil surface change in degraded drylands: shrub encroachment and effects of fire and vegetation removal in a desert grassland

Woody plant encroachment, a worldwide phenomenon, is a major driver of land degradation in desert grasslands. Woody plant encroachment by shrub functional types ultimately leads to the formation of a patchy landscape with fertile shrub patches interspaced with nutrient-depleted bare soil patches. This is considered to be an irreversible process of land and soil degradation. Recent studies have indicated that in the early stages of shrub encroachment, when there is sufficient herbaceous connectivity, fires (prescribed or natural) might provide some reversibility to the shrub encroachment process by negatively affecting shrub demography and homogenizing soil resources across patches within weeks to months after burning. A comprehensive understanding of longer term changes in microtopography and spatial patterning of soil properties following fire in shrub-encroached grasslands is desirable. Here, we investigate the changes in microtopography with LiDAR (light detection and ranging), vegetation recovery, and spatial pattering of soil properties in replicated burned, clipped, and control areas in a shrub-grass transition zone in the northern Chihuahuan Desert four years after prescribed fire or clipping. Results indicate a greater homogeneity in soil, microtopography, and vegetation patterning on burned relative to clipped and control treatments. Findings provide further evidence that disturbance by prescribed fire may allow for reversal of the shrub encroachment process, if the event occurs in the early stages of the vegetation shift. Improved understanding of longer-term effects of fire and associated changes in soil patterning can inform the use and role of fire in the context of changing disturbance regimes and climate.

Chihuahuan Desert↗