Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth in Space”

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

North Dakota and Landsat

The State of North Dakota once did not figure prominently in the Nation’s economy. The sparsely populated State supported food production, and hunters and anglers were drawn to its lakes, rivers, and wide-open spaces, but its economy was overshadowed by that of other States. However, the State and its prairie expanses recently rocketed from an economic afterthought to a national energy leader with the soaring production of oil and natural gas in the Bakken oil patch. The Bakken development has been transformative for North Dakota’s landscapes in myriad ways. It has boosted economic output, drawn thousands of new residents to cities like Williston and Watford City, and led to a proliferation of oil and gas pads. In the past two decades, North Dakota experienced other major changes, such as the expansion of the depressional wetlands of the Prairie Pothole Region on the eastern side of the State. These critical breeding areas for waterfowl, which stretch across Minnesota, South Dakota, North Dakota, and Canada, are home to more than 50 percent of North America’s migratory birds. Changes from oil and gas production, urban development, and wetland resurgence can all be tracked over time using the unparalleled Earth observation record of the U.S. Geological Survey Landsat data archive. Its 50-year record of repeat imagery also aids in the monitoring, cataloging, and management of cropland, invasive insect species, and natural or human-made disaster recovery. Here are just a few examples of the benefits offered to North Dakota by the Landsat Program.

North Dakota↗

Tsunamis generated by subaerial mass flows

Tsunamis generated in lakes and reservoirs by subaerial mass flows pose distinctive problems for hazards assessment because the domain of interest is commonly the "near field," beyond the zone of complex splashing but close enough to the source that wave propagation effects are not predominant. Scaling analysis of the equations governing water wave propagation shows that near-field wave amplitude and wavelength should depend on certain measures of mass flow dynamics and volume. The scaling analysis motivates a successful collapse (in dimensionless space) of data from two distinct sets of experiments with solid block "wave makers." To first order, wave amplitude/water depth is a simple function of the ratio of dimensionless wave maker travel time to dimensionless wave maker volume per unit width. Wave amplitude data from previous laboratory investigations with both rigid and deformable wave makers follow the same trend in dimensionless parameter space as our own data. The characteristic wavelength/water depth for all our experiments is simply proportional to dimensionless wave maker travel time, which is itself given approximately by a simple function of wave maker length/water depth. Wave maker shape and rigidity do not otherwise influence wave features. Application of the amplitude scaling relation to several historical events yields "predicted" near-field wave amplitudes in reasonable agreement with measurements and observations. Together, the scaling relations for near-field amplitude, wavelength, and submerged travel time provide key inputs necessary for computational wave propagation and hazards assessment.

Journal of Geophysical Research B: Solid Earth↗

Forecasting, detecting, and tracking volcanic eruptions from space

Satellite monitoring of volcanic activity typically includes four primary observations: (1) deformation and surface change, (2) gas emissions, (3) thermal anomalies, and (4) ash plumes. These phenomena are imaged by remote sensing data that span the electromagnetic spectrum, from microwave to ultraviolet energy and including visible and infrared wavelengths. The primary uses of satellite data in volcanology are forecasting, detecting, and tracking eruptive activity. Eruptions are often preceded by a number of indicators that are detectable from space, including surface deformation, subtle increases in surface temperature, and elevated gas emissions. The first indications of eruption, especially at remote volcanoes, are often identified in satellite data by strong thermal anomalies and/or the presence of ash and gas in the atmosphere, the recognition of which can be automated for rapid eruption detection. Once an eruption is in progress, space-based imagery of all types can track activity over time, providing information on the emplacement of volcanic deposits, the presence and character of ash plumes, and potential changes in the character of the eruption, all of which aid hazards assessment. Activity at Agung volcano, Indonesia, during 2017–2019, offers an excellent example of the importance of remote sensing datasets for forecasting, detecting, and tracking eruptions. Challenges to exploiting current and future satellite data include ensuring regular acquisitions over active volcanoes and developing tools for automated analysis of the massive volume of imagery for volcano-related signals.

Remote Sensing in Earth Systems Sciences↗

Planning for the Global Earth Observation System of Systems (GEOSS)

The Group on Earth Observations was established to promote comprehensive, coordinated, and sustained Earth observations. Its mandate is to implement the Global Earth Observation System of Systems (GEOSS) in accord with the GEOSS 10-Year Implementation Plan and Reference Document. During the months over which the GEOSS Implementation Plan was developed, many issues surfaced and were addressed. This article discusses several of the more interesting or challenging of those issues-e.g. fitting in with existing organizations and securing stable funding - some of which have yet to be resolved fully as of this writing. Despite the relatively short period over which the Implementation Plan had to be developed, there is a good chance that the work undertaken will be influential for decades to come. ?? 2005 Elsevier Ltd. All rights reserved.

Space Policy↗

Cassini/VIMS observations of the moon

In this paper, we present preliminary scientific results obtained from the analysis of VIMS (Visible and Infrared Mapping Spectrometer) lunar images and spectra. These data were obtained during the Cassini Earth flyby in August 1999. Spectral ratios have been produced in order to derive lunar mineralogical maps. Some spectra observed at the north-east lunar limb, show few unusual absorption features located at 0.357, 0.430 and 0.452 ??m, the origin of which is presently unknown. ?? 2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Advances in Space Research↗

The 1954 Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes: A triggered normal faulting sequence

In 1954, four earthquakes of M > 6.0 occurred within a 30 km radius in a period of six months. The Rainbow Mountain-Fairview Peak-Dixie Valley earthquakes are among the largest to have been recorded geodetically in the Basin and Range province. The Fairview Peak earthquake (M = 7.2, December 12, 1954) followed two events in the Rainbow Mountains ( M = 6.2, July 6, and M = 6.5, August 24, 1954) by 6 months. Four minutes later the Dixie Valley fault ruptured ( M = 6.7, December 12, 1954). The changes in static stresses caused by the events are calculated using the Coulomb-Navier failure criterion and assuming uniform slip on rectangular dislocations embedded in an elastic half-space. Coulomb stress changes are resolved on optimally oriented faults and on each of the faults that ruptured in the chain of events. These calculations show that each earthquake in the Rainbow Mountain-Fairview Peak-Dixie Valley sequence was preceded by a static stress change that encouraged failure. The magnitude of the stress increases transferred from one earthquake to another ranged from 0.01 MPa (0.1 bar) to over 0.1 MPa (1 bar). Stresses were reduced by up to 0.1 MPa over most of the Rainbow Mountain-Fairview Peak area as a result of the earthquake sequence.

Journal of Geophysical Research B: Solid Earth↗

Digital elevation models: Terminology and definitions

Digital elevation models (DEMs) provide fundamental depictions of the three-dimensional shape of the Earth’s surface and are useful to a wide range of disciplines. Ideally, DEMs record the interface between the atmosphere and the lithosphere using a discrete two-dimensional grid, with complexities introduced by the intervening hydrosphere, cryosphere, biosphere, and anthroposphere. The treatment of DEM surfaces, affected by these intervening spheres, depends on their intended use, and the characteristics of the sensors that were used to create them. DEM is a general term, and more specific terms such as digital surface model (DSM) or digital terrain model (DTM) record the treatment of the intermediate surfaces. Several global DEMs generated with optical (visible and near-infrared) sensors and synthetic aperture radar (SAR), as well as single/multi-beam sonars and products of satellite altimetry, share the common characteristic of a georectified, gridded storage structure. Nevertheless, not all DEMs share the same vertical datum, not all use the same convention for the area on the ground represented by each pixel in the DEM, and some of them have variable data spacings depending on the latitude. This paper highlights the importance of knowing, understanding and reflecting on the sensor and DEM characteristics and consolidates terminology and definitions of key concepts to facilitate a common understanding among the growing community of DEM users, who do not necessarily share the same background.

Remote Sensing↗

Elemental abundance anomalies in the late Cenomanian extinction interval: A search for the source(s)

Elemental abundances have been measured by neutron activation methods across the Cenomanian-Turonian (late Cretaceous) extinction interval in samples collected from sixteen sites in the Western Interior Basin of North America and from twelve widely separated locations around the globe, including six ODP/DSDP sites. In most Western Interior Basin sites, in Colombia, and in western Europe (weaker), two closely spaced elemental abundance peaks occur in the upper Cenomanian (∼ 92 m.y.), spanning the ammonite zones of Sciponoceras gracile through Neocardioceras juddii . Elements with anomalously high concentrations include Sc, Ti, V, Cr, Mn, Co, Ni, Ir, Pt and Au. The lower peak coincides with the disappearance (extinction) of the foraminifer Rotalipora cushmani . In North American sections R. greenhornensis also disappears at or just below this horizon, but in Europe it disappears considerably earlier than R. cushmani . A series of molluscan extinction and speciation or migration events also begins near the stratigraphic level of the lower elemental abundance peak.

Earth and Planetary Science Letters↗

Petrology and tectonics of Phanerozoic continent formation: From island arcs to accretion and continental arc magmatism

Mesozoic continental arcs in the North American Cordillera were examined here to establish a baseline model for Phanerozoic continent formation. We combine new trace-element data on lower crustal xenoliths from the Mesozoic Sierra Nevada Batholith with an extensive grid-based geochemical map of the Peninsular Ranges Batholith, the southern equivalent of the Sierras. Collectively, these observations give a three-dimensional view of the crust, which permits the petrogenesis and tectonics of Phanerozoic crust formation to be linked in space and time. Subduction of the Farallon plate beneath North America during the Triassic to early Cretaceous was characterized by trench retreat and slab rollback because old and cold oceanic lithosphere was being subducted. This generated an extensional subduction zone, which created fringing island arcs just off the Paleozoic continental margin. However, as the age of the Farallon plate at the time of subduction decreased, the extensional environment waned, allowing the fringing island arc to accrete onto the continental margin. With continued subduction, a continental arc was born and a progressively more compressional environment developed as the age of subducting slab continued to young. Refinement into a felsic crust occurred after accretion, that is, during the continental arc stage, wherein a thickened crustal and lithospheric column permitted a longer differentiation column. New basaltic arc magmas underplate and intrude the accreted terrane, suture, and former continental margin. Interaction of these basaltic magmas with pre-existing crust and lithospheric mantle created garnet pyroxenitic mafic cumulates by fractional crystallization at depth as well as gabbroic and garnet pyroxenitic restites at shallower levels by melting of pre-existing lower crust. The complementary felsic plutons formed by these deep-seated differentiation processes rose into the upper crust, stitching together the accreted terrane, suture and former continental margin. The mafic cumulates and restites, owing to their high densities, eventually foundered into the mantle, leaving behind a more felsic crust. Our grid-based sampling allows us to estimate an unbiased average upper crustal composition for the Peninsular Ranges Batholith. Major and trace-element compositions are very similar to global continental crust averaged over space and time, but in detail, the Peninsular Ranges are slightly lower in compatible to mildly incompatible elements, MgO, Mg#, V, Sc, Co, and Cr. The compositional similarities suggest a strong arc component in global continental crust, but the slight discrepancies suggest that additional crust formation processes are also important in continent formation as a whole. Finally, the delaminated Sierran garnet pyroxenites have some of the lowest U/Pb ratios ever measured for silicate rocks. Such material, if recycled and stored in the deep mantle, would generate a reservoir with very unradiogenic Pb, providing one solution to the global Pb isotope paradox.

Earth and Planetary Science Letters↗

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Estimating vegetation biomass and cover across large plots in shrub and grass dominated drylands using terrestrial lidar and machine learning

Terrestrial laser scanning (TLS) has been shown to enable an efficient, precise, and non-destructive inventory of vegetation structure at ranges up to hundreds of meters. We developed a method that leverages TLS collections with machine learning techniques to model and map canopy cover and biomass of several classes of short-stature vegetation across large plots. We collected high-definition TLS scans of 26 1-ha plots in desert grasslands and big sagebrush shrublands in southwest Idaho, USA. We used the Random Forests machine learning algorithm to develop decision tree models predicting the biomass and canopy cover of several vegetation classes from statistical descriptors of the aboveground heights of TLS points. Manual measurements of vegetation characteristics collected within each plot served as training and validation data. Models based on five or fewer TLS descriptors of vegetation heights were developed to predict the canopy cover fraction of shrubs (R 2 = 0.77, RMSE = 7%), annual grasses (R 2 = 0.70, RMSE = 21%), perennial grasses (R 2 = 0.36, RMSE = 12%), forbs (R 2 = 0.52, RMSE = 6%), bare earth or litter (R 2 = 0.49, RMSE = 19%), and the biomass of shrubs (R 2 = 0.71, RMSE = 175 g) and herbaceous vegetation (R 2 = 0.61, RMSE = 99 g) (all values reported are out-of-bag). Our models explained much of the variability between predictions and manual measurements, and yet we expect that future applications could produce even better results by reducing some of the methodological sources of error that we encountered. Our work demonstrates how TLS can be used efficiently to extend manual measurement of vegetation characteristics from small to large plots in grasslands and shrublands, with potential application to other similarly structured ecosystems. Our method shows that vegetation structural characteristics can be modeled without classifying and delineating individual plants, a challenging and time-consuming step common in previous methods applying TLS to vegetation inventory. Improving application of TLS to studies of shrub-steppe ecosystems will serve immediate management needs by enhancing vegetation inventories, environmental modeling studies, and the ability to train broader datasets collected from air and space.

Idaho↗

Spectroscopic remote sensing for material identification, vegetation characterization, and mapping

Identifying materials by measuring and analyzing their reflectance spectra has been an important procedure in analytical chemistry for decades. Airborne and space-based imaging spectrometers allow materials to be mapped across the landscape. With many existing airborne sensors and new satellite-borne sensors planned for the future, robust methods are needed to fully exploit the information content of hyperspectral remote sensing data. A method of identifying and mapping materials using spectral feature analyses of reflectance data in an expert-system framework called MICA (Material Identification and Characterization Algorithm) is described. MICA is a module of the PRISM (Processing Routines in IDL for Spectroscopic Measurements) software, available to the public from the U.S. Geological Survey (USGS) at http://pubs.usgs.gov/of/2011/1155/. The core concepts of MICA include continuum removal and linear regression to compare key diagnostic absorption features in reference laboratory/field spectra and the spectra being analyzed. The reference spectra, diagnostic features, and threshold constraints are defined within a user-developed MICA command file (MCF). Building on several decades of experience in mineral mapping, a broadly-applicable MCF was developed to detect a set of minerals frequently occurring on the Earth's surface and applied to map minerals in the country-wide coverage of the 2007 Afghanistan HyMap data set. MICA has also been applied to detect sub-pixel oil contamination in marshes impacted by the Deepwater Horizon incident by discriminating the C-H absorption features in oil residues from background vegetation. These two recent examples demonstrate the utility of a spectroscopic approach to remote sensing for identifying and mapping the distributions of materials in imaging spectrometer data.

Conference Paper↗

Spatial and temporal variations in landscape evolution: historic and longer-term sediment flux through global catchments

Sediment generation and transport through terrestrial catchments influence soil distribution, geochemical cycling of particulate and dissolved loads, and the character of the stratigraphic record of Earth history. To assess the spatiotemporal variation in landscape evolution, we compare global compilations of stream gauge–derived () and cosmogenic radionuclide (CRN)–derived (predominantly 10 Be; ) denudation of catchments (mm/yr) and sediment load of rivers (Mt/yr). Stream gauges measure suspended sediment loads of rivers during several to tens of years, whereas CRNs provide catchment-integrated denudation rates at 10 2 –10 5 -yr time scales. Stream gauge–derived and CRN-derived sediment loads in close proximity to one another (<500 km) exhibit broad similarity ( stream gauge samples; CRN samples). Nearly two-thirds of CRN-derived sediment loads exceed historic loads measured at the same locations (). Excessive longer-term sediment loads likely are a result of longer-term recurrence of large-magnitude sediment-transport events. Nearly 80% of sediment loads measured at approximately the same locations exhibit stream gauge loads that are within an order of magnitude of CRN loads, likely as a result of the buffering capacity of large flood plains. Catchments in which space for deposition exceeds sediment supply have greater buffering capacity. Superior locations in which to evaluate anthropogenic influences on landscape evolution might be buffered catchments, in which temporary storage of sediment in flood plains can provide stream gauge–based sediment loads and denudation rates that are applicable over longer periods than the durations of gauge measurements. The buffering capacity of catchments also has implications for interpreting the stratigraphic record; delayed sediment transfer might complicate the stratigraphic record of external forcings and catchment modification.

The Journal of Geology↗

Mantle to surface degassing of alkalic magmas at Erebus volcano, Antarctica

Continental intraplate volcanoes, such as Erebus volcano, Antarctica, are associated with extensional tectonics, mantle upwelling and high heat flow. Typically, erupted magmas are alkaline and rich in volatiles (especially CO 2 ), inherited from low degrees of partial melting of mantle sources. We examine the degassing of the magmatic system at Erebus volcano using melt inclusion data and high temporal resolution open-path Fourier transform infrared (FTIR) spectroscopic measurements of gas emissions from the active lava lake. Remarkably different gas signatures are associated with passive and explosive gas emissions, representative of volatile contents and redox conditions that reveal contrasting shallow and deep degassing sources. We show that this unexpected degassing signature provides a unique probe for magma differentiation and transfer of CO 2 -rich oxidised fluids from the mantle to the surface, and evaluate how these processes operate in time and space. Extensive crystallisation driven by CO 2 fluxing is responsible for isobaric fractionation of parental basanite magmas close to their source depth. Magma deeper than 4 kbar equilibrates under vapour-buffered conditions. At shallower depths, CO 2 -rich fluids accumulate and are then released either via convection-driven, open-system gas loss or as closed-system slugs that ascend and result in Strombolian eruptions in the lava lake. The open-system gases have a reduced state (below the QFM buffer) whereas the closed-system gases preserve their deep oxidised signatures (close to the NNO buffer).

Earth and Planetary Science Letters↗

Magnetic map of the Irish Hills and surrounding areas, San Luis Obispo County, central California

A magnetic map of the Irish Hills and surrounding areas was created as part of a cooperative research and development agreement with the Pacific Gas and Electric Company and is intended to promote further understanding of the areal geology and structure by serving as a basis for geophysical interpretations and by supporting geological mapping, mineral and water resource investigations, and other topical studies. Local spatial variations in the Earth's magnetic field (evident as anomalies on magnetic maps) reflect the distribution of magnetic minerals, primarily magnetite, in the underlying rocks. In many cases the volume content of magnetic minerals can be related to rock type, and abrupt spatial changes in the amount of magnetic minerals can be related to either lithologic or structural boundaries. Magnetic susceptibility measurements from the area indicate that bodies of serpentinite and other mafic and ultramafic rocks tend to produce the most intense magnetic anomalies, but such generalizations must be applied with caution because some sedimentary units also can produce measurable magnetic anomalies. Remanent magnetization does not appear to be a significant source for magnetic anomalies because it is an order of magnitude less than the induced magnetization. The map is a mosaic of three separate surveys collected by (1) fixed-wing aircraft at a nominal height of 305 m, (2) by boat with the sensor at sea level, and (3) by helicopter. The helicopter survey was flown by New-Sense Geophysics in October 2009 along flight lines spaced 150-m apart and at a nominal terrain clearance of 50 to 100 m. Tie lines were flown 1,500-m apart. Data were adjusted for lag error and diurnal field variations. Further processing included microleveling using the tie lines and subtraction of the reference field defined by International Geomagnetic Reference Field (IGRF) 2005 extrapolated to August 1, 2008.

California↗

ECCOE Landsat quarterly Calibration and Validation report—Quarter 2, 2023

Executive Summary The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote-sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level. This report provides observed geometric and radiometric analysis results for Landsats 7–8 for quarter 2 (April–June) of 2023. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov . One specific activity that the ECCOE Landsat Cal/Val Team closely monitored was a Landsat 8 Thermal Infrared Sensor (TIRS) Scene Select Mechanism (SSM) excursion anomaly. On April 21, 2023, a TIRS SSM excursion error flag was indicated in telemetry during a calibration activity when the SSM encoder was powered on and the mirror was between the nadir position and the deep space position. An initial recovery plan indicated the SSM was moving erratically, so the instrument was put into a safe state for additional troubleshooting. A second recovery plan was developed and successfully executed on April 23, 2023. Additional information about the Landsat 8 TIRS SSM excursion anomaly is available at https://www.usgs.gov/landsat-missions/news/landsat-8-level-1-product-processing-resumes .

Open-File Report↗

Landsat data availability from the EROS Data Center and status of future plans

The Department of Interior's EROS Data Center, managed by the U.S. Geological Survey, was established in 1972, in Sioux Falls, South Dakota, to serve as a principal dissemination facility for Landsat and other remotely Sensed data. Through the middle of 1977, the Center has supplied approximately 1.7 million copies of images from the more than 5 million images of the Earth's surface archived at the Center. Landsat accounted for half of these images plus approximately 5,800 computer-compatible tapes of Landsat data were also supplied to users. New methods for processing data products to make them more useful are being developed, and new accession aids for determining data availability are being placed in operation. The Center also provides assistance and training to resource specialists and land managers in the use of Landsat and other remotely sensed data. A Data Analysis Laboratory is operated at the Center to provide both digital and analog multispectral/multitemporal image analysis capabilities in support of the training and assistance programs. In addition to conventionally processed data products, radiometrically enhanced Landsat imagery are now available from the Center in limited quantities. In mid-1978, the Center will convert to an all-digital processing system for Landsat data that will provide improved products for user analysis in production quantities. The Department of Interior and NASA are currently studying concepts that use communication satellites to relay Landsat data between U.S. ground stations, Goddard Space Flight Center and the EROS Data Center which would improve the timeliness of data availability. The Data Center also works closely with the remote sensing programs and Landsat data receiving and processing facilities being developed in foreign countries.

Conference Paper↗

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth↗