Recent crustal subsidence at Yellowstone Caldera,Wyoming

Following a period of net uplift at an average rate of 15±1 mm/year from 1923 to 1984, the east-central floor of Yellowstone Caldera stopped rising during 1984–1985 and then subsided 25±7 mm during 1985–1986 and an additional 35±7 mm during 1986–1987. The average horizontal strain rates in the northeast part of the caldera for the period from 1984 to 1987 were: ₁ = 0.10 ± 0.09 strain/year oriented N33° E±9° and ₂ = 0.20 ± 0.09 strain/year oriented N57° W±9° (extension reckoned positive). A best-fit elastic model of the 1985–1987 vertical and horizontal displacements in the eastern part of the caldera suggests deflation of a horizontal tabular body located 10±5 km beneath Le Hardys Rapids, i.e., within a deep hydrothermal system or within an underlying body of partly molten rhyolite. Two end-member models each explain most aspects of historical unrest at Yellowstone, including the recent reversal from uplift to subsidence. Both involve crystallization of an amount of rhyolitic magma that is compatible with the thermal energy requirements of Yellowstone's vigorous hydrothermal system. In the first model, injection of basalt near the base of the rhyolitic system is the primary cause of uplift. Higher in the magmatic system, rhyolite crystallizes and releases all of its magmatic volatiles into the shallow hydrothermal system. Uplift stops and subsidence starts whenever the supply rate of basalt is less than the subsidence rate produced by crystallization of rhyolite and associated fluid loss. In the second model, uplift is caused primarily by pressurization of the deep hydrothermal system by magmatic gas and brine that are released during crystallization of rhyolite and them trapped at lithostatic pressure beneath an impermeable self-sealed zone. Subsidence occurs during episodic hydrofracturing and injection of pore fluid from the deep lithostatic-pressure zone into a shallow hydrostatic-pressure zone. Heat input from basaltic intrusions is required to maintain Yellowstone's silicic magmatic system and shallow hydrothermal system over time scales longer than about 10⁵ years, but for the historical time period crystallization of rhyolite can account for most aspects of unrest at Yellowstone, including seismicity, uplift, subsidence, and hydrothermal activity.     

Sensitivity of glaciation to initial snow cover,CO2, snow albedo,and oceanic roughness in the NCAR CCM

We present results from numerical experiments made with a GCM, the NCAR CCM1, that were designed to estimate the annual balance between snow-fall accumulation and ablation for geographically important land regions for a variety of conditions. We also attempt to assess the reliability of these results by investigating model sensitivity to changes in prescribed physical parameters. Experiments were run with an initial imposition of 1 m of (midwinter) snowcover over all northern hemisphere land points. Over Alaska, western Canada, Siberia, and the Tibetan Plateau the model tended to retain this snow cover through the summer and in some cases increase its depth as well. We define these regions as glaciation sensitive and note some correspondence between them and source regions for the Pleistocene ice sheets. An experiment with greatly reduced CO₂ (100 ppm) showed a tendency towards spontaneous glaciation, i.e., the model remained snow-covered throughout the summer over the same geographic regions noted above. With 200 ppm CO₂ (roughly equal to values at the last glacial maximum), snow cover over these regions did not quite survive the summer on a consistent basis. Combining 200 ppm CO₂ and 1 m of initial northern hemisphere snow cover yielded glaciation-sensitive conditions, agreeing remarkably well with locations undergoing glaciation during the Pleistocene. To assess the reliability of these results, we have determined minimal model uncertainty by varying two of the empirical coefficients in the model within physically plausible ranges. In one case surface roughness of all ocean gridpoints was reduced by an order of magnitude, leading to local 10% reductions in precipitation (snowfall), a change hard to distinguish from inherent model variability. In the other case, the fraction of a land grid square assumed to be occupied by snow cover for albedo purposes was varied from one-half to unity. Large changes occurred in the degree of summer melting, and in some cases the sign of the net balance changed as fractional snow cover was changed. We conclude that the model may be able to reveal regions sensitive to glaciation, but that it cannot yield a reliable quantitative computation of the magnitude of the net snow accumulation that can be implicitly or explicitly integrated through time.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dilmenil     

Magmatism in the eastern Aleutian Arc: temporal characteristic of igneous activity on Akutan Island

Lavas from Akutan Island, located in the eastern Aleutian arc at the transition between continental and oceanic crust, show a gradual change in their petrologic and chemical characteristics over the last 4 million years. The oldest lavas exposed on the island, the Hot Springs Bay Volcanics (HSBV), range from magnesian basalt to dacite (45%–62% SiO₂). The most mafic basalts contain salitic clinopyroxene, Cr- and Al-rich spinel, and pargasitic amphibole suggesting that they were derived from relatively hydrous magmas at greater pressures than lavas from the younger Akutan Volcanics (AKV) and the modern volcano (MOD). AKV lavas also range between basalt and dacite (46%–63% SiO₂), but contain no hydrous phenocrysts and seem to have fractionated within a shallow level magma chamber. Lavas from the modern volcano are andesitic (52%–57% SiO₂) and have a mineral assemblage similar to that of AKV lavas of similar composition. With the exception of clinopyroxene and spinel in the most mafic lavas, the compositions of plagioclase (An_92?45), olivine (Fo_88?51), orthopyroxene (En_69?56), and titanomagnetite (15%–21% TiO₂) phenocrysts found in these lavas are within the range observed in lavas from other Aleutian volcanoes. Variations in the major element chemistry of the older lavas can be reproduced by fractional crystallization of the observed mineral assemblages, however closed system crystal fractionation models are inadequate to explain the trace element variations. During the last 4 million years, La/Yb ratios have decreased (6.5–3.3 for HSBV lavas and 2.9–1.9 for MOD lavas) whereas Ba/La ratios appear to have increased slightly (37–43 for HSBV and AKV, and 41–45 of MOD). The lower La/Yb ratios of MOD lavas correspond with lower total abundances of the REE and slightly higher Sr and Pb isotopic ratios. The increased⁸⁷Sr/⁸⁶Sr ratios and Pb isotopic ratios in the MOD lavas, the less enriched LREE, and the higher Ba/La ratios may result from partial melting of an arc source which has experienced previous melting events but has continued to be contaminated by a component from the subducting slab. It may also indicate a change in the degree of partial melting of the underlying mantle, which corresponds to a different percentage of a slab derived component being incorporated into the overlying mantle.     

The depth vs. depth plot: An aid to selecting piezometer installation depths in boreholes

Dick, R.C., 1976. The depth vs. depth plot: an aid to selecting piezometer installation depths in boreholes. Eng. geol., 10: 37–42.The location of the piezometric surface within a soil or rock mass is often the critical factor in stability analysis and in evaluations of reservoir watertightness. Sealed piezometers are now almost universally used to measure water pressure within the mass. The location of these piezometers must be decided immediately upon completion of each borehole. A common aid to making this decision is a plot of borehole depth and borehole water level as functions of time. This paper introduces the reader to another way to evaluate the same data, with the borehole water depth plotted as a function of borehole depth. Anomalies can be detected more easily, and piezometers installed within these zones. Examples are given.     

A maximal background spectrum of gravitational radiation

A maximal spectrum of gravitational radiation from sources outside our galaxy is calculated. The sources are galaxies, quasars and events that occur in the early history of the universe. The major contribution is from galaxies whose effect extends over the frequency region 10^–810⁺⁴Hz, peaking at 10^–110 Hz, with a spectral flux of 10 erg cm^–2, s^–1. The main processes of gravitational radiation in the galaxies are stellar collapse into a black hole and dying binary systems. In the region 10^–410⁴ Hz the background spectrum is well above the detection levels of currently proposed detectors. FromMinimal considerations of this spectrum it is determined that the density of gravitational radiation is 10^–39g cm^–3. This background spectrum is sensitive to galactic evolution and especially sensitive to the upper mass limits and mass distribution of stars in galactic models. Therefore, the spectrum could provide information about galactic evolution complementary to that obtained by electromagnetic investigations.     

Observations of cosmic gamma-ray bursts

The occasional occurrence of brief but intense bursts of cosmic gamma-rays was disclosed through a systematic search of data acquired from the Vela satellites. Confirmation of the nature of the events and additional detail of their characteristics has subsequently been provided by 15 other groups of experimenters with instruments on 13 spacecraft. Thirty-nine such events have been identified from data spanning a period of four and a half years. The record of intensity as a function of time varies considerably for different events, with total durations ranging from 0.1 to 60s. Time-integrated flux density ranges from 10^–6 to 10^–3 erg cm^–2. Spectral measurements have been accomplished by several groups of experimenters, showing a broad maximum in the energy distribution at about 150 keV. The distribution of source directions implies either near galactic or extragalactic locations. The existing data are not sufficient to distinguish between the various models proposed to explain the phenomenon; no model is completely consistent with all observed characteristics.Work performed under the auspices of the United States Energy Research and Development Administration.Paper presented at the COSPAR Symposium on Fast Transients in X-and Gamma-Rays, held at Varna, Bulgaria, 29–31 May, 1975.     

Transport and emplacement of crater and basin deposits

Material is ejected from impact craters in ballastic trajectories; it impacts first near the crater rim and then at progressively greater ranges. Ejecta from craters smaller than approximately 1 km is laid predominantly on top of the surrounding surface. With increasing crater size, however, more and more surrounding surface will be penetrated by secondary cratering action and these preexisting materials will be mixed with primary crater ejecta. Ejecta from large craters and especially basin forming events not only excavate preexisting, local materials, but also are capable of moving large amounts of material away from the crater. Thus mixing and lateral transport give rise to continuous deposits that contain materials from within and outside the primary crater. As a consequence ejecta of basins and large highland craters have eroded and mixed highland materials throughout geologic time and deposited them in depressions inside and between older crater structures.Because lunar mare surfaces contain few large craters, the mare regolith is built up by successive layers of predominantly primary ejecta. In contrast, the lunar highlands are dominated by the effects of large scale craters formed early in lunar history. These effects lead to thick fragmental deposits which are a mixture of primary crater material and local components. These deposits may also properly be named regolith though the term has been traditionally applied only to the relatively thin fine grained surficial deposit on mare and highland terranes generated during the past few billion year. We believe that the surficial highland regolith - generated over long periods of time - rests on massive fragmental units that have been produced during the early lunar history.