USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics Research”

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,405 records · Page 78Linked to original sources
USGS2012

Seismicity and detection/location threshold in the southern Great Basin seismic network

A spatially varying model of the detection/location capabilities of the Southern Great Basin seismic network (SGBSN) has been derived that is based on simple empirical relations and statistics. This permits use of almost all the catalog data gathered; instead of ignoring data that are below the threshold of completeness, a spatially varying threshold model is developed so that subregions having lower completeness levels than the network as a whole can be outlined and the completeness level of each sub-region determined. Such a model is required to unambiguously identify regions that are aseismic due to natural processes rather than to limited detection and/or location capabilities. Accounting for spatial variations in detection/location threshold is also important for studies in which magnitude-frequency distributions are interpreted in terms of source scaling properties. The characteristics of the spatial distribution of earthquakes, where earthquake clusters and aseismic regions locate, appear to be stable at all magnitude levels so that inferences about where strain is being accommodated will be the same whether numbers of earthquakes or strain estimated from seismic moments are examined. For the southern Great Basin region these principal characteristics include clusters at the northern end of the Furnace Creek fault and in the Pahranagat Shear Zone, and a relatively large number of earthquakes in the northern and southeastern portions of the Nevada Test Site. These clusters cover regions much larger than the surface projections of any of the mapped faults. The extent to which seismicity is induced by nuclear testing is unclear. The predominantly aseismic regions include the area west of the Death Valley/Furnace Creek fault system and an almost complete absence of events at Yucca Mountain. Finally, a considerable number of isolated events in the SGBSN catalog cannot be correlated with mapped faults.

Journal of Geophysical Research↗
USGS2012

Late Cretaceous paleomagnetism of the Tucson Mountains: Implications for vertical axis rotations in south central Arizona

The Tucson Mountains of southern Arizona are the site of an Upper Cretaceous caldera from which the rhyolitic Cat Mountain Tuff was erupted at about 72 Ma. Two magnetic units within the Cat Mountain Tuff are distinguished by paleomagnetic data in both the northern and southern Tucson Mountains. The resurgent Amole pluton (≃72 Ma) in the northern Tucson Mountains was emplaced soon after eruption of the Cat Mountain Tuff but cooled and was magnetized after northeastward tilting (50°–85°) of the adjacent caldera-fill sequence. Petrologic and paleomagnetic data indicate that the lower magnetic unit of the Cat Mountain Tuff caps the Silver Bell Mountains to the northwest. A previous paleomagnetic investigation ( N = 34) indicates that the Silver Bell Mountains have been rotated clockwise 30°±16° (95% confidence level) about a vertical axis relative to cratonic North America. A similar paleomagnetic study of Upper Cretaceous volcanic, volcaniclastic, and intrusive units in the Tucson Mountains ( N = 26) indicates that these rocks have been rotated 7°±14° clockwise relative to stable North America. A direct comparison of paleomagnetic directions for the lower unit of the Cat Mountain Tuff shows a 17°±10° clockwise rotation between the Silver Bell Mountains and the Tucson Mountains which supports the relative accuracy of the absolute rotations determined for these two mountain ranges. Preliminary paleomagnetic directions for middle Tertiary units from the Silver Bell and Tucson Mountains are consistent with clockwise rotation having occurred prior to deposition of these rocks. Clockwise rotation of crustal blocks in southern Arizona likely was associated with strike-slip movement on major northwest trending faults in the region, and this movement may have been related to oblique subduction of oceanic plates along the western continental margin in Late Cretaceous and early Tertiary time. The available paleomagnetic data indicate that rocks in southern Arizona have not remained unrotated with respect to North America since Late Cretaceous time and that vertical axis rotations may have played an important role in the region during Laramide deformation.

Journal of Geophysical Research↗
USGS2012

Comparison of Vibroseis and explosive source methods for deep crustal seismic reflection profiling in the Basin and Range province

Direct comparison of low-fold, high-energy explosive and high-fold, lower-energy Vibroseis methods for acquiring deep crustal seismic reflection data in the Basin and Range Province suggests that the high-fold common midpoint (CMP) method there does not provide the best possible image of lower crustal structure. During the recent acquisition of a Vibroseis profile in the Basin and Range Province we fired single deep shot holes to obtain a coincident single-fold explosive section. Within the upper crust (upper 3 s) the explosive source and Vibroseis records are nearly equivalent. For record times below 3 s, however, comparison of the explosive source gathers and the coincident final 60-fold Vibroseis section demonstrates that low-fold explosive profiling provides a higher-quality image of the midcrust to lower crust (3–10 s). The higher record quality of the explosive sources results primarily from the larger seismic energy levels produced by the explosives, making them less sensitive to common noise sources. Whereas deeper than 4–5 s the Vibroseis energy levels on individual source efforts fall to that of ambient noise levels, the explosions provide signal-generated energy exceeding ambient noise levels down to 18–19 s. Although individual reflections can be correlated on explosion and Vibroseis shot gathers, reflection events on the 60-fold Vibroseis stack do not correlate to those on the single-fold explosion profile, suggesting that the high-fold CMP method in our study did not maintain the integrity of the weak lower crustal reflected arrivals. Reasons why the high-fold CMP method apparently failed include complex, even time-varying, statics, nonhyperbolic moveout at long offsets, and the difficulty in resolving stacking velocities with data having low signal-to-noise ratios. Reflections on the explosion section are longer and imply a greater degree of layering than one would infer from the lower-energy Vibroseis section.

Journal of Geophysical Research↗
USGS2012

Visible and near-infrared (0.4-2.5 μm) reflectance spectra of playa evaporite minerals

Visible and near-infrared (VNIR; 0.4–2.4 μm) reflectance spectra were recorded for 35 saline minerals that represent the wide range of mineral and brine chemical compositions found in playa evaporite settings. The spectra show that many of the saline minerals exhibit diagnostic near-infrared absorption bands, chiefly attributable to vibrations of hydrogen-bonded structural water molecules. VNIR reflectance spectra can be used to detect minor hydrate phases present in mixtures dominated by anhydrous halite or thenardite, and therefore will be useful in combination with X ray diffraction data for characterizing natural saline mineral assemblages. In addition, VNIR reflectance spectra are sensitive to differences in sample hydration state and should facilitate in situ studies of minerals that occur as fragile, transitory dehydration products in natural salt crusts. The use of spectral reflectance measurements in playa studies should aid in mapping evaporite mineral distributions and may provide insight into the geochemical and hydrological controls on playa mineral and brine development.

Journal of Geophysical Research↗
USGS2012

Chemical bonding in the outer core: high-pressure electronic structures of oxygen and sulfur in metallic iron

From its density the outer core is believed to be an alloy of iron and a light element such as sulfur or oxygen. The nature of the light element in the core is an important constraint for theories of the Earth's formation. In this paper the electronic structures of oxygen and sulfur impurities in metallic iron are investigated to determine if pressure, temperature, and composition-induced changes in bonding might affect phase equilibria along the Fe-FeS and Fe-FeO binaries. The electronic structure of sulfur in metallic iron is consistent with the miscibility between Fe and FeS liquids. Volume compression strengthens the Fe-S bond, and it is expected that at sufficiently high pressure, sulfur can substitute for Fe and give solid solution behavior between Fe and FeS. In contrast, the electronic structure of oxygen in metallic iron shows that oxygen cannot act as a substitutional impurity (replacing Fe). This explains the observed miscibility gap on the Fe-FeO binary at 1 atm pressure. Volume compression does not greatly change the electronic structure if oxygen substitutes for iron in bcc and fcc iron. Iron-oxygen bonding does occur, however, if oxygen occupies interstitial sites. Insofar as the molar volume of FeO incorporated as interstitial oxygen in metallic iron is smaller than that of pure FeO, the incorporation of oxygen into metallic iron may be favored under the pressures of the Earth's core.

Journal of Geophysical Research↗
USGS2012

Imaging the Juan de Fuca plate beneath southern Oregon using teleseismic P wave residuals

The Juan de Fuca plate subducts under the western margin of the North American plate in southern British Columbia, Washington, Oregon, and northern California. Benioff zone seismicity delineates shallow parts of the plate in Washington and northern California, whereas in Oregon there is a lack of seismicity. This study images the Juan de Fuca plate in southern Oregon using seismic tomography. We inverted P wave travel time residuals from a 366-km-long seismic array operated in southern Oregon in 1982. The southeast striking array extended from the Coast ranges to the Modoc Plateau and crossed the High Cascades at Crater Lake, Oregon. We imaged three features under the array: one high-velocity zone and two low-velocity zones. The high-velocity zone is 3–4% faster than the surrounding upper mantle. It dips steeply at 65° to the east beneath the Cascade Range and extends down to at least 200 km. We propose that this high-velocity feature is subducted Juan de Fuca plate. We also imaged two low-velocity zones, both of which are 3–4% slower than the surrounding earth structure. The southeastern low-velocity zone may be caused by partially molten crust underlying the Crater Lake volcano region.

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

Rift flank uplifts and Hinterland Basins: Comparison of the Transantarctic Mountains with the Great Escarpment of southern Africa

Uplifted rift margins are a common feature of continents and oceans. Two variants of rift flank morphologies have been recognized: One in which the topography warps down from an inland high toward the continental margin, and one where the tropographic peak lies close to the continental margin. The Great Escarpment of southern Africa and the Transantarctic Mountains are examples of the first and the second variants of rift flanks, respectively. Both rift flanks are bordered on their landward side by broad continental basins: the Kalahari and the Wilkes hinterland basins. If these basins are interpreted as flexural “outer lows” that deepen in unison with the uplift of the rift flanks, the lithosphere on the uplifted side is very rigid in both cases (elastic thickness T e of 100 ± 20 km for southern Africa and 110 ± 20 km for East Antarctica). We suggest that the variation in rift flank morphology is caused by the isostatic response to uplift forces of elastic plates sharing different boundary conditions. We model the uplift of the Transantarctic Mountains as an upward deflection of an elastic plate which is broken at the front of the Transantarctic Mountains, and we model the uplift of the Great Escarpment as an upward deflection of a continuous elastic plate that is modified by the downward load of sediments on the continental margin. Although the Transantarctic Mountain uplift is young (60–0 Ma) and the southern African uplift is old (<100 Ma), the different isostatic responses of the two margins are not a function of age, because most loading (sedimentation) and unloading (erosion) took place shortly after rifting. Detailed modeling of topography, gravity, geological markers, and the locations of depocenters suggests that lithospheric rigidity decreases under the Transantarctic Mountains, whereas in southern Africa the decrease occurs not under the Great Escarpment but far seaward under the continental shelf and slope. If the distribution of lithospheric rigidity is indicative of the thermal regime of the lithosphere, then uplifted rift flanks are not always underlain by a thermal anomaly. This and other geological evidence indicate that a single mechanism cannot explain the uplift of both the Antarctic and the African margins.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Variations in the reflectivity of the moho transition zone beneath the Midcontinent Rift System of North America: results from true amplitude analysis of GLIMPCE data

True amplitude processing of The Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection data from the Midcontinent Rift System of North America shows large differences in the reflectivity of the Moho transition zone beneath the axial rift, beneath the rift flanks, and outside of the rift. The Moho reflection from the axial rift has a discontinuous, diffractive character marginally stronger (several decibels) than an otherwise transparent lower crust and upper mantle. Beneath the axial rift, Moho is interpreted to be a synrift igneous feature. Beneath the rift flanks, the reflectivity of the Moho transition is generally well developed with two identifiable boundaries, although in places it is weakly reflective to nonreflective, similar to Moho outside the rift. The two boundaries are interpreted as the base of essentially intact, although stretched, prerift Archean crust (upper boundary) and new synrift Moho 1-2 s (6-7 km) deeper (lower boundary). Beneath the rift flanks, the layered reflection Moho transition results from the preexisting crustal composition and fabric modified by synrift igneous processes and extensional tectonic/metamorphic processes. The geologic evidence for extensive basaltic magmatism in the rift is the basis for interpreting the Moho signature as a Keweenawan structure that has been preserved for 1.1 b.y. Extension and magmatism appear to enhance reflectivity in the lower crust and Moho transition zone only where stretching factors are moderate (rift flanks) and not where they are extreme (axial rift). This leads to the prediction that the reflectivity across analogous volcanic passive continental margins should be greatest beneath the moderately stretched continental shelves and should decrease towards the ocean-continent boundary.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Strain accumulation across the Wasatch Fault near Ogden, Utah

Deformation of a 70 by 40‐km trilateration network spanning the north trending Wasatch fault near Ogden, Utah, has been monitored from 1972 through 1990. All but nine of the 200 measurements are consistent with deformation that is linear in time. We presume that those nine observations are contaminated by some blunder in making the measurements and that deformation is linear in time. The strain rate over the 1972–1990 interval across a 40‐km‐wide zone west of the Wasatch fault (hanging wall block) is adequately approximated by a 0.04 ± 0.01 μstrain/yr N85°E±5° uniaxial extension. Adequate coverage of the footwall block is available only over the 1981–1990 interval, and the strain rate across a 30‐km‐wide zone there is approximated by a 0.07 ± 0.03 μstrain/yr N20°E±12° uniaxial extension. The observed extension of the hanging wall block in the direction normal to fault strike suggests that the Wasatch fault may be listric; a planar fault would produce a zone of contraction normal to the fault at the surface of the hanging wall block. The component of extension parallel to fault strike observed on the footwall block is unexplained.

Utah↗
USGS2012

Thin, low‐velocity crust beneath the southern Yukon‐Tanana Terrane, east central Alaska: Results from Trans‐Alaska crustal transect refraction/wide‐angle reflection data

A seismic refraction/wide‐angle reflection survey for the Trans‐Alaska Crustal Transect program reveals a thin, reflective crust beneath the southern Yukon‐Tanana terrane (YTT) in east central Alaska. These data are the first detailed refraction survey of the southern YTT and compose a 130‐km‐long reversed profile along the Alaska and Richardson highways. Results from this study indicate that low‐velocity (≤ 6.4 km/s) rocks extend to approximately 27 km in depth. Based on these low velocities and an average Poisson's ratio of 0.23 determined for depths of ≤27 km, an overall silicic composition is interpreted for this portion of the crust beneath the Yukon‐Tanana terrane. From approximately 8 to 27 km depth the crust exhibits an increase in reflectivity. This middle to lower crustal reflectivity is modeled as alternating high‐ and low‐velocity lamellae with an average velocity of 6.1 km/s at 10 km depth to an average velocity of 6.4 km/s at 27 km depth. Beneath these reflective, low‐velocity rocks a 3‐ to 5‐km‐thick, 7.0 km/s basal crustal layer produces a prominent reflection that extends to offsets of up to 280 km. The crust‐mantle boundary, modeled at an average depth of 30 km, produces a variable PmP reflection, which may indicate lateral heterogeneity of this boundary, and a weak and emergent Pn refraction with a velocity of 8.2 km/s. We interpret the crustal section as follows: the low‐velocity rocks of the southern YTT extend from the surface to depths of approximately 10 km; underthrust Mesozoic flysch of the Kahiltna terrane, rocks of the Gravina arc, and basement of the Wrangellia(?) terrane extend from 10 to 27 km depth; a 3‐ to 5‐km‐thick layer of mantle‐derived mafic rocks, relic oceanic crust, or Wrangellia(?) terrane lower crust extends from 27 to approximately 30 km depth; a tectonically young Moho beneath the southern YTT is found at an average depth of 30 km; and it is underlain by a mantle that may be relatively cool and/or olivine rich. In this interpretation, the Yukon‐Tanana terrane is a thin‐skinned terrane. Our results indicate that tectonic, and possibly magmatic, underplating has played a significant role in crustal growth for central Alaska.

Alaska↗
USGS2012

Seismicity and geometry of a 110‐km‐long blind thrust fault 1. The 1985 Kettleman Hills, California, earthquake

The August 4, 1985, Kettleman Hills earthquake was the third in a sequence of moderate shocks to occur beneath the northern half of a 110‐km‐long fold chain bounding the eastern California Coast Ranges. The 1982 M W =5.4 New Idria, 1983 M W =6.5 Coalinga, and 1985 M W =6.1 Kettleman Hills events define a southward progression of seismic activity beneath the fold. We use teleseismic waveforms, geodetic modeling, hypocenters relocated in a three‐dimensional velocity model, and subsurface structural data to investigate the Kettleman Hills earthquake. The main shock results from motion on a shallowly dipping thrust fault buried at ∼10 km depth. Aftershocks and coseismic fault slip extend 20 km along the fold axis, nearly the full extent of the Kettleman Hills North Dome anticline. Aftershocks occur primarily several kilometers in front of the fault tip and in the core of the anticline. The main shock and several foreshocks occurred at a 2‐km right step in the Quaternary fold axis, which also corresponds to the southern end of the 1983 Coalinga and northern end of the 1985 Kettleman Hills aftershock zones. From this we infer that the step in the fold is caused by an offset or tear in the underlying fault. The scalar seismic moment is 1.6 × 10 18 N m, consistent with the geodetic deformation, and the duration of rupture is 16 s, 3–4 times greater than for the average earthquake with this scalar moment. The slow rate of moment release provides an explanation for the low level of ground shaking and low local magnitude reported for the event. The peak of the geodetic uplift is located 5 km perpendicular to the Quaternary fold axis. We argue that the fault is propagating northeast into the undeformed San Joaquin Valley sediments and that the overlying fold is growing at about 0.5 mm/yr.

California↗
USGS2012

Application of continuum models to deformation of the Aleutian Island Arc

Continuum models were constructed to describe large‐scale deformation of the Aleutian Island Arc over the past 5 m.y. These models consider the island arc as a continuum in the horizontal plane with the velocity boundary condition at the Pacific edge stated as a fraction of Pacific plate convergence transferred to the arc. First, a simple model of uniformly distributed strain is formulated to illustrate the mechanics of continuous deformation. Lineaments along the arc massif rotated about a vertical axis are matched by small‐element rotation calculated from the model. However, this model does not predict across‐arc variations in deformation and produces an unrealistic amount of crustal thickening after 5 m.y. A physically more meaningful model of deformation is the thin viscous sheet model based on averages of stress and rheology throughout the lithosphere. The amount of motion transferred from the Pacific plate to the arc is constrained by the rotated lineaments, while the effective stress‐strain exponent ( n ) and the ability the lithosphere has to sustain crustal thickness contrasts (the Argand number) are independent variables. Primarily, bathymetry, earthquake focal mechanisms, and styles of faulting are used to evaluate the models. The preferred model is one where the amount of motion transferred from the Pacific plate is greater in an arc‐parallel direction than in an arc‐normal direction, producing stresses consistent with strike‐slip faulting at the far western end of the arc and tensional stresses consistent with transverse normal faulting elsewhere in the arc massif. This model agrees with observations of slip vectors by Ekström and Engdahl (1989), who conclude that a portion of the arc‐parallel component of relative plate motion is taken up in the overriding plate. This model implies that compressive stress transferred to the arc is small in comparison to along‐arc shear stress and that stresses conducive to strikeslip faulting are prevalent throughout the arc.

Alaska↗
USGS2012

Shear stress and bed roughness estimates for combined wave and current flows over a rippled bed

High‐quality bottom boundary layer measurements and bottom photographs were obtained over a sand substrate during a 10‐day deployment of the GEOPROBE tripod at an inner shelf (35‐m water depth) location off northern California. The seafloor surrounding the tripod was composed of well‐sorted medium‐grained (mean diameter, 0.25 mm) sand which was formed into symmetrical wave ripples with heights of 3–4 cm and wavelengths of 22–30 cm. Mean velocity profiles in the region from 23 cm to 102 cm above the rippled bed were highly logarithmic ( R > 0.95) approximately 30% of the time. Nineteen profiles exhibiting R > 0.997 were analyzed to obtain the shear velocity ( U * c ) and roughness length (Z 0 c ) for the mean current. The near‐bottom flow field was composed of quasi‐steady currents (up to 12 cm s −1 at z = 102 cm) and wave‐induced, oscillatory currents (up to 14 cm s −1 ). The data‐derived estimates of U * c and Z 0 c were 0.3–0.93 cm s −1 and 0.82–1.5 cm, respectively. The mean shear estimates are 50–100% larger than those predicted using a drag coefficient ( C D ) of 3 × 10 −3 that is typical for rough boundaries, and the roughness lengths are up to an order of magnitude larger than the maximum expected values based on the observed wave‐rippled bottom. These results indicate the importance of the combined flow turbulent interaction in producing a large apparent Z 0 c . However, comparison of the shear and roughness estimates derived from the velocity profile analysis to predictions made by the combined flow model of Grant and Madsen (1979) show that direct application of the wave‐ripple roughness equation of Grant and Madsen (1982) yields large overestimates of z 0 c and U * c . Selecting the physical roughness length k b (=30 z 0 ) that produced the best agreement with the data resulted in z 0 values ranging from 0.03 to 0.43 cm. Moreover, a direct correlation exists between these physical roughness estimates and the angle (θ cr ) formed by the mean current and the trend of the wave ripple crests. A simple linear relationship between k b and θ cr is suggested by our limited data set.

California↗
USGS2012

A model for the generation of two-dimensional surf beat

A finite difference model predicting group-forced long waves in the nearshore is constructed with two interacting parts: an incident wave model providing time-varying radiation stress gradients across the nearshore, and a long-wave model which solves the equations of motion for the forcing imposed by the incident waves. Both shallow water group-bound long waves and long waves generated by a time-varying breakpoint are simulated. Model-generated time series are used to calculate the cross correlation between wave groups and long waves through the surf zone. The cross-correlation signal first observed by Tucker (1950) is well predicted. For the first time, this signal is decomposed into the contributions from the two mechanisms of leaky mode forcing. Results show that the cross-correlation signal can be explained by bound long waves which are amplified, though strongly modified, through the surf zone before reflection from the shoreline. The breakpoint-forced long waves are added to the bound long waves at a phase of pi/2 and are a secondary contribution owing to their relatively small size.

Journal of Geophysical Research C: Oceans↗
USGS2012

Use of weather types to disaggregate general circulation model predictions

General circulation models (GCMs) simulate climatic conditions with a grid cell resolution on the order of 100,000 km 2 . This resolution is inadequate to assess the effects of climatic change on water resources at a regional scale. A method has been developed that uses weather-type analysis as a tool to spatially disaggregate GCM predictions to make them useful for water resource studies. The method has been applied to the Delaware River basin to predict the effects of doubling atmospheric carbon dioxide on precipitation patterns in the region. An application of the technique to the Delaware River basin indicates that future climatic conditions will show minimal changes in weather-type frequency, implying that air circulation patterns will remain unchanged. Results of this study indicate that changes in regional precipitation patterns under a doubling of atmospheric carbon dioxide will be a result of within-type changes in weather characteristics.

Journal of Geophysical Research↗
USGS2012

Sensitivity of ground - water recharge estimates to climate variability and change, Columbia Plateau, Washington

The sensitivity of groundwater recharge estimates was investigated for the semiarid Ellensburg basin, located on the Columbia Plateau, Washington, to historic and projected climatic regimes. Recharge was estimated for predevelopment and current (1980s) land use conditions using a daily energy-soil-water balance model. A synthetic daily weather generator was used to simulate lengthy sequences with parameters estimated from subsets of the historical record that were unusually wet and unusually dry. Comparison of recharge estimates corresponding to relatively wet and dry periods showed that recharge for predevelopment land use varies considerably within the range of climatic conditions observed in the 87-year historical observation period. Recharge variations for present land use conditions were less sensitive to the same range of historical climatic conditions because of irrigation. The estimated recharge based on the 87-year historical climatology was compared with adjustments to the historical precipitation and temperature records for the same record to reflect CO 2 -doubling climates as projected by general circulation models (GCMs). Two GCM scenarios were considered: an average of conditions for three different GCMs with CO 2 doubling, and a most severe &ldquo;maximum&rdquo; case. For the average GCM scenario, predevelopment recharge increased, and current recharge decreased. Also considered was the sensitivity of recharge to the variability of climate within the historical and adjusted historical records. Predevelopment and current recharge were less and more sensitive, respectively, to the climate variability for the average GCM scenario as compared to the variability within the historical record. For the maximum GCM scenario, recharge for both predevelopment and current land use decreased, and the sensitivity to the CO 2 -related climate change was larger than sensitivity to the variability in the historical and adjusted historical climate records.

Journal of Geophysical Research D: Atmospheres↗
USGS2012

Reconciliation of stress and structural histories of the Tharsis region of Mars

We present a new compilation of the structural and stratigraphic evolution of the Tharsis region of Mars that incorporates recent advances in understanding its stratigraphy, and we introduce a lithospheric deformation model that can account for the observations. The first period in the formation of Tharsis occurred in Late Noachian/Early Hesperian time with the deposition of volcanic plains materials throughout the surrounding highlands (e.g., Lunae Planum) and on the Tharsis rise (which includes the giant volcanoes and surrounding, elevated lava fields). Extensive radial normal faulting occurred on the rise, locally extending outward at Valles Marineris and Tempe Terra, and concentric wrinkle ridges formed along the edge of the rise. This regional deformation appears to have been modulated by a global compressional stress field due to rapid planetary cooling and contraction. The second period occurred during the Late Hesperian/Amazonian with Tharsis volcanism centered on the rise and radial extensional deformation that extended from the center of the rise for thousands of kilometers. We propose a model in which the lithosphere beneath Tharsis consists of a thin elastic crustal cap on the rise that is mechanically detached from the strong upper mantle by a volcanically thickened, hot, weak lower crust. These layers merge into a single cooler, strong lithospheric layer around the edges of the rise. This model is capable of generating large extensional hoop stresses throughout much of the western hemisphere, in agreement with observations. The tectonic interpretation of the stresses predicted by this model requires the reconciliation of extensional strain within narrow grabens and compressional strain within wrinkle ridges with (1) processes in the deeper lithosphere, (2) the sparsity of strike-slip faults, and (3) other global or locally important stress fields. Stresses predicted by global models affect the entire thickness of the lithosphere, and they can be reconciled with narrow, closely spaced grabens that accommodate large amounts of extensional strain in the upper few kilometers of the lithosphere if the grabens are underlain and kinematically linked with dikes or other tension cracks, such as hydrofractures. Deeper levels of the lithosphere can accommodate this strain by elastic expansion if grabens are spaced far apart (many tens to hundreds of kilometers). Mechanical considerations suggest that deformation on faults beneath wrinkle ridges could extend through a significant thickness of the brittle crust. A number of factors, including stresses generated by the addition of overburden, intrusion of dikes, weakness of geologic materials under extension, and the laterally constrained nature of a single lithospheric plate, may have inhibited the formation of strike-slip faults on Mars. Stresses generated from the removal of overburden could have augmented planetwide wrinkle ridge formation during the Late Noachian/Early Hesperian and in Kasei Valles and western Chryse Planitia during the Early Amazonian. The nonuniform distribution of tectonic features around Tharsis can be understood in terms of the concentration of regional stresses and strain near weaker volcanotectonic centers.

Journal of Geophysical Research↗