Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 937 records · Page 52Linked to original sources
USGS2004

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↗
USGS2004

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↗
USGS2004

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↗
USGS2004

Precise relocation of earthquakes following the 15 June 1991 eruption of Mount Pinatubo (Philippines)

The 15 June 1991 climactic eruption of Mount Pinatubo (Philippines) was followed by intense seismicity that remained at a high level for several months. We located 10,839 events recorded between 1 July and mid-December 1991. In contrast to the preeruptive seismicity which was focused in two groups below the summit area, posteruptive events were widely distributed below and around the volcano. The classification of the events indicates the presence of several large multiplets, and the application of relative relocation techniques to the similar events by calculating high-precision delays between traces outlines a number of clear seismogenic structures. We used different methods to confirm the validity of our results; these tests indicate that reliable features can be detected with a small monitoring network. While the main cluster of activity can be attributed to an intrusive process starting from below the 15 June crater, the volcanic origin of the seismic activity in the other areas is more difficult to establish. Away from the summit, relocations define streaks or planes which are oriented predominantly southwest-northeast, with in several cases the presence of northwest-southeast conjugate structures. Most of the composite focal mechanisms that we could determine indicate predominantly strike-slip, right-lateral faulting. Our results indicate that most of the seismicity that occurred after the 15 June eruption is related to the east-west regional compressional stress field related to the subduction. We suggest that the regional stress field induces seismicity along new or preexisting faults in the medium surrounding the volcano where the stress field was locally disturbed by the volcanic eruption.

Mount Pinatubo↗
USGS2004

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↗
USGS2004

Effects of the fungicides mancozeb and chlorothalonil on fluxes of CO2, N2O, and CH4 in a fertilized Colorado grassland soil

[1] Management of agricultural soil plays an important role in present and future atmospheric concentrations of the greenhouse gases carbon dioxide (CO 2 ), nitrous oxide (N 2 O), and methane (CH 4 ). Pesticides are used as management tools in crop production, but little is known about their effects on soil‐atmosphere exchange of CO 2 , N 2 O, and CH 4 . Field studies described in this paper determined the effect of two commonly used fungicides, mancozeb and chlorothalonil, on trace gas exchange. Separate experimental plots, 1 m 2 , were established in nitrogen fertilized no‐tilled native grassland and tilled soils with and without fungicide application. Two studies were conducted. The first study was initiated in June 1999 and lasted for 1 year with monthly flux measurements from tilled and no‐till soils. The second study commenced in August 2001 with twelve weekly measurements from tilled soils only. From both studies mancozeb suppressed emissions of CO 2 and N 2 O in the tilled soil by an average of 28% and 47%, respectively. This suppression corresponded with efficacy periods of 14–29 and 56–77 days, respectively. From the no‐till soils mancozeb decreased CO 2 and N 2 O emissions by 33% and 80% for periods of 29 and 94 days, respectively. Mancozeb inhibited CH 4 consumption in the first study by 46% and 71% in the tilled and no‐till soil for periods of 8 and 29 days, respectively, but had no effect in the second study. From both studies chlorothalonil initially suppressed CO 2 and N 2 O emissions and enhanced CH 4 uptake in the tilled soil by an average of 37%, 40%, and 115%, respectively. These effects corresponded with efficacy periods of 14–29, 21–56, and 1–14 days, respectively. In the no‐till soil chlorothalonil inhibited CO 2 and N 2 O emissions and enhanced CH 4 uptake by 29%, 48%, and 86% for periods of 29, 56, and 56 days, respectively. Following the initial period of suppression, chlorothalonil subsequently enhanced N 2 O emissions in the tilled soil by an average of 51% and in the no‐till soil by 81% before returning to near background levels. The beginning of increased N 2 O emissions from the chlorothalonil‐amended plots corresponded with a maximum soil concentration of the chlorothalonil degradate, 4‐hydroxy‐2, 5, 6‐trichloroisophthalonitrile. The site specific global warming potential (GWP) resulting from the fluxes of CO 2 , N 2 O, and CH 4 from all soils was determined to decrease by an average 26% and 21% as a result of a single application of mancozeb or chlorothalonil, respectively. The decrease in CO 2 emissions in the fungicide‐amended plots potentially could result in the conservation of as much as 1200 and 2400 kg C ha −1 yr −1 organic carbon in the tilled and no‐till plots, respectively. Therefore it is feasible that application of certain fungicides to agricultural soil might lead to enhanced soil carbon sequestration and thus have additional positive effects on atmospheric CO 2 concentrations.

Journal of Geophysical Research D: Atmospheres↗
USGS2004

Effects of the herbicides prosulfuron and metolachlor on fluxes of CO2,N2O, and CH4 in a fertilized Colorado grassland soil

The effect that pesticides have on trace gas production and consumption in agricultural soils is often overlooked. Independent field and laboratory experiments were used to measure the effects that the commonly used herbicides prosulfuron and metolachlor have on trace gas fluxes (CO 2 , N 2 O, and CH 4 ) from fertilized soil of the Colorado shortgrass steppe. Separate sample plots (1 m 2 ) on tilled and no-till soil at the sites included the following treatments: 1) a control without fertilizer or herbicide, 2) a fertilized (NH 4 NO 3 equivalent to 244 kg ha −1 ) control without herbicide, 3) and fertilized plots amended with an herbicide (prosulfuron equivalent to 0.46 kg ha −1 57% by weight active ingredient or metolachlor equivalent to 5.7 L ha −1 , 82.4% by weight active ingredient). During an initial study of one year duration, measurement of gas exchange revealed that prosulfuron-amendment stimulated N 2 O emissions and CH 4 consumption by as much as 1600% and 1300% during a single measurement, respectively. During a second set of flux measurements beginning in August 2001, more frequent weekly measurements were made during a twelve week period. From this second study an increased N 2 O efflux and CH 4 uptake occurred after a 7-week lag period that persisted for about 5 weeks. These changes in gas flux amounted to an overall increase of 41% and 30% for N 2 O emission and CH 4 consumption, respectively. The co-occurrence of stimulated N 2 O and CH 4 fluxes suggests a similar cause that is related to prosulfuron degradation. Evidence suggested that prosulfuron degradation stimulated microbial activity responsible for trace gas flux. Ultimately, prosulfuron-amendment led to an ∼50% reduction in the global warming potential from N 2 O and CH 4 fluxes at this field site, which is equivalent to a reduction of the global warming potential of 0.18 mols CO 2 m −2 d −1 from these gases. Metolachlor application did not significantly affect the trace gas fluxes measured. These results demonstrate the potential impact that pesticides have on trace gas fluxes from agricultural soils, which could mean that the effects of other agricultural practices have been over or under estimated.

Colorado↗
USGS2003

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↗
USGS2003

Selection of the Mars Exploration Rover landing sites

The selection of Meridiani Planum and Gusev crater as the Mars Exploration Rover landing sites took over 2 years, involved broad participation of the science community via four open workshops, and narrowed an initial ∼155 potential sites (80–300 × 30 km) to four finalists based on science and safety. Engineering constraints important to the selection included (1) latitude (10°N–15°S) for maximum solar power, (2) elevation (less than −1.3 km) for sufficient atmosphere to slow the lander, (3) low horizontal winds, shear, and turbulence in the last few kilometers to minimize horizontal velocity, (4) low 10‐m‐scale slopes to reduce airbag spin‐up and bounce, (5) moderate rock abundance to reduce abrasion or strokeout of the airbags, and (6) a radar‐reflective, load‐bearing, and trafficable surface safe for landing and roving that is not dominated by fine‐grained dust. The evaluation of sites utilized existing as well as targeted orbital information acquired from the Mars Global Surveyor and Mars Odyssey. Three of the final four landing sites show strong evidence for surface processes involving water and appear capable of addressing the science objectives of the missions, which are to determine the aqueous, climatic, and geologic history of sites on Mars where conditions may have been favorable to the preservation of evidence of possible prebiotic or biotic processes. The evaluation of science criteria placed Meridiani and Gusev as the highest‐priority sites. The evaluation of the three most critical safety criteria (10‐m‐scale slopes, rocks, and winds) and landing simulation results indicated that Meridiani and Elysium Planitia are the safest sites, followed by Gusev and Isidis Planitia.

Journal of Geophysical Research E: Planets↗
USGS2003

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↗
USGS2003

Meter-scale slopes of candidate MER landing sites from point photoclinometry

Photoclinometry was used to analyze the small-scale roughness of areas that fall within the proposed Mars Exploration Rover (MER) 2003 landing ellipses. The landing ellipses presented in this study were those in Athabasca Valles, Elysium Planitia, Eos Chasma, Gusev Crater, Isidis Planitia, Melas Chasma, and Meridiani Planum. We were able to constrain surface slopes on length scales comparable to the image resolution (1.5 to 12 m/pixel). The MER 2003 mission has various engineering constraints that each candidate landing ellipse must satisfy. These constraints indicate that the statistical slope values at 5 m baselines are an important criterion. We used our technique to constrain maximum surface slopes across large swaths of each image, and built up slope statistics for the images in each landing ellipse. We are confident that all MER 2003 landing site ellipses in this study, with the exception of the Melas Chasma ellipse, are within the small-scale roughness constraints. Our results have provided input into the landing hazard assessment process. In addition to evaluating the safety of the landing sites, our mapping of small-scale roughnesses can also be used to better define and map morphologic units. The morphology of a surface is characterized by the slope distribution and magnitude of slopes. In looking at how slopes are distributed, we can better define landforms and determine the boundaries of morphologic units.

Journal of Geophysical Research E: Planets↗
USGS2003

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↗
USGS2003

High-resolution topomapping of candidate MER landing sites with Mars Orbiter Camera narrow-angle images

We analyzed narrow‐angle Mars Orbiter Camera (MOC‐NA) images to produce high‐resolution digital elevation models (DEMs) in order to provide topographic and slope information needed to assess the safety of candidate landing sites for the Mars Exploration Rovers (MER) and to assess the accuracy of our results by a variety of tests. The mapping techniques developed also support geoscientific studies and can be used with all present and planned Mars‐orbiting scanner cameras. Photogrammetric analysis of MOC stereopairs yields DEMs with 3‐pixel (typically 10 m) horizontal resolution, vertical precision consistent with ∼0.22 pixel matching errors (typically a few meters), and slope errors of 1–3°. These DEMs are controlled to the Mars Orbiter Laser Altimeter (MOLA) global data set and consistent with it at the limits of resolution. Photoclinometry yields DEMs with single‐pixel (typically ∼3 m) horizontal resolution and submeter vertical precision. Where the surface albedo is uniform, the dominant error is 10–20% relative uncertainty in the amplitude of topography and slopes after “calibrating” photoclinometry against a stereo DEM to account for the influence of atmospheric haze. We mapped portions of seven candidate MER sites and the Mars Pathfinder site. Safety of the final four sites (Elysium, Gusev, Isidis, and Meridiani) was assessed by mission engineers by simulating landings on our DEMs of “hazard units” mapped in the sites, with results weighted by the probability of landing on those units; summary slope statistics show that most hazard units are smooth, with only small areas of etched terrain in Gusev crater posing a slope hazard.

Journal of Geophysical Research E: Planets↗
USGS2004

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 ± 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)↗
USGS2004

Frequency and initiation of debris flows in Grand Canyon, Arizona

Debris flows from 740 tributaries transport sediment into the Colorado River in Grand Canyon, Arizona, creating rapids that control its longitudinal profile. Debris flows mostly occur when runoff triggers failures in colluvium by a process termed “the fire hose effect.” Debris flows originate from a limited number of geologic strata, almost exclusively shales or other clay‐rich, fine‐grained formations. Observations from 1984 through 2003 provide a 20 year record of all debris flows that reached the Colorado River in Grand Canyon, and repeat photography provides a 100 year record of debris flows from 147 tributaries. Observed frequencies are 5.1 events/year from 1984 to 2003, and historic frequencies are 5.0 events/year from 1890 to 1983. Logistic regression is used to model historic frequencies based on drainage basin parameters observed to control debris flow initiation and transport. From 5 to 7 of the 16 parameters evaluated are statistically significant, including drainage area, basin relief, and the height of and gradient below debris flow source areas, variables which reflect transport distance and potential energy. The aspect of the river channel, which at least partially reflects storm movement within the canyon, is also significant. Model results are used to calculate the probability of debris flow occurrence at the river over a century for all 740 tributaries. Owing to the variability of underlying geomorphic controls, the distribution of this probability is not uniform among tributaries of the Colorado River in Grand Canyon.

Arizona↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 1. Theory and computation

To establish a theoretical basis for predicting and interpreting the behavior of rapid mass movements on Earth's surface, we develop and test a new computational model for gravity-driven motion of granular avalanches across irregular, three-dimensional (3-D) terrain. The principles embodied in the model are simple and few: continuum mass and momentum conservation and intergranular stress generation governed by Coulomb friction. However, significant challenges result from the necessity of satisfying these principles when deforming avalanches interact with steep and highly variable 3-D terrain. We address these challenges in four ways. (1) We formulate depth-averaged governing equations that are referenced to a rectangular Cartesian coordinate system (with z vertical) and that account explicitly for the effect of nonzero vertical accelerations on depth-averaged mass and momentum fluxes and stress states. (2) We compute fluxes of mass and momentum across vertical cell boundaries using a high-resolution finite volume method and Roe-type Riemann solver. Our algorithm incorporates flux difference splitting, an entropy correction for the flux, and eigenvector decomposition to embed the effects of driving and resisting forces in Riemann solutions. (3) We use a finite element method and avalanche displacements predicted by Riemann solutions to compute Coulomb stresses conjugate to the displacements in 3-D stress space. (4) We test the model output against analytical solutions, a sand cone conceptual experiment, and (in a companion paper) data from detailed laboratory experiments. Model results illustrate a complex interplay of basal traction and internal stress, and they successfully predict not only the gross behavior but also many details of avalanche motion from initiation to deposition.

Journal of Geophysical Research Earth Surface↗