Uncertainties in field-line tracing in the magnetosphere. Part II: the complete internal geomagnetic field

The discussion in the preceding paper is restricted to the uncertainties in magnetic-field-iine tracing in the magnetosphere resulting from published standard errors in the spherical harmonic coefficients that define the axisymmetric part of the internal geomagnetic field (i.e. g_n⁰ ± g_n⁰). Numerical estimates of these uncertainties based on an analytic equation for axisymmetric field lines are in excellent agreement with independent computational estimates based on stepwise numerical integration along magnetic field lines. This comparison confirms the accuracy of the computer program used in the present paper to estimate the uncertainties in magnetic-field-line tracing that arise from published standard errors in the full set of spherical harmonic coefficients, which define the complete (non-axisymmetric) internal geomagnetic field (i.e. g_n^m ± g_n^m and h_n^m ± h_n^m). An algorithm is formulated that greatly reduces the computing time required to estimate these uncertainties in magnetic-field-line tracing. The validity of this algorithm is checked numerically for both the axisymmetric part of the internal geomagnetic field in the general case (1 n 10) and the complete internal geomagnetic field in a restrictive case (0 m n, 1 n 3). On this basis it is assumed that the algorithm can be used with confidence in those cases for which the computing time would otherwise be prohibitively long. For the complete internal geomagnetic field, the maximum characteristic uncertainty in the geocentric distance of a field line that crosses the geomagnetic equator at a nominal dipolar distance of 2 R_E is typically 100 km. The corresponding characteristic uncertainty for a field line that crosses the geomagnetic equator at a nominal dipolar distance of 6 R_E is typically 500 km. Histograms and scatter plots showing the characteristic uncertainties associated with magnetic-field-line tracing in the magnetosphere are presented for a range of illustrative examples. Finally, estimates are given for the maximum uncertainties in the locations of the conjugate points of selected geophysical observatories. Numerical estimates of the uncertainties in magnetic-field-line tracing in the magnetosphere, including the associated uncertainties in thelocations of the conjugate points of geophysical observatories, should be regarded as first approximations in the sense that these estimates are only as accurate as the published standard errors in the full set of spherical harmomic coefficients. As in the preceding paper, howerver, all computational techniques developed in this paper can be used to derive more realistic estimates of the uncertainties in magnetic-field-line tracing in the magnetosphere, following further progress in the determination of more accurate standard errors in the spherical harmonic coefficients.Also Visiting Reader in Physics, University of Sussex, Palmer, Brighton, BN1 9QH, UK     

Sulfur isotope fractionation during bacterial sulfate reduction in organic-rich sediments

Isotope fractionation during sulfate reduction by natural populations of sulfate-reducing bacteria was investigated in the cyanobacterial microbial mats of Solar Lake, Sinai and the sediments of Logten Lagoon sulfuretum, Denmark. Fractionation was measured at different sediment depths, sulfate concentrations, and incubation temperatures. Rates of sulfate reduction varied between 0.1 and 37 micromoles cm-3 d-1, with the highest rates among the highest ever reported from natural sediments. The depletion of 34S during dissimilatory sulfate reduction ranged from 16% to 42%, with the largest 34S-depletions associated with the lowest rates of sulfate reduction and the lowest 34S-depletions with the highest rates. However, at high sulfate reduction rates (>10 micromoles cm-3 d-1) the lowest fractionation was 20% independent of the rates. Overall, there was a similarity between the fractionation obtained by the natural populations of sulfate reducers and previous measurements from pure cultures. This was somewhat surprising given the extremely high rates of sulfate reduction in the experiments. Our results are explained if we conclude that the fractionation was mainly controlled by the specific rate of sulfate reduction (mass cell-1 time-1) and not by the absolute rate (mass volume-1 time-1). Sedimentary sulfides (mainly FeS2) were on average 40% depleted in 34S compared to seawater sulfate. This amount of depletion was more than could be explained by the isotopic fractionations that we measured during bacterial sulfate reduction. Therefore, additional processes contributing to the fractionation of sulfur isotopes in the sediments are indicated. From both Solar Lake and Logten Lagoon we were able to enrich cultures of elemental sulfur-disproportionating bacteria. We suggest that isotope fractionation accompanying elemental sulfur disproportionation contributes to the 34S depletion of sedimentary sulfides at our study sites.     

Investigation by backscatter radar of artificial irregularities produced in ionospheric plasma heating experiments

The artificial ionospheric turbulence occurs in the ionosphere illuminated by high power HF radio waves. There are a lot of irregularities stretched along the geomagnetic field in this region. The investigation of the artificially disturbed ionospheric region is based on the reception of back scattered signals (BSS) which permits the basic parameters of this region to be estimated and its inhomogeneous structure to be described.Experiments were carried out using ‘Sura’ heating facility in the frequency range of ? = 4.7–9 MHz (ordinary mode) with the effective radiated power P_eff = 50–70 MW beamed vertically upwards. The most important dynamic and statistical BSS characteristics (the built-up time, the relaxation and autocorrelation times, the BSS amplitude spatial correlation function and power spectra) were measured using probe waves in the frequency range of that made it possible to obtain the basic parameters of the artificial irregularities. The model representation of a disturbed region in a form of a periodic structure gives a possibity to evaluate the scale of the structure, the whole size of the disturbance and its power and to calculate the main BSS characteristics.     

The mass accommodation coefficient of ammonia on water

The mass accommodation coefficient of ammonia gas on water has been determined by measuring the absorption rate of 50–200 ppm NH₃ in one atm of air or helium into a liquid jet of 97 µm diameter as function of the exposed jet length, and comparing the results with numerical simulations which treat as the only free parameter. The model considers in detail transport of NH₃ by molecular diffusion, penetration of the gas/water interface, hydrolysis in the acidified water, and transport of the solutes from the surface into the jet. A correction is applied for the time evolution of the jet surface speed, using literature data on the fluid mechanics of liquid jets. The result of nine sets of independent measurements is

Stable boundary layers: observations,models and variability part II: Data analysis and averaging effects

As argued in Part I (Derbyshire, 1995), variability is a key issue in stable boundary layers, and differences in variability between observations and idealized models may imply sytematic biases. Here we discuss how data analysis can be geared to allow for variability and thus consistency with models. Instrumental errors, smoothing methods and vertical discretization are considered. We then show how statistical averaging broadly improves the agreement of heterogeneous results in Part I with the Brost-Wyngaard closure. Recommendations are made for the information needed to analyze apparent differences between homogeneous and heterogeneous stable boundary layers.Part of UK Meteorological Office Atmospheric Processes Research Division     

Mechanisms of shrubland expansion: land use,climate or CO2?

Encroachment of trees and shrubs into grasslands and the thicketization of savannas has occurred worldwide over the past century. These changes in vegetation structure are potentially relevant to climatic change as they may be indicative of historical shifts in climate and as they may influence biophysical aspects of land surface-atmosphere interactions and alter carbon and nitrogen cycles. Traditional explanations offered to account for the historic displacement of grasses by woody plants in many arid and semi-arid ecosystems have centered around changes in climatic, livestock grazing and fire regimes. More recently, it has been suggested that the increase in atmospheric CO₂ since the industrial revolution has been the driving force. In this paper we evaluate the CO₂ enrichment hypotheses and argue that historic, positive correlations between woody plant expansion and atmospheric CO₂ are not cause and effect.Please direct all correspondence to the senior author.     

A fast response,open path,infrared hygrometer,using a semiconductor source

Investigations were carried out as to the feasibility of using a semiconductor source in the design of a new rapid response, open-path hygrometer. A single-beam instrument was constructed employing an infrared light emitting diode (LED) as a source instead of the usual high energy, wideband filament. The spectral emission envelope encompassed the 1.87 m water absorption band. Electronic modulation and thermoelectrical cooling of the diode eliminated the conventional chopper wheel and stabilized the peak wavelength emission. Path length was 200 mm. Over a water vapour concentration range of 0–16 g m^–3, absorption varied by 2% in a linear fashion. At 10 Hz, the noise level was 0.1 g m^–3 rms. Hygrometer resolution and stability are constrained by the detector noise level, the small source emission in the absorption band and low frequency drift in the optical filter. Despite these problems, the new instrument showed comparable performance characteristics to a commercial Lyman- hygrometer. Latent heat fluxes measured with both instruments and a Kaijo-Denki, 3-D sonic anemometer agreed to within 4% over a range 0–350 W m^–2. Further improvements in performance can be anticipated with advances in detector and LED technology.     

Geology of the volcanic-hosted Brockman rare-metals deposit,Halls Creek Mobile Zone,northwest Australia. l. Volcanic environment,geochronology and petrography of the Brockman volcanics

Summary Rare-metals mineralization of the Brockman deposit (Halls Creek Mobile Zone, NW Australia) is hosted in a fluorite-bearing, rhyolitic volcaniclastic unit informally termed the Niobium Tuff. The Tuff, more correctly described as a tuffaceous volcaniclastic deposit, is the lowermost unit of a sequence of trachyte-to-rhyolite lavas, trachyandesite subvolcanic rocks, and volcaniclastic units of the Brockman volcanics located within the Halls Creek Group, a thick, early Proterozoic volcano-sedimentary sequence. High precision SHRIMP ion-microprobe zircon dating of the Niobium Tuff gives an eruption age of 1870 ± 4 Ma. Regional geochronological constraints indicate the Niobium Tuff was deposited about 15 Ma before major orogenic activity affected the area. Despite folding, faulting and low-grade metamorphism, the Brockman volcanics show excellent preservation of primary volcanological features, including pillow-lavas and pillow-breccias, that suggest a dominantly subaqueous, below-wave-base emplacement environment. The style of eruption products and magma volume constraints suggest the trachyte-dominated volcanics were erupted from a small shield volcanic complex probably in a rift-related basin in a shallow-marine setting. The tectonic setting was intraplate but differs from most Tertiary to recent trachyte volcanic complexes which are largely subaerial, are built on relatively thick continental crust, and show no post-eruptive orogenic history. Brockman-style rare-metal deposits are characterized by preservation of subaqueous volcanics beneath a thick sedimentary sequence, eruption of early incompatible-element enriched products followed by less differentiated magmas, and fine-grained mineralogy influenced by alteration processes. Prospects exist for discovery of analogous deposits, particularly in early Proterozoic mobile belts and Tertiary intraplate shield volcanic provinces.

---

Geologie der an Vulkanite gebundenen Seltene-Metalle-Lagerstätte Brockman, Halls Creek Mobile Zone, Nordwest-Australien. I. Vulkanologie, Geochronologie und Petrographie der Brockman-Vulkanite

Zusammenfassung Selten-Metall-Vererzung der Brockman-Lagerstätten (Halls Creek Mobile Zone, NW Australien) sitzt in einer Fluorit-führenden, rhyolitischen vulkanoklastischen Einheit auf, die informell als der Niob-Tuff bezeichnet wird. Es handelt sich hier um ein vulkanoklastisches Tuff-Sediment; dieses ist die unterste Einheit einer Abfolge vor Trachyt bis Rhyolit-Laven, trachyandesitischen Subvulkaniten and vulkanoklastischen Einheiten der Brockman Vulkanite innerhalb der Halls Creek Gruppe, einer mächtigen frühproterozoischen vulkanosedimentdren Abfolge. Prüzisionsdatierungen des Niob Tuffs mit der SHRIMP Ionen-Mikrosonde ergeben ein Eruptionsalter von 1870 +-4 Millionen Jahren. Regionale geochronologische Zusammenhänge zeigen, daß der Niob Tuff etwa 15 Millionen Jahre vor einer größeren Orogenese, die das Gebiet betroffen hat, abgelagert wurde. Trotz Faltung, Bruchtektonik and niedriggradiger Metamorphose zeigen die Brockman-Vulkanite einen hervorragenden Erhaltungszustand primärer vulkanologischer Erscheinungen. Diese urnfassen auch Kissenlaven und KissenBreckzien, die eine vorwiegend subaquatische Ablagerung in ruhigem Wasser erkennen lassen. Die Art der Eruptionsprodukte and das Magmavolumen zeigen, daß die Trachytdominierten Vulkanite von einem kleinen Schildvulkan stammer, wahrscheinlich in einem Becken in einer Rift-Situation im seichten marinen Milieu. Die tektonische Situation war intraplate, aber unterscheidet sich von den meisten tertiären bis rezenten trachytischen Vulkan-Komplexen, die hauptsächlich subaerisch sind, auf einer relativ mächtigen kontinentalen Kruste aufsitzen, and keine post-eruptive orogene Entwicklung zeigen. Seltene-Metalle-Lagerstatten des Brockman-Typs rind durch die Erhaltung subaquatischer-Vulkanite unterhalb einer machtigen sedimentdren Abfolge gekennzeichnet; welters durch frühe Eruptionsprodukte, die an inkompatiblen Elementen angereichert sind, auf die dann weniger differenzierte Magmen folgten, und schließlich durch einen feinkörnigen Mineralbestand, der vor Umwandlungsprozessen betroffen war. Es besteht die Möglichkeit der Entdeckung analoger Lagerstätten, besonders in frühproterozoischen mobilen Gürteln und in tertidren intraplate Schildvulkan Provinzen.

---

---

With 5 Figures

---

Deceased     

Textures exhibited by feldspars in the Giant Mottled Anorthosite (GMA) of the bastard unit in the Upper Critical Zone,Western Bushveld Complex

Summary The Bastard Unit, overlying the Merensky Unit, is the uppermost and most complete of the cyclic units in the Upper Critical Zone and its upper contact is taken as the boundary between the Critical and Main Zones (as recommended by SACS, 1980). The Unit is different from underlying units for three reasons: it does not have a well-defined top contact (with an overlying unit); it is much thicker than the Merensky and Footwall Units; and it has an unusually thick sequence of mottled anorthosite at its top. This distinctive composite anorthositic sequence is termed the Giant Mottled Anorthosite (GMA) which is predominantly made up of poikilitic anorthosite. The GMA can broadly be subdivided into three distinct parts; the Lower Giant Mottled Anorthosite (LGMA), Giant Mottled Middling (GMM) and the Upper Giant Mottled Anorthosite (UGMA).Petrographic examination of 82 samples taken in the GMA from 11 borehole profiles around the western limb of the Complex revealed that the LGMA and GMM are essentially adcumulates while the UGMA is orthocumulate in character. By its very nature the dominant phase throughout the GMA is cumulus plagioclase feldspar ( 85%) with the balance being made up of intercumulus orthopyroxene and clinopyroxene and minor biofite — all of which constitute the darker mottles. Intercumulus inverted pigeonite occurs in the upper part of the GMM, throughout the UGMA and for c. 30 m into the Main Zone. K-feldspar also joins the paragenesis within the UGMA and occurs both as intercumulus orthoclase and antiperthite. Cumulus plagioclase grains, showing varying degrees of complex oscillatory-zonal growth are common in the upper part of the GMA. These grains are interpreted as records of growth episodes during earlier influxes of magma which were responsible for the development of the underlying cyclic units, but escaped incorporation into the cumulate pile.

---

Feldspat-Texturen im Giant Mottled Anorthosite der Bastard-Einheit in der oberen kritischen Zone des westlichen Bushveld-Komplexes

Zusammenfassung Die Bastard-Einheit, die die Merensky-Einheit überlagert, ist die oberste und am besten entwickelte der zyklischen Einheiten in der oberen kritischen Zone; ihr oberer Kontakt gilt als die Grenze zwischen der Kritischen und der Main-Zone (entsprechend den Empfehlungen von SACS, 1980). Die Einheit unterscheidet sich aus drei Gründen von den Einheiten im liegenden: sie hat keinen gut definierten oberen Kontakt (mit einer darüberliegenden Einheit), sie ist wesentlich mächtiger als die Merensky und die Footwall-Einheiten, und sie hat eine ungewöhnliche mächtige Abfolge von Mottled Anorthosite als obersten Teil der Abfolge. Diese wohldefinierte anorthositische Abfolge wird als Giant Mottled Anorthosite (GMA) bezeichnet, und besteht hauptsächlich aus poikilitischem Anorthosit. Die GMA kann in 3 Teile, den unteren (LGMA), den mittleren (GMM) un den oberen (UGMA) unterteilt werden.Die petrographische Untersuchung von 82 Proben in der GMA aus 11 Bohrlöchern im Westteil des Bushveld-Komplexes zeigt, daß LGMA und GMM im wesentlichen Adkumulate sind, während die UGMA ein Orthokumulat darstellt. Die Hauptphase in der gesamten GMA ist Kumulus-Plagioklas ( 85 %); der Rest besteht aus Intercumulus-Orthopyroxen, Klinopyroxen, und geringen Mengen von Biotit. Diese sind die Komponenten der dunkleren mottles. Intercumulus invertierter Pigeonit kommt im oberen Teil der GMM, in der gesamten UGMA und in den unteren 30 Metern der Main-Zone vor. Innerhalb der UGMA tritt K-Feldspat hinzu, dieser kommt als Intercumulus Orthoklas und Antiperthit vor. Körner von Cumulusplagioklas mit verschiedenen Intensitäten komplexer Zonierung sind im oberen Teil der GMA verbreitet. Diese Körner werden als Hinweise auf Wachstumsepisoden während früherer Magmenzufuhr-Phasen interpretiert, die für die Entwicklung der zyklischen Einheiten im Liegenden verantwortlich waren, aber selbst nicht Teil der Cumulatabfolge wurden.

---

---

With 5 Figures     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.