SEARCH · USGS Science

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46Linked to original sources

USGS2003

Experimental and textural constraints on mafic enclave formation in volcanic rocks

We have used experiments and textural analysis to investigate the process of enclave formation during magma mixing at Southwest Trident volcano, Alaska. Andesite enclaves are present throughout the four dacite lava flows produced by the eruption, and resemble mafic enclaves commonly found in other volcanic rocks. Our experiments replicate the pressure–temperature–time path taken by enclave-forming andesite magma as it is engulfed in dacite during magma mixing. Pressure and temperature information for the andesite and dacite are from [Coombs et al., Contrib. Mineral. Petrol. 140 (2000) 99–118]. The andesite was annealed at 1000°C, and then cooled to 890°C at rates of 110°C h −1 , 10°C h −1 , and 2°C h −1 . Once cooled to 890°C, andesite was held at this lower temperature from times ranging from 1 to 40 h. The andesite that was cooled at the slower rates of 2°C h −1 and 10°C h −1 most resembles enclave groundmass texturally and compositionally. Based on simple thermal calculations, these rates are more consistent with cooling of the andesite groundmass below an andesite–dacite interface than with cooling of enclave-sized spheres. If enclaves do crystallize as spheres, post-crystallization disaggregation must occur. Calculations using the MELTS algorithm [Ghiorso and Sack, Contrib. Mineral. Petrol. 119 (1995) 197–212] show that for incoming andesite to become less dense than the dacite ∼34 volume % of its groundmass must crystallize to undergo ∼18 volume % vesiculation; these values are similar to those determined for Southwest Trident enclaves. Thus such crystallization may lead to ‘flotation’ of enclaves and be a viable mechanism for enclave formation and dispersal. The residual melt in the cooling experiments did not evolve to rhyolitic compositions such as seen in natural enclaves due to a lack of a decompression step in the experiments. Decompression experiments on Southwest Trident dacite suggest an average ascent rate for the eruption of ∼2–3 MPa h −1 . An andesite experiment that was cooled and then decompressed at this rate contains melt that matches that of the natural enclaves. It is apparent that decompression (ascent)-induced crystallization occurs in enclaves, but not in the form of microlites as happens in the dacite host, due either to insufficient residence time at chamber temperatures or to the pre-existing microphenocrysts which act as sites for new growth.

Alaska↗

USGS2006

Chapter 13 Petrogenesis of the Campanian Ignimbrite: Implications for crystal-melt separation and open-system processes from major and trace elements and Th isotopic data

The Campanian Ignimbrite is a large-volume trachytic to phonolitic ignimbrite that was deposited at ≈39.3 ka and represents one of a number of highly explosive volcanic events that have occurred in the region near Naples, Italy. Thermodynamic modeling using the MELTS algorithm reveals that major element variations are dominated by crystal-liquid separation at 0.15 GPa. Initial dissolved H 2 O content in the parental melt is ∼3 wt.% and the magmatic system fugacity of oxygen was buffered along QFM+1. Significantly, MELTS results also indicate that the liquid line of descent is marked by a large change in the proportion of melt (from 0.46 to 0.09) at ∼884°C, which leads to a discontinuity in melt composition (i.e., a compositional gap) and different thermodynamic and transport properties of melt and magma across the gap. Crystallization of alkali feldspar and plagioclase dominates the phase assemblage at this pseudo-invariant point temperature of ∼884°C. Evaluation of the variations in the trace elements Zr, Nb, Th, U, Rb, Sm, and Sr using a mass balance equation that accounts for changing bulk mineral-melt partition coefficients as crystallization occurs indicates that crystal-liquid separation and open-system processes were important. Th isotope data yield an apparent isochron that is ∼20 kyr younger than the age of the deposit, and age-corrected Th isotope data indicate that the magma body was an open system at the time of eruption. Because open-system behavior can profoundly change isotopic and elemental characteristics of a magma body, these Th results illustrate that it is critical to understand the contribution that open-system processes make to magmatic systems prior to assigning relevance to age or timescale information derived from such systems. Fluid-magma interaction has been proposed as a mechanism to change isotopic and elemental characteristics of magma bodies, but an evaluation of the mass and thermal constraints on such a process suggests large-scale interaction is unlikely. In the case of the magma body associated with the Campanian Ignimbrite, the most likely source of the open-system signatures is assimilation of partial melts of compositionally heterogeneous basement composed of cumulates and intrusive equivalents of volcanic activity that has characterized the Campanian region for over 300 kyr.

Developments in Volcanology↗

USGS2008

Modeling potential habitats for alien species Dreissena polymorpha in continental USA

The effective measure to minimize the damage of invasive species is to block the potential invasive species to enter into suitable areas. 1864 occurrence points with GPS coordinates and 34 environmental variables from Daymet datasets were gathered, and 4 modeling methods, i.e., Logistic Regression (LR), Classification and Regression Trees (CART), Genetic Algorithm for Rule-Set Prediction (GARP), and maximum entropy method (Maxent), were introduced to generate potential geographic distributions for invasive species Dreissena polymorpha in Continental USA. Then 3 statistical criteria of the area under the Receiver Operating Characteristic curve (AUC), Pearson correlation (COR) and Kappa value were calculated to evaluate the performance of the models, followed by analyses on major contribution variables. Results showed that in terms of the 3 statistical criteria, the prediction results of the 4 ecological niche models were either excellent or outstanding, in which Maxent outperformed the others in 3 aspects of predicting current distribution habitats, selecting major contribution factors, and quantifying the influence of environmental variables on habitats. Distance to water, elevation, frequency of precipitation and solar radiation were 4 environmental forcing factors. The method suggested in the paper can have some reference meaning for modeling habitats of alien species in China and provide a direction to prevent Mytilopsis sallei on the Chinese coast line.

Acta Ecologica Sinica↗

USGS2022

The future of coastal monitoring through satellite remote sensing

Satellite remote sensing is transforming coastal science from a “data-poor” field into a “data-rich” field. Sandy beaches are dynamic landscapes that change in response to long-term pressures, short-term pulses, and anthropogenic interventions. Until recently, the rate and breadth of beach change have outpaced our ability to monitor those changes, due to the spatiotemporal limitations of our observational capacity. Over the past several decades, only a handful of beaches worldwide have been regularly monitored with accurate yet expensive in situ surveys. The long-term coastal-change data of these few well-monitored beaches have led to in-depth understanding of many site-specific coastal processes. However, because the best-monitored beaches are not representative of all beaches, much remains unknown about the processes and fate of the other >99% of unmonitored beaches worldwide. The fleet of Earth-observing satellites has enabled multiscale monitoring of beaches, for the very first time, by providing imagery with global coverage and up to daily frequency. The long-standing and ever-expanding archive of satellite imagery will enable coastal scientists to investigate coastal change at sites vulnerable to future sea-level rise, that is, (almost) everywhere. In the past decade, our capability to observe coastal change from space has grown substantially with computing and algorithmic power. Yet, further advances are needed in automating monitoring using machine learning, deep learning, and computer vision to fully leverage this massive treasure trove of data. Extensive monitoring and investigation of the causes and effects of coastal change at the requisite spatiotemporal scales will provide coastal managers with additional, valuable information to evaluate problems and solutions, addressing the potential for widespread beach loss due to accelerated sea-level rise, development, and reduced sediment supply. Monitoring from Earth-observing satellites is currently the only means of providing seamless data with high spatiotemporal resolution at the global scale of the impending impacts of climate change on coastal systems.

Cambridge Prisms: Coastal Futures↗

USGS2022

Assessing reproducibility in sedimentary macroscopic charcoal count data

Current understanding of global late Quaternary fire history is largely drawn from sedimentary charcoal data. Since publication, CharAnalysis increasingly has been relied upon as a robust method for analyzing these data. However, several underlying assumptions of the algorithm have not been tested. This study uses replicated charcoal count data to examine the assumption of Poisson distribution and reproducibility of peak detection. Results show <10% of the replicate counts are Poisson distributed, a maximum peak replication rate of 60%, and, for >90% of the data, intra-level count differences were larger than the threshold used to identify significance in inter-level differences. A pronounced “edge effect” was observed at the beginning and end of the records, cautioning against validation of results based on sections corresponding to the historical period. The proximal cause for low reproducibility is likely a lack of spatial randomness of charcoal particles at the scale of a core diameter. Until and unless decomposition methods can be developed that accommodate the observed limitations inherent in particle count data, best practices for interpreting charcoal records may be to rely on qualitative interpretations based on smoothed influx values and minimum particle count values in the hundreds.

Quaternary Research↗

USGS2011

Tsunami modelling with adaptively refined finite volume methods

Numerical modelling of transoceanic tsunami propagation, together with the detailed modelling of inundation of small-scale coastal regions, poses a number of algorithmic challenges. The depth-averaged shallow water equations can be used to reduce this to a time-dependent problem in two space dimensions, but even so it is crucial to use adaptive mesh refinement in order to efficiently handle the vast differences in spatial scales. This must be done in a 'wellbalanced' manner that accurately captures very small perturbations to the steady state of the ocean at rest. Inundation can be modelled by allowing cells to dynamically change from dry to wet, but this must also be done carefully near refinement boundaries. We discuss these issues in the context of Riemann-solver-based finite volume methods for tsunami modelling. Several examples are presented using the GeoClaw software, and sample codes are available to accompany the paper. The techniques discussed also apply to a variety of other geophysical flows. ?? 2011 Cambridge University Press.

Acta Numerica↗

USGS2015

Photometric properties of Vesta

The Dawn spacecraft orbited Asteroid (4) Vesta for a year, and returned disk-resolved images and spectra covering visible and near-infrared wavelengths at scales as high as 20 m/pix. The visible geometric albedo of Vesta is ~ 0.36. The disk-integrated phase function of Vesta in the visible wavelengths derived from Dawn approach data, previous ground-based observations, and Rosetta OSIRIS observations is consistent with an IAU H-G phase law with H=3.2 mag and G=0.28. Hapke's modeling yields a disk-averaged single-scattering albedo of 0.50, an asymmetry factor of -0.25, and a roughness parameter of ~20 deg at 700 nm wavelength. Vesta's surface displays the largest albedo variations observed so far on asteroids, ranging from ~0.10 to ~0.76 in geometric albedo in the visible wavelengths. The phase function of Vesta displays obvious systematic variations with respect to wavelength, with steeper slopes within the 1- and 2-micron pyroxene bands, consistent with previous ground-based observations and laboratory measurement of HED meteorites showing deeper bands at higher phase angles. The relatively high albedo of Vesta suggests significant contribution of multiple scattering. The non-linear effect of multiple scattering and the possible systematic variations of phase function with albedo across the surface of Vesta may invalidate the traditional algorithm of applying photometric correction on airless planetary surfaces.

Proceedings of the International Astronomical Unio↗

USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗

USGS2023

Wide-area debris field and seabed characterization of a deep ocean dump site surveyed by autonomous underwater vehicles

Disposal of industrial and hazardous waste in the deep ocean was a pervasive global practice near coastlines in the 20th century. Uncertainty in the quantity, location, and contents of dumped materials at historical disposal sites underscores ongoing risks to marine ecosystems and human health. This study presents analysis of a 150-km2 wide-area sidescan sonar survey conducted in March 2021 with two autonomous underwater vehicles (AUVs) at an offshore dumpsite in San Pedro Basin, California. Previous camera surveys located 60 barrels and other debris; sediment analysis in the region showed varying concentrations of the insecticidal chemical Dichlorodiphenyl-trichloroethane (DDT), for which an estimated 350-700 metric tons was discarded in the San Pedro Basin between 1947 and 1961. A lack of primary historical documents from Montrose Chemical or California Salvage Company that specify the DDT acid waste disposal method has contributed to the ambiguity surrounding whether the dumping method was via bulk discharge or containerized units. Barrels and debris visually detected during previous surveys by Remotely Operated Vehicles (ROVs) Jason and SuBastian and AUV Sentry served as the basis for ground truth testing and object classification using size and acoustic intensity characteristics. Target detection and classification algorithms were developed using image and signal processing techniques that resulted in the identification of over 74,000 debris targets within the survey region. Statistical, spectral, and machine learning techniques are used to characterize seabed variability and for bottom-type classification. These analytical techniques combined with the survey capabilities of AUVs provide a framework for efficient mapping and characterization of uncharted deep-water disposal sites.

California↗

USGS2026

Microplastics undergo fragmentation during pressurized membrane filtration

Low-micrometer microplastics (<10 μm) have been detected in drinking water, driving growing interest in using pressure-driven membrane filtration to remove these particles and ensure drinking water safety. However, little attention is paid to the fate of microplastics concentrated in the filtration byproduct. In this study, using well-defined polystyrene (PS) and poly(methyl methacrylate) microspheres as model particles, we observed microplastic fragmentation during nanofiltration and the subsequent release of smaller fragments. After operating for 3 h at 20 bar, 67.9 ± 8.0% of the PS spheres in the concentrate, based on particles counted in selected fields of view, were transformed into fractured particles. The characteristic Raman band signals of microplastics were significantly weakened after fragmentation, leading to detection challenges. To address this, a data processing algorithm was developed to identify PS fragments as small as 1 μm. Preliminary statistical analysis reveals that within the tested pressure range, operation time has a stronger influence on fragmentation than pressure magnitude alone and that fragmentation is further governed by the intrinsic mechanical properties of the tested model polymers. This work provides direct evidence that pressure-driven membrane processes induce microplastic fragmentation and highlights the environmental risks associated with the discharge of fragmented microplastics into the concentrate.

Environmental Science & Technology Letters↗

USGS2004

Turbulence effects on volatilization rates of liquids and solutes

Volatilization rates of neat liquids (benzene, toluene, fluorobenzene, bromobenzene, ethylbenzene, m-xylene, o-xylene, o-dichlorobenzene, and 1-methylnaphthalene) and of solutes (phenol, m-cresol, benzene, toluene, ethylbenzene, o-xylene, and ethylene dibromide) from dilute water solutions have been measured in the laboratory over a wide range of air speeds and water-stirring rates. The overall transfer coefficients (KL) for individual solutes are independent of whether they are in single- or multi-solute solutions. The gas-film transfer coefficients (kG) for solutes in the two-film model, which have hitherto been estimated by extrapolation from reference coefficients, can now be determined directly from the volatilization rates of neatliquids through anew algorithm. The associated liquid-film transfer coefficients (KL) can then be obtained from measured KL and kG values and solute Henry law constants (H). This approach provides a novel means for checking the precision of any kL and kG estimation methods for ultimate prediction of KL. The improved kG estimation enables accurate K L predictions for low-volatility (i.e., low-H) solutes where K L and kGH are essentially equal. In addition, the prediction of KL values for high-volatility (i.e., high-H) solutes, where KL ??? kL, is also improved by using appropriate reference kL values.

Environmental Science & Technology↗

USGS2010

Measurement of dissolved organic matter fluorescense in aquatic environments: An interlaboratory comparison

The fluorescent properties of dissolved organic matter (DOM) are often studied in order to infer DOM characteristics in aquatic environments, including source, quantity, composition, and behavior. While a potentially powerful technique, a single widely implemented standard method for correcting and presenting fluorescence measurements is lacking, leading to difficulties when comparing data collected by different research groups. This paper reports on a large-scale interlaboratory comparison in which natural samples and well-characterized fluorophores were analyzed in 20 laboratories in the U.S., Europe, and Australia. Shortcomings were evident in several areas, including data quality-assurance, the accuracy of spectral correction factors used to correct EEMs, and the treatment of optically dense samples. Data corrected by participants according to individual laboratory procedures were more variable than when corrected under a standard protocol. Wavelength dependency in measurement precision and accuracy were observed within and between instruments, even in corrected data. In an effort to reduce future occurrences of similar problems, algorithms for correcting and calibrating EEMs are described in detail, and MATLAB scripts for implementing the study's protocol are provided. Combined with the recent expansion of spectral fluorescence standards, this approach will serve to increase the intercomparability of DOM fluorescence studies.

Environmental Science & Technology↗

USGS2002

A statistical evaluation of non-ergodic variogram estimators

Geostatistics is a set of statistical techniques that is increasingly used to characterize spatial dependence in spatially referenced ecological data. A common feature of geostatistics is predicting values at unsampled locations from nearby samples using the kriging algorithm. Modeling spatial dependence in sampled data is necessary before kriging and is usually accomplished with the variogram and its traditional estimator. Other types of estimators, known as non-ergodic estimators, have been used in ecological applications. Non-ergodic estimators were originally suggested as a method of choice when sampled data are preferentially located and exhibit a skewed frequency distribution. Preferentially located samples can occur, for example, when areas with high values are sampled more intensely than other areas. In earlier studies the visual appearance of variograms from traditional and non-ergodic estimators were compared. Here we evaluate the estimators' relative performance in prediction. We also show algebraically that a non-ergodic version of the variogram is equivalent to the traditional variogram estimator. Simulations, designed to investigate the effects of data skewness and preferential sampling on variogram estimation and kriging, showed the traditional variogram estimator outperforms the non-ergodic estimators under these conditions. We also analyzed data on carabid beetle abundance, which exhibited large-scale spatial variability (trend) and a skewed frequency distribution. Detrending data followed by robust estimation of the residual variogram is demonstrated to be a successful alternative to the non-ergodic approach.

Environmental and Ecological Statistics↗

USGS2003

Preliminary investigation of submerged aquatic vegetation mapping using hyperspectral remote sensing

The use of airborne hyperspectral remote sensing imagery for automated mapping of submerged aquatic vegetation (SAV) in the tidal Potomac River was investigated for near to realtime resource assessment and monitoring. Airborne hyperspectral imagery and field spectrometer measurements were obtained in October of 2000. A spectral library database containing selected ground-based and airborne sensor spectra was developed for use in image processing. The spectral library is used to automate the processing of hyperspectral imagery for potential real-time material identification and mapping. Field based spectra were compared to the airborne imagery using the database to identify and map two species of SAV (Myriophyllum spicatum and Vallisneria americana). Overall accuracy of the vegetation maps derived from hyperspectral imagery was determined by comparison to a product that combined aerial photography and field based sampling at the end of the SAV growing season. The algorithms and databases developed in this study will be useful with the current and forthcoming space-based hyperspectral remote sensing systems.

Environmental Monitoring and Assessment↗

USGS2003

A hydrologic network supporting spatially referenced regression modeling in the Chesapeake Bay watershed

The U.S. Geological Survey has developed a methodology for statistically relating nutrient sources and land-surface characteristics to nutrient loads of streams. The methodology is referred to as SPAtially Referenced Regressions On Watershed attributes (SPARROW), and relates measured stream nutrient loads to nutrient sources using nonlinear statistical regression models. A spatially detailed digital hydrologic network of stream reaches, stream-reach characteristics such as mean streamflow, water velocity, reach length, and travel time, and their associated watersheds supports the regression models. This network serves as the primary framework for spatially referencing potential nutrient source information such as atmospheric deposition, septic systems, point-sources, land use, land cover, and agricultural sources and land-surface characteristics such as land use, land cover, average-annual precipitation and temperature, slope, and soil permeability. In the Chesapeake Bay watershed that covers parts of Delaware, Maryland, Pennsylvania, New York, Virginia, West Virginia, and Washington D.C., SPARROW was used to generate models estimating loads of total nitrogen and total phosphorus representing 1987 and 1992 land-surface conditions. The 1987 models used a hydrologic network derived from an enhanced version of the U.S. Environmental Protection Agency's digital River Reach File, and course resolution Digital Elevation Models (DEMs). A new hydrologic network was created to support the 1992 models by generating stream reaches representing surface-water pathways defined by flow direction and flow accumulation algorithms from higher resolution DEMs. On a reach-by-reach basis, stream reach characteristics essential to the modeling were transferred to the newly generated pathways or reaches from the enhanced River Reach File used to support the 1987 models. To complete the new network, watersheds for each reach were generated using the direction of surface-water flow derived from the DEMs. This network improves upon existing digital stream data by increasing the level of spatial detail and providing consistency between the reach locations and topography. The hydrologic network also aids in illustrating the spatial patterns of predicted nutrient loads and sources contributed locally to each stream, and the percentages of nutrient load that reach Chesapeake Bay.

Environmental Monitoring and Assessment↗

USGS2004

Sampling issues affecting accuracy of likelihood-based classification using genetical data

We demonstrate the effectiveness of a genetic algorithm for discovering multi-locus combinations that provide accurate individual assignment decisions and estimates of mixture composition based on likelihood classification. Using simulated data representing different levels of inter-population differentiation (F st ~ 0.01 and 0.10), genetic diversities (four or eight alleles per locus), and population sizes (20, 40, 100 individuals in baseline populations), we show that subsets of loci can be identified that provide comparable levels of accuracy in classification decisions relative to entire multi-locus data sets, where 5, 10, or 20 loci were considered. Microsatellite data sets from hatchery strains of lake trout, Salvelinus namaycush, representing a comparable range of inter-population levels of differentiation in allele frequencies confirmed simulation results. For both simulated and empirical data sets, assignment accuracy was achieved using fewer loci (e.g., three or four loci out of eight for empirical lake trout studies). Simulation results were used to investigate properties of the 'leave-one-out' (L1O) method for estimating assignment error rates. Accuracy of population assignments based on L1O methods should be viewed with caution under certain conditions, particularly when baseline population sample sizes are low (<50).

Environmental Biology of Fishes↗

USGS2011

Workshop discusses community models for coastal sediment transport

Numerical models of coastal sediment transport are increasingly used to address problems ranging from remediation of contaminated sediments, to siting of sewage outfalls and disposal sites, to evaluating impacts of coastal development. They are also used as a test bed for sediment-transport algorithms, to provide realistic settings for biological and geochemical models, and for a variety of other research, both fundamental and applied. However, there are few full-featured, publicly available coastal sediment-transport models, and fewer still that are well tested and have been widely applied. This was the motivation for a workshop in Woods Hole, Massachusetts, on June 22–23, 2000, that explored the establishment of community models for coastal sediment-transport processes.

Eos Science News↗

USGS1999

Dynamic 3D simulations of earthquakes on en echelon faults

One of the mysteries of earthquake mechanics is why earthquakes stop. This process determines the difference between small and devastating ruptures. One possibility is that fault geometry controls earthquake size. We test this hypothesis using a numerical algorithm that simulates spontaneous rupture propagation in a three-dimensional medium and apply our knowledge to two California fault zones. We find that the size difference between the 1934 and 1966 Parkfield, California, earthquakes may be the product of a stepover at the southern end of the 1934 earthquake and show how the 1992 Landers, California, earthquake followed physically reasonable expectations when it jumped across en echelon faults to become a large event. If there are no linking structures, such as transfer faults, then strike-slip earthquakes are unlikely to propagate through stepover s >5 km wide.

Geophysical Research Letters↗