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.

SourceResearch connectionFrom yearThrough year

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

USGS1967

Volcanic hazards at Mount Rainier, Washington

Mount Rainier is a large stratovolcano of andesitic rock in the Cascade Range of western Washington. Although the volcano as it now stands was almost completely formed before the last major glaciation, geologic formations record a variety of events that have occurred at the volcano in postglacial time. Repetition of some of these events today without warning would result in property damage and loss of life on a catastrophic scale. It is appropriate, therefore, to examine the extent, frequency, and apparent origin of these phenomena and to attempt to predict the effects on man of similar events in the future. The present report was prompted by a contrast that we noted during a study of surficial geologic deposits in Mount Rainier National Park, between the present tranquil landscape adjacent to the volcano and the violent events that shaped parts of that same landscape in the recent past. Natural catastrophes that have geologic causes - such as eruptions, landslides, earthquakes, and floods - all too often are disastrous primarily because man has not understood and made allowance for the geologic environment he occupies. Assessment of the potential hazards of a volcanic environment is especially difficult, for prediction of the time and kind of volcanic activity is still an imperfect art, even at active volcanoes whose behavior has been closely observed for many years. Qualified predictions, however, can be used to plan ways in which hazards to life and property can be minimized. The prediction of eruptions is handicapped because volcanism results from conditions far beneath the surface of the earth, where the causative factors cannot be seen and, for the most part, cannot be measured. Consequently, long-range predictions at Mount Rainier can be based only on the past behavior of the volcano, as revealed by study of the deposits that resulted from previous eruptions. Predictions of this sort, of course, cannot be specific as to time and locale of future events, and clearly are valid only if the past behavior is, as we believe, a reliable guide. The purpose of this report is to infer the events recorded by certain postglacial deposits at Mount Rainier and to suggest what bearing similar events in the future might have on land use within and near the park. In addition, table 2 (page 22) gives possible warning signs of an impending eruption. We want to increase man's understanding of a possibly hazardous geologic environment around Mount Rainier volcano, yet we do not wish to imply for certain that the hazards described are either immediate or inevitable. However, we do believe that hazards exist, that some caution is warranted, and that some major hazards can be avoided by judicious planning. Most of the events with which we are concerned are sporadic phenomena that have resulted directly or indirectly from volcanic eruptions. Although no eruptions (other than steam emission) of the volcano in historic time are unequivocally known (Hopson and others, 1962), pyroclastic (air-laid) deposits of pumice and rock debris attest to repeated, widely spaced eruptions during the 10,000 years or so of postglacial time. In addition, the constituents of some debris flows indicate an origin during eruptions of molten rock; other debris flows, because of their large size and constituents, are believed to have been caused by steam explosions. Some debris flows, however, are not related to volcanism at all.

Washington↗

USGS1987

Basin-ring spacing on the Moon, Mercury, and Mars

Radial spacing between concentric rings of impact basins that lack central peaks is statistically similar and nonrandom on the Moon, Mercury, and Mars, both inside and outside the main ring. One spacing interval, (2.0 ?? 0.3)0.5D, or an integer multiple of it, dominates most basin rings. Three analytical approaches yield similar results from 296 remapped or newly mapped rings of 67 multi-ringed basins: least-squares of rank-grouped rings, least-squares of rank and ring diameter for each basin, and averaged ratios of adjacent rings. Analysis of 106 rings of 53 two-ring basins by the first and third methods yields an integer multiple (2 ??) of 2.00.5D. There are two exceptions: (1) Rings adjacent to the main ring of multi-ring basins are consistently spaced at a slightly, but significantly, larger interval, (2.1 ?? 0.3)0.5D; (2) The 88 rings of 44 protobasins (large peak-plus-inner-ring craters) are spaced at an entirely different interval (3.3 ?? 0.6)0.5D. The statistically constant and target-invariant spacing of so many rings suggests that this characteristic may constrain formational models of impact basins on the terrestrial planets. The key elements of such a constraint include: (1) ring positions may not have been located by the same process(es) that formed ring topography; (2) ring location and emplacement of ring topography need not be coeval; (3) ring location, but not necessarily the mode of ring emplacement, reflects one process that operated at the time of impact; and (4) the process yields similarly-disposed topographic features that are spatially discrete at 20.5D intervals, or some multiple, rather than continuous. These four elements suggest that some type of wave mechanism dominates the location, but not necessarily the formation, of basin rings. The waves may be standing, rather than travelling. The ring topography itself may be emplaced at impact by this and/or other mechanisms and may reflect additional, including post-impact, influences. ?? 1987 D. Reidel Publishing Company.

Earth, Moon and Planets↗

USGS2020

RestoreNet: An emerging restoration network reveals controls on seeding success across dryland ecosystems

Drylands are Earth's largest terrestrial biome and support one‐third of the global population. However, they are also highly vulnerable to land degradation. Despite widespread demand for dryland restoration and rehabilitation, little information is available to help land managers effectively re‐establish native perennial vegetation across drylands. RestoreNet is an emerging dryland restoration network that systematically tests revegetation techniques across environmental gradients. Using the RestoreNet experimental framework, we tested the effectiveness of restoration treatments (i.e. ConMod nurse plant structures, mulch, pits) that increase soil moisture and seed mixes with different climatic niches to achieve revegetation goals. Across sites, seedling recruitment was consistently influenced by treatment and seed mix type. Pit and mulch treatments increased total seedling density, with pits promoting the highest seeded species recruitment while limiting non‐native species establishment. Seeding increased total seedling density regardless of seed mix type, but cooler‐adapted seed mixes promoted greater seeded species density and resulted in lower density of unseeded (non‐native) species relative to warmer‐adapted mixes. Seedling recruitment was also controlled by the temporal and environmental context of restoration with the positive effect of high precipitation greatest in the weeks immediately following seeding. Above‐average precipitation during the study period across most of the sites may partially explain why the highest seeded species recruitment occurred in pit treatments and seed mixes with cooler, wetter niche requirements. Synthesis and applications . Results from the dryland restoration network, RestoreNet help to better understand variation in seeding and restoration treatment success across space and time in drylands. Relationships between restoration practices and environmental conditions in our study suggest the importance of anticipatory restoration strategies that forecast seasonal and sub‐seasonal weather conditions and select plant species with climate niche requirements appropriate for current and future climate conditions. This information is critical to land managers tasked with improving ecosystem conditions across degraded dryland regions.

Arizona↗

USGS2012

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California↗

USGS2022

Seismometer records of ground tilt induced by debris flows

A change in surface loading causes the Earth’s surface to deform. Mass movements, such as debris flows, can cause a tilt large enough to be recorded by nearby instruments, but the signal is strongly dependent on the mass loading and subsurface parameters. Specifically designed sensors for such measurements (tiltmeters) are cumbersome to install. Alternatively, broadband seismometers record translational motion and also tilt signals, often at periods of tens to hundreds of seconds. Their horizontal components are thereby the most sensitive to tilt. In this study, we show how to obtain tilt caused by the passing by of debris flows from seismic measurements recorded within tens of meters of the flow and investigate the usefulness of this signal for flow characterization. We investigate the problem on three scales (1) large‐scale laboratory experiments at the U.S. Geological Survey debris‐flow flume, where broadband seismometers and tiltmeters were installed for six 8 – 10    m 3 "> 8 – 10 m 3 experiments, (2) the Illgraben torrent in Switzerland, one of the most active mass wasting sites in the European Alps, where a broadband seismometer placed within a few meters of the channel recorded 15 debris‐flow events with volumes up to 10 5    m 3 "> 10 5 m 3 ⁠ , and (3) Volcán de Fuego, Guatemala, where a broadband seismometer recorded two lahars. We investigate how the tilt signals compare to debris‐flow parameters such as mean normal stresses, usually measured by expensive force plates, and debris‐flow height. We model the elastic ground deformation as the response of an elastic half‐space to a moving surface load. In addition, we use the model with some simplifications to determine the maximum debris‐flow heights of Volcán de Fuego events, where no force plate measurements are available. Finally, we address how and under what assumptions the relatively affordable and straightforward tilt measurements may be utilized to infer debris‐flow parameters, as opposed to force plates and other complicated instrument setups.

Oregon↗

USGS2007

Space geodetic observation of expansion of the San Gabriel Valley, California, aquifer system, during heavy rainfall in winter 2004-2005

[1] Starting early in 2005, the positions of GPS stations in the San Gabriel valley region of southern California showed statistically significant departures from their previous behavior. Station LONG moved up by about 47 mm, and nearby stations moved away from LONG by about 10 mm. These changes began during an extremely rainy season in southern California and coincided with a 16-m increase in water level at a nearby well in Baldwin Park and a regional uplift detected by interferometric synthetic aperture radar. No equivalent signals were seen in GPS station position time series elsewhere in southern California. Our preferred explanation, supported by the timing and by a hydrologic simulation, is deformation due to recharging of aquifers after near-record rainfall in 2004–2005. We cannot rule out an aseismic slip event, but we consider such an event unlikely because it requires slip on multiple faults and predicts other signals that are not observed.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Origin of the pulse-like signature of shallow long-period volcano seismicity

Short-duration, pulse-like long-period (LP) events are a characteristic type of seismicity accompanying eruptive activity at Mount Etna in Italy in 2004 and 2008 and at Turrialba Volcano in Costa Rica and Ubinas Volcano in Peru in 2009. We use the discrete wave number method to compute the free surface response in the near field of a rectangular tensile crack embedded in a homogeneous elastic half space and to gain insights into the origin of the LP pulses. Two source models are considered, including (1) a vertical fluid-driven crack and (2) a unilateral tensile rupture growing at a fixed sub-Rayleigh velocity with constant opening on a vertical crack. We apply cross correlation to the synthetics and data to demonstrate that a fluid-driven crack provides a natural explanation for these data with realistic source sizes and fluid properties. Our modeling points to shallow sources (<1 km depth), whose signatures are representative of the Rayleigh pulse sampled at epicentral distances >∼1 km. While a slow-rupture failure provides another potential model for these events, the synthetics and resulting fits to the data are not optimal in this model compared to a fluid-driven source. We infer that pulse-like LP signatures are parts of the continuum of responses produced by shallow fluid-driven sources in volcanoes.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Medium spatial resolution mapping of global land cover and land cover change across multiple decades from Landsat

Land cover maps are essential for characterizing the biophysical properties of the Earth’s land areas. Because land cover information synthesizes a rich array of information related to both the ecological condition of land areas and their exploitation by humans, they are widely used for basic and applied research that requires information related to land surface properties (e.g., terrestrial carbon models, water balance models, weather, and climate models) and are core inputs to models and analyses used by natural resource scientists and land managers. As the Earth’s global population has grown over the last several decades rates of land cover change have increased dramatically, with enormous impacts on ecosystem services (e.g., biodiversity, water supply, carbon sequestration, etc.). Hence, accurate information related to land cover is essential for both managing natural resources and for understanding society’s ecological, biophysical, and resource management footprint. To address the need for high-quality land cover information we are using the global record of Landsat observations to compile annual maps of global land cover from 2001 to 2020 at 30 m spatial resolution. To create these maps we use features derived from time series of Landsat imagery in combination with ancillary geospatial data and a large database of training sites to classify land cover at annual time step. The algorithm that we apply uses temporal segmentation to identify periods with stable land cover that are separated by breakpoints in the time series. Here we provide an overview of the methods and data sets we are using to create global maps of land cover. We describe the algorithms used to create these maps and the core land cover data sets that we are creating through this effort, and we summarize our approach to accuracy assessment. We also present a synthesis of early results and discuss the strengths and weaknesses of our early map products and the challenges that we have encountered in creating global land cover data sets from Landsat. Initial accuracy assessment for North America shows good overall accuracy (77.0 ± 2.0% correctly classified) and 79.8% agreement with the European Space Agency (ESA) WorldCover product. The land cover mapping results we report provide the foundation for robust, repeatable, and accurate mapping of global land cover and land cover change across multiple decades at 30 m spatial resolution from Landsat.

Frontiers in Remote Sensing↗

USGS2013

Archive of post-Hurricane Isabel coastal oblique aerial photographs collected during U.S. Geological Survey Field Activity 03CCH01 from Ocean City, Maryland, to Fort Caswell, North Carolina and Inland from Waynesboro to Redwood, Virginia, September 21 - 23, 2003

On September 21 - 23, 2003, the United States Geological Survey (USGS) conducted an oblique aerial photographic survey along the Atlantic coast from Ocean City, Md., to Fort Caswell, N.C., and inland oblique aerial photographic survey from Waynesboro to Redwood, Va., aboard a Navajo Piper twin-engine airplane. The coastal survey was conducted at an altitude of 500 feet (ft) and approximately 1,000 ft offshore. For the inland photos, the aircraft tried to stay approximately 500 ft above the terrain. These coastal photos were used to document coastal changes like beach erosion and overwash caused by Hurricane Isabel, while the inland photos looked for potential landslides caused by heavy rains. The photos may also be used as baseline data for future coastal change analysis. The USGS and the National Aeronautics and Space Administration (NASA) surveyed the impact zone of Hurricane Isabel to better understand the changes in vulnerability of the Nation’s coasts to extreme storms (Morgan, 2009). This report serves as an archive of photographs collected during the September 21 - 23, 2003, post-Hurricane Isabel coastal and inland oblique aerial survey along with associated survey maps, KML files, navigation files, digital Field Activity Collection System (FACS) logs, and Federal Geographic Data Committee (FGDC) metadata. Refer to the Acronyms page for expansions of all acronyms and abbreviations used in this report. The USGS St. Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 03CCH01 tells us the data were collected in 2003 for the Coastal Change Hazards (CCH) study and the data were collected during the first field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the ID number. The photographs provided here are Joint Photographic Experts Group (JPEG) scanned images of the analog 35 millimeter (mm) color positive slides. The photograph locations are estimates of the location of the plane (see the Navigation page). The metadata values for photo creation time, GPS latitude, GPS longitude, GPS position (latitude and longitude), keywords, credit, artist, caption, copyright, and contact were added to each photograph's EXIF header using EXIFtool (Subino and others, 2012). Photographs can be opened directly with any JPEG-compatible image viewer by clicking on a thumbnail on the contact sheet, or, when viewing the Google Earth KML file, by clicking on the marker and then clicking on either the thumbnail or the link below the thumbnail. Nathaniel Plant (USGS - St. Petersburg, Fla.), and Ann Marie Ascough (formerly contracted at the USGS - St. Petersburg, Fla.) helped with the creation of KML files. To view the photos and survey maps, proceed to the Photos and Maps page.

Maryland, North Carolina, Virginia↗

USGS2012

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↗

USGS2022

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

USGS2023

Debris-flow process controls on steepland morphology in the San Gabriel Mountains, California

Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris-flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady-state fluvial channel longitudinal profiles are concave-up and exhibit a power-law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope-area signature in log-log space. Here, we use a one-dimensional landform evolution model that incorporates debris-flow erosion to reproduce the relationship between this curved slope-area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris-flow erosion, valley-bottom slope, and debris-flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris-flow incision law. Finally, modeled fluvial incision outpaces debris-flow erosion at drainage areas less than those at which valleys morphologically transition from near-invariant slopes to concave profiles. This result emphasizes the critical role of debris-flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

California↗

USGS2007

Postseismic relaxation associated with transient creep rheology

[1] Perfettini and Avouac (2004) postulated that both the aftershock rate (assumed proportional to the local stressing rate) and the postseismic relaxation are driven by the loading imposed by postseismic slip on the brittle creep fault zone (BCFZ), the downdip extension of the fault zone below the coseismic rupture. I explore the consequences of that hypothesis for a long, strike-slip fault in the case where the BCFZ rheology is compatible with ordinary transient creep (creep strain proportional to log e (1 + t / τ 2 )). Because the important relaxation occurs near the bottom of the coseismic rupture, I calculate the postearthquake response with a model in which the BCFZ is represented by a viscoelastic half-space below the coseismic rupture. I find that both the predicted postseismic relaxation and the cumulative number of aftershocks can be approximated by the same temporal dependence N MO ( t ) = a MO (1−(1 + t / τ ) 1− p )/( p − 1), where t is the time after the earthquake and a MO , τ , and p are the constants chosen to fit either data set. Notice that d N MO ( t )/d t = ( a MO / τ )/(1 + t / τ ) p is the modified Omori law used to describe the rate of aftershock occurrence. Thus, the modified Omori law can be understood as a consequence of the Perfettini–Avouac hypothesis (aftershocks driven by slip on the BCFZ) and a BCFZ rheology compatible with ordinary transient creep. Moreover, the temporal dependence N MO ( t ) has been shown to fit postseismic surface deformation following at least 9 earthquakes. I also show that the conventional, one-dimensional, spring-block model of a BFCZ with a rheology compatible with ordinary transient creep leads to the same temporal dependence ( N MO ( t )).

Journal of Geophysical Research B: Solid Earth↗

USGS2007

Hydrothermal fluid flow and deformation in large calderas: Inferences from numerical simulations

[1] Inflation and deflation of large calderas is traditionally interpreted as being induced by volume change of a discrete source embedded in an elastic or viscoelastic half-space, though it has also been suggested that hydrothermal fluids may play a role. To test the latter hypothesis, we carry out numerical simulations of hydrothermal fluid flow and poroelastic deformation in calderas by coupling two numerical codes: (1) TOUGH2 [Pruess et al., 1999], which simulates flow in porous or fractured media, and (2) BIOT2 [Hsieh, 1996], which simulates fluid flow and deformation in a linearly elastic porous medium. In the simulations, high-temperature water (350°C) is injected at variable rates into a cylinder (radius 50 km, height 3–5 km). A sensitivity analysis indicates that small differences in the values of permeability and its anisotropy, the depth and rate of hydrothermal injection, and the values of the shear modulus may lead to significant variations in the magnitude, rate, and geometry of ground surface displacement, or uplift. Some of the simulated uplift rates are similar to observed uplift rates in large calderas, suggesting that the injection of aqueous fluids into the shallow crust may explain some of the deformation observed in calderas.

Journal of Geophysical Research B: Solid Earth↗

USGS2023

Estimating geomagnetically induced currents in southern Brazil using 3-D Earth resistivity model

Geomagnetically induced currents (GICs) result from the interaction of the time variation of ground magnetic field during a geomagnetic disturbance with the Earth's deep electrical resistivity structure. In this study, we simulate induced GICs in a hypothetical representation of a low-latitude power transmission network located mainly over the large Paleozoic Paraná basin (PB) in southern Brazil. Two intense geomagnetic storms in June and December 2015 are chosen and geoelectric fields are calculated by convolving a three-dimensional (3-D) Earth resistivity model with recorded geomagnetic variations. The dB / dt proxy often used to characterize GIC activity fails during the June storm mainly due to the relationship of the instantaneous geoelectric field to previous magnetic field values. Precise resistances of network components are unknown, so assumptions are made for calculating GIC flows from the derived geoelectric field. The largest GICs are modeled in regions of low conductance in the 3-D resistivity model, concentrated in an isolated substation at the northern edge of the network and in a cluster of substations in its central part where the east-west (E-W) oriented transmission lines coincide with the orientation of the instantaneous geoelectric field. The maximum magnitude of the modeled GIC was obtained during the main phase of the June storm, modeled at a northern substation, while the lowest magnitudes were found over prominent crustal anomalies along the PB axis and bordering the continental margin. The simulation results will be used to prospect the optimal substations for installation of GIC monitoring equipment.

Space Weather↗

USGS2012

Monitoring glacier surface seismicity in time and space using Rayleigh waves

Sliding glaciers and brittle ice failure generate seismic body and surface wave energy characteristic to the source mechanism. Here we analyze continuous seismic recordings from an array of nine short-period passive seismometers located on Bench Glacier, Alaska (USA) (61.033°N, 145.687°W). We focus on the arrival-time and amplitude information of the dominant Rayleigh wave phase. Over a 46-hour period we detect thousands of events using a cross-correlation based event identification method. Travel-time inversion of a subset of events (7% of the total) defines an active crevasse, propagating more than 200 meters in three hours. From the Rayleigh wave amplitudes, we estimate the amount of volumetric opening along the crevasse as well as an average bulk attenuation ( = 42) for the ice in this part of the glacier. With the remaining icequake signals we establish a diurnal periodicity in seismicity, indicating that surface run-off and subglacial water pressure changes likely control the triggering of these surface events. Furthermore, we find that these events are too weak (i.e., too noisy) to locate individually. However, stacking individual events increases the signal-to-noise ratio of the waveforms, implying that these periodic sources are effectively stationary during the recording period.

Journal of Geophysical Research F: Earth Surface↗

USGS1998

On the influence of biomass burning on the seasonal CO2 signal as observed at monitoring stations

We investigated the role of biomass burning in simulating the seasonal signal in both prognostic and diagnostic analyses. The prognostic analysis involved the High-Resolution Biosphere Model, a prognostic terrestrial biosphere model, and the coupled vegetation fire module, which together produce a prognostic data set of biomass burning. The diagnostic analysis involved the Simple Diagnostic Biosphere Model (SDBM) and the Hao and Liu [1994] diagnostic data set of biomass burning, which have been scaled to global 2 and 4 Pg C yr −1 , respectively. The monthly carbon exchange fields between the atmosphere and the biosphere with a spatial resolution of 0.5° × 0.5°, the seasonal atmosphere-ocean exchange fields, and the emissions from fossil fuels have been coupled to the three-dimensional atmospheric transport model TM2. We have chosen eight monitoring stations of the National Oceanic and Atmospheric Administration network to compare the predicted seasonal atmospheric CO 2 signals with those deduced from atmosphere-biosphere carbon exchange fluxes without any contribution from biomass burning. The prognostic analysis and the diagnostic analysis with global burning emissions of 4 Pg C yr −1 agree with respect to the change in the amplitude of the seasonal CO 2 concentration introduced through biomass burning. We find that the seasonal CO 2 signal at stations in higher northern latitudes (north of 30°N) is marginally influenced by biomass burning. For stations in tropical regions an increase in the CO 2 amplitude of more than 1 ppmv (up to 50% with respect to the observed trough to peak amplitude) has been calculated. Biomass burning at stations farther south accounts for an increase in the CO 2 amplitude of up to 59% (0.6 ppmv). A change in the phase of the seasonal CO 2 signal at tropical and southern stations has been shown to be strongly influenced by the onset of biomass burning in southern tropical Africa and America. Comparing simulated and observed seasonal CO 2 signals, we find higher discrepancies at southern tropical stations if biomass burning emissions are included. This is caused by the additional increase in the amplitude in the prognostic analysis and a phase shift in a diagnostic analysis. In contrast, at the northern tropical stations biomass burning tends to improve the estimates of the seasonal CO 2 signal in the prognostic analysis because of strengthening of the amplitude. Since the SDBM predicts the seasonal CO 2 signal reasonably well for the northern hemisphere tropical stations, no general improvement of the fit occurs if biomass burning emissions are considered.

Global Biogeochemical Cycles↗

USGS2002

Fault structure and mechanics of the Hayward Fault, California from double-difference earthquake locations

The relationship between small-magnitude seismicity and large-scale crustal faulting along the Hayward Fault, California, is investigated using a double-difference (DD) earthquake location algorithm. We used the DD method to determine high-resolution hypocenter locations of the seismicity that occurred between 1967 and 1998. The DD technique incorporates catalog travel time data and relative P and S wave arrival time measurements from waveform cross correlation to solve for the hypocentral separation between events. The relocated seismicity reveals a narrow, near-vertical fault zone at most locations. This zone follows the Hayward Fault along its northern half and then diverges from it to the east near San Leandro, forming the Mission trend. The relocated seismicity is consistent with the idea that slip from the Calaveras Fault is transferred over the Mission trend onto the northern Hayward Fault. The Mission trend is not clearly associated with any mapped active fault as it continues to the south and joins the Calaveras Fault at Calaveras Reservoir. In some locations, discrete structures adjacent to the main trace are seen, features that were previously hidden in the uncertainty of the network locations. The fine structure of the seismicity suggests that the fault surface on the northern Hayward Fault is curved or that the events occur on several substructures. Near San Leandro, where the more westerly striking trend of the Mission seismicity intersects with the surface trace of the (aseismic) southern Hayward Fault, the seismicity remains diffuse after relocation, with strong variation in focal mechanisms between adjacent events indicating a highly fractured zone of deformation. The seismicity is highly organized in space, especially on the northern Hayward Fault, where it forms horizontal, slip-parallel streaks of hypocenters of only a few tens of meters width, bounded by areas almost absent of seismic activity. During the interval from 1984 to 1998, when digital waveforms are available, we find that fewer than 6.5% of the earthquakes can be classified as repeating earthquakes, events that rupture the same fault patch more than one time. These most commonly are located in the shallow creeping part of the fault, or within the streaks at greater depth. The slow repeat rate of 2–3 times within the 15-year observation period for events with magnitudes around M = 1.5 is indicative of a low slip rate or a high stress drop. The absence of microearthquakes over large, contiguous areas of the northern Hayward Fault plane in the depth interval from ∼5 to 10 km and the concentrations of seismicity at these depths suggest that the aseismic regions are either locked or retarded and are storing strain energy for release in future large-magnitude earthquakes.

California↗