Present limitations of accurate satellite Doppler positioning for tectonics. An example: Djibouti

One of the possibilities of the Doppler positioning from satellite is to provide geodetic measurements continuous in time without the need for reference stations. If measurements of sufficient accuracy can be achieved they may be used to study local surface displacements in relation to tectonic activity. A Doppler receiver of the MEDOC network is located near Djibouti in the Ghoubhat-Asal rift region which corresponds to the accreting plate boundary between the Arabian and African plates. In November 1978, a seismic and volcanic crisis occured in this area. Surface geodetic measurements and levellings performed in 1973 and 1978–79, just after the crisis, reveal a 60–80 cm sinking of the graben floor and a lateral extension of about 2 meters. Here we analyse Doppler measurements for the period January 1977 to November 1980. Point positions are computed for 7 to 10 day intervals using a precise ephemeris, and a moving window analysis is applied to the data. An apparent 2 meter uplift preceding the November 1978 seismic crisis is detected at Djibouti, whereas no similar phenomenon is observed at the two closest stations, Pretoria and Uccle-Brussels. However, field observations rule out a tectonic origin for this uplift. In Djibouti, the correlation between the apparent vertical station position and the electron density in the ionospheric F-layer reveals that a bias may be induced by the third order term of the ionospheric refractive index not previously taken into account, or more probably by the ray curvature through the ionosphere. This bias is particularly strong for our data set, from a station located close to the magnetic equator, and related to a period of rapid increase in the solar activity. Although our analysis fails to detect any tectonic displacement related to the rift seismic activity, it shows that similar tectonic studies by the Doppler method will be possible once the residual ionospheric errors are removed, for example by use of higher radio frequencies.     

Compositional and thermal convection in magma chambers

Magma chambers cool and crystallize at a rate determined by the heat flux from the chamber. The heat is lost predominantly through the roof, whereas crystallization takes place mainly at the floor. Both processes provide destabilizing buoyancy fluxes which drive highly unsteady, chaotic convection in the magma. Even at the lowest cooling rates the thermal Rayleigh number Ra is found to be extremely large for both mafic and granitic magmas. Since the compositional and thermal buoyancy fluxes are directly related it can be shown that the compositional Rayleigh number Rs (and therefore a total Rayleigh number) is very much greater than Ra. In the case of basaltic melt crystallizing olivine Rs is up to 10⁶ times greater than Ra. However compositional and thermal buoyancy fluxes are roughly equal. Therefore thermal and compositional density gradients contribute equally to convection velocities in the interior of the magma. Effects of thermal buoyancy generated by latent heat release at the floor are included.The latent heat boundary layer at the floor of a basaltic chamber is shown to be of the order of 1 m thick with very low thermal gradients whereas the compositional boundary layer is about 1 cm thick with large compositional gradients. As a consequence, the variation in the degree of supercooling in front of the crystal-liquid interface is dominated by compositional effects. The habit and composition of the growing crystals is also controlled by the nature of the compositional boundary layer. Elongate crystals are predicted to form when the thickness of the compositional boundary layer is small compared with the crystal size (as in laboratory experiments with aqueous solutions). In contrast, equant crystals form when the boundary layer is thicker than the crystals (as in most magma chambers). Instability of the boundary layer in the latter case gives rise to zoning within crystals. Diffusion of compatible trace elements through the boundary layer can also explain an inverse correlation, observed in layered intrusions, between Ni concentration in olivine and the proportion of Ni-bearing phases in the crystallizing assemblage.     

A skeletal assemblage classification system for non-tropical carbonate deposits based on New Zealand Cenozoic limestones

Cenozoic limestones in New Zealand are mainly skeletal grainstones and packstones formed under non-tropical climatic conditions in open marine shelf or ramp environments. Following petrographic analysis of the nature and abundance of the skeletal components in nearly 500 samples of these limestones, a complete linkage cluster analysis identified seven major skeletal assemblages that may be regarded as subdivisions of the single foramol skeletal association defined by Lees and Buller (1972) for temperate-region carbonate deposits. The seven assemblages are given contracted names, as follows: (a) BARNAMOL = barnacle/bivalve-dominated; (b) BIMOL = bivalve-dominated; (c) BRYOMOL = bryozoan/bivalve-dominated; (d) ECHINOFOR = echinoderm/benthic foraminiferal-dominated; (e) NANNOFOR = nannofossil/planktonic foraminiferal-dominated; (f) RHODALGAL = calcareous red algal-dominated; and (g) RHODECHFOR = calcareous red algal/echinoderm/benthic foraminiferal-dominated. A composite triangular classification diagram has been devised for naming the skeletal assemblage of an unknown sample on the basis of its three main skeletal components. The diagram successfully characterises more than 85% of the New Zealand Cenozoic limestone samples and also appears to be applicable for the skeletal assemblage designation of many overseas examples of non-tropical carbonate deposits. Limitations relate mainly to locally common skeletal types (e.g. serpulids, brachiopods) that are presently not incorporated into the New Zealand-based scheme. The general ecological preferences of the main skeletal contributors in each of the seven skeletal assemblages form a basis for relating the assemblages to broad shelf habitats. Consequently, as well as the benefits of providing a more consistent skeletal assemblage terminology for comparative studies between different workers, the scheme can assist with the paleoenvironmental interpretation of non-tropical skeletal carbonate facies.     

Compositions of group IVB iron meteorites and their parent melt

The concentrations of P, V, Cr, Fe, Co, Ni, Cu, Ga, Ge, As, Mo, Ru, Rh, Pd, W, Re, Os, Ir, Pt, and Au in the group IVB iron meteorites Cape of Good Hope, Hoba, Skookum, Santa Clara, Tawallah Valley, Tlacotepec, and Warburton Range have been measured by laser ablation inductively coupled plasma mass spectrometry. The data were fitted to a model of fractional crystallization of the IVB parent body core, from which the composition of the parent melt and metal/melt distribution coefficients for each element in the system were determined, for a chosen value of D(Ni). Relative to Ni and chondritic abundances, the parent melt was enriched in refractory siderophiles, with greatest enrichment of 5× chondritic in the most refractory elements, and was strongly volatile-depleted, down to 0.00014× chondritic in Ge. Comparison to an equilibrium condensation sequence from a gas of solar composition indicates that no single temperature satisfactorily explains the volatility trend in the IVB parent melt; a small (<1%) complement of ultrarefractory components added to metal that is volatile-depleted but otherwise has nearly chondritic abundances (for Fe, Co and Ni) best explains the volatility trend. In addition to this volatility processing, which probably occurred in a nebular setting, there was substantial oxidation of the metal in the IVB parent body, leading to loss of Fe and other moderately siderophile elements such as Cr, Ga, and W, and producing the high Ni contents that are observed in the IVB irons. By assuming that the entire IVB parent body underwent a similar chemical history as its core, the composition of the silicate that is complementary to the IVB parent melt was also estimated, and appears to be similar to that of the angrite parent.     

First seismic imaging results of tectonically complex structures at shallow depths beneath the northwest Canterbury Plains,New Zealand

Seismogenic structures underlie many regions of the vast Canterbury Plains on the South Island of New Zealand. Most of these structures are hidden beneath a layer of rapidly deposited Late Pleistocene sediments, the youth and thickness of which make the general application of conventional paleoseismological studies impractical. In an attempt to improve our understanding of potentially active structures in this region, we have acquired, processed and interpreted shallow seismic reflection data across the northwest Canterbury Plains. To separate the useful reflected signals from unusually high amplitude ambient and source-generated noise, we subjected the data to a specially tailored processing scheme that included time- and space-variant spectral balancing, custom static corrections and mutes, F-X deconvolution, DMO corrections and finite-difference migration. The final stacked and migrated seismic sections supply high-resolution images of the basement and overlying layered Cretaceous- to Quaternary-age supracrustal rocks that have been complexly faulted and folded. At one location, the uppermost Late Pleistocene layers appear to have been gently buckled.     

Radon transfer from groundwater used in showers to indoor air

Estimation of Rn transfer from water to indoor air based on multi-day measurements may underestimate alpha exposure that occurs at short time scales in confined spaces, such as from showering, in houses with high Rn activities in the water supply. In order to examine one such incremental increase in exposure, variations in Rn in water and indoor air in 18 houses with private wells in western North Carolina (USA) were investigated. Radon in well water ranged from 158 to 811 Bq L⁻¹ (median 239 Bq L⁻¹). After 20-min showers in bathrooms with closed doors, peak Rn in air increases (above background) ranged from 71 to 4420 Bq m⁻³ (median 1170 Bq m⁻³). Calculated transfer coefficients at the scale of a 40-min closed bathroom (20-min shower plus 20 min post-shower) are described by a lognormal distribution whose geometric mean exceeds the widely-used ∼10⁻⁴ whole-house transfer coefficient by about one order of magnitude. As short-lived decay products grow from shower-derived Rn, short-term alpha energy exposure occurs in bathrooms in addition to the exposure caused by Rn mixed throughout the volume of the house. Due to the increasing ratio of Rn decay products to Rn, alpha energy exposure is greatest several minutes after the shower is turned off. For a 7.2-min shower with 10 min of additional exposure before opening the door, a geometric mean 5.6% increase in exposure over the ∼10⁻⁴ whole-house transfer coefficient derived from longer measurement periods was estimated. In addition to Rn activity in water, short-term shower exposure to Rn progeny depends on exposure time, ventilation, attachment and deposition, among other variable factors that characterize individual houses and residents.     

A study of the diurnal variation in the thermosphere as derived by satellite drag

An analysis was performed on 29,574 densities derived from the drag of 10 satellites to determine simultaneously the parameters of the solar-activity effect in the thermosphere on the one hand, and the amplitude and shape of the diurnal-variation curve on the other. This paper reports on the study of the diurnal variation only.Although a considerable amount of smoothing is inherent in the drag method, it seemed useful to see whether we could detect any change in the shape of the diurnal-variation curve with height, latitude and solar activity. None was detected: the curve remains remarkably stable, symmetric, with a maximum at 14 hr 20 min L.S.T. and a minimum at 2 hr 20 min L.S.T. A systematic variation of the temperature range with height is observed when static models are used to derive it.