Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

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 433 records · Page 24Linked to original sources

Posteruption suspended sediment transport at Mount St. Helens: Decadal‐scale relationships with landscape adjustments and river discharges

Widespread landscape disturbance by the cataclysmic 1980 eruption of Mount St. Helens abruptly increased sediment supply in surrounding watersheds. The magnitude and duration of the redistribution of sediment deposited by the eruption as well as decades‐ to centuries‐old sediment remobilized from storage have varied chiefly with the style of disturbance. Posteruption suspended sediment transport has been greater and more persistent from zones of channel disturbance than from zones of hillslope disturbance. Despite the severe landscape disturbances caused by the eruption, relationships between discharge magnitudes and frequencies and suspended sediment transport have been remarkably consistent. Discharges smaller than mean annual flows generally have transported <5%, but locally ∼15%, of the annual suspended sediment loads, and infrequent ( p < 0.01), large floods have transported as much as 50% of the annual suspended sediment loads in a single day. However, moderate‐magnitude discharges (those greater than mean annual flows but less than 2‐year floods) have transported the greatest amounts of sediment from all disturbance zones. Such discharges have transported, on average, 60% to ∼95% of the annual suspended sediment loads, usually within cumulative periods of 1–3 weeks each year. Although small‐magnitude and large‐magnitude discharges have locally and episodically transported considerable amounts of suspended sediment, there has not been any notable change in the overall nature of the effective discharges; moderate‐magnitude flows have been the predominant discharges responsible for transporting the majority of suspended sediment during 20 years of posteruption landscape adjustment.

Washington↗

Stress field variations in the Swiss Alps and the northern Alpine foreland derived from inversion of fault plane solutions

This study is devoted to a systematic analysis of the state of stress of the central European Alps and northern Alpine foreland in Switzerland based on focal mechanisms of 138 earthquakes with magnitudes between 1 and 5. The most robust feature of the results is that the azimuth of the minimum compressive stress, S 3 , is generally well constrained for all data subsets and always lies in the NE quadrant. However, within this quadrant, the orientation of S 3 changes systematically both along the structural strike of the Alpine chain and across it. The variation in stress along the mountain belt from NE to SW involves a progressive, counterclockwise rotation of S 3 and is most clear in the foreland, where it amounts to 45°–50°. This pattern of rotation is compatible with the disturbance to the stress field expected from the indentation of the Adriatic Block into the central European Plate, possibly together with buoyancy forces arising from the strongly arcuate structure of the Moho to the immediate west of our study area. Across the Alps, the variation in azimuth of S 3 is defined by a progressive, counterclockwise rotation of about 45° from the foreland in the north across the Helvetic domain to the Penninic nappes in the south and is accompanied by a change from a slight predominance of strike-slip mechanisms in the foreland to a strong predominance of normal faulting in the high parts of the Alps. The observed rotation can be explained by the perturbation of the large-scale regional stress by a local uniaxial deviatoric tension with a magnitude similar to that of the regional differential stress and with an orientation perpendicular to the strike of the Alpine belt. The tensile nature and orientation of this stress is consistent with the “spreading” stress expected from lateral density changes due to a crustal root beneath the Alps.

European Alps↗

Crystallographic controls on the frictional behavior of dry and water-saturated sheet structure minerals

We compare the frictional strengths of 17 sheet structure mineral powders, measured under dry and water-saturated conditions, to identify the factors that cause many of them to be relatively weak. The dry coefficient of friction μ ranges upward from 0.2 for graphite, leveling off at 0.8 for margarite, clintonite, gibbsite, kaolinite, and lizardite. The values of μ (dry) correlate directly with calculated (001) interlayer bond strengths of the minerals. This correlation occurs because shear becomes localized along boundary and Riedel shears and the platy minerals in them rotate into alignment with the shear planes. For those gouges with μ (dry) < 0.8, shear occurs by breaking the interlayer bonds to form new cleavage surfaces. Where μ (dry) = 0.8, consistent with Byerlee's law, the interlayer bonds are sufficiently strong that other frictional processes dominate. The transition in dry friction mechanisms corresponds to calculated surface energies of 2–3 J/m 2 . Adding water causes μ to decrease for every mineral tested except graphite. If the minerals are separated into groups with similar crystal structures, μ (wet) increases with increasing interlayer bond strength within each group. This relationship also holds for the swelling clay montmorillonite, whose water-saturated strength is consistent with the strengths of nonswelling clays of similar crystal structure. Water in the saturated gouges forms thin, structured films between the plate surfaces. The polar water molecules are bonded to the plate surfaces in proportion to the mineral's surface energy, and μ (wet) reflects the stresses required to shear through the water films.

Journal of Geophysical Research B: Solid Earth↗

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↗

Strain accumulation across the Coast Ranges at the latitude of San Francisco, 1994-2000

A 66-monument geodetic array spanning the Coast Ranges near San Francisco has been surveyed more than eight times by GIPS between late 1993 and early 2001. The measured horizontal velocities of the monuments are well represented by uniform, right-lateral, simple shear parallel to N29°W. (The local strike of the San Andreas Fault is ∼N34°W.) The observed areal dilatation rate of 6.9 ± 10.0 nstrain yr −1 (quoted uncertainty is one standard deviation and extension is reckoned positive) is not significantly different from zero, which implies that the observed strain accumulation could be released by strike-slip faulting alone. Our results are consistent with the slip rates assigned by the Working Group on California Earthquake Probabilities [2003] to the principal faults (San Gregorio, San Andreas, Hayward-Rodgers Creek, Calaveras-Concord-Green Valley, and Greenville Faults) cutting across the GPS array. The vector sum of those slip rates is is 39.8 ± 2.6 mm yr −1 N29.8°W ± 2.8°, whereas the motion across the GPS array (breadth 120 km) inferred from the uniform strain rate approximation is 38.7 ± 1.2 mm yr −1 N29.0°W ± 0.9° right-lateral shear and 0.4 ± 0.9 mm yr −1 N61°E ± 0.9° extension. We interpret the near coincidence of these rates and the absence of significant accumulation of areal dilatation to imply that right-lateral slip on the principal faults can release the accumulating strain; major strain release on reverse faults subparallel to the San Andreas Fault within the Coast Ranges is not required.

California↗

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗

Interseismic strain and rotation rates in the northeast Mojave domain, eastern California

The northeast Mojave domain, a type locality for bookshelf faulting, is a region of east striking, left-lateral faults in the northeast corner of the Mojave block, a block otherwise dominated by ∼N40°W striking, right-lateral faults. Paleomagnetic evidence suggests that blocks within the domain have rotated clockwise about a vertical axis as much as 60° since 12.8 Ma [Schermer et al., 1996] . In 1994, and again in 2002, the U.S. Geological Survey surveyed an array of 14 geodetic monuments distributed across the northeast Mojave domain. The 2002 survey results were adjusted to remove the coseismic offsets imposed by the nearby Hector Mine earthquake (16 October 1999, M w = 7.1). The adjusted deformation across the array appears to be uniform and can be approximated by the principal strain rates ε 1 = 28.9 ± 9.1 N77.2°W ± 4.8° and ε 2 = −48.2 ± 8.9 N12.8°E ± 4.8° nstrain yr −1 ; extension reckoned positive, and quoted uncertainties are standard deviations. That strain accumulation could be released by slip on faults striking N32°W but not by bookshelf faulting on the east striking faults alone. The vertical axis rotation rate of the northeast Mojave domain as a whole relative to fixed North America is 71.0 ± 6.4 nrad yr −1 (4.07° ± 0.37° Myr −1 ) clockwise, about twice the maximum tensor shear strain rate. The observed rotation rate acting over 12.8 Myr would produce a clockwise rotation of 52.1° ± 4.7°, exclusive of possible coseismic rotations. That rotation is in rough agreement with the paleomagnetic rotation accumulated in the individual fault blocks within the northeast Mojave domain since 12.8 Ma.

California↗

Contemporary tectonic deformation of the Basin and Range province, western United States: 10 years of observation with the Global Positioning System

[1] We have estimated patterns and rates of crustal movement across 800 km of the Basin and Range at ∼39° north latitude with Global Positioning System surveys in 1992, 1996, 1998, and 2002. The total rate of motion tangent to the small circle around the Pacific‐North America pole of rotation is 10.4 ± 1.0 mm/yr, and motion normal to this small circle is 3.9 ± 0.9 mm/yr compared to the east end of our network. On the Colorado Plateau the east end of our network moves by ∼1–2 mm/yr westerly with respect to North America. Transitions in strain rates delimit six major tectonic domains within the province. These deformation zones coincide with areas of modern seismicity and are, from east to west, (1) east‐west extension in the Wasatch Fault zone, (2) low rate east‐west extension centered near the Nevada‐Utah border, (3) low rate east‐west contraction between 114.7°W and 117.9°W, (4) extension normal to and strike‐slip motion across the N10°E striking Central Nevada Seismic Zone, (5) right lateral simple shear oriented N13°W inside the Walker Lane Belt, and (6) shear plus extension near the Sierra Nevada frontal faults. Concentration of shear and dilatational deformation across the three westernmost zones suggests that the Walker Lane Belt lithosphere is rheologically weak. However, we show that linear gradients in viscosity and gravitational potential energy can also effectively concentrate deformation. In the Basin and Range, gradients in gravitational potential are spatially anticorrelated with dilatational strain rates, consistent with the presence of horizontal variations in viscosity of the lithosphere.

Arizona, California, Nevada, Utah↗

Magmatic effects of the Cobb hot spot on the Juan de Fuca Ridge

The interaction of the Juan de Fuca Ridge with the Cobb hot spot has had a considerable influence on the magmatism of the Axial Segment of the ridge, the second-order segment that overlies the hot spot. In addition to the construction of the large volcanic edifice of Axial Seamount, the Axial Segment has shallow bathymetry and a prevalence of constructional volcanic features along its 100-km length, suggesting that hot spot-derived magmas supplement and oversupply the ridge. Lavas are generally more primitive at Axial Seamount and more evolved in the Axial Segment rift zones, suggesting that fractional crystallization is enhanced with increasing distance from the hot spot because of a reduced magma supply and more rapid cooling. Although the Cobb hot spot is not an isotopically enriched plume, it produces lavas with some distinct geochemical characteristics relative to normal mid-ocean ridge basalt, such as enrichments in alkalis and highly incompatible trace elements, that can be used as tracers to identify the presence and prevalence of the hot spot influence along the ridge. These characteristics are most prominent at Axial Seamount and decline in gradients along the Axial Segment. The physical model that can best explain the geochemical observations is a scenario in which hot spot and mid-ocean ridge basalt (MORB) magmas mix to varying degrees, with the proportions controlled by the depth to the MORB source. Modeling of two-component mixing suggests that MORB is the dominant component in most Axial Segment basalts. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Noise in two-color electronic distance meter measurements revisited

Frequent, high-precision geodetic data have temporally correlated errors. Temporal correlations directly affect both the estimate of rate and its standard error; the rate of deformation is a key product from geodetic measurements made in tectonically active areas. Various models of temporally correlated errors are developed and these provide relations between the power spectral density and the data covariance matrix. These relations are applied to two-color electronic distance meter (EDM) measurements made frequently in California over the past 15-20 years. Previous analysis indicated that these data have significant random walk error. Analysis using the noise models developed here indicates that the random walk model is valid for about 30% of the data. A second 30% of the data can be better modeled with power law noise with a spectral index between 1 and 2, while another 30% of the data can be modeled with a combination of band-pass-filtered plus random walk noise. The remaining 10% of the data can be best modeled as a combination of band-pass-filtered plus power law noise. This band-pass-filtered noise is a product of an annual cycle that leaks into adjacent frequency bands. For time spans of more than 1 year these more complex noise models indicate that the precision in rate estimates is better than that inferred by just the simpler, random walk model of noise.

Journal of Geophysical Research B: Solid Earth↗

Estimating hydrodynamic roughness in a wave-dominated environment with a high-resolution acoustic Doppler profiler

Hydrodynamic roughness is a critical parameter for characterizing bottom drag in boundary layers, and it varies both spatially and temporally due to variation in grain size, bedforms, and saltating sediment. In this paper we investigate temporal variability in hydrodynamic roughness using velocity profiles in the bottom boundary layer measured with a high-resolution acoustic Doppler profiler (PCADP). The data were collected on the ebb-tidal delta off Grays Harbor, Washington, in a mean water depth of 9 m. Significant wave height ranged from 0.5 to 3 m. Bottom roughness has rarely been determined from hydrodynamic measurements under conditions such as these, where energetic waves and medium-to-fine sand produce small bedforms. Friction velocity due to current u*c and apparent bottom roughness z0a were determined from the PCADP burst mean velocity profiles using the law of the wall. Bottom roughness kB was estimated by applying the Grant-Madsen model for wave-current interaction iteratively until the model u*c converged with values determined from the data. The resulting kB values ranged over 3 orders of magnitude (10-1 to 10-4 m) and varied inversely with wave orbital diameter. This range of kB influences predicted bottom shear stress considerably, suggesting that the use of time-varying bottom roughness could significantly improve the accuracy of sediment transport models. Bedform height was estimated from kB and is consistent with both ripple heights predicted by empirical models and bedforms in sonar images collected during the experiment. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

Suspended sediment and hydrodynamics above mildly sloped long wave ripples

We investigate the spatial and temporal distribution of suspended sediment and the associated hydrodynamics over mildly sloped long wave ripples on the inner shelf. These bedforms had wavelengths of approximately 1 m and heights of approximately 5 cm, in a mean water depth of 4 m. The vertical and temporal structures of the suspended sediment concentration (SSC) are consistent with the entrainment of sediment on the offshore flank of the ripple, and rapid vertical mixing at the time of flow reversal, followed by advection onshore by the onshore fluid motion. This work confirms that the mechanism for sediment suspension above low-amplitude, long wave ripples is similar to the vortex formation process expected over steeper vortex ripples. Numerical simulations of the flow using the Dune2d model indicate that a separated rotational flow structure is generated at the flank of the ripple on the seaward side of the ripple crest, near the time of flow reversal. The simulations indicate that only one vortex is formed during each wave period, in agreement with the field observations. This asymmetry is due mainly to the presence of an offshore mean near-bed current of approximately 6 to 8 cm/s. The SSC is calculated by the model and compared to the field observations. A hydraulic bed roughness of 10 to 15 median grain diameters ( d 50 ) was used in order to match the model prediction to the observed SSC approximately 1 cm above the seabed (cab). However, the modeled SSC and turbulent kinetic energy were significantly lower than the field observations at elevations exceeding approximately 2 cab.

North Carolina↗

Mars Exploration Rover Athena Panoramic Camera (Pancam) investigation

The Panoramic Camera (Pancam) investigation is part of the Athena science payload launched to Mars in 2003 on NASA's twin Mars Exploration Rover (MER) missions. The scientific goals of the Pancam investigation are to assess the high‐resolution morphology, topography, and geologic context of each MER landing site, to obtain color images to constrain the mineralogic, photometric, and physical properties of surface materials, and to determine dust and aerosol opacity and physical properties from direct imaging of the Sun and sky. Pancam also provides mission support measurements for the rovers, including Sun‐finding for rover navigation, hazard identification and digital terrain modeling to help guide long‐term rover traverse decisions, high‐resolution imaging to help guide the selection of in situ sampling targets, and acquisition of education and public outreach products. The Pancam optical, mechanical, and electronics design were optimized to achieve these science and mission support goals. Pancam is a multispectral, stereoscopic, panoramic imaging system consisting of two digital cameras mounted on a mast 1.5 m above the Martian surface. The mast allows Pancam to image the full 360° in azimuth and ±90° in elevation. Each Pancam camera utilizes a 1024 × 1024 active imaging area frame transfer CCD detector array. The Pancam optics have an effective focal length of 43 mm and a focal ratio of f /20, yielding an instantaneous field of view of 0.27 mrad/pixel and a field of view of 16° × 16°. Each rover's two Pancam “eyes” are separated by 30 cm and have a 1° toe‐in to provide adequate stereo parallax. Each eye also includes a small eight position filter wheel to allow surface mineralogic studies, multispectral sky imaging, and direct Sun imaging in the 400–1100 nm wavelength region. Pancam was designed and calibrated to operate within specifications on Mars at temperatures from −55° to +5°C. An onboard calibration target and fiducial marks provide the capability to validate the radiometric and geometric calibration on Mars.

Journal of Geophysical Research E: Planets↗

Athena Microscopic Imager investigation

The Athena science payload on the Mars Exploration Rovers (MER) includes the Microscopic Imager (MI). The MI is a fixed‐focus camera mounted on the end of an extendable instrument arm, the Instrument Deployment Device (IDD). The MI was designed to acquire images at a spatial resolution of 30 microns/pixel over a broad spectral range (400–700 nm). The MI uses the same electronics design as the other MER cameras but has optics that yield a field of view of 31 × 31 mm across a 1024 × 1024 pixel CCD image. The MI acquires images using only solar or skylight illumination of the target surface. A contact sensor is used to place the MI slightly closer to the target surface than its best focus distance (about 66 mm), allowing concave surfaces to be imaged in good focus. Coarse focusing (∼2 mm precision) is achieved by moving the IDD away from a rock target after the contact sensor has been activated. The MI optics are protected from the Martian environment by a retractable dust cover. The dust cover includes a Kapton window that is tinted orange to restrict the spectral bandpass to 500–700 nm, allowing color information to be obtained by taking images with the dust cover open and closed. MI data will be used to place other MER instrument data in context and to aid in petrologic and geologic interpretations of rocks and soils on Mars.

Journal of Geophysical Research E: Planets↗

Athena Mars rover science investigation

Each Mars Exploration Rover carries an integrated suite of scientific instruments and tools called the Athena science payload. The primary objective of the Athena science investigation is to explore two sites on the Martian surface where water may once have been present, and to assess past environmental conditions at those sites and their suitability for life. The remote sensing portion of the payload uses a mast called the Pancam Mast Assembly (PMA) that provides pointing for two instruments: the Panoramic Camera (Pancam), and the Miniature Thermal Emission Spectrometer (Mini‐TES). Pancam provides high‐resolution, color, stereo imaging, while Mini‐TES provides spectral cubes at mid‐infrared wavelengths. For in‐situ study, a five degree‐of‐freedom arm called the Instrument Deployment Device (IDD) carries four more tools: a Microscopic Imager (MI) for close‐up imaging, an Alpha Particle X‐Ray Spectrometer (APXS) for elemental chemistry, a Mössbauer Spectrometer (MB) for the mineralogy of Fe‐bearing materials, and a Rock Abrasion Tool (RAT) for removing dusty and weathered surfaces and exposing fresh rock underneath. The payload also includes magnets that allow the instruments to study the composition of magnetic Martian materials. All of the Athena instruments have undergone extensive calibration, both individually and using a set of geologic reference materials that are being measured with all the instruments. Using a MER‐like rover and payload in a number of field settings, we have devised operations processes that will enable us to use the MER rovers to formulate and test scientific hypotheses concerning past environmental conditions and habitability at the landing sites.

Journal of Geophysical Research E: Planets↗

Luminescence dating of the Wabar meteorite craters, Saudi Arabia

Luminescence dating has been used to find the age of meteorite impact craters at Wabar (Al Hadida) in Saudi Arabia. The luminescence characteristics of the shocked material were determined. Using a variety of luminescence dating techniques applied to impactite formed by the meteorite, and to the underlying sand, the age is found to be 290 &plusmn; 38 years. A comparison is made with two possible historically recorded ages. An impact as young as this has implications for the assessment of hazards from the impact on Earth of small meteorites.

Al Hadida, Ar-Rub' Al-Khali (Empty Quarter)↗

Photomosaics of the cathodoluminescence of 60 sections of meteorites and lunar samples

Cathodoluminescence (CL) petrography provides a means of observing petrographic and compositional properties of geological samples not readily observable by other techniques. We report the low-magnification CL images of 60 sections of extraterrestrial materials. The images we report include ordinary chondrites (including type 3 ordinary chondrites and gas-rich regolith breccias), enstatite chondrites, CO chondrites and a CM chondrite, eucrites and a howardite, lunar highland regolith breccias, and lunar soils. The CL images show how primitive materials respond to parent body metamorphism, how the metamorphic history of EL chondrites differs from that of EH chondrites, how dark matrix and light clasts of regolith breccias relate to each other, how metamorphism affects eucrites, the texture of lunar regolith breccias and the distribution of crystallized lunar spherules ("lunar chondrules"), and how regolith working affects the mineral properties of lunar soils. More particularly, we argue that such images are a rich source of new information on the nature and history of these materials and that our efforts to date are a small fraction of what can be done. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗