USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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 379 records · Page 21Linked to original sources
USGS2000

Crustal deformation associated with glacial fluctuations in the eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993–1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 km 3 from the surge reservoir area to the terminus region. We examined the background seismicity (M L ≥ 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M L ≤ 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, σ 3 , could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1–6 m/yr in the ablation region is 1–12 mm/yr. For our reference model with a viscosity of 5×10 19 Pa s for depths between ≈ 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in σv between 0 and 10 km was ≈ 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (≈ 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M S = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Estimation of wave phase speed and nearshore bathymetry from video imagery

A new remote sensing technique based on video image processing has been developed for the estimation of nearshore bathymetry. The shoreward propagation of waves is measured using pixel intensity time series collected at a cross-shore array of locations using remotely operated video cameras. The incident band is identified, and the cross-spectral matrix is calculated for this band. The cross-shore component of wavenumber is found as the gradient in phase of the first complex empirical orthogonal function of this matrix. Water depth is then inferred from linear wave theory's dispersion relationship. Full bathymetry maps may be measured by collecting data in a large array composed of both cross-shore and longshore lines. Data are collected hourly throughout the day, and a stable, daily estimate of bathymetry is calculated from the median of the hourly estimates. The technique was tested using 30 days of hourly data collected at the SandyDuck experiment in Duck, North Carolina, in October 1997. Errors calculated as the difference between estimated depth and ground truth data show a mean bias of −35 cm (rms error = 91 cm). Expressed as a fraction of the true water depth, the mean percent error was 13% (rms error = 34%). Excluding the region of known wave nonlinearities over the bar crest, the accuracy of the technique improved, and the mean (rms) error was −20 cm (75 cm). Additionally, under low-amplitude swells (wave height H ≤ 1 m), the performance of the technique across the entire profile improved to 6% (29%) of the true water depth with a mean (rms) error of −12 cm (71 cm).

North Carolina↗
USGS2001

Sources of global warming in upper ocean temperature during El Niño

Global average sea surface temperature (SST) from 40°S to 60°N fluctuates ±0.3°C on interannual period scales, with global warming (cooling) during El Niño (La Niña). About 90% of the global warming during El Niño occurs in the tropical global ocean from 20°S to 20°N, half because of large SST anomalies in the tropical Pacific associated with El Niño and the other half because of warm SST anomalies occurring over ∼80% of the tropical global ocean. From examination of National Centers for Environmental Prediction [ Kalnay et al. , 1996] and Comprehensive Ocean-Atmosphere Data Set [ Woodruff et al. , 1993] reanalyses, tropical global warming during El Niño is associated with higher troposphere moisture content and cloud cover, with reduced trade wind intensity occurring during the onset phase of El Niño. During this onset phase the tropical global average diabatic heat storage tendency in the layer above the main pycnocline is 1–3 W m −2 above normal. Its principal source is a reduction in the poleward Ekman heat flux out of the tropical ocean of 2–5 W m −2 . Subsequently, peak tropical global warming during El Niño is dissipated by an increase in the flux of latent heat to the troposphere of 2–5 W m −2 , with reduced shortwave and longwave radiative fluxes in response to increased cloud cover tending to cancel each other. In the extratropical global ocean the reduction in poleward Ekman heat flux out of the tropics during the onset of El Niño tends to be balanced by reduction in the flux of latent heat to the troposphere. Thus global warming and cooling during Earth's internal mode of interannual climate variability arise from fluctuations in the global hydrological balance, not the global radiation balance. Since it occurs in the absence of extraterrestrial and anthropogenic forcing, global warming on decadal, interdecadal, and centennial period scales may also occur in association with Earth's internal modes of climate variability on those scales.

Journal of Geophysical Research C: Oceans↗
USGS2000

Time-averaged fluxes of lead and fallout radionuclides to sediments in Florida Bay

Recent, unmixed sediments from mud banks of central Florida Bay were dated using 210 Pb/ 226 Ra, and chronologies were verified by comparing sediment lead temporal records with Pb/Ca ratios in annual layers of coral ( Montastrea annularis ) located on the ocean side of the Florida Keys. Dates of sediment lead peaks (1978±2) accord with prior observations of a 6 year lag between the occurrence of maximum atmospheric lead in 1972 and peak coral lead in 1978. Smaller lags of 1–2 years occur between the maximum atmospheric radionuclide fallout and peaks in sediment temporal records of 137 Cs and Pu. Such lags are consequences of system time averaging (STA) in which atmospherically delivered particle-associated constituents accumulate and mix in a (sedimentary?) reservoir before transferring to permanent sediments and coral. STA model calculations, using time-dependent atmospheric inputs, produced optimized profiles in excellent accord with measured sediment 137 Cs, Pu, lead, and coral lead distributions. Derived residence times of these particle tracers (16±1, 15.7±0.7, 19±3, and 16±2 years, respectively) are comparable despite differences in sampling locations, in accumulating media, and in element loading histories and geochemical properties. For a 16 year weighted mean residence time, STA generates the observed 6 year lead peak lag. Evidently, significant levels of nondegradable, particle-associated contaminants can persist in Florida Bay for many decades following elimination of external inputs. Present results, in combination with STA model analysis of previously reported radionuclide profiles, suggest that decade-scale time averaging may occur widely in recent coastal marine sedimentary environments.

Florida↗
USGS1999

Spatially averaged flow over a wavy boundary revisited

Vertical profiles of streamwise velocity measured over bed forms are commonly used to deduce boundary shear stress for the purpose of estimating sediment transport. These profiles may be derived locally or from some sort of spatial average. Arguments for using the latter procedure are based on the assumption that spatial averaging of the momentum equation effectively removes local accelerations from the problem. Using analogies based on steady, uniform flows, it has been argued that the spatially averaged velocity profiles are approximately logarithmic and can be used to infer values of boundary shear stress. This technique of using logarithmic profiles is investigated using detailed laboratory measurements of flow structure and boundary shear stress over fixed two-dimensional bed forms. Spatial averages over the length of the bed form of mean velocity measurements at constant distances from the mean bed elevation yield vertical profiles that are highly logarithmic even though the effect of the bottom topography is observed throughout the water column. However, logarithmic fits of these averaged profiles do not yield accurate estimates of the measured total boundary shear stress.

Journal of Geophysical Research C: Oceans↗
USGS1999

Precipitation structure in the Sierra Nevada of California during winter

Influences of upper air characteristics along the coast of California upon winter-time (November-April) precipitation in the Sierra Nevada are investigated. Precipitation events in the Sierra Nevada region occur mostly during wintertime, irrespective of station location (leeside or windside) and elevation. Most precipitation episodes in the region are associated with moist southwesterly winds (coming from the southwest direction) and also tend to occur when the 700-mbar temperature at the upwind direction is close to −2°C. This favored wind direction and temperature signify the importance of both moisture transport and orographic lifting in augmenting precipitation in the region. By utilizing the observed dependency of the precipitation upon the upper air conditions, a linear model is formulated to quantify the precipitation observed at different sites as a function of moisture transport. The skill of the model increases with timescale of aggregation, reaching more than 50% variance explained at an aggregation period of 5–7 days. This indicates that upstream air moisture transport can be used to estimate the precipitation totals in the Sierra Nevada region.

Journal of Geophysical Research D: Atmospheres↗
USGS1999

Potential impacts on Colorado Rocky Mountain weather due to land use changes on the adjacent Great Plains

Evidence from both meteorological stations and vegetational successional studies suggests that summer temperatures are decreasing in the mountain-plain system in northeast Colorado, particularly since the early 1980s. These trends are coincident with large changes in regional land cover. Trends in global, Northern Hemisphere and continental surface temperatures over the same period are insignificant. These observations suggest that changes in the climate of this mountain-plain system may be, in some part, a result of localized forcing mechanisms. In this study the effects of land use change on the northern Colorado plains, where large regions of grasslands have been transformed into both dry and irrigated agricultural lands, on regional weather is examined in an effort to understand this local deviation from larger-scale trends. We find with high-resolution numerical simulations of a 3-day summer period using a regional atmospheric-land surface model that replacing grasslands with irrigated and dry farmland can have impacts on regional weather and therefore climate which are not limited to regions of direct forcing. Higher elevations remote from regions of land use change are affected as well. Specifically, cases with altered landcover had cooler, moister boundary layers, and diminished low-level upslope winds over portions of the plains. At higher elevations, temperatures also were lower as was low-level convergence. Precipitation and cloud cover were substantially affected in mountain regions. We advance the hypothesis that observed land use changes may have already had a role in explaining part of the observed climate record in the northern Colorado mountain-plain system.

Journal of Geophysical Research D: Atmospheres↗
USGS1999

Nevada STORMS project: Measurement of mercury emissions from naturally enriched surfaces

Diffuse anthropogenic and naturally mercury-enriched areas represent long-lived sources of elemental mercury to the atmosphere. The Nevada Study and Tests of the Release of Mercury From Soils (STORMS) project focused on the measurement of mercury emissions from a naturally enriched area. During the project, concurrent measurements of mercury fluxes from naturally mercury-enriched substrate were made September 1–4, 1997, using four micrometeorological methods and seven field flux chambers. Ambient air mercury concentrations ranged from 2 to nearly 200 ng m −3 indicating that the field site is a source of atmospheric mercury. The mean daytime mercury fluxes, during conditions of no precipitation, measured with field chambers were 50 to 360 ng m −2 h −1 , and with the micrometeorological methods were 230 to 600 ng m −2 h −1 . This wide range in mercury emission rates reflects differences in method experimental designs and local source strengths. Mercury fluxes measured by many field chambers were significantly different ( p <0.05) but linearly correlated. This indicates that field chambers responded similarly to environmental conditions, but differences in experimental design and site heterogeneity had a significant influence on the magnitude of mercury fluxes. Data developed during the field study demonstrated that field flux chambers are ideal for assessment of the physicochemical processes driving mercury flux and development of an understanding of the magnitude of the influence of individual factors on flux. In general, mean mercury fluxes measured with micrometeorological methods during daytime periods were nearly 3 times higher than mean fluxes measured with field flux chambers. Micrometeorological methods allow for derivation of a representative mercury flux occurring from an unconstrained system and provide an assessment of the actual magnitude and variability of fluxes occurring from an area.

Journal of Geophysical Research D: Atmospheres↗
USGS2000

Chamber measurement of surface-atmosphere trace gas exchange: Numerical evaluation of dependence on soil interfacial layer, and source/sink products

We employed a three-dimensional finite difference gas diffusion model to simulate the performance of chambers used to measure surface-atmosphere tace gas exchange. We found that systematic errors often result from conventional chamber design and deployment protocols, as well as key assumptions behind the estimation of trace gas exchange rates from observed concentration data. Specifically, our simulationshowed that (1) when a chamber significantly alters atmospheric mixing processes operating near the soil surface, it also nearly instantaneously enhances or suppresses the postdeployment gas exchange rate, (2) any change resulting in greater soil gas diffusivity, or greater partitioning of the diffusing gas to solid or liquid soil fractions, increases the potential for chamber-induced measurement error, and (3) all such errors are independent of the magnitude, kinetics, and/or distribution of trace gas sources, but greater for trace gas sinks with the same initial absolute flux. Finally, and most importantly, we found that our results apply to steady state as well as non-steady-state chambers, because the slow rate of gas diffusion in soil inhibits recovery of the former from their initial non-steady-state condition. Over a range of representative conditions, the error in steady state chamber estimates of the trace gas flux varied from -30 to +32%, while estimates computed by linear regression from non-steadystate chamber concentrations were 2 to 31% too small. Although such errors are relatively small in comparison to the temporal and spatial variability characteristic of trace gas exchange, they bias the summary statistics for each experiment as well as larger scale trace gas flux estimates based on them.

Journal of Geophysical Research D: Atmospheres↗
USGS2001

The failure of earthquake failure models

In this study I show that simple heuristic models and numerical calculations suggest that an entire class of commonly invoked models of earthquake failure processes cannot explain triggering of seismicity by transient or "dynamic" stress changes, such as stress changes associated with passing seismic waves. The models of this class have the common feature that the physical property characterizing failure increases at an accelerating rate when a fault is loaded (stressed) at a constant rate. Examples include models that invoke rate state friction or subcritical crack growth, in which the properties characterizing failure are slip or crack length, respectively. Failure occurs when the rate at which these grow accelerates to values exceeding some critical threshold. These accelerating failure models do not predict the finite durations of dynamically triggered earthquake sequences (e.g., at aftershock or remote distances). Some of the failure models belonging to this class have been used to explain static stress triggering of aftershocks. This may imply that the physical processes underlying dynamic triggering differs or that currently applied models of static triggering require modification. If the former is the case, we might appeal to physical mechanisms relying on oscillatory deformations such as compaction of saturated fault gouge leading to pore pressure increase, or cyclic fatigue. However, if dynamic and static triggering mechanisms differ, one still needs to ask why static triggering models that neglect these dynamic mechanisms appear to explain many observations. If the static and dynamic triggering mechanisms are the same, perhaps assumptions about accelerating failure and/or that triggering advances the failure times of a population of inevitable earthquakes are incorrect.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Permeability reduction in granite under hydrothermal conditions

The formation of impermeable fault seals between earthquake events is a feature of many models of earthquake generation, suggesting that earthquake recurrence may depend in part on the rate of permeability reduction of fault zone materials under hydrothermal conditions. In this study, permeability measurements were conducted on intact, fractured, and gouge-bearing Westerly granite at an effective pressure of 50 MPa and at temperatures from 150° to 500°C, simulating conditions in the earthquake-generating portions of fault zones. Pore fluids were cycled back and forth under a 2 MPa pressure differential for periods of up to 40 days. Permeability of the granite decreased with time t , following the exponential relation k = c (10 − rt ). For intact samples run between 250° and 500°C the time constant for permeability decrease r was proportional to temperature and ranged between 0.001 and 0.1 days −1 (i.e., between 0.4 and 40 decades year −1 loss of permeability). Values of r for the lower-temperature experiments differed little from the 250°C runs. In contrast, prefractured samples showed higher rates of permeability decrease at a given temperature. The surfaces of the fractured samples showed evidence of dissolution and mineral growth that increased in abundance with both temperature and time. The experimentally grown mineral assemblages varied with temperature and were consistent with a rock-dominated hydrothermal system. As such mineral deposits progressively seal the fractured samples, their rates of permeability decrease approach the rates for intact rocks at the same temperature. These results place constraints on models of precipitation sealing and suggest that fault rocks may seal at a rate consistent with earthquake recurrence intervals of typical fault zones.

California↗
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↗
USGS2001

Strain accumulation and rotation in the Eastern California Shear Zone

Although the Eastern California Shear Zone (ECSZ) (strike ∼N25°W) does not quite coincide with a small circle drawn about the Pacific-North America pole of rotation, trilateration and GPS measurements demonstrate that the motion within the zone corresponds to right-lateral simple shear across a vertical plane (strike N33°W±5°) roughly parallel to the tangent to that local small circle (strike ∼N40°W). If the simple shear is released by slip on faults subparallel to the shear zone, the accumulated rotation is also released, leaving no secular rotation. South of the Garlock fault the principal faults (e.g., Calico-Blackwater fault) strike ∼N40°W, close enough to the strike of the vertical plane across which maximum right-lateral shear accumulates to almost wholly accommodate that accumulation of both strain and rotation by right-lateral slip. North of the Garlock fault dip slip as well as strike slip on the principal faults (strike ∼N20°W) is required to accommodate the simple shear accumulation. In both cases the accumulated rotation is released with the shear strain. The Garlock fault, which transects the ECSZ, is not offset by north-northwest striking faults nor, despite geological evidence for long-term left-lateral slip, does it appear at the present time to be accumulating left-lateral simple shear strain across the fault due to slip at depth. Rather the motion is explained by right-lateral simple shear across the orthogonal ECSZ. Left-lateral slip on the Garlock fault will release the shear strain accumulating there but would augment the accumulating rotation, resulting in a secular clockwise rotation rate ∼80 nrad yr −1 (4.6° Myr −1 ).

California↗
USGS2002

Complex earthquake rupture and local tsunamis

In contrast to far-field tsunami amplitudes that are fairly well predicted by the seismic moment of subduction zone earthquakes, there exists significant variation in the scaling of local tsunami amplitude with respect to seismic moment. From a global catalog of tsunami runup observations this variability is greatest for the most frequently occurring tsunamigenic subduction zone earthquakes in the magnitude range of 7 < M w < 8.5. Variability in local tsunami runup scaling can be ascribed to tsunami source parameters that are independent of seismic moment: variations in the water depth in the source region, the combination of higher slip and lower shear modulus at shallow depth, and rupture complexity in the form of heterogeneous slip distribution patterns. The focus of this study is on the effect that rupture complexity has on the local tsunami wave field. A wide range of slip distribution patterns are generated using a stochastic, self-affine source model that is consistent with the falloff of far-field seismic displacement spectra at high frequencies. The synthetic slip distributions generated by the stochastic source model are discretized and the vertical displacement fields from point source elastic dislocation expressions are superimposed to compute the coseismic vertical displacement field. For shallow subduction zone earthquakes it is demonstrated that self-affine irregularities of the slip distribution result in significant variations in local tsunami amplitude. The effects of rupture complexity are less pronounced for earthquakes at greater depth or along faults with steep dip angles. For a test region along the Pacific coast of central Mexico, peak nearshore tsunami amplitude is calculated for a large number ( N = 100) of synthetic slip distribution patterns, all with identical seismic moment ( M w = 8.1). Analysis of the results indicates that for earthquakes of a fixed location, geometry, and seismic moment, peak nearshore tsunami amplitude can vary by a factor of 3 or more. These results indicate that there is substantially more variation in the local tsunami wave field derived from the inherent complexity subduction zone earthquakes than predicted by a simple elastic dislocation model. Probabilistic methods that take into account variability in earthquake rupture processes are likely to yield more accurate assessments of tsunami hazards.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗
USGS2000

Geochemical evidence for a magmatic CO 2 degassing event at Mammoth Mountain, California, September-December 1997

Recent time series soil CO 2 concentration data from monitoring stations in the vicinity of Mammoth Mountain, California, reveal strong evidence for a magmatic degassing event during the fall of 1997 lasting more than 2 months. Two sensors at Horseshoe Lake first recorded the episode on September 23, 1997, followed 10 days later by a sensor on the north flank of Mammoth Mountain. Direct degassing from shallow intruding magma seems an implausible cause of the degassing event, since the gas released at Horseshoe Lake continued to be cold and barren of other magmatic gases, except for He. We suggest that an increase in compressional strain on the area south of Mammoth Mountain driven by movement of major fault blocks in Long Valley caldera may have triggered an episode of increased degassing by squeezing additional accumulated CO 2 from a shallow gas reservoir to the surface along faults and other structures where it could be detected by the CO 2 monitoring network. Recharge of the gas reservoir by CO 2 emanating from the deep intrusions that probably triggered deep long-period earthquakes may also have contributed to the degassing event. The nature of CO 2 discharge at the soil-air interface is influenced by the porous character of High Sierra soils and by meteorological processes. Solar insolation is the primary source of energy for the Earth atmosphere and plays a significant role in most diurnal processes at the Earth surface. Data from this study suggest that external forcing due largely to local orographic winds influences the fine structure of the recorded CO 2 signals.

California↗
USGS2000

Seismovolcanic signals at Deception Island volcano, Antarctica: Wave field analysis and source modeling

The seismovolcanic signals associated with the volcanic activity of Deception Island (Antarctica), recorded during three Antarctic summers (1994–1995, 1995–1996 and 1996–1997), are analyzed using a dense small-aperture (500 m) seismic array. The visual and spectral classification of the seismic events shows the existence of long-period and hybrid isolated seismic events, and of low-frequency, quasi-monochromatic and spasmodic continuous tremors. All spectra have the highest amplitudes in the frequency band between 1 and 4 Hz, while hybrids and spasmodic tremors have also significant amplitudes in the high-frequency band (4–10 Hz). The array analysis indicates that almost all the well-correlated low-frequency signals share similar array parameters (slowness and back azimuth) and have the same source area, close to the array site. The polarization analysis shows that phases at high-frequency are mostly composed of P waves, and those phases dominated by low frequencies can be interpreted as surface waves. No clear shear waves are evidenced. From the energy evaluation, we have found that the reduced displacement values for surface and body waves are confined in a narrow interval. Volcano-tectonic seismicity is located close to the array, at a depth shallower than 1 km. The wave-field properties of the seismovolcanic signals allow us to assume a unique source model, a shallow resonating fluid-filled crack system at a depth of some hundreds of meters. All of the seismic activity is interpreted as the response of a reasonably stable stationary geothermal process. The differences observed in the back azimuth between low and high frequencies are a near-field effect. A few episodes of the degassification process in an open conduit were observed and modeled with a simple organ pipe.

Journal of Geophysical Research: Solid Earth↗