Search USGS⌕ Search

SEARCH · Search USGS

Results for “The R 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 721 records · Page 40Linked to original sources
USGS2006

Effects of tributary debris on the longitudinal profile of the Colorado River in Grand Canyon

The Colorado River in Grand Canyon has long been known as a "rapids-and-pools" river, with the rapids owing their existence primarily to tributary debris flows. The debris flows deposit subaerial debris fans that constrict the channel laterally and, when they enter the river, raise the bed elevation. The rapids are short-wavelength (???0.1 to ???1 km), small-amplitude (??????5 m) convexities in the river's longitudinal profile, arising from the shallow gradient in the upstream pool and the steep gradient through the rapid itself. Analysis of the entire longitudinal profile through Grand Canyon reveals two long-wavelength (???100 km), large-amplitude (15-30 m) river profile convexities: the eastern canyon convexity between river mile (RM) 30 and RM 80 and the western canyon convexity between RM 150 and RM 250. Convexities of intermediate scale are also identified in the longitudinal profile. These longer-wavelength, larger-amplitude convexities have strong spatial correlations with high rates of debris flow occurrence, high densities of Holocene debris fans, the largest debris fans along the river, and alluvial thicknesses of 10 m or more. River profile convexities are unstable and require an active and powerful geologic process to maintain them, in this case the abundant, frequent, and voluminous Holocene debris flow activity in Grand Canyon. At all wavelengths the most likely cause for these river profile convexities is Holocene aggradation of the riverbed beneath them, driven by the coarse particles of tributary debris flows. Large enough debris flows will slow river flow for kilometers upstream, causing it to drop much of its suspended load. Integrated over time and all of the tributary point source contributions, this process will build short-wavelength convexities into long-wavelength convexities. For most if not all of the Holocene the Colorado River has been dissipating most of its energy in the rapids and expending the remainder in transporting fine sediment through Grand Canyon, with little or no regional incision of bedrock.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

February 2003 marine atmospheric conditions and the bora over the northern Adriatic

A winter oceanographic field experiment provided an opportunity to examine the atmospheric marine conditions over the northern Adriatic. Mean February winds are from a northeasterly direction over most of the Adriatic and a more northerly direction along the western coast. Wind speeds are fastest in jets over the NE coast during bora events and weakest in the mid-northwestern Adriatic. Diurnal air temperature cycles are smallest on the NE coast and largest in the midwestern Adriatic. The maximum sea-air difference is +10??C on the eastern coast and near zero on the midwestern Adriatic. Boras are northeasterly (from) wind events that sweep off Croatia and Slovenia, bringing slightly colder and drier air over the northern Adriatic. The main bora season is December to March. Winter 2002-2003 was normal for bora events. Synoptic-scale temporal variations are correlated over the northern Adriatic. Fastest Bora winds and highest wind stress over the northern Adriatic is concentrated in four topographically controlled jets. The strongest is the Senj Jet, while the Trieste Jet extends across the entire northern Adriatic. Between each two jets is a weak wind zone. The greatest mean net heat loss is in bora jets in the NE Adriatic, where it was -438 W m-2 and is weakest in the midwestern northern Adriatic, where it was near zero. Wind stress is concentrated over the NE half of Adriatic in four bora jets, while wind stress is weak in the NW Adriatic. There is significant variation in wind stress mean and standard deviation structure over the northern Adriatic with each bora event. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2006

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2006

Aqueous processes at Gusev crater inferred from physical properties of rocks and soils along the Spirit traverse

Gusev crater was selected as the landing site for Spirit on the basis of morphological evidence of long-lasting water activity, including possibly fluvial and lacustrine episodes. From the Columbia Memorial Station to the Columbia Hills, Spirit's traverse provides a journey back in time, from relatively recent volcanic plains showing little evidence for aqueous processes up to the older hills, where rock and soil composition are drastically different. For the first 156 sols, the only evidence of water action was weathering rinds, vein fillings, and soil crust cementation by salts. The trenches of Sols 112-145 marked the first significant findings of increased concentrations of sulfur and magnesium varying in parallel, suggesting that they be paired as magnesium-sulfate. Spirit's arrival at West Spur coincided with a shift in rock and soil composition with observations hinting at substantial amounts of water in Gusev's past. We used the Microscopic Imager data up to Sol 431 to analyze rock and soil properties and infer plausible types and magnitude of aqueous processes through time. We show the role played early by topography and structure. The morphology, texture, and deep alteration shown by the rocks in West Spur and the Columbia Hills Formation (CHF) suggest conditions that are not met in present-day Mars and required a wetter environment, which could have included transport of sulfur, chlorine, and bromine in water, vapor in volcanic gases, hydrothermal circulation, or saturation in a briny fluid containing the same elements. Changing conditions that might have affected flow circulation are suggested by different textural and morphological characteristics between the rocks in the CHF and those of the plains, with higher porosity proxy, higher void ratio, and higher water storage potential in the CHF. Soils were used to assess aqueous processes and water pathways in the top layers of modern soils. We conclude that infiltration might have become more difficult with time.

Journal of Geophysical Research E: Planets↗
USGS2006

Acoustic bed velocity and bed load dynamics in a large sand bed river

Development of a practical technology for rapid quantification of bed load transport in large rivers would represent a revolutionary advance for sediment monitoring and the investigation of fluvial dynamics. Measurement of bed load motion with acoustic Doppler current profiles (ADCPs) has emerged as a promising approach for evaluating bed load transport. However, a better understanding of how ADCP data relate to conditions near the stream bed is necessary to make the method practical for quantitative applications. In this paper, we discuss the response of ADCP bed velocity measurements, defined as the near-bed sediment velocity detected by the instrument's bottom-tracking feature, to changing sediment-transporting conditions in the lower Missouri River. Bed velocity represents a weighted average of backscatter from moving bed load particles and spectral reflections from the immobile bed. The ratio of bed velocity to mean bed load particle velocity depends on the concentration of the particles moving in the bed load layer, the bed load layer thickness, and the backscatter strength from a unit area of moving particles relative to the echo strength from a unit area of unobstructed bed. A model based on existing bed load transport theory predicted measured bed velocities from hydraulic and grain size measurements with reasonable success. Bed velocities become more variable and increase more rapidly with shear stress when the transport stage, defined as the ratio of skin friction to the critical shear stress for particle entrainment, exceeds a threshold of about 17. This transition in bed velocity response appears to be associated with the appearance of longer, flatter bed forms at high transport stages.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Modeling the evolution of channel shape: Balancing computational efficiency with hydraulic fidelity

The cross-sectional shape of a natural river channel controls the capacity of the system to carry water off a landscape, to convey sediment derived from hillslopes, and to erode its bed and banks. Numerical models that describe the response of a landscape to changes in climate or tectonics therefore require formulations that can accommodate evolution of channel cross-sectional geometry. However, fully two-dimensional (2-D) flow models are too computationally expensive to implement in large-scale landscape evolution models, while available simple empirical relationships between width and discharge do not adequately capture the dynamics of channel adjustment. We have developed a simplified 2-D numerical model of channel evolution in a cohesive, detachment-limited substrate subject to steady, unidirectional flow. Erosion is assumed to be proportional to boundary shear stress, which is calculated using an approximation of the flow field in which log-velocity profiles are assumed to apply along vectors that are perpendicular to the local channel bed. Model predictions of the velocity structure, peak boundary shear stress, and equilibrium channel shape compare well with predictions of a more sophisticated but more computationally demanding ray-isovel model. For example, the mean velocities computed by the two models are consistent to within ???3%, and the predicted peak shear stress is consistent to within ???7%. Furthermore, the shear stress distributions predicted by our model compare favorably with available laboratory measurements for prescribed channel shapes. A modification to our simplified code in which the flow includes a high-velocity core allows the model to be extended to estimate shear stress distributions in channels with large width-to-depth ratios. Our model is efficient enough to incorporate into large-scale landscape evolution codes and can be used to examine how channels adjust both cross-sectional shape and slope in response to tectonic and climatic forcing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

THEMIS high-resolution digital terrain: Topographic and thermophysical mapping of Gusev Crater, Mars

We discuss a new technique to generate high-resolution digital terrain models (DTMs) and to quantitatively derive and map slope-corrected thermophysical properties such as albedo, thermal inertia, and surface temperatures. This investigation is a continuation of work started by Kirk et al. (2005), who empirically deconvolved Thermal Emission Imaging System (THEMIS) visible and thermal infrared data of this area, isolating topographic information that produced an accurate DTM. Surface temperatures change as a function of many variables such as slope, albedo, thermal inertia, time, season, and atmospheric opacity. We constrain each of these variables to construct a DTM and maps of slope-corrected albedo, slope- and albedo-corrected thermal inertia, and surface temperatures across the scene for any time of day or year and at any atmospheric opacity. DTMs greatly facilitate analyses of the Martian surface, and the MOLA global data set is not finely scaled enough (128 pixels per degree, ???0.5 km per pixel near the equator) to be combined with newer data sets (e.g., High Resolution Imaging Science Experiment, Context Camera, and Compact Reconnaissance Imaging Spectrometer for Mars at ???0.25, ???6, and ???20 m per pixel, respectively), so new techniques to derive high-resolution DTMs are always being explored. This paper discusses our technique of combining a set of THEMIS visible and thermal infrared observations such that albedo and thermal inertia variations within the scene are eliminated and only topographic variations remain. This enables us to produce a high-resolution DTM via photoclinometry techniques that are largely free of albedo-induced errors. With this DTM, THEMIS observations, and a subsurface thermal diffusion model, we generate slope-corrected maps of albedo, thermal inertia, and surface temperatures. In addition to greater accuracy, these products allow thermophysical properties to be directly compared with topography.

Journal of Geophysical Research E: Planets↗
USGS2009

Compositions of modern dust and surface sediments in the Desert Southwest, United States

Modern dusts across southwestern United States deserts are compositionally similar to dust-rich Av soil horizons (depths of 0-0.5 cm and 1-4 cm at 35 sites) for common crustal elements but distinctly different for some trace elements. Chemical compositions and magnetic properties of the soil samples are similar among sites relative to dust sources, geographic areas, and lithologic substrates. Exceptions are Li, U, and W, enriched in Owens Valley, California, and Mg and Sr, enriched in soils formed on calcareous fan gravel in southeast Nevada. The Av horizons are dominated by dust and reflect limited mixing with substrate sediments. Modern dust samples are also similar across the region, except that Owens Valley dusts are higher in Mg, Ba, and Li and dusts both there and at sites to the north on volcanic substrates are higher in Sb and W. Thus, dust and Av horizons consist of contributions from many different sources that are well mixed before deposition. Modern dusts contain significantly greater amounts of As, Cd, Cr, Cu, Ni, Pb, and Sb than do Av horizons, which record dust additions over hundreds to thousands of years. These results suggest that modern dust compositions are influenced by anthropogenic sources and emissions from Owens (dry) Lake after its artificial desiccation in 1926. Both modern dusts and Av horizons are enriched in As, Ba, Cu, Li, Sb, Th, U, and W relative to average crustal composition, which we interpret to indicate that the geologic sources of dust in the southwestern United States are geochemically distinctive.

Journal of Geophysical Research F: Earth Surface↗
USGS2010

The 2007 M5.4 Alum Rock, California, earthquake: Implications for future earthquakes on the central and southern Calaveras Fault

[1] The similarity of seismograms recorded by two seismic stations demonstrate that the 31 October 2007 moment magnitude M 5.4 Alum Rock earthquake is a repeat of a 1955 M L 5.5 earthquake. Both occurred on Oppenheimer et al.'s (1990) Zone V “stuck patch” on the central Calaveras fault, providing new support for their model of Calaveras fault earthquake activity. We suggest that Zone V fails only in a family of recurring M ∼ 5.4–5.5 earthquakes. The 1955 and 2007 earthquakes are the penultimate and ultimate Zone V events. Earthquakes in 1891 and 1864 are possible earlier Zone V events. The next Zone V event is not expected in the next few decades, assuming a time‐dependent recurrence model: the mean forecast date is 2064 (2035–2104, 95% confidence range). We further suggest that Zones I, II, III, and IV fail in recurring M ∼ 5.1–5.3, M ∼ 5.6–5.8, M ∼ 6.1–6.3, and M ∼ 4.9–5.0 earthquakes, respectively. If our earthquake recurrence model is correct, the next Zone I event is overdue and could occur anytime, and M 5–6 earthquakes should not occur on Zones II, III, and IV before 2014, 2012, and 2026, respectively. We cannot rule out the possibility that Zone VI, which lies at the southern end of the Mission Seismic Trend, where the southern Hayward and central Calaveras faults appear to connect at depth, fails aseismically or in large events on the southern Hayward fault, such as last occurred in 1868, or in large events on the adjoining northern Calaveras fault segment.

California↗
USGS2019

Azimuthal seismic anisotropy of 70 Ma Pacific‐plate upper mantle

Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2D corner flow and plate‐driven shear near mid‐ocean ridges. We measure 6.0 ± 0.3% anisotropy at the Moho in 70 Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2D mantle flow near a mid‐ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Paleomagnetic constraints on the geometry and timing of deformation at Yucca Mountain, Nevada

Paleomagnetic data from Miocene ash flow sheets indicate that southern Yucca Mountain has undergone about 30° of clockwise vertical-axis rotation since emplacement of the Tiva Canyon Member of the Paintbrush Tuff at about 13 Ma. Declinations of remanent magnetization from 32 sites in the reversely magnetized Tiva Canyon Member display a systematic southward increase over the 25 km north–south extent of Yucca Mountain. Probable errors in estimating structural attitudes are a major source of uncertainty in determining the amount of rotation between any two sites. Analysis of tilt-corrected directions of remanent magnetization from the Tiva Canyon sites indicates that structural attitudes determined from outcrop patterns on geologic maps are generally more accurate than field measurements of attitudes of compaction foliation or contacts. Attitudes determined from map patterns apparently minimize effects of (1) initial dips, (2) poorly defined compaction foliation, and (3) tilting that occurred after compaction but before blocking of remanence; however, such attitudes cannot accurately represent tilting related to unmapped local structures. Rotations implied by data from three older ash flow sheets, the Topopah Spring Member of the Paintbrush Tuff (seven sites), and the Prow Pass (four sites) and Bullfrog (three sites) Members of the Crater Flat Tuff are generally consistent with the amount and sense of rotation documented by data from the Tiva Canyon Member. Uncertainties in the declinations of tilt-corrected site mean directions from these older units preclude evaluating whether older units have undergone more rotation than the Tiva Canyon Member.

Journal of Geophysical Research↗
USGS2012

Magnetic susceptibility and relation to initial 87Sr/86Sr for granitoids of the central Sierra Nevada, California

Measurement of the magnetic susceptibility of more than 6000 samples of granitic rock from the Mariposa 1° by 2° quadrangle, which crosses the central part of the Sierra Nevada batholith between 37° and 38°N latitude, shows that magnetic susceptibility values are above 10 −2 SI units in the east and central parts of the batholith and drop abruptly to less than 10 −3 SI units in the western foothills. In a narrow transitional zone, intermediate values (10 −3 to 10 −2 ) prevail. Magnetic susceptibility appears to decrease slightly westward within the zones of both high and low values. Magnetic susceptibility in plutonic rocks is chiefly a function of the abundance of magnetite, which depends, in turn, on the total iron content of the rocks and their oxidation ratio. Lower magnetic susceptibilities of felsic members of Sierran intrusive suites and of some felsic rock units relative to adjacent mafic rock units commonly reflect differences in total iron content, but the differences of magnetic susceptibility that define the regional pattern generally are much larger and are determined chiefly by the oxidation ratios of the rocks. The relatively unaltered condition of the samples and restriction of δ 18 O to the range of +7 to +10.3‰ indicate that neither hydrothermal fluids nor subsolidus alteration were important in modifying oxidation ratios. Correlations of magnetic susceptibility with initial 87 Sr/ 86 Sr suggest that oxidation ratios have been inherited from the source regions for the magmas from which the rocks crystallized. Reduction of Fe 3+ to Fe 2+ by organic carbon or other reducing substances may also have affected magnetic susceptibility.

Journal of Geophysical Research↗
USGS2012

Relation of peralkaline magmatism to heterogeneous extension during the Middle Miocene, southeastern Nevada

Volcanism migrated southward in the northern Basin and Range province in the Oligocene and early Miocene to produce voluminous calcalkaline silicic ash flow tuffs. Alkaline volcanism became dominant by middle Miocene (17–14 Ma) as smaller volumes of rhyolite-trachyte-basalt suites were erupted from the relatively small Kane Springs Wash caldera complex including the Narrow Canyon, Boulder Canyon, and Kane Springs Wash calderas in southeastern Nevada. Only minor extension affected the Kane Wash area before the end of calcalkaline activity, but extension expressed by rate of progressive stratal tilt peaked (15–13.5 Ma) with peralkaline magmatism (14.7–14.4 Ma). Variations in distribution, degree, style, and timing of deformation demonstrate heterogeneous extension in the Kane Wash area. Only minor extension and tilting persisted post-middle Miocene (<12 Ma). All major eruptive sources overlap domains of rapid extension. Most of the eruptive volumes from the two oldest calderas of the complex apparently pooled within their calderas, creating outflow deficits. Denudation faulting associated with magmatic tumescence may have followed preexisting active extensional fault systems to unload magma chambers, thus triggering eruptions into structural depressions. Evolution of alkaline magmas is demonstrated by progressive increases in peralkalinity and high field strength elements such as Zr, Y, and Nb. Nd, Pb, and Sr isotopic compositions provide evidence that significantly less crustal interaction affected middle Miocene peralkaline magmas than pre-middle Miocene calcalkaline magmas. ε Nd values are −5 to −7 for peralkaline magmas and −7 to −11 for calcalkaline magmas; 208 Pb/ 204 Pb ratios are 38.2–38.6 for peralkaline magmas and 38.5–38.9 for calcalkaline magmas. Regional cooling, short duration of magmatism, small volumes of magma, and local extension caused less crustal interaction in peralkaline Kane Wash magmas than in earlier magmas. North of the Kane Wash area, older more voluminous calcalkaline magmas intruded hotter crust for a longer period and thus interacted with the crust to a greater degree in spite of synvolcanic extension.

Journal of Geophysical Research↗
USGS1997

Timing of detachment faulting in the Bullfrog Hills and Bare Mountain area, southwest Nevada: Inferences from 40Ar/39Ar, K-Ar, U-Pb, and fission track thermochronology

Crustal extension in the Bullfrog Hills and Bare Mountain area of southwest Nevada is associated with movement along a regional detachment fault. Normal faulting in the upper plate and rapid cooling (denudation) of the lower plate were coeval with Miocene silicic volcanism and with west-northwest transport along the detachment fault. A west-northwest progression of tilting along upper plate normal faults is indicated by ages of the volcanic rocks in relation to angular unconformities. Near the breakaway, tilting in the upper plate occurred between 12.7 and 11.6 Ma, continued less strongly past 10.7 Ma, and was over by 8.2 Ma. Ten to 20 km west of the breakaway, tilting occurred between 10.7 and 10.33 Ma, continued less strongly after 10.33 Ma, and was over by 8.1 Ma. The cooling histories of the lower plate metamorphic rocks were determined by thermochronologic dating methods: K-Ar and 40 Ar/ 39 Ar on muscovite, biotite, and hornblende, 40 Ar/ 39 Ar on K-feldspar, U-Pb on apatite, zircon, and sphene, and fission track on apatite, zircon, and sphene. Lower plate rocks 10 km west of the breakaway cooled slowly from Early Cretaceous lower-amphibolite facies conditions through 350±50° to 300±50°C between 57 and 38 Ma, then cooled rapidly from 205±50° to 120±50°C between 12.6±1.6 and 11.1±1.9 Ma. Lower plate rocks 20 km west of the breakaway cooled slowly from Early Cretaceous upper-amphibolite facies conditions through 500±50°C at 78–67 Ma, passed through 350±50° to 300±50°C between 16.3±0.4 and 10.5±0.3 Ma, then cooled rapidly from 285±50° to 120±50°C between 10.2 and 8.6 Ma. Upper plate tilting and rapid cooling (denudation) of the lower plate occurred simultaneously in the respective areas. The early slow-cooling part of the lower plate thermal histories was probably related to erosion at the Earth's surface, which stripped off about 9 km of material in 50 to 100 m.y. The results indicate an initial fault dip ≥30° and a 12 mm yr −1 west-northwest migration of the locus of rapid tilting in the upper plate.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Variable rates of late Quaternary strike slip on the San Jacinto fault zone, southern California.

3 strike slip displacements of strata with known approximate ages have been measured at 2 locations on the San Jacinto fault zone. Minimum horizontal offset between 5.7 and 8.6km in no more than 0.73Myr NE of Anza indicates 8-12 mm/yr average slip rate since late Pleistocene time. Horizontal slip of 1.7m has been calculated for the youngest sediment of Lake Cahuilla since its deposition 271- 510 yr BP. The corresponding slip rate is 2.8-5.0 mm/yr. Right lateral offset of 10.9m measured on a buried stream channel older than 5060 yr BP but younger than 6820 yr BP yields average slip rates for the intermediate time periods, 400 to 6000 yr BP of 1-2 mm/yr. The rates of slip suggest a relatively quiescent period from about 4000 BC to about 1600 AD.-from Author

Journal of Geophysical Research↗
USGS2012

Tectonic setting for ophiolite obduction in Oman

The Samail ophiolite is part of an elongate belt in the Middle East that forms an integral part of the Alpine mountain chains that make up the northern boundary of the Arabian-African plate. The Samail ophiolite represents a portion of the Tethyan ocean crust formed at a spreading center of Middle Cretaceous age (Cenomanian). During the Cretaceous spreading of the Tethyan Sea, Gondwana Land continued its dispersal, and the Arabian-African plate drifted northward about 10°. These events combined with the opposite rotation of Eurasia and Africa initiated the closing of the Tethyan during the Late Cretaceous. At the early stages of closure, downwarping of the Arabian continental margin combined with the compressional forces of closure from the Eurasian plate initiated obduction of the Tethyan oceanic crust along preexisting transform faults, and still hot oceanic crust was detached along oblique northeast dipping thrust faults. Amphibolites developed at the base of the detached hot peridotite as it was thrust southward over oceanic volcanic and sedimentary rocks. Plate configurations combined with palinspastic reconstructions show that subduction and attendant large-scale island arc volcanism did not commence until after the Tethyan sea began to close and after the Samail ophiolite was emplaced southward across the Arabian continental margin. The Samail ophiolite nappe now rests upon a melange consisting mainly of pelagic sediments, volcanics, and detached fragments of the basal amphibolites which in turn rest on autochthonous shelf carbonates of the Arabian platform. Laterites and conglomerates with reworked laterites on the eroded upper surface of the ophiolite indicate a period of emergence prior to the deposition of shallow water Maestrichtian carbonates. Following emplacement (Eocene) of the Samail ophiolite, the Tethyan oceanic crust began northward subduction, and active arc volcanism started just north of the present Jaz Murian depression in Iran.

Journal of Geophysical Research Solid Earth↗
USGS2012

An assessment of the accuracy of the geodetic measurements that define the southern California uplift

Examination of the charge that the geodetic measurements which define the southern California uplift are seriously flawed by height-dependent systematic errors indicates that this charge is unfounded. Our review of those factors designed to minimize the effects of rod error shows that systematic errors attributable to rod miscalibration are trivial and would tend to cancel during the particularly pertinent period 1955–1965. Comparisons between the results of measurements in which the correlation between topography and signal is generally poor reveal large aseismic tilts in a number of places within and around the margins of the uplift. Especially significant in this context are the results of preuplift and postuplift levelings over routes characterized by diverse length, topography, and atmospheric conditions that produce closely matching, temporally equivalent heights for a representative bench mark within the uplift. Statistical analyses of the short-wavelength components of uplift signal and elevation along a frequently repeated survey line indicate variable correlations between the two. However, the spatial and temporal patterns of these correlations are inconsistent with their attribution to height-dependent systematic errors. Statistically significant short-wavelength correlations along this survey line are reasonably explained in part as due to real movement associated with differentially subsiding bench marks, and there is a strong likelihood that this movement tends to dominate many of the correlations.

Journal of Geophysical Research Solid Earth↗
USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗