USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Observer”

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 1,243 records · Page 69Linked to original sources

Grain-size-induced weakening of H2O ices I and II and associated anisotropic recrystallization

Grain-size-dependent flow mechanisms tend to be favored over dislocation creep at low differential stresses and can potentially influence the rheology of low-stress, low-strain rate environments such as those of planetary interiors. We experimentally investigated the effect of reduced grain size on the solid-state flow of water ice I, a principal component of the asthenospheres of many icy moons of the outer solar system, using techniques new to studies of this deformation regime. We fabricated fully dense ice samples of approximate grain size 2±1 μm by transforming “standard” ice I samples of 250±50 μm grain size to the higher-pressure phase ice II, deforming them in the ice II field, and then rapidly releasing the pressure deep into the ice I stability field. At T ≤200 K, slow growth and rapid nucleation of ice I combine to produce a fine grain size. Constant-strain rate deformation tests conducted on these samples show that deformation rates are less stress sensitive than for standard ice and that the fine-grained material is markedly weaker than standard ice, particularly during the transient approach to steady state deformation. Scanning electron microscope examination of the deformed fine-grained ice samples revealed an unusual microstructure dominated by platelike grains that grew normal to the compression direction, with c axes preferentially oriented parallel to compression. In samples tested at T ≥220 K the elongation of the grains is so pronounced that the samples appear finely banded, with aspect ratios of grains approaching 50:1. The anisotropic growth of these crystallographically oriented neoblasts likely contributes to progressive work hardening observed during the transient stage of deformation. We have also documented remarkably similar microstructural development and weak mechanical behavior in fine-grained ice samples partially transformed and deformed in the ice II field.

Journal of Geophysical Research B: Solid Earth↗

Synthetic aperture radar interferometry coherence analysis over Katmai volcano group, Alaska

The feasibility of measuring volcanic deformation or monitoring deformation of active volcanoes using space-borne synthetic aperture radar (SAR) interferometry depends on the ability to maintain phase coherence over appropriate time intervals. Using ERS 1 C band (λ=5.66 cm) SAR imagery, we studied the seasonal and temporal changes of the interferometric SAR coherence for fresh lava, weathered lava, tephra with weak water reworking, tephra with strong water reworking, and fluvial deposits representing the range of typical volcanic surface materials in the Katmai volcano group, Alaska. For interferograms based on two passes with 35 days separation taken during the same summer season, we found that coherence increases after early June, reaches a peak between the middle of July and the middle of September, and finally decreases until the middle of November when coherence is completely lost for all five sites. Fresh lava has the highest coherence, followed by either weathered lava or fluvial deposits. These surfaces maintain relatively high levels of coherence for periods up to the length of the summer season. Coherence degrades more rapidly with time for surfaces covered with tephra. For images taken in different summers, only the lavas maintained coherence well enough to provide useful interferometric images, but we found only a small reduction in coherence after the first year for surfaces with lava. Measurement of volcanic deformation is possible using summer images spaced a few years apart, as long as the surface is dominated by lavas. Our studies suggest that in order to make volcanic monitoring feasible along the Aleutian arc or other regions with similar climatic conditions, observation intervals of the satellite with C band SAR should be at least every month from July through September, every week during the late spring/early summer or late fall, and every 2–3 days during the winter.

Journal of Geophysical Research B: Solid Earth↗

Strain accumulation in southern California, 1973-1980

Frequent surveys of seven trilateration networks in southern California over the interval 1973–1980 suggest that a regional increment in strain may have occurred in 1978–1979. Prior to 1978 and after late 1979 the strain accumulation has been predominantly a uniaxial north-south compression. This secular trend was interrupted sometime in 1978–1979 by an increment in both north-south and east-west extension in five of the seven networks. The onset of this change appears to have occurred first in the networks farthest south. The changes occurred without any unusual seismicity within the networks, but the overall seismicity in southern California was unusually low prior to and has been unusually high since the occurrence. The average principal strain rates for the seven networks in the 1973–1980 interval are 0.17 μstrain/yr north-south contraction and 0.08 μstrain/yr east-west extension. Although the observed increment in strain could be related to unidentified systematic error in the measuring system, a careful review of the measurements and comparisons with three other measuring systems reveal no appreciable cumulative systematic error.

Journal of Geophysical Research Solid Earth↗

Transport of atrazine and dicamba through silt and loam soils

The objectives of this research were to determine the role of preferential flow paths in the transport of atrazine (2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine) and dicamba (3-6-dichloro-2-methoxybenzoic acid) through silt and loam soils overlying the High Plains aquifer in Nebraska. In a previous study, 3 of 6 study areas demonstrated high percentages of macropores; those three areas were used in this study for analysis of chemical transport. As a subsequent part of the study, 12 intact soil cores (30-cm diameter by 40-cm height), were excavated sequentially, two from each of the following depths: 0-40cm and 40-80cm. These cores were used to study preferential flow characteristics using dye staining and to determine hydraulic properties. Two undisturbed experimental field plots, each with a 3-m2 surface area, were installed in three study areas in Nebraska. Each was instrumented with suction lysimeters and tensiometers at depths of 10cm to 80cm in 10-cm increments. Additionally, each plot was planted with corn (Zea mays). A neutron probe access tube was installed in each plot to determine soil w ater content at 15-cm intervals. All plots were enclosed w ith a raised frame (of 8-cm height) to prevent surface runoff. All suction lysimeters were purged monthly for three months and were sampled immediately prior to pre-plant herbicide application to obtain background chemical concentrations. Atrazine and dicamba moved rapidly through the soil, but only after a heavy rainfall event, probably owing to the presence of preferential flow paths and lack of microbial degradation in these soil areas. Staining of laboratory cores showed a positive correlation between the percent area stained by depth and the subsequent breakthrough of Br- in the laboratory and leaching of field-applied herbicides owing to large rainfall events. Suction lysimeter samples in the field showed increases in concentrations of herbicides at depths where laboratory data indicated greater percentages of what appeared to be preferential flow paths. Concentrations of atrazine and dicamba exceeding 0.30 and 0.05µg m1-1 were observed at depths of 10-30cm and 50-70cm after two months following heavy rainfall events. It appears from the laboratory experiment that preferential flow paths were a significant factor in transport of atrazine and dicamba.

Nebraska↗

Spatial and seasonal responses of precipitation in the Ganges and Brahmaputra river basins to ENSO and Indian Ocean dipole modes: Implications for flooding and drought

We evaluated the spatial and seasonal responses of precipitation in the Ganges and Brahmaputra basins as modulated by the El Niño Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) modes using Global Precipitation Climatology Centre (GPCC) full data reanalysis of monthly global land-surface precipitation data from 1901 to 2010 with a spatial resolution of 0.5° × 0.5°. The GPCC monthly total precipitation climatology targeting the period 1951–2000 was used to compute gridded monthly anomalies for the entire time period. The gridded monthly anomalies were averaged for the years influenced by combinations of climate modes. Occurrences of El Niño alone significantly reduce (88% of the long-term average (LTA)) precipitation during the monsoon months in the western and southeastern Ganges Basin. In contrast, occurrences of La Niña and co-occurrences of La Niña and negative IOD events significantly enhance (110 and 109% of LTA in the Ganges and Brahmaputra Basin, respectively) precipitation across both basins. When El Niño co-occurs with positive IOD events, the impacts of El Niño on the basins' precipitation diminishes. When there is no active ENSO or IOD events (occurring in 41 out of 110 years), precipitation remains below average (95% of LTA) in the agriculturally intensive areas of Haryana, Uttar Pradesh, Rajasthan, Madhya Pradesh, and Western Nepal in the Ganges Basin, whereas precipitation remains average to above average (104% of LTA) across the Brahmaputra Basin. This pattern implies that a regular water deficit is likely, especially in the Ganges Basin, with implications for the agriculture sector due to its reliance on consistent rainfall for successful production. Historically, major droughts occurred during El Niño and co-occurrences of El Niño and positive IOD events, while major flooding occurred during La Niña and co-occurrences of La Niña and negative IOD events in the basins. This observational analysis will facilitate well-informed decision making in minimizing natural hazard risks and climate impacts on agriculture, and supports development of strategies ensuring optimized use of water resources in best management practice under a changing climate.

Brahmaputra Basin, Ganges Basin↗

Volcanoes: observations and impact

Volcanoes are critical geologic hazards that challenge our ability to make long-term forecasts of their eruptive behaviors. They also have direct and indirect impacts on human lives and society. As is the case with many geologic phenomena, the time scales over which volcanoes evolve greatly exceed that of a human lifetime. On the other hand, the time scale over which a volcano can move from inactivity to eruption can be rather short: months, weeks, days, and even hours. Thus, scientific study and monitoring of volcanoes is essential to mitigate risk. There are thousands of volcanoes on Earth, and it is impractical to study and implement ground-based monitoring at them all. Fortunately, there are other effective means for volcano monitoring , including increasing capabilities for satellite-based technologies.

Book chapter↗

Introduction to the fifth Mars Polar Science special issue: key questions, needed observations, and recommended investigations

The Fifth International Conference on Mars Polar Science and Exploration – which was held from September 12–16, 2011, at the Pike’s Waterfront Lodge in Fairbanks, Alaska – is the latest in a continuing series of meetings that are intended to promote the exchange of knowledge and ideas between planetary and terrestrial scientists interested in Mars polar and climate research (http://www.lpi.usra.edu/meetings/polar2011/polar20113rd.html). The conference was sponsored by the Lunar and Planetary Institute, National Aeronautics and Space Administration, NASA’s Mars Program Office, University of Alaska Fairbanks, International Association of Cryospheric Sciences and the Centre for Research in Earth and Space Sciences at York University.

Icarus↗

Differentiation of volcanic ash-fall and water-borne detrital layers in the Eocene Senakin coal bed, Tanjung Formation, Indonesia

The Sangsang deposit of the Eocene Senakin coal bed, Tanjung Formation, southeastern Kalimantan, Indonesia, contains 11 layers, which are thin (<5 cm) and high in ash (> 70%). These layers are characterized by their pelitic macroscopic texture. Examination of eight of the layers by scanning-electron microscopy, energy-dispersive X-ray, and X-ray diffraction analyses show that they are composed primarily of fairly well-crystallized kaolinite, much of which is vermicular. Accessory minerals include abundant Ti oxide, rare-earth element-rich Ca and A1 phosphates, quartz that luminescences in the blue color range, and euhedral to subhedral pyroxene, hornblende, zircon, and sanidine. Although this mineral suite is suggestive of volcanic ash-fall material, only the four pelitic layers in the middle of the bed are thought to be solely derived from volcanic ash-falls on the basis of diagnostic minerals, replaced glass shards, and lithostratigraphic relationships observed in core and outcrop. The three uppermost pelitic layers contain octahedral chromites, some quartz grains that luminesce in teh orange color range, and some quartz grains that contain two-phase fluid inclusions. These layers are interpreted to be derived from a combination of volcanic ash-fall material and hydrologic transport of volcaniclastic sediment. In contrast, the lowermost pelitic layer, which contains large, rounded FeMg-rich chromites, is thought to have been dominantly deposited by water. The source of the volcanic ash-fall material may have been middle Tertiary volcanism related to plate tectonic activity between Kalimantan and Sulawesi. The volcanic ash was deposited in sufficient amounts to be preserved as layers within the coal only in the northern portions of the Senakin region: the southern coal beds in the region do not contain pelitic layers.

southeastern Kalimantan↗

Voyager 2 at Neptune: Imaging science results

Voyager 2 images of Neptune reveal a windy planet characterized by bright clouds of methane ice suspended in an exceptionally clear atmosphere above a lower deck of hydrogen sulfide or ammonia ices. Neptune's atmosphere is dominated by a large anticyclonic storm system that has been named the Great Dark Spot (GDS). About the same size as Earth in extent, the GDS bears both many similarities and some differences to the Great Red Spot of Jupiter. Neptune's zonal wind profile is remarkably similar to that of Uranus. Neptune has three major rings at radii of 42,000, 53,000, and 63,000 kilometers. The outer ring contains three higher density arc-like segments that were apparently responsible for most of the ground-based occultation events observed during the current decade. Like the rings of Uranus, the Neptune rings are composed of very dark material; unlike that of Uranus, the Neptune system is very dusty. Six new regular satellites were found, with dark surfaces and radii ranging from 200 to 25 kilometers. All lie inside the orbit of Triton and the inner four are located within the ring system. Triton is seen to be a differentiated body, with a radius of 1350 kilometers and a density of 2.1 grams per cubic centimeter; it exhibits clear evidence of early episodes of surface melting. A now rigid crust of what is probably water ice is overlain with a brilliant coating of nitrogen frost, slightly darkened and reddened with organic polymer material. Streaks of organic polymer suggest seasonal winds strong enough to move particles of micrometer size or larger, once they become airborne. At least two active plumes were seen, carrying dark material 8 kilometers above the surface before being transported downstream by high level winds. The plumes may be driven by solar heating and the subsequent violent vaporization of subsurface nitrogen.

Science↗

Flood lavas on Earth, Io and Mars

Flood lavas are major geological features on all the major rocky planetary bodies. They provide important insight into the dynamics and chemistry of the interior of these bodies. On the Earth, they appear to be associated with major and mass extinction events. It is therefore not surprising that there has been significant research on flood lavas in recent years. Initial models suggested eruption durations of days and volumetric fluxes of order 10 7 m 3 s −1 with flows moving as turbulent floods. However, our understanding of how lava flows can be emplaced under an insulating crust was revolutionized by the observations of actively inflating pahoehoe flows in Hawaii. These new ideas led to the hypothesis that flood lavas were emplaced over many years with eruption rates of the order of 10 4 m 3 s −1 . The field evidence indicates that flood lava flows in the Columbia River Basalts, Deccan Traps, Etendeka lavas, and the Kerguelen Plateau were emplaced as inflated pahoehoe sheet flows. This was reinforced by the observation of active lava flows of ≥100 km length on Io being formed as tube-fed flows fed by moderate eruption rates (10 2 –10 3 m 3 s −1 ). More recently it has been found that some flood lavas are also emplaced in a more rapid manner. New high-resolution images from Mars revealed ‘platy–ridged’ flood lava flows, named after the large rafted plates and ridges formed by compression of the flow top. A search for appropriate terrestrial analogues found an excellent example in Iceland: the 1783–1784 Laki Flow Field. The brecciated Laki flow top consists of pieces of pahoehoe, not aa clinker, leading us to call this ‘rubbly pahoehoe’. Similar flows have been found in the Columbia River Basalts and the Kerguelen Plateau. We hypothesize that these flows form with a thick, insulating, but mobile crust, which is disrupted when surges in the erupted flux are too large to maintain the normal pahoehoe mode of emplacement. Flood lavas emplaced in this manner could have (intermittently) reached effusion rates of the order of 10 6 m 3 s −1 .

Journal of the Geological Society↗

Mapping crustal heterogeneity using Lg propagation efficiency throughout the Middle East, Mediterranean, Southern Europe and Northern Africa

In this paper we describe a technique for mapping the lateral variation of Lg characteristics such as Lg blockage, efficient Lg propagation, and regions of very high attenuation in the Middle East, North Africa, Europe and the Mediterranean regions. Lg is used in a variety of seismological applications from magnitude estimation to identification of nuclear explosions for monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These applications can give significantly biased results if the Lg phase is reduced or blocked by discontinuous structure or thin crust. Mapping these structures using quantitative techniques for determining Lg amplitude attenuation can break down when the phase is below background noise. In such cases Lg blockage and inefficient propagation zones are often mapped out by hand. With our approach, we attempt to visually simplify this information by imaging crustal structure anomalies that significantly diminish the amplitude of Lg. The visualization of such anomalies is achieved by defining a grid of cells that covers the entire region of interest. We trace Lg rays for each event/ station pair, which is simply the great circle path, and attribute to each cell a value equal to the maximum value of the Lg/P-coda amplitude ratio for all paths traversing that particular cell. The resulting map, from this empirical approach, is easily interpreted in terms of crustal structure and can successfully image small blockage features often missed by analysis of raypaths alone. This map can then be used to screen out events with blocked Lg prior to performing Q tomography, and to avoid using Lg-based methods of event identification for the CTBT in regions where they cannot work. For this study we applied our technique to one of the most tectonically complex regions on the earth. Nearly 9000 earthquake/station raypaths, traversing the vast region comprised of the Middle East, Mediterranean, Southern Europe and Northern Africa, have been analyzed. We measured the amplitude of Lg relative to the P-coda and mapped the lateral variation of Lg propagation efficiency. With the relatively dense coverage provided by the numerous crossing paths we are able to map out the pattern of crustal heterogeneity that gives rise to the observed character of Lg propagation. We observe that the propagation characteristics of Lg within the region of interest are very complicated but are readily correlated with the different tectonic environments within the region. For example, clear strong Lg arrivals are observed for paths crossing the stable continental interiors of Northern Africa and the Arabian Shield. In contrast, weakened to absent Lg is observed for paths crossing much of the Middle East, and Lg is absent for paths traversing the Mediterranean. Regions that block Lg transmission within the Middle East are very localized and include the Caspian Sea, the Iranian Plateau and the Red Sea. Resolution is variable throughout the region and strongly depends on the distribution of seismicity and recording stations. Lg propagation is best resolved within the Middle East where regions of crustal heterogeneity on the order of 100 km are imaged (e.g., South Caspian Sea and Red Sea). Crustal heterogeneity is resolvable but is poorest in seismically quiescent Northern Africa.

Pure and Applied Geophysics↗

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

Factors controlling landslide frequency-area distributions

A power‐law relation for the frequency–area distribution (FAD) of medium and large landslides (e.g. tens to millions of square meters) has been observed by numerous authors. But the FAD of small landslides diverges from the power‐law distribution, with a rollover point below which frequencies decrease for smaller landslides. Some studies conclude that this divergence is an artifact of unmapped small landslides due to lack of spatial or temporal resolution; others posit that it is caused by the change in the underlying failure process. An explanation for this dilemma is essential both to evaluate the factors controlling FADs of landslides and power‐law scaling, which is a crucial factor regarding both landscape evolution and landslide hazard assessment. This study examines the FADs of 45 earthquake‐induced landslide inventories from around the world in the context of the proposed explanations. We show that each inventory probably involves some combination of the proposed explanations, though not all explanations contribute to each case. We propose an alternative explanation to understand the reason for the divergence from a power‐law. We suggest that the geometry of a landslide at the time of mapping reflects not just one single movement but many, including the propagation of numerous smaller landslides before and after the main failure. Because only the resulting combination of these landslides can be observed due to a lack of temporal resolution, many smaller landslides are not taken into account in the inventory. This reveals that the divergence from the power‐law is not necessarily attributed to the incompleteness of an inventory. This conceptual model will need to be validated by ongoing observation and analysis. Also, we show that because of the subjectivity of mapping procedures, the total number of landslides and total landslide areas in inventories differ significantly, and therefore the shapes of FADs also differ considerably.

Earth Surface Processes and Landforms↗

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth↗

Stress interaction between subduction earthquakes and forearc strike-slip faults: Modeling and application to the northern Caribbean plate boundary

Strike-slip faults in the forearc region of a subduction zone often present significant seismic hazard because of their proximity to population centers. We explore the interaction between thrust events on the subduction interface and strike-slip faults within the forearc region using three-dimensional models of static Coulomb stress change. Model results reveal that subduction earthquakes with slip vectors subparallel to the trench axis enhance the Coulomb stress on strike-slip faults adjacent to the trench but reduce the stress on faults farther back in the forearc region. In contrast, subduction events with slip vectors perpendicular to the trench axis enhance the Coulomb stress on strike-slip faults farther back in the forearc, while reducing the stress adjacent to the trench. A significant contribution to Coulomb stress increase on strike-slip faults in the back region of the forearc comes from "unclamping" of the fault, i.e., reduction in normal stress due to thrust motion on the subduction interface. We argue that although Coulomb stress changes from individual subduction earthquakes are ephemeral, their cumulative effects on the pattern of lithosphere deformation in the forearc region are significant. We use the Coulomb stress models to explain the contrasting deformation pattern between two adjacent segments of the Caribbean subduction zone. Subduction earthquakes with slip vectors nearly perpendicular to the Caribbean trench axis is dominant in the Hispaniola segment, where the strike-slip faults are more than 60 km inland from the trench. In contrast, subduction slip motion is nearly parallel to the Caribbean trench axis along the Puerto Rico segment, where the strike-slip fault is less than 15 km from the trench. This observed jump from a strike-slip fault close to the trench axis in the Puerto Rico segment to the inland faults in Hispaniola is explained by different distributions of Coulomb stress in the forearc region of the two segments, as a result of the change from the nearly trench parallel slip on the Puerto Rico subduction interface to the more perpendicular subduction slip beneath Hispaniola. The observations and modeling suggest that subduction-induced strike-slip seismic hazard to Puerto Rico may be smaller than previously assumed but the hazard to Hispaniola remains high. Copyright 2004 by the American Geophysical Union.

Puerto Rico↗

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.

Journal of Geophysical Research F: Earth Surface↗

Permeability and porosity of the Illinois UPH 3 drillhole granite and a comparison with other deep drillhole rocks

Permeability, porosity, and volumetric strain measurements were conducted on granite cores obtained at depths of 0.7 to 1.6 km from the Illinois UPH 3 drillhole at effective confining pressures from 5 to 100 MPa. Initial permeabilities were in the range of 10-17 to 10-19 m2 and dropped rapidly with applied pressure to values between 10-20 and 10-24 m2 at 100 MPa, typical of other deep granite core samples. These values are several decades lower than equivalent weathered surface granites at comparable effective confining pressures, where weathering products in cracks and pores inhibit crack closure with applied pressure. Permeabilities of the Illinois cores were inversely related to sample depth, suggesting that stress relief and thermal microfractures induced during core retrieval dominated the fluid flow. Thus these samples provide an upper bound on in situ matrix permeability values. A comparison of core permeability from UPH 3 and other deep drillholes shows that stress relief damage can often dominate laboratory permeability measurements. We conclude that it may be difficult to make meaningful estimates of in situ permeability based on either borehole samples (possible damage during retrieval) or surface-derived analogs (altered by weathering). Volumetric strain determined from porosity measurements was compared with differential strain analysis (DSA) data reported by other investigators on samples from the same depths in the drillhole. Our strain measurements (0.002 to 0.005 at 100 MPa) were nearly twice as large as the DSA values, probably because of the crack-enhancing effects of fluids present in our samples that are absent in the dry DSA cores, as well as other time-dependent deformation effects. This difference in observed strain magnitudes between the two measurement methods may be an important consideration if strain and/or porosity data from deep core samples are used in models of stress, fluid circulation, and excess fluid pressure generation in the midcrust.

Journal of Geophysical Research B: Solid Earth↗

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗