Potential environmental issues of CO2 storage in deep saline aquifers: Geochemical results from the Frio-I Brine Pilot test,Texas, USA

Sedimentary basins in general, and deep saline aquifers in particular, are being investigated as possible repositories for large volumes of anthropogenic CO₂ that must be sequestered to mitigate global warming and related climate changes. To investigate the potential for the long-term storage of CO₂ in such aquifers, 1600 t of CO₂ were injected at 1500 m depth into a 24-m-thick “C” sandstone unit of the Frio Formation, a regional aquifer in the US Gulf Coast. Fluid samples obtained before CO₂ injection from the injection well and an observation well 30 m updip showed a Na–Ca–Cl type brine with ∼93,000 mg/L TDS at saturation with CH₄ at reservoir conditions; gas analyses showed that CH₄ comprised ∼95% of dissolved gas, but CO₂ was low at 0.3%. Following CO₂ breakthrough, 51 h after injection, samples showed sharp drops in pH (6.5–5.7), pronounced increases in alkalinity (100–3000 mg/L as HCO₃) and in Fe (30–1100 mg/L), a slug of very high DOC values, and significant shifts in the isotopic compositions of H₂O, DIC, and CH₄. These data, coupled with geochemical modeling, indicate corrosion of pipe and well casing as well as rapid dissolution of minerals, especially calcite and iron oxyhydroxides, both caused by lowered pH (initially ∼3.0 at subsurface conditions) of the brine in contact with supercritical CO₂.     

Low induction number, ground conductivity meters: A correction procedure in the absence of magnetic effects

Ground conductivity meters, comprising a variety of coil–coil configurations, are intended to operate within the limits provided by a low induction number (LIN), electromagnetic condition. They are now routinely used across a wide range of application areas and the measured apparent conductivity data may be spatially assembled and examined/correlated alongside information obtained from many other earth science, environmental, soil and land use assessments. The theoretical behaviour of the common systems is examined in relation to both the prevailing level of subsurface conductivity and the instrument elevation. It is demonstrated that, given the inherent high level of accuracy of modern instruments, the prevailing LIN condition may require operation in environments restricted to very low (< 12 mS/m) conductivities. Beyond this limit, non-linear departures from the apparent conductivity that would be associated with a LIN condition occur and are a function of the coil configuration, the instrument height and the prevailing conductivity. Using both theory and experimental data, it is demonstrated that this has the potential to provide biased and spatially distorted measurements. A simple correction procedure that can be applied to the measured data obtained from any of the LIN instruments is developed. The correction procedure would, in the limit of a uniform subsurface, return the same (correct) conductivity, irrespective of the ground conductivity meter used, the prevailing conductivity or the measurement height.     

Toward micro-scale spatial modeling of gentrification

 A simple preliminary model of gentrification is presented. The model is based on an irregular cellular automaton architecture drawing on the concept of proximal space, which is well suited to the spatial externalities present in housing markets at the local scale. The rent gap hypothesis on which the model's cell transition rules are based is discussed. The model's transition rules are described in detail. Practical difficulties in configuring and initializing the model are described and its typical behavior reported. Prospects for further development of the model are discussed. The current model structure, while inadequate, is well suited to further elaboration and the incorporation of other interesting and relevant effects. Received: 1 November 2001 / Accepted: 12 April 2002 Town centres data were made available by kind permission of the UK Department of the Environment, Transport and the Regions.     

A sedimentary record of the rise and fall of the metal industry in Bergslagen,south central Sweden

Detailed stratigraphic analyses of sediments deposited in Lake Botjärnen, a small boreal forest lake in the shield terrain of central Sweden, clearly reflect progressively increasing human impact on terrestrial and aquatic ecosystems following settlement and establishment of an iron industry in the 17th century. Rising frequencies of pollen and spores from light-demanding plants provide evidence of extensive forest clearance for charcoal and timber production, which peaked in the early 20th century. An associated increase in catchment erosion is reflected by changing carbon and nitrogen elemental content and carbon–isotope composition of sediment organic matter and by increased magnetic susceptibility of the sediments. Records of air-borne pollutants (lead, zinc and sulphur) can be correlated to the development of local and regional mining and metal industry as inferred from historical accounts. Rapid recession of the iron industry led to re-forestation and recovery of the aquatic nutrient status to pre-industrial conditions over the past 100 years. The chronology of the sediment succession, which is based on ²¹⁰Pb and ¹³⁷Cs radionuclide data in combination with radiocarbon dating, is confirmed by historical lead pollution trends established for the region.     

Avoiding numerical stability problems of long duration DGPS/INS Kalman filters

A current pursuit of the geodetic community is the optimal integration of differential GPS (DGPS) and inertial navigation system (INS) data streams for precise and efficient position and gravity vector surveying. Therein a complete INS and multiple-antenna GPS receiver payload, mounted on a moving platform, is used in conjunction with a network of ground-fixed single antenna GPS receivers. This paper presents a complete, GPS-based, external updating measurement model for the applicable Kalman filter. The model utilizes four external observation types for every GPS satellite in-view: DGPS range differences, single phase differences, and single phase-rate differences; as well as the mobile, multipleantenna GPS receiver's measurement of theerrors in the INS's estimate of the phase difference between any two vehicle-borne GPS antennae. Although not widely conveyed in the geodetic world, the inertial navigation community has long known that traditional Kalman filter covariance propagation recurrences are inherently unstable when such highly accurate external updates are repeatedly applied (every 1 second) over long time durations. A hybrid square root covariance/U — D covariance factorization approach is a numerically stable alternative and is reviewed herein. The hybrid makeup of the algorithm is necessitated by the correlated nature of the fourth type of GPS external measurement listed above (each vehicle-borne GPS antenna formstwo baselines). Such measurement correlations require a functional transformation of the overall external updating model to permit the multiple updates (simultaneously available at each updating epoch) to be sequentially (and efficiently) processed. An appropriate transformation is given. Stable covariance propagation relationships are presented and the transformed Kalman gain is also furnished and its use in the determination of the externally updated error states is discussed. Specific DGPS/INS instabilities produced by the traditional recurrences are displayed. The stable alternative method requires about 25% more CPU time than the traditional Kalman recurrences. With the ever-increasing computational speeds of microprocessors, this added CPU time is of no real concern.     

SOME ANALYSES OF RECENT CHINESE PROVINCIAL DATA

A complete set of Chinese provincial data has recently been published and may be employed to study the regional variations in China. Great regional inequalities are found in industrial outputs, provisions of university places and distribution of qualified doctors.     

Marymia: an Archean, amphibolite facies-hosted, orogenic lode-gold deposit overprinted by Palaeoproterozoic orogenesis and base metal mineralisation, Western Australia

The Marymia gold deposit, comprising two orebodies, Keillor 1 and Keillor 2, is at the northern end of the Plutonic Well greenstone belt in the Marymia Inlier, in the southern Capricorn Orogen, just north of the Yilgarn craton. The Marymia Inlier is a discrete fault-bounded Archean gneiss-granitoid-greenstone domain surrounded by sedimentary basins that were formed and variably metamorphosed and deformed during several Palaeoproterozoic orogenic cycles. The greenstone sequence at Marymia is stratigraphically and geochemically similar to greenstone sequences in the Yilgarn craton, but was subjected to further deformation and metamorphism in the Palaeoproterozoic. Late Archean deformation (D1-D2) was ductile to brittle-ductile in style, whereas Palaeoproterozoic deformation was predominantly brittle. Equilibrium mineral assemblages indicate that peak amphibolite-facies metamorphism (540-575 °C, <3 kb) was overprinted by greenschist-facies metamorphism (300-360 °C). Petrographic textures indicate that prograde metamorphism was coeval with D1-D2, with peak metamorphism early to syn D2. Gold mineralisation at Marymia is hosted in metamorphosed tholeiitic basalts and banded iron formation. On a gross scale, the distribution of gold is controlled by D2 folds and shear zones. Lithological contacts with strong rheological or chemical contrasts provide local controls. Gold-related alteration comprises subtle millimetre- to centimetre-wide zones of silicification with variable amounts of quartz, hornblende, biotite, K-feldspar, plagioclase, calcite/siderite, scheelite, titanite, epidote, sulfide and telluride minerals. Quartz veins are generally narrow and discontinuous with low total volume of quartz. Gold is sited in the wall rock, at vein salvedges or within stringers of wall rock within veins. There are two distinct opaque-mineral assemblages: pyrite-pyrrhotite-chalcopyrite-galena and hessite-petzite-altaite-Bi-telluride-galena. Ore samples are variably enriched in Ag, Te, Pb, W, Cu, S and Fe reflecting heterogeneity of the ore mineralogy. Structural timing and temperature of formation of alteration and ore minerals support deposition of gold during late peak amphibolite-facies metamorphism from neutral to alkaline (pH=5-6), moderately oxidising (log P_O2,-21-22) and CO₂-bearing (X_CO2 Ƹ.2) fluids. The total sulfur content of the fluid is estimated at 1mDS. Lead isotope compositions support derivation of lead from within the local greenstone sequence. Gold lodes were deformed by faults and shear zones in the Palaeoproterozoic, with only limited remobilisation. Subeconomic, carbonate vein- and breccia-hosted base metal mineralisation is locally hosted within Palaeoproterozoic fault zones, which clearly cut gold lodes. Base-metal-related alteration is characterised by intense carbonatisation, chloritisation, and albitisation of the mafic host rocks. Mineral assemblages are consistent with formation at greenschist facies conditions. Lead isotope compositions support crystallisation at ca. 1.7 Ga from lead that is similar in composition to earlier gold-related galena.     

Phase Relations and Melting of Anhydrous K-bearing Eclogite from 1200 to 1600{degrees}C and 3 to 5 GPa

To investigate eclogite melting under mantle conditions, wehave performed a series of piston-cylinder experiments usinga homogeneous synthetic starting material (GA2) that is representativeof altered mid-ocean ridge basalt. Experiments were conductedat pressures of 3·0, 4·0 and 5·0 GPa andover a temperature range of 1200–1600°C. The subsolidusmineralogy of GA2 consists of garnet and clinopyroxene withminor quartz–coesite, rutile and feldspar. Solidus temperaturesare located at 1230°C at 3·0 GPa and 1300°C at5·0 GPa, giving a steep solidus slope of 30–40°C/GPa.Melting intervals are in excess of 200°C and increase withpressure up to 5·0 GPa. At 3·0 GPa feldspar, rutileand quartz are residual phases up to 40°C above the solidus,whereas at higher pressures feldspar and rutile are rapidlymelted out above the solidus. Garnet and clinopyroxene are theonly residual phases once melt fractions exceed 20% and garnetis the sole liquidus phase over the investigated pressure range.With increasing melt fraction garnet and clinopyroxene becomeprogressively more Mg-rich, whereas coexisting melts vary fromK-rich dacites at low degrees of melting to basaltic andesitesat high melt fractions. Increasing pressure tends to increasethe jadeite and Ca-eskolaite components in clinopyroxene andenhance the modal proportion of garnet at low melt fractions,which effects a marked reduction in the Al₂O₃ and Na₂O contentof the melt with pressure. In contrast, the TiO₂ and K₂O contentsof the low-degree melts increase with increasing pressure; thusNa₂O and K₂O behave in a contrasted manner as a function ofpressure. Altered oceanic basalt is an important component ofcrust returned to the mantle via plate subduction, so GA2 maybe representative of one of many different mafic lithologiespresent in the upper mantle. During upwelling of heterogeneousmantle domains, these mafic rock-types may undergo extensivemelting at great depths, because of their low solidus temperaturescompared with mantle peridotite. Melt batches may be highlyvariable in composition depending on the composition and degreeof melting of the source, the depth of melting, and the degreeof magma mixing. Some of the eclogite-derived melts may alsoreact with and refertilize surrounding peridotite, which itselfmay partially melt with further upwelling. Such complex magma-genesisconditions may partly explain the wide spectrum of primitivemagma compositions found within oceanic basalt suites. KEY WORDS: eclogite; experimental petrology; mafic magmatism; mantle melting; oceanic basalts