Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Oceans”

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 541 records · Page 30Linked to original sources
USGS2005

Modeling effects of bank friction and woody bank vegetation on channel flow and boundary shear stress in the Rio Puerco, New Mexico

[1] We have applied a physically based model for steady, horizontally uniform flow to calculate reach-averaged velocity and boundary shear-stress distributions in a natural stream with woody vegetation on the channel banks. The model calculates explicitly the form drag on woody plant stems and includes the effects of vegetation on the boundary shear stress, velocity, and turbulence fields. Average channel shapes, bed gradients, and shrub characteristics were determined for four long, internally similar reaches covering an 81 -km segment of the lower Rio Puerco, in north-central New Mexico. Hydraulic geometries for each of three flow events were determined from reach-average channel shapes, slopes, and measured silt-line heights above the bed. Model results show that friction on the lateral boundaries reduced the boundary shear stress in the center of the Rio Puerco channel by as much as 20%. In reaches with moderate to dense bank shrubs, perimeter-averaged boundary shear stresses were reduced by almost 40% during near bankfull flows. Model-calculated discharges for all three flow events indicate there was a loss of about 40% of the upstream discharge through the 81-km river segment. Infiltration rates estimated from discharge losses during the quasi-steady peak flows and areas of the bed and banks suggest rates of flow loss into the silty sand of the upper banks are about 8 times greater than rates of loss into the clayey silt-covered bed and lower banks.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Acoustic bed velocity and bed load dynamics in a large sand bed river

Development of a practical technology for rapid quantification of bed load transport in large rivers would represent a revolutionary advance for sediment monitoring and the investigation of fluvial dynamics. Measurement of bed load motion with acoustic Doppler current profiles (ADCPs) has emerged as a promising approach for evaluating bed load transport. However, a better understanding of how ADCP data relate to conditions near the stream bed is necessary to make the method practical for quantitative applications. In this paper, we discuss the response of ADCP bed velocity measurements, defined as the near-bed sediment velocity detected by the instrument's bottom-tracking feature, to changing sediment-transporting conditions in the lower Missouri River. Bed velocity represents a weighted average of backscatter from moving bed load particles and spectral reflections from the immobile bed. The ratio of bed velocity to mean bed load particle velocity depends on the concentration of the particles moving in the bed load layer, the bed load layer thickness, and the backscatter strength from a unit area of moving particles relative to the echo strength from a unit area of unobstructed bed. A model based on existing bed load transport theory predicted measured bed velocities from hydraulic and grain size measurements with reasonable success. Bed velocities become more variable and increase more rapidly with shear stress when the transport stage, defined as the ratio of skin friction to the critical shear stress for particle entrainment, exceeds a threshold of about 17. This transition in bed velocity response appears to be associated with the appearance of longer, flatter bed forms at high transport stages.

Journal of Geophysical Research F: Earth Surface↗
USGS2007

Suspended sediment transport in an ephemeral stream following wildfire

We examine the impacts of a stand-clearing wildfire on the characteristics and magnitude of suspended sediment transport in ephemeral streams draining the burn area. We report the results of a monitoring program that includes 2 years of data prior to the Cerro Grande fire in New Mexico, and 3 years of postfire data. Suspended sediment concentration (SSC) increased by about 2 orders of magnitude following the fire, and the proportion of silt and clay increased from 50% to 80%. For a given flow event, SSC is highest at the flood bore and decreases monotonically with time, a pattern evident in every flood sampled both before and after the fire. We propose that the accumulation of flow and wash load at the flow front is an inherent characteristic of ephemeral stream flows, due to amplified momentum losses at the flood bore. We present a new model for computing suspended sediment transport in ephemeral streams (in the presence or absence of wildfire) by relating SSC to the time following the arrival of the flood bore, rather than to instantaneous discharge. Using this model and a rainfall history, we estimate that in the 3 years following the fire, floods transported in suspension a mass equivalent to about 3 mm of landscape lowering across the burn area, 20% of this following a single rainstorm. We test the model by computing fine sediment delivery to a small reservoir in an adjacent watershed, where we have a detailed record of postfire sedimentation based on repeat surveys. Systematic discrepancies between modeled and measured sedimentation rates in the reservoir suggest rapid reductions in suspended sediment delivery in the first several years after the fire.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Effect of tree line advance on carbon storage in NW Alaska

We investigated the size, distribution, and temporal dynamics of ecosystem carbon (C) pools in an area of recent tree line advance, northwest Alaska. Repeat aerial photographs show forest cover increased ???10% in our study area since 1949. We sampled C pools of four principal ecosystem types, tussock tundra, shrub tundra, woodland, and forest, all located on a 600-800 year old river terrace. Significant differences between ecosystem C pools, both above ground and below ground existed. Tundra sites store >22.2 kg C/m2, shrub tundra sites and woodland sites store 9.7 kg C/m2 and 14.3 kg C/m2, respectively, and forest sites store 14.4 kg C/m2. Landscape variation of total ecosystem C was primarily due to organic soil C and was secondarily due to C stored in trees. Soil C/N profiles of shrub tundra sites and woodland sites showed similarities with forest site soils at surface and tundra site soils at depth. We hypothesize that tundra systems transformed to forest systems in this area under a progression of permafrost degradation and enhanced drainage. On the basis of C pool estimates for the different ecosystem types, conversion of tundra sites to forest may have resulted in a net loss of > 7.8 kg C/m2, since aboveground C gains were more than offset by belowground C losses to decomposition in the tundra sites. Tree line advance therefore might not increase C storage in high-latitude ecosystems and thus might not, as previously suggested, act as a negative feedback to warming. Key to this hypothesis and to its projection to future climate response is the fate of soil carbon upon warming and permafrost drainage. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Geomorphic control of landscape carbon accumulation

We use the CREEP process-response model to simulate soil organic carbon accumulation in an undisturbed prairie site in Iowa. Our primary objectives are to identify spatial patterns of carbon accumulation, and explore the effect of erosion on basin-scale C accumulation. Our results point to two general findings. First, redistribution of soil carbon by erosion results in a net increase in basin-wide carbon storage relative to a noneroding environment. Landscape-average mean residence times are increased in an eroding landscape owing to the burial/preservation of otherwise labile C. Second, field observations taken along a slope transect may overlook significant intraslope variations in carbon accumulation. Spatial patterns of modeled deep C accumulation are complex. While surface carbon with its relatively short equilibration time is predictable from surface properties, deep carbon is strongly influenced by the landscape's geomorphic and climatic history, resulting in wide spatial variability. Convergence and divergence associated with upland swales and interfluves result in bimodal carbon distributions in upper and mid slopes; variability in carbon storage within modeled mid slopes was as high as simulated differences between erosional shoulders and depositional valley bottoms. The bimodality of mid-slope C variability in the model suggests that a three-dimensional sampling strategy is preferable over the traditional two-dimensional analog or "catena" approach. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Phosphate oxygen isotope ratios as a tracer for sources and cycling of phosphate in North San Francisco Bay, California

[1] A seasonal analysis assesing variations in the oxygen isotopic composition of dissolved inorganic phosphate (DIP) was conducted in the San Francisco Bay estuarine system, California. Isotopic fractionation of oxygen in DIP (exchange of oxygen between phosphate and environmental water) at surface water temperatures occurs only as a result of enzyme‐mediated, biological reactions. Accordingly, if phospate demand is low relative to input and phosphate is not heavily cycled in the ecosystem, the oxygen isotopic composition of DIP (δ 18 O p ) will reflect the isotopic composition of the source of phosphate to the system. Such is the case for the North San Francisco Bay, an anthropogenically impacted estuary with high surface water phosphate concentrations. Variability in the δ 18 O p in the bay is primarily controlled by mixing of water masses with different δ 18 O p signatures. The δ 18 O p values range from 11.4‰ at the Sacramento River to 20.1‰ at the Golden Gate. Deviations from the two‐component mixing model for the North Bay reflect additional, local sources of phosphate to the estuary that vary seasonally. Most notably, deviations from the mixing model occur at the confluence of a major river into the bay during periods of high river discharge and near wastewater treatment outlets. These data suggest that δ 18 O p can be an effective tool for identifying P point sources and understanding phosphate dynamics in estuarine systems.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Modeling physical and biogeochemical controls over carbon accumulation in a boreal forest soil

Boreal soils are important to the global C cycle owing to large C stocks, repeated disturbance from fire, and the potential for permafrost thaw to expose previously stable, buried C. To evaluate the primary mechanisms responsible for both short- and long-term C accumulation in boreal soils, we developed a multi-isotope (12,14C) Soil C model with dynamic soil layers that develop through time as soil organic matter burns and reaccumulates. We then evaluated the mechanisms that control organic matter turnover in boreal regions including carbon input rates, substrate recalcitrance, soil moisture and temperature, and the presence of historical permafrost to assess the importance of these factors in boreal C accumulation. Results indicate that total C accumulation is controlled by the rate of carbon input, decomposition rates, and the presence of historical permafrost. However, unlike more temperate ecosystems, one of the key mechanisms involved in C preservation in boreal soils examined here is the cooling of subsurface soil layers as soil depth increases rather than increasing recalcitrance in subsurface soils. The propagation of the 14C bomb spike into soils also illustrates the importance of historical permafrost and twentieth century warming in contemporary boreal soil respiration fluxes. Both 14C and total C simulation data also strongly suggest that boreal SOM need not be recalcitrant to accumulate; the strong role of soil temperature controls on boreal C accumulation at our modeling test site in Manitoba, Canada, indicates that carbon in the deep organic soil horizons is probably relatively labile and thus subject to perturbations that result from changing climatic conditions in the future. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Role of microbial iron reduction in the dissolution of iron hydroxysulfate minerals

Iron-hydroxysulfate minerals can be important hosts for metals such as lead, mercury, copper, zinc, silver, chromium, arsenic, and selenium and for radionuclides such as 226Ra. These mineral-bound contaminants are considered immobilized under oxic conditions. However, when anoxic conditions develop, the activities of sulfate- or iron-reducing bacteria could result in mineral dissolution, releasing these bound contaminants. Reduction of structural sulfate in the iron-hydroxysulfate mineral jarosite by sulfate-reducing bacteria has previously been demonstrated. The primary objective of this work was to evaluate the potential for anaerobic dissolution of the iron-hydroxysulfate minerals jarosite and schwertmannite at neutral PH by iron-reducing bacteria. Mineral dissolution was tested using a long-term cultivar, Geobacter metallireducens strain GS-15, and a fresh isolate Geobacter sp. strain ENN1, previously undescribed. ENN1 was isolated from the discharge site of Shadle Mine, in the southern anthracite coalfield of Pennsylvania, where schwertmannite was the predominant iron-hydroxysulfate mineral. When jarosite from Elizabeth Mine (Vermont) was provided as the sole terminal electron acceptor, resting cells of both G. metallireducens and ENN1 were able to reduce structural Fe(III), releasing Fe+2, SO4-2, and K+ ions. A lithified jarosite sample from Utah was more resistant to microbial attack, but slow release of Fe+2 was observed. Neither bacterium released Fe+2 from poorly crystalline synthetic schwertmannite. Our results indicate that exposure of jarosite to iron-reducing conditions at neutral pH is likely to promote the mobility of hazardous constituents and should therefore be considered in evaluating waste disposal and/or reclamation options involving jarosite-bearing materials.

Journal of Geophysical Research G: Biogeosciences↗
USGS2006

Effects of permafrost melting on CO2 and CH4 exchange of a poorly drained black spruce lowland

[1] Permafrost melting is occurring in areas of the boreal forest region where large amounts of carbon (C) are stored in organic soils. We measured soil respiration, net CO 2 flux, and net CH 4 flux during May–September 2003 and March 2004 in a black spruce lowland in interior Alaska to better understand how permafrost thaw in poorly drained landscapes affects land‐atmosphere CO 2 and CH 4 exchange. Sites included peat soils underlain by permafrost at ∼0.4 m depth (permafrost plateau, PP), four thermokarst wetlands (TW) having no permafrost in the upper 2.2 m, and peat soils bordering the thermokarst wetlands having permafrost at ∼0.5 m depth (thermokarst edges, TE). Soil respiration rates were not significantly different among the sites, and 5‐cm soil temperature explained 50–91% of the seasonal variability in soil respiration within the sites. Groundcover vegetation photosynthesis (calculated as net CO 2 minus soil respiration) was significantly different among the sites (TW > TE > PP), which can be partly attributed to the difference in photosynthetically active radiation reaching the ground at each site type. Methane emission rates were 15 to 28 times greater from TW than from TE and PP. We modeled annual soil respiration and groundcover vegetation photosynthesis using soil temperature and radiation data, and CH 4 flux by linear interpolation. We estimated all sites as net C gas sources to the atmosphere (not including tree CO 2 uptake at PP and TE), although the ranges in estimates when accounting for errors were large enough that TE and TW may have been net C sinks.

Journal of Geophysical Research G: Biogeosciences↗
USGS2007

Crustal deformation of the Yellowstone-Snake River Plain volcano-tectonic system-Campaign and continuous GPS observations, 1987-2004

The Yellowstone-Snake River Plain tectonomagmatic province resulted from Late Tertiary volcanism in western North America, producing three large, caldera-forming eruptions at the Yellowstone Plateau in the last 2 Myr. To understand the kinematics and geodynamics of this volcanic system, the University of Utah conducted seven GPS campaigns at 140 sites between 1987 and 2003 and installed a network of 15 permanent stations. GPS deployments focused on the Yellowstone caldera, the Hebgen Lake and Teton faults, and the eastern Snake River Plain. The GPS data revealed periods of uplift and subsidence of the Yellowstone caldera at rates up to 15 mm/yr. From 1987 to 1995, the caldera subsided and contracted, implying volume loss. From 1995 to 2000, deformation shifted to inflation and extension northwest of the caldera. From 2000 to 2003, uplift continued to the northwest while caldera subsidence was renewed. The GPS observations also revealed extension across the Hebgen Lake fault and fault-normal contraction across the Teton fault. Deformation rates of the Yellowstone caldera and Hebgen Lake fault were converted to equivalent total moment rates, which exceeded historic seismic moment release and late Quaternary fault slip-derived moment release by an order of magnitude. The Yellowstone caldera deformation trends were superimposed on regional southwest extension of the Yellowstone Plateau at up to 4.3 ± 0.2 mm/yr, while the eastern Snake River Plain moved southwest as a slower rate at 2.1 ± 0.2 mm/yr. This southwest extension of the Yellowstone-Snake River Plain system merged into east-west extension of the Basin-Range province. Copyright 2007 by the American Geophysical Union.

Idaho, Oregon, Utah, Wyoming↗
USGS2007

Cosmogenic 10Be and 36Cl geochronology of offset alluvial fans along the northern Death Valley fault zone: Implications for transient strain in the eastern California shear zone

The northern Death Valley fault zone (NDVFZ) has long been recognized as a major right-lateral strike-slip fault in the eastern California shear zone (ECSZ). However, its geologic slip rate has been difficult to determine. Using high-resolution digital topographic imagery and terrestrial cosmogenic nuclide dating, we present the first geochronologically determined slip rate for the NDVFZ. Our study focuses on the Red Wall Canyon alluvial fan, which exposes clean dextral offsets of seven channels. Analysis of airborne laser swath mapping data indicates ???297 ?? 9 m of right-lateral displacement on the fault system since the late Pleistocene. In situ terrestrial cosmogenic 10Be and 36C1 geochronology was used to date the Red Wall Canyon fan and a second, correlative fan also cut by the fault. Beryllium 10 dates from large cobbles and boulders provide a maximum age of 70 +22/-20 ka for the offset landforms. The minimum age of the alluvial fan deposits based on 36Cl depth profiles is 63 ?? 8 ka. Combining the offset measurement with the cosmogenic 10Be date yields a geologic fault slip rate of 4.2 +1.9/-1.1 mm yr-1, whereas the 36Cl data indicate 4.7 +0.9/-0.6 mm yr-1 of slip. Summing these slip rates with known rates on the Owens Valley, Hunter Mountain, and Stateline faults at similar latitudes suggests a total geologic slip rate across the northern ECSZ of ???8.5 to 10 mm yr-1. This rate is commensurate with the overall geodetic rate and implies that the apparent discrepancy between geologic and geodetic data observed in the Mojave section of the ECSZ does not extend north of the Garlock fault. Although the overall geodetic rates are similar, the best estimates based on geology predict higher strain rates in the eastern part of the ECSZ than to the west, whereas the observed geodetic strain is relatively constant. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

M ≥ 7.0 earthquake recurrence on the San Andreas fault from a stress renewal model

Forecasting M ≥ 7.0 San Andreas fault earthquakes requires an assessment of their expected frequency. I used a three-dimensional finite element model of California to calculate volumetric static stress drops from scenario M ≥ 7.0 earthquakes on three San Andreas fault sections. The ratio of stress drop to tectonic stressing rate derived from geodetic displacements yielded recovery times at points throughout the model volume. Under a renewal model, stress recovery times on ruptured fault planes can be a proxy for earthquake recurrence. I show curves of magnitude versus stress recovery time for three San Andreas fault sections. When stress recovery times were converted to expected M ≥ 7.0 earthquake frequencies, they fit Gutenberg-Richter relationships well matched to observed regional rates of M ≤ 6.0 earthquakes. Thus a stress-balanced model permits large earthquake Gutenberg-Richter behavior on an individual fault segment, though it does not require it. Modeled slip magnitudes and their expected frequencies were consistent with those observed at the Wrightwood paleoseismic site if strict time predictability does not apply to the San Andreas fault.

Journal of Geophysical Research B: Solid Earth↗
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↗
USGS2007

Earthquake-by-earthquake fold growth above the Puente Hills blind thrust fault, Los Angeles, California: Implications for fold kinematics and seismic hazard

Boreholes and high-resolution seismic reflection data collected across the forelimb growth triangle above the central segment of the Puente Hills thrust fault (PHT) beneath Los Angeles, California, provide a detailed record of incremental fold growth during large earthquakes on this major blind thrust fault. These data document fold growth within a discrete kink band that narrows upward from ???460 m at the base of the Quaternary section (200-250 m depth) to <150 m at 2.5 m depth, with most growth during the most recent folding event occurring within a zone only ???60 m wide. These observations, coupled with evidence from petroleum industry seismic reflection data, demonstrate that most (>82% at 250 m depth) folding and uplift occur within discrete kink bands, thereby enabling us to develop a paleoseismic history of the underlying blind thrust fault. The borehole data reveal that the youngest part of the growth triangle in the uppermost 20 m comprises three stratigraphically discrete growth intervals marked by southward thickening sedimentary strata that are separated by intervals in which sediments do not change thickness across the site. We interpret the intervals of growth as occurring after the formation of now-buried paleofold scarps during three large PHT earthquakes in the past 8 kyr. The intervening intervals of no growth record periods of structural quiescence and deposition at the regional, near-horizontal stream gradient at the study site. Minimum uplift in each of the scarp-forming events, which occurred at 0.2-2.2 ka (event Y), 3.0-6.3 ka (event X), and 6.6-8.1 ka (event W), ranged from ???1.1 to ???1.6 m, indicating minimum thrust displacements of ???2.5 to 4.5 m. Such large displacements are consistent with the occurrence of large-magnitude earthquakes (Mw > 7). Cumulative, minimum uplift in the past three events was 3.3 to 4.7 m, suggesting cumulative thrust displacement of ???7 to 10.5 m. These values yield a minimum Holocene slip rate for the PHT of ???0.9 to 1.6 mm/yr. The borehole and seismic reflection data demonstrate that dip within the kink band is acquired incrementally, such that older strata that have been deformed by more earthquakes dip more steeply than younger strata. Specifically, strata dip 0.4?? at 4 m depth, 0.7?? at 20 m depth, 8?? at 90 m, 16?? at 110 m, and 17?? at 200 m. Moreover, structural restorations of the borehole data show that the locus of active folding (the anticlinal active axial surface) does not extend to the surface in exactly the same location from earthquake to earthquake. Rather, that the axial surfaces migrate from earthquake to earthquake, reflecting a component of fold growth by kink band migration. The incremental acquisition of bed dip in the growth triangle may reflect some combination of fold growth by limb rotation in addition to kink band migration, possibly through a component of trishear or shear fault bend folding. Alternatively, the component of limb rotation may result from curved hinge fault bend folding, and/or the mechanical response of loosely consolidated granular sediments in the shallow subsurface to folding at depth. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Mechanical properties of sand, silt, and clay containing tetrahydrofuran hydrate

The mechanical behavior of hydrate-bearing sediments subjected to large strains has relevance for the stability of the seafloor and submarine slopes, drilling and coring operations, and the analysis of certain small-strain properties of these sediments (for example, seismic velocities). This study reports on the results of comprehensive axial compression triaxial tests conducted at up to 1 MPa confining pressure on sand, crushed silt, precipitated silt, and clay specimens with closely controlled concentrations of synthetic hydrate. The results show that the stress-strain behavior of hydrate-bearing sediments is a complex function of particle size, confining pressure, and hydrate concentration. The mechanical properties of hydrate-bearing sediments at low hydrate concentration (probably < 40% of pore space) appear to be determined by stress-dependent soil stiffness and strength. At high hydrate concentrations (>50% of pore space), the behavior becomes more independent of stress because the hydrates control both stiffness and strength and possibly the dilative tendency of sediments by effectively increasing interparticle coordination, cementing particles together, and filling the pore space. The cementation contribution to the shear strength of hydrate-bearing sediments decreases with increasing specific surface of soil minerals. The lower the effective confining stress, the greater the impact of hydrate formation on normalized strength.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Hydrogeologic controls imposed by mechanical stratigraphy in layered rocks of the Chateauguay River Basin, a U.S.-Canada transborder aquifer

[1] The Ch&acirc;teauguay River Basin delineates a transborder watershed with roughly half of its surface area located in northern New York State and half in southern Qu&eacute;bec Province, Canada. As part of a multidisciplinary study designed to characterize the hydrogeologic properties of this basin, geophysical logs were obtained in 12 wells strategically located to penetrate the four major sedimentary rock formations that constitute the regional aquifers. The layered rocks were classified according to their elastic properties into three primary units: soft sandstone, hard sandstone, and dolostone. Downhole measurements were analyzed to identify fracture patterns associated with each unit and to evaluate their role in controlling groundwater flow. Fracture networks are composed of orthogonal sets of laterally extensive, subhorizontal bedding plane partings and bed-delimited, subvertical joints with spacings that are consistent with rock mechanics principles and stress models. The vertical distribution of transmissive zones is confined to a few select bedding plane fractures, with soft sandstone having the fewest (one per 70-m depth) and hard sandstone the most (five per 70-m depth). Bed-normal permeability is examined using a probabilistic model that considers the lengths of flow paths winding along joints and bedding plane fractures. Soft sandstone has the smallest bed-normal permeability primarily because of its wide, geomechanically undersaturated joint spacing. Results indicate that the three formations have similar values of bulk transmissivity, within roughly an order of magnitude, but that each rock unit has its own unique system of groundwater flow paths that constitute that transmissivity.

Chateauguay River Basin↗
USGS2007

Waveform inversion of volcano-seismic signals for an extended source

[1] We propose a method to investigate the dimensions and oscillation characteristics of the source of volcano-seismic signals based on waveform inversion for an extended source. An extended source is realized by a set of point sources distributed on a grid surrounding the centroid of the source in accordance with the source geometry and orientation. The source-time functions for all point sources are estimated simultaneously by waveform inversion carried out in the frequency domain. We apply a smoothing constraint to suppress short-scale noisy fluctuations of source-time functions between adjacent sources. The strength of the smoothing constraint we select is that which minimizes the Akaike Bayesian Information Criterion (ABIC). We perform a series of numerical tests to investigate the capability of our method to recover the dimensions of the source and reconstruct its oscillation characteristics. First, we use synthesized waveforms radiated by a kinematic source model that mimics the radiation from an oscillating crack. Our results demonstrate almost complete recovery of the input source dimensions and source-time function of each point source, but also point to a weaker resolution of the higher modes of crack oscillation. Second, we use synthetic waveforms generated by the acoustic resonance of a fluid-filled crack, and consider two sets of waveforms dominated by the modes with wavelengths 2 L /3 and 2 W /3, or L and 2 L /5, where W and L are the crack width and length, respectively. Results from these tests indicate that the oscillating signature of the 2 L /3 and 2 W /3 modes are successfully reconstructed. The oscillating signature of the L mode is also well recovered, in contrast to results obtained for a point source for which the moment tensor description is inadequate. However, the oscillating signature of the 2 L /5 mode is poorly recovered owing to weaker resolution of short-scale crack wall motions. The triggering excitations of the oscillating cracks are successfully reconstructed.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Postseismic relaxation and aftershocks

[1] Perfettini et al. (2005) suggested that the temporal dependence of surface displacements u ( t ) measured in the epicentral area following an earthquake is related to N ( t ), the cumulative number of aftershocks, by the equation u ( t ) = a + bt + cN ( t ) + d (1 − e − αt ), where a , b , c , d , and α are constants chosen to fit the data and t is the postearthquake time. N ( t ) appears in the expression for u ( t ) because both the aftershocks and a portion of u ( t ) are thought to be driven by the same source, postseismic fault creep at subseismogenic depths on the downdip extension of the coseismic rupture. We show that this equation with the actually observed N ( t ) fits the postseismic displacements recorded on several baselines following each of five earthquakes: 1999 M7.6 Chi-Chi (Taiwan), 1999 M7.1 Hector Mine (southern California), 2002 M7.9 Denali (central Alaska), 2003 M6.5 San Simeon (central California), and 2004 M6.0 Parkfield (central California) earthquakes. Although there are plausible physical interpretations for each of the terms in the expression for u ( t ), the large number of adjustable constants ( a , b , c , d , and α ) involved in fitting the rather simple postseismic displacements diminishes the significance of the fit. Because the observed N ( t ) is well fit by the modified Omori's law, fault creep at depth presumably exhibits the same temporal dependence. That dependence could be explained if the rheology of the fault downdip from the coseismic rupture is consistent with ordinary transient creep. Montesi (2004) demonstrated that power law creep across a shear zone at depth would also produce that temporal signal.

Journal of Geophysical Research B: Solid Earth↗