USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics Research”

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
USGS2002

Concordant paleolatitudes for Neoproterozoic ophiolitic rocks of the Trinity Complex, Klamath Mountains, California

New paleomagnetic results from the eastern Klamath Mountains of northern California show that Neoproterozoic rocks of the Trinity ophiolitic complex and overlying Middle Devonian volcanic rocks are latitudinally concordant with cratonal North America. Combining paleomagnetic data with regional geologic and faunal evidence suggests that the Trinity Complex and related terranes of the eastern Klamath plate were linked in some fashion to the North American craton throughout that time, but that distance between them may have varied considerably. A possible model that is consistent with our paleomagnetic results and the geologic evidence is that the Trinity Complex formed and migrated parallel to paleolatitude in the basin between Laurasia and Australia–East Antarctica as the Rodinian supercontinent began to break up. It then continued to move parallel to paleolatitude at least through Middle Devonian time. Although the eastern Klamath plate served as a nucleus against which more western components of the Klamath Mountains province amalgamated, the Klamath superterrane was not accreted to North America until Early Cretaceous time.

California↗
USGS2002

Episodic thermal perturbations associated with groundwater flow: An example from Kilauea Volcano, Hawaii

Temperature measurements in deep drill holes on volcano summits or upper flanks allow a quantitative analysis of groundwater induced heat transport within the edifice. We present a new temperature-depth profile from a deep well on the summit of Kilauea Volcano, Hawaii, and analyze it in conjunction with a temperature profile measured 26 years earlier. We propose two groundwater flow models to interpret the complex temperature profiles. The first is a modified confined lateral flow model (CLFM) with a continuous flux of hydrothermal fluid. In the second, transient flow model (TFM), slow conductive cooling follows a brief, advective heating event. We carry out numerical simulations to examine the timescales associated with each of the models. Results for both models are sensitive to the initial conditions, and with realistic initial conditions it takes between 750 and 1000 simulation years for either model to match the measured temperature profiles. With somewhat hotter initial conditions, results are consistent with onset of a hydrothermal plume ∼550 years ago, coincident with initiation of caldera subsidence. We show that the TFM is consistent with other data from hydrothermal systems and laboratory experiments and perhaps is more appropriate for this highly dynamic environment. The TFM implies that volcano-hydrothermal systems may be dominated by episodic events and that thermal perturbations may persist for several thousand years after hydrothermal flow has ceased.

Hawaii↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Fire dynamics and implications for nitrogen cycling in boreal forests

We used a dynamic, long-term mass balance approach to track cumulative carbon (C) and nitrogen (N) losses to fire in boreal Manitoba over the 6500 years since deglaciation. Estimated C losses to decomposition and fire, combined with measurements of N pools in mature and burned forest floors, suggest that loss of N by combustion has likely resulted in a long-term loss that exceeds the amount of N stored in soil today by 2 to 3 times. These estimates imply that biological N fixation rates could be as high as 5 to 10 times atmospheric deposition rates in boreal regions. At the site scale, the amount of N lost is due to N content of fuels, which varies by stand type and fire severity, which in turn vary with climate and fire dynamics. The interplay of fire frequency, fire severity, and N partitioning during regrowth are important for understanding rates and sustainability of nutrient and carbon cycling over millenia and over broad regions.

Manitoba↗
USGS2001

Moss and lichen cover mapping at local and regional scales in the boreal forest ecosystem of central Canada

Mosses and lichens are important components of boreal landscapes [Vitt et al., 1994; Bubier et al., 1997]. They affect plant productivity and belowground carbon sequestration and alter the surface runoff and energy balance. We report the use of multiresolution satellite data to map moss and lichens over the BOREAS region at a 10 m, 30 m, and 1 km scales. Our moss and lichen classification at the 10 m scale is based on ground observations of associations among soil drainage classes, overstory composition, and cover type among four broad classes of ground cover (feather, sphagnum, and brown mosses and lichens). For our 30 m map, we used field observations of ground cover-overstory associations to map mosses and lichens in the BOREAS southern study area (SSA). To scale up to a 1 km (AVHRR) moss map of the BOREAS region, we used the TM SSA mosaics plus regional field data to identify AVHRR overstory-ground cover associations. We found that: 1) ground cover, overstory composition and density are highly correlated, permitting inference of moss and lichen cover from satellite-based land cover classifications; 2) our 1 km moss map reveals that mosses dominate the boreal landscape of central Canada, thereby a significant factor for water, energy, and carbon modeling; 3) TM and AVHRR moss cover maps are comparable; 4) satellite data resolution is important; particularly in detecting the smaller wetland features, lakes, and upland jack pine sites; and 5) distinct regional patterns of moss and lichen cover correspond to latitudinal and elevational gradients. Copyright 2001 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope composition of waters in the Great Basin, United States 1. Air-mass trajectories

Isentropic trajectories, calculated using the NOAA/Climate Monitoring and Diagnostics Laboratory's isentropic transport model, were used to determine air-parcel origins and the influence of air mass trajectories on the isotopic composition of precipitation events that occurred between October 1991 and September 1993 at Cedar City, Utah, and Winnemucca, Nevada. Examination of trajectories that trace the position of air parcels backward in time for 10 days indicated five distinct regions of water vapor origin: (1) Gulf of Alaska and North Pacific, (2) central Pacific, (3) tropical Pacific, (4) Gulf of Mexico, and (5) continental land mass. Deuterium (??D) and oxygen-18 (??18O) analyses were made of precipitation representing 99% of all Cedar City events. Similar analyses were made on precipitation representing 66% of the precipitation falling at Winnemucca during the same period. The average isotopic composition of precipitation derived from each water vapor source was determined. More than half of the precipitation that fell at both sites during the study period originated in the tropical Pacific and traveled northeast to the Great Basin; only a small proportion traversed the Sierra Nevada. The isotopic composition of precipitation is determined by air-mass origin and its track to the collection station, mechanism of droplet formation, reequilibration within clouds, and evaporation during its passage from cloud to ground. The Rayleigh distillation model can explain the changes in isotopic composition of precipitation as an air mass is cooled pseudo-adiabatically during uplift. However, the complicated processes that take place in the rapidly convecting environment of cumulonimbus and other clouds that are common in the Great Basin, especially in summer, require modification of this model because raindrops that form in the lower portion of those clouds undergo isotopic change as they are elevated to upper levels of the clouds from where they eventually drop to the ground.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope compositions of waters in the Great Basin, United States 2. Modern precipitation

Precipitation was collected between 1991 and 1997 at 41 locations within and adjacent to parts of the Great Basin lying in California, Oregon, Nevada, and Utah. These samples were analyzed for their deuterium (??D) and oxygen-18 (??18O) contents. Separate collections were made of summer and winter season precipitation at stations ranging in elevation from -65 m to 3246 m. The ??D per mil values of stations that were closely spaced but at different elevations showed an average ??D decrease of approximately 10???/km rise in elevation. Data for all samples representing winter precipitation, when plotted on a ??D versus ??18O plot, fall close to the Meteoric Water Line (??D = 8 ??18O + 10); samples representing summer precipitation define a line of slightly lower slope due to evaporation of the raindrops during their passage from cloud to ground. Comparison of our 1991-1997 ??D data with those from the same three stations reported by an earlier study in the southeastern California shows seasonal differences ranging from 0 per mil to 19??? (average: 15) and annual differences ranging from 0 to 13 per mil (average: 2), illustrating the degree of annual and seasonal variability in this region. When contoured, the ??D values display gradients indicating a north to northwest decrease in deuterium, with values ranging from -60 to -125??? in winter precipitation and from -40 to -110??? in summer precipitation. These gradient trends can be explained by the predominance of air mass trajectories originating in the tropical Pacific, the Gulf of California, and (in summer) the Gulf of Mexico.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope compositions of waters in the Great Basin, United States 3. Comparison of groundwaters with modern precipitation

Groundwater samples from wells and springs, scattered over most of the Great Basin province, were collected and analyzed for their isotopic makeup. They were augmented by previously published isotopic data on groundwaters from southeast California and by several hundred unpublished isotopic analyses. The ratio of 2H (deuterium, D) to 1H, in water samples from valleys in parts of California, Idaho, Nevada, Oregon, and Utah, are here compared with the winter, summer, and annual isotopic compositions of precipitation falling in or near the sampled areas. The main goal of this study was to identify basins where the groundwaters have isotopic compositions that are "lighter" (depleted in the heavier isotope, D) relative to modern winter precipitation. Where these basins do not adjoin substantially higher terrain, we consider those light groundwaters to be of Pleistocene age and thus more than 10,000 years old. Where the groundwater is 10 to 19??? lighter than local winter precipitation, we consider it to be possibly an indication of Pleistocene water; where the ??D makeup is >20??? lighter, we consider it to be probably Pleistocene water. More than 80 sites underlain by waters of possible or probable Pleistocene age were identified.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Modeling soil thermal and carbon dynamics of a fire chronosequence in interior Alaska

In this study, the dynamics of soil thermal, hydrologic, and ecosystem processes were coupled to project how the carbon budgets of boreal forests will respond to changes in atmospheric CO 2 , climate, and fire disturbance. The ability of the model to simulate gross primary production and ecosystem respiration was verified for a mature black spruce ecosystem in Canada, the age-dependent pattern of the simulated vegetation carbon was verified with inventory data on aboveground growth of Alaskan black spruce forests, and the model was applied to a postfire chronosequence in interior Alaska. The comparison between the simulated soil temperature and field-based estimates during the growing season (May to September) of 1997 revealed that the model was able to accurately simulate monthly temperatures at 10 cm ( R > 0.93) for control and burned stands of the fire chronosequence. Similarly, the simulated and field-based estimates of soil respiration for control and burned stands were correlated ( R = 0.84 and 0.74 for control and burned stands, respectively). The simulated and observed decadal to century-scale dynamics of soil temperature and carbon dynamics, which are represented by mean monthly values of these variables during the growing season, were correlated among stands ( R = 0.93 and 0.71 for soil temperature at 20- and 10-cm depths, R = 0.95 and 0.91 for soil respiration and soil carbon, respectively). Sensitivity analyses indicate that along with differences in fire and climate history a number of other factors influence the response of carbon dynamics to fire disturbance. These factors include nitrogen fixation, the growth of moss, changes in the depth of the organic layer, soil drainage, and fire severity.

Alaska↗
USGS2002

Valley floor climate observations from the McMurdo dry valleys, Antarctica, 1986-2000

Climate observations from the McMurdo dry valleys, East Antarctica are presented from a network of seven valley floor automatic meteorological stations during the period 1986 to 2000. Mean annual temperatures ranged from -14.8??C to -30.0??C, depending on the site and period of measurement. Mean annual relative humidity is generally highest near the coast. Mean annual wind speed increases with proximity to the polar plateau. Site-to-site variation in mean annual solar flux and PAR is due to exposure of each station and changes over time are likely related to changes in cloudiness. During the nonsummer months, strong katabatic winds are frequent at some sites and infrequent at others, creating large variation in mean annual temperature owing to the warming effect of the winds. Katabatic wind exposure appears to be controlled to a large degree by the presence of colder air in the region that collects at low points and keeps the warm less dense katabatic flow from the ground. The strong influence of katabatic winds makes prediction of relative mean annual temperature based on geographical position (elevation and distance from the coast) alone, not possible. During the summer months, onshore winds dominate and warm as they progress through the valleys creating a strong linear relationship (r2 = 0.992) of increasing potential temperature with distance from the coast (0.09??C km-1). In contrast to mean annual temperature, summer temperature lends itself quite well to model predictions, and is used to construct a statistical model for predicting summer dry valley temperatures at unmonitored sites. Copyright 2002 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Thermal infrared spectroscopy of experimentally shocked anorthosite and pyroxenite: Implications for remote sensing of Mars

The feldspar and pyroxene mineralogies on Mars revealed by the Thermal Emission Spectrometer (TES) on Mars Global Surveyor likely record a variety of shock effects, as suggested by petrologic analyses of the Martian meteorites and the abundance of impact craters on the planet's surface. To study the effects of shock pressures on thermal infrared spectra of these minerals, we performed shock recovery experiments on orthopyroxenite and anorthosite samples from the Stillwater Complex (Montana) over peak pressures from 17 to 63 GPa. We acquired emissivity and hemispherical reflectance spectra (350–1400 cm −1 ; ∼7–29 μm) of both coherent chips and fine-grained powders of shocked and unshocked samples. These spectra are more directly comparable to remotely sensed data of Mars (e.g., TES) than previously acquired absorption or transmission spectra of shocked minerals. The spectra of experimentally shocked feldspar show systematic changes with increasing pressure due to depolymerization of the silica tetrahedra. For the spectra of chips, this includes the disappearance of small bands in the 500–650 cm −1 region and a strong band at 1115 cm −1 , and changes in positions of a strong band near 940 cm −1 and the Christiansen feature near 1250 cm −1 . Spectra of the shocked powders show the gradual disappearance of a transparency feature near 830 cm −1 . Fewer changes are observed in the pyroxene spectra at pressures as high as 63 GPa. Spectra of experimentally shocked minerals will help identify more precisely the mineralogy of rocks and soils not only from TES but also from Mars instruments such as miniTES and THEMIS.

Journal of Geophysical Research E: Planets↗
USGS2002

Geology and origin of Europa's "Mitten" feature (Murias Chaos)

The “Mitten” (provisionally named Murias Chaos by the International Astronomical Union) is a region of elevated chaos-like terrain in the leading hemisphere of Europa. Its origin had been explained under the currently debated theories of melting through a thin lithosphere or convection within a thick one. Galileo observations reveal several characteristics that suggest that the Mitten is distinct from typical chaos terrain and point to a different formational process. Photoclinometric elevation estimates suggest that the Mitten is slightly elevated with respect to the surrounding terrain; geologic relations indicate that it must have raised significantly from the plains in its past, resembling disrupted domes on Europa's trailing hemisphere. Moreover, the Mitten material appears to have extruded onto the plains and flowed for tens of kilometers. The area subsequently subsided as a result of isostatic adjustment, viscous relaxation, and/or plains loading. Using plate flexure models, we estimated the elastic lithosphere in the area to be several kilometers thick. We propose that the Mitten originated by the ascent and extrusion of a large thermal diapir. Thermal-mechanical modeling shows that a Mitten-sized plume would remain sufficiently warm and buoyant to pierce through the crust and flow unconfined on the surface. Such a diapir probably had an initial radius between 5 and 8 km and an initial depth of 20–40 km, consistent with a thick-lithosphere model. In this scenario the Mitten appears to represent the surface expression of the rare ascent of a large diapir, in contrast to lenticulae and chaos terrain, which may form by isolated and clustered small diapirs, respectively.

Journal of Geophysical Research E: Planets↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress-induced, time-dependent fracture closure at hydrothermal conditions

Time-dependent closure of fractures in quartz was measured in situ at 22–530°C temperature and 0.1-150 MPa water pressure. Unlike previous crack healing and rock permeability studies, in this study, fracture aperture is monitored directly and continuously using a windowed pressure vessel, a long-working-distance microscope, and reflected-light interferometry. Thus the fracture volume and geometry can be measured as a function of time, temperature, and water pressure. Relatively uniform closure occurs rapidly at temperatures and pressures where quartz becomes significantly soluble in water. During closure the aperture is reduced by as much as 80% in a few hours. We infer that this closure results from the dissolution of small particles or asperities that prop the fracture open. The driving force for closure via dissolution of the prop is the sum of three chemical potential terms: (1) the dissolution potential, proportional to the logarithm of the degree of undersaturation of the solution; (2) the coarsening potential, proportional to the radius of curvature of the prop; and (3) the pressure solution potential, proportional to the effective normal stress at the contact between propping particles and the fracture wall. Our observations suggest that closure is controlled by a pressure solution-like process. The aperture of dilatant fractures and microcracks in the Earth that are similar to those in our experiments, such as ones generated from thermal stressing or brittle failure during earthquake rupture and slip, will decrease rapidly with time, especially if the macroscopic stress is nonhydrostatic.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

A fault constitutive relation accounting for thermal pressurization of pore fluid

The heat generated in a slip zone during an earthquake can raise fluid pressure and thereby reduce frictional resistance to slip. The amount of fluid pressure rise depends on the associated fluid flow. The heat generated at a given time produces fluid pressure that decreases inversely with the square root of hydraulic diffusivity times the elapsed time. If the slip velocity function is crack-like, there is a prompt fluid pressure rise at the onset of slip, followed by a slower increase. The stress drop associated with the prompt fluid pressure rise increases with rupture propagation distance. The threshold propagation distance at which thermally induced stress drop starts to dominate over frictionally induced stress drop is proportional to hydraulic diffusivity. If hydraulic diffusivity is 0.02 m 2 /s, estimated from borehole samples of fault zone material, the threshold propagation distance is 300 m. The stress wave in an earthquake will induce an unknown amount of dilatancy and will increase hydraulic diffusivity, both of which will lessen the fluid pressure effect. Nevertheless, if hydraulic diffusivity is no more than two orders of magnitude larger than the laboratory value, then stress drop is complete in large earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications

Using P and S arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations, this paper presents a three‐dimensional crustal and upper mantle velocity structure of southwestern China (21°–34°N, 97°–105°E). Southwestern China lies in the transition zone between the uplifted Tibetan plateau to the west and the Yangtze continental platform to the east. In the upper crust a positive velocity anomaly exists in the Sichuan Basin, whereas a large‐scale negative velocity anomaly exists in the western Sichuan Plateau, consistent with the upper crustal structure under the southern Tibetan plateau. The boundary between these two anomaly zones is the Longmen Shan Fault. The negative velocity anomalies at 50‐km depth in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle. The Red River Fault is the boundary between the positive and negative velocity anomalies at 50‐km depth. The overall features of the crustal and the upper mantle structures in southwestern China are a low average velocity, large crustal thickness variations, the existence of a high‐conductivity layer in the crust or/and upper mantle, and a high heat flow value. All these features are closely related to the collision between the Indian and the Asian plates.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Basin-centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350??C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes. Published in 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗