Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research – Earth Surface”

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
USGS2016

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

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

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The competition between thermal contraction and differentiation in the stress history of the Moon

The scarcity of both extension and compression features on the Moon strongly constrains the history of the lunar radius—to variations of less than ±1 km over the past 3.8 Gyr. This limit has traditionally been interpreted as requiring a delicate balance between thermal contraction of the near‐surface and expansion of a substantial cold interior region. Recent theories of lunar origin (e.g., giant impact), in contrast, favor a “hot” initial state. We propose that a reconciliation may be possible by taking account of the volume change Δ V / V | d due to differentiation. We calculate STP densities based on simplified normative mineralogies for a suite of estimates of the bulk lunar composition, of primary lunar basalt, and of the residuum left when the maximum amount of the latter is extracted from the former. Typically Δ V / V | d ≃ 2 to 5%—an expansion equivalent to heating by ∼10 3 K. Provided the timing of differentiation is correct, one might offset the cooling of a magma ocean as much as 630 km deep by differentiation of the remainder of the Moon (which need not start much below the solidus temperature). A large but not impossible amount of gabbroic melt production is implied: ∼100 times the volume of mare basalts known to have been extruded. We do not address the detailed genetic relationship of this melt to the basalts observed on the lunar surface but point out that it need not have reached the surface directly or even have entered the crust in order for the expansion to have occurred. To assess the timing of melt formation, we investigate a simple conductive lunar thermal model which takes account of both Δ V / V | d and thermal contraction. Our initial state is characterized by a central temperature T c and a depth Z 0 above which the material (derived from the magma ocean) is already at the solidus and is not suceptible to volume changes upon further differentiation. We find a range of models satisfying the limits on radius increase and decrease. The hottest has T c = 1210 K, Z 0 = 400 km; without Δ V / V | d , we would need a larger or colder (or both) core, e.g., T c ≲ 700 K for Z 0 = 200–400 km, in agreement with previous investigators. Our modeling thus lends credence to the idea that the Moon could have been initially ≳50% molten (with the remainder relatively close to the solidus) and yet experienced little volume change over the last 3.8 Gyr.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗
USGS2004

Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults

We argue that key features of thrust earthquake triggering, inhibition, and clustering can be explained by Coulomb stress changes, which we illustrate by a suite of representative models and by detailed examples. Whereas slip on surface-cutting thrust faults drops the stress in most of the adjacent crust, slip on blind thrust faults increases the stress on some nearby zones, particularly above the source fault. Blind thrusts can thus trigger slip on secondary faults at shallow depth and typically produce broadly distributed aftershocks. Short thrust ruptures are particularly efficient at triggering earthquakes of similar size on adjacent thrust faults. We calculate that during a progressive thrust sequence in central California the 1983 M w = 6.7 Coalinga earthquake brought the subsequent 1983 M w = 6.0 Nuñez and 1985 M w = 6.0 Kettleman Hills ruptures 10 bars and 1 bar closer to Coulomb failure. The idealized stress change calculations also reconcile the distribution of seismicity accompanying large subduction events, in agreement with findings of prior investigations. Subduction zone ruptures are calculated to promote normal faulting events in the outer rise and to promote thrust-faulting events on the periphery of the seismic rupture and its downdip extension. These features are evident in aftershocks of the 1957 M w = 9.1 Aleutian and other large subduction earthquakes. We further examine stress changes on the rupture surface imparted by the 1960 M w = 9.5 and 1995 M w = 8.1 Chile earthquakes, for which detailed slip models are available. Calculated Coulomb stress increases of 2–20 bars correspond closely to sites of aftershocks and postseismic slip, whereas aftershocks are absent where the stress drops by more than 10 bars. We also argue that slip on major strike-slip systems modulates the stress acting on nearby thrust and strike-slip faults. We calculate that the 1857 M w = 7.9 Fort Tejon earthquake on the San Andreas fault and subsequent interseismic slip brought the Coalinga fault ∼1 bar closer to failure but inhibited failure elsewhere on the Coast Ranges thrust faults. The 1857 earthquake also promoted failure on the White Wolf reverse fault by 8 bars, which ruptured in the 1952 M w = 7.3 Kern County shock but inhibited slip on the left-lateral Garlock fault, which has not ruptured since 1857. We thus contend that stress transfer exerts a control on the seismicity of thrust faults across a broad spectrum of spatial and temporal scales.

California↗
USGS2002

Imaging the mantle beneath Iceland using integrated seismological techniques

Using a combination of body wave and surface wave data sets to reveal the mantle plume and plume head, this study presents a tomographic image of the mantle structure beneath Iceland to 400 km depth. Data comes primarily from the PASSCAL-HOTSPOT deployment of 30 broadband instruments over a period of 2 years, and is supplemented by data from the SIL and ICEMELT networks. Three sets of relative teleseismic body wave arrival times are generated through cross correlation: S and SKS arrivals at 0.03–0.1 Hz, and P and PKIKP arrivals at 0.03–0.1 and 0.8–2.0 Hz. Prior to inversion the crustal portion of the travel time anomalies is removed using the crustal model ICECRTb. This step has a significant effect on the mantle velocity variations imaged down to a depth of ∼250 km. Inversion of relative arrival times only provides information on lateral velocity variations. Surface waves are therefore used to provide absolute velocity information for the uppermost mantle beneath Iceland. The average wave number for the Love wave fundamental mode at 0.020 and 0.024 Hz is measured and used to invert for the average S velocity. Combination of the body wave and surface wave information reveals a predominantly horizontal low-velocity anomaly extending from the Moho down to ∼250 km depth, interpreted as a plume head. Below the plume head a near-cylindrical low-velocity anomaly with a radius of ∼100 km and peak V P and V S anomalies of −2% and −4%, respectively, extends down to the maximum depth of resolution at 400 km. Within the plume head, in the uppermost mantle above the core of the plume, there is a relatively high velocity with a maximum V P and V S anomaly of +2%. This high-velocity anomaly may be the result of the extreme degree of melt extraction necessary to generate the thick (46 km) crust in central Iceland. Comparison of the plume volumetric flux implied by our images, the crustal generation rate, and the degree of melting suggested by rare earth element inversions, suggests that (1) mantle material must be flowing horizontally away from the plume core faster than the overlying lithosphere and (2) the bulk of the plume material does not participate in melting beneath Iceland.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Three-dimensional velocity structure of Siletzia and other accreted terranes in the Cascadia forearc of Washington

Eocene mafic crust with high seismic velocities underlies much of the Oregon and Washington forearc and acts as a backstop for accretion of marine sedimentary rocks from the obliquely subducting Juan de Fuca slab. Arc-parallel migration of relatively strong blocks of this terrane, known as Siletzia, focuses upper crustal deformation along block boundaries, which are potential sources of earthquakes. In a three-dimensional velocity model of coastal Washington, we have combined surface geology, well data, and travel times from earthquakes and controlled source seismic experiments to resolve the major boundaries of the Siletz terrane with the adjacent accreted sedimentary prism and volcanic arc. In southern Washington and northern Oregon the Siletz terrane appears to be a thick block (???20 km) that extends west of the coastline and makes a high-angle contact with the offshore accreted sedimentary prism. On its east flank the high-velocity Siletz terrane boundary coincides with an en echelon zone of seismicity in the arc. In northern Washington the western edge of Siletzia makes a lower-angled, fault-bound contact with the accretionary prism. In addition, alternating, east-west trending uplifts and downwarps of the Siletz terrane centered on the antiformal Olympic Mountains may reflect focusing of north-south compression in the northern part of the Siletz terrane. This compressional strain may result from northward transport and clockwise rotation of the Siletz terrane into the relatively fixed Canadian Coast Mountains restraining bend along the coast.

Oregon, Washington↗
USGS2023

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗
USGS2012

Limitations of quantitative analysis of deep crustal seismic reflection data: Examples from GLIMPCE

Amplitude preservation in seismic reflection data can be obtained by a relative true amplitude (RTA) processing technique in which the relative strength of reflection amplitudes is preserved vertically as well as horizontally, after compensating for amplitude distortion by near-surface effects and propagation effects. Quantitative analysis of relative true amplitudes of the Great Lakes International Multidisciplinary Program on Crustal Evolution seismic data is hampered by large uncertainties in estimates of the water bottom reflection coefficient and the vertical amplitude correction and by inadequate noise suppression. Processing techniques such as deconvolution, F-K filtering, and migration significantly change the overall shape of amplitude curves and hence calculation of reflection coefficients and average reflectance. Thus lithological interpretation of deep crustal seismic data based on the absolute value of estimated reflection strength alone is meaningless. The relative strength of individual events, however, is preserved on curves generated at different stages in the processing. We suggest that qualitative comparisons of relative strength, if used carefully, provide a meaningful measure of variations in reflectivity. Simple theoretical models indicate that peg-leg multiples rather than water bottom multiples are the most severe source of noise contamination. These multiples are extremely difficult to remove when the water bottom reflection coefficient is large (>0.6), a condition that exists beneath parts of Lake Superior and most of Lake Huron.

Great Lakes↗
USGS2012

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗
USGS2012

Three-dimensional P velocity structure of the crust below Newberry Volcano, Oregon

Utilizing teleseismic P residuals, we have detected a column of high P velocity material extending from within 10 km of the surface below the summit of Newberry Volcano, Oregon, to midcrustal depths near 25 km. We interpret this column to be the expression of a swarm of predominantly subsolidus gabbroic sills and dikes which were intruded as the volcano was built. The high P velocities observed below the volcano severely limit the size of magma chambers which could presently exist in the crust below Newberry Volcano. Those possible include a few percent of partial melt distributed through large volumes of a mafic intrusion zone in the midcrust; a few, smaller, higher melt fraction zones in the midcrust with dimensions less than 6 km and whose aggregate volume is only a few percent of enclosing volumes of 200 km 3 ; small magma bodies with dimensions of a few kilometers located within the upper 10 km of the crust; or a mafic, crystal-rich magma of arbitrary dimensions located in the upper few km. The low P velocities detected in the upper 4 km beneath the center of the summit caldera may be partially caused by a magma chamber in the second of these catagories.

Journal of Geophysical Research Solid Earth↗
USGS2009

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

The Augustine magmatic system as revealed by seismic tomography and relocated earthquake hypocenters from 1994 through 2009

We incorporate 14 years of earthquake data from the Alaska Volcano Observatory with data from a 1975 controlled‐source seismic experiment to obtain the three‐dimensional P and S wave velocity structure and the first high‐precision earthquake locations at Augustine Volcano to be calculated in a fully three‐dimensional velocity model. Velocity tomography shows two main features beneath Augustine: a narrow, high‐velocity column beneath the summit, extending from ∼2 km depth to the surface, and elevated velocities on the south flank. Our relocation results allow a thorough analysis of the spatio‐temoral patterns of seismicity and the relationship to the magmatic and eruptive activity. Background seismicity is centered beneath the summit at an average depth of 0.6 km above sea level. In the weeks leading to the January 2006 eruption of Augustine, seismicity focused on a NW‐SE line along the trend of an inflating dike. A series of drumbeat earthquakes occurred in the early weeks of the eruption, indicating further magma transport through the same dike system. During the six months following the onset of the eruption, the otherwise quiescent region 1 to 5 km below sea level centered beneath the summit became seismically active with two groups of earthquakes, differentiated by frequency content. The deep longer‐period earthquakes occurred during the eruption and are interpreted as resulting from the movement of magma toward the summit, and the post‐eruptive shorter‐period earthquakes may be due to the relaxation of an emptied magma tube. The seismicity subsequently returned to its normal background rates and patterns.

Alaska↗
USGS1999

Overview of the Mars Pathfinder Mission: Launch through landing, surface operations, data sets, and science results

Mars Pathfinder successfully landed at Ares Vallis on July 4, 1997, deployed and navigated a small rover about 100 m clockwise around the lander, and collected data from three science instruments and ten technology experiments. The mission operated for three months and returned 2.3 Gbits of data, including over 16,500 lander and 550 rover images, 16 chemical analyses of rocks and soil, and 8.5 million individual temperature, pressure and wind measurements. Path‐finder is the best known location on Mars, having been clearly identified with respect to other features on the surface by correlating five prominent horizon features and two small craters in lander images with those in high‐resolution orbiter images and in inertial space from two‐way ranging and Doppler tracking. Tracking of the lander has fixed the spin pole of Mars, determined the precession rate since Viking 20 years ago, and indicates a polar moment of inertia, which constrains a central metallic core to be between 1300 and ∼2000 km in radius. Dark rocks appear to be high in silica and geochemically similar to anorogenic andesites; lighter rocks are richer in sulfur and lower in silica, consistent with being coated with various amounts of dust. Rover and lander images show rocks with a variety of morphologies, fabrics and textures, suggesting a variety of rock types are present. Rounded pebbles and cobbles on the surface as well as rounded bumps and pits on some rocks indicate these rocks may be conglomerates (although other explanations are also possible), which almost definitely require liquid water to form and a warmer and wetter past. Air‐borne dust is composed of composite silicate particles with a small fraction of a highly magnetic mineral, interpreted to be most likely maghemite; explanations suggest iron was dissolved from crustal materials during an active hydrologic cycle with maghemite freeze dried onto silicate dust grains. Remote sensing data at a scale of a kilometer or greater and an Earth analog correctly predicted a rocky plain safe for landing and roving with a variety of rocks deposited by catstrophic floods, which are relatively dust free. The surface appears to have changed little since it formed billions of years ago, with the exception that eolian activity may have deflated the surface by ∼3–7 cm, sculpted wind tails, collected sand into dunes, and eroded ventifacts (fluted and grooved rocks). Pathfinder found a dusty lower atmosphere, early morning water ice clouds, and morning near‐surface air temperatures that changed abruptly with time and height. Small scale vortices, interpreted to be dust devils, were observed repeatedly in the afternoon by the meteorology instruments and have been imaged.

Journal of Geophysical Research E: Planets↗
USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗
USGS2012

The sea slope problem revisited

Discrepancies in sea surface topography based on comparisons between the results of steric leveling and repeated geodetic levelings have identified what is known as the ‘sea slope problem.’ This problem is actually twofold: (1) the sea surface relief based on steric leveling differs significantly from that based on geodetic leveling along several generally north-south coasts, and (2) successively propagated levelings between several widely separated tide stations indicate that the stationary sea slope seemingly has been changing with time, whereas differenced sea level means between these stations indicate that the sea surface relief has remained virtually invariant during the same intervals. Reexamination of the three reported discrepancies between geodetically and sterically determined sea slopes indicates that the Australian example is based on leveling of a quality inappropriate to the comparison. The discrepancy developed along the Atlantic coast of the United States is limited to the reach between Portsmouth–Hampton Roads, Virginia, and Key West, Florida, where the accuracy of steric leveling may be especially vulnerable owing to the dynamic effects of the Gulf Stream. Reconsideration of the example along the Pacific coast of the United States indicates that the various discrepancies are due chiefly to intrasurvey movement and resultant distortion of geodetically defined height differences between tide stations. Agreement between the results of steric and geodetic leveling along tectonically inactive north-south coasts devoid of strong boundary currents is generally good. This observation supports the conclusions (1) that any directionally dependent systematic error in geodetic leveling is measurably insignificant and (2) that where allowance is made for the possible effects of major boundary currents or intrasurvey movement during levelings between tide stations, the sea slope problem tends to vanish.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Sedimentary deposits in the northern lowland plains, Mars

The lowland plains on Mars have surfaces marked by large polygonal fracture patterns. It was recently proposed that the fracture patterns were developed on sedimentary deposits from outflow channels. We support this hypothesis because of the following observations. (1) Polygonal fracture patterns tend to occur in low areas on Mars that apparently received influx of sediments; the area of northern plains, in which the fractures occur, also forms a regional low. (2) The fracture patterns are concentrated near reentrants of the plains projecting into the southern highlands. Not only are these reentrants associated with outflow channels, but channels and polygonally fractured deposits have similar ages. (3) Sinuous ridges at the mouths of Martian outflow channels resemble ridges in Antarctic ice streams and ice shelves; the similarity suggests that the Martian material was transported from the southern highlands toward the northern plains, as the Antarctic material is transported from the continent toward the adjacent ocean. Overall, it appears that a massive transfer of material took place midway in Martian history, either through the outflow channels or elsewhere along the northern highland scarp, and that this material probably formed the deposits now characterized by polygonal fracture patterns.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength

The great 27 February 2010 M w 8.8 earthquake off the coast of southern Chile ruptured a ∼600 km length of subduction zone. In this paper, we make two independent estimates of shear stress in the crust in the region of the Chile earthquake. First, we use a coseismic slip model constrained by geodetic observations from interferometric synthetic aperture radar (InSAR) and GPS to derive a spatially variable estimate of the change in static shear stress along the ruptured fault. Second, we use a static force balance model to constrain the crustal shear stress required to simultaneously support observed fore‐arc topography and the stress orientation indicated by the earthquake focal mechanism. This includes the derivation of a semianalytic solution for the stress field exerted by surface and Moho topography loading the crust. We find that the deviatoric stress exerted by topography is minimized in the limit when the crust is considered an incompressible elastic solid, with a Poisson ratio of 0.5, and is independent of Young's modulus. This places a strict lower bound on the critical stress state maintained by the crust supporting plastically deformed accretionary wedge topography. We estimate the coseismic shear stress change from the Maule event ranged from −6 MPa (stress increase) to 17 MPa (stress drop), with a maximum depth‐averaged crustal shear‐stress drop of 4 MPa. We separately estimate that the plate‐driving forces acting in the region, regardless of their exact mechanism, must contribute at least 27 MPa trench‐perpendicular compression and 15 MPa trench‐parallel compression. This corresponds to a depth‐averaged shear stress of at least 7 MPa. The comparable magnitude of these two independent shear stress estimates is consistent with the interpretation that the section of the megathrust fault ruptured in the Maule earthquake is weak, with the seismic cycle relieving much of the total sustained shear stress in the crust.

Maule↗