Samovar: a thermomechanical code for modeling of geodynamic processes in the lithosphere—application to basin evolution

We present a new 2D finite difference code, Samovar, for high-resolution numerical modeling of complex geodynamic processes. Examples are collision of lithospheric plates (including mountain building and subduction) and lithosphere extension (including formation of sedimentary basins, regions of extended crust, and rift zones). The code models deformation of the lithosphere with viscoelastoplastic rheology, including erosion/sedimentation processes and formation of shear zones in areas of high stresses. It also models steady-state and transient conductive and advective thermal processes including partial melting and magma transport in the lithosphere. The thermal and mechanical parts of the code are tested for a series of physical problems with analytical solutions. We apply the code to geodynamic modeling by examining numerically the processes of lithosphere extension and basin formation. The results are directly applicable to the Basin and Range province, western USA, and demonstrate the roles of crust–mantle coupling, preexisting weakness zones, and erosion rate on the evolutionary trends of extending continental regions. Modeling of basin evolution indicates a critical role of syn-rift sedimentation on the basin depth and a governing role of Peierls deformation in cold lithospheric mantle. While the former may increase basin depth by 50%, the latter limits the depth of rift basins by preventing faulting in the subcrustal lithosphere.     

Possibilities for water conservation in Bulgaria

The drought during the last 10 years in Bulgaria caused the interruption of water supply to more than 70% of the population of the country. Investigations on recent water demand showed that the role of water losses and irrational use so considerable that existing water resources became inqdequate. Bringing new water to existing supply systems is connected with considerable investments, time, and controversy. Water conservation could solve some water supply problems by providing more reasonable water consumption and additional time for development of new resources. In this paper the possibilities for reducing water use in Bulgaria and the need for developing a national conservation program are discussed. An analysis of the variables which effect water use is performed. The reasons for the water losses and using more than the necessary water quantities are discussed. Developing of water-use efficiency standards for plumbing fixtures and conservation-oriented water rates and measures are proposed.     

Morphotectonics, Volcanism and Hydrothermal Activity on the East Pacific Rise between 21°12′ S and 22°40′ S

The morphotectonic setting of the East Pacific Rise (EPR) between21°12 and 22°40 S and its recent and past hydrothermalactivity were the focus of the Russian R/V Geolog Fersmans expeditionin 1987–1988.The EPR axial zone in the study area is comprised of three segmentsseparated by overlapping spreading centers (OSCs) near 21°44 and22°08 S. The northern segment is the shallowest of three and hasa distinct massive axial ridge, trapeziodal in cross-section, toppedby a very wide flat summit surface and cut by a well-developedcentral graben. These features testify to intense magmatism and to avoluminous crustal magmatic chamber underlying the whole segment.Fine-scale segmentation is most clearly revealed in the structure ofthe central graben within which several 4th-order segments can bedistinguished. This scale of segmentation is also reflected on flanks of theaxis by variations in the character and intensity of faulting.According to structural and petrologic data, the magmatism is mostintense in the central part of the segment which is probably locateddirectly over a magmatic diapir supplying the melt to the whole segment.Magma migration at the subcrustal level from the center towards the ends ofthe segment with discrete injection into the crustal magmatic chamber ispresumed.The central segment is broken into two morphologically distinct partsseparated by a deval. In the subsided northern part, the wide summit of theaxial ridge is cut by a well-developed, intensely fractured axialgraben. In the southern part, the axial ridge is relatively elevated, butnarrow with an ephemeral graben along its crest. The character and intensityof faulting on the axial flanks are also considerably different in thenorthern and southern parts of the segment. Thus, the magmatic supply tothese two parts is thought to originate from two different sources. If so,then at present the magma chamber underlying the southern part of thesegment is probably at the stage of replenishment, while in the north it isat the stage of deep cooling.The southern segment is structurally similar to the central one. Howeverthere is considerably less intensive magmatic activity in this region,especially south of 22°30 S where the axial ridge is narrow, andtriangular in cross-section.Both OSCs studied are marked by abrupt narrowing and sharp subsidence ofthe tips of axial ridges within the northern limbs. The southern OSC limbsare morphologically similar to normal sections of axial ridges. In bothcases the flanks are structurally and morphologically disrupted adjacent tothe OSCs and oblique structures can be traced far southward of the OSCflanks. Due to the spatial position of oblique structures on the the flanksit is presumed that the OSC near 22°07 S is migrating northward.The 21°44 S OSC zone has apparently undergone small spatialoscillations. In spite of the small amplitude of lateral displacement, thiszone is marked by prominent bathymetric anomalies.Numerous massive sulfide deposits were discovered atop the axial ridgealong the entire length of the uplifted and hydrothermally active northernsegment. Ore metal concentrations in near-bottom waters are maximumover the southern part of the northern segment, while maximum concentrationsof the same metals in surficial sediments are confined to the central partof the same segment. We surmise that there has been a recentalong-axis shift of the zone of maximum hydrothermal activity fromthe middle of the segment to its present position in the southern part ofthe segment. Considering sedimentation rates, the age of this shift can beapproximately estimated to be 5 to 10 thousand years before the present.The relatively Mg-enriched basalts of the middle part of thenorthern segment represent a tike of a more primitive pattern, while therelatively Fe-rich rocks of its southern part probably reflect alarge degree of fractionation at shallow crustal levels. Considering thistrend, in addition to morphotectonic data we presume that subaxial magmaflow from the middle to the southern part of the segment is responsible forthe along-axis shift of hydrothermal activity.In the central segment of the study area, massive sulfides have only beendiscovered south of the 21°55 S deval, where the axial ridgeshoals and where the existence of a subjacent magma chamber is presumed.The very weak manifestations of recent volcanism within the southernsegment explain the absence of hydrothermal activity and sulfide depositswithin this segment.     

Chemistry of fluids from a natural analogue for a geological CO2 storage site (Montmiral, France): Lessons for CO2–water–rock interaction assessment and monitoring

Chemical and isotope studies of natural CO₂ accumulations aid in assessing the chemical effects of CO₂ on rock and thus provide a potential for understanding the long-term geochemical processes involved in CO₂ geological storage. Several natural CO₂ accumulations were discovered during gas and oil exploration in France’s carbogaseous peri-Alpine province (south-eastern France) in the 1960s. One of these, the Montmiral accumulation at a depth of more than 2400 m, is currently being exploited. The chemical composition of the water collected at the wellhead has changed in time and the final salinity exceeds 75 g/L. These changes in time can be explained by assuming that the fraction of the reservoir brine in the recovered brine–CO₂–H₂O mixture varies, resulting in variable proportions of H₂O and brine in the sampled water. The proportions can be estimated in selected samples due to the availability of gas and water flowrate data. These data enabled the reconstruction of the chemical and isotope composition of the brine. The proportions of H₂O and brine can also be estimated from isotope (δ²H, δ¹⁸O) composition of collected water and δ¹⁸O of the sulfates or CO₂. The reconstituted brine has a salinity of more than 85 g/L and, according to its Br⁻ content and isotope (δ²H, δ¹⁸O, δ³⁴S) composition, originates from an evaporated Triassic seawater that underwent dilution by meteoric water. The reconstitution of the brine’s chemical composition enabled an evaluation of the CO₂–water–rock interactions based on: (1) mineral saturation indices; and (2) comparison with initial evaporated Triassic seawater. Dissolution of K- and SO₄-containing minerals such as K-feldspar and anhydrite, and precipitation of Ca and Mg containing minerals that are able to trap CO₂ (carbonates) are highlighted. The changes in concentration of these elements in the brine, which are attributed to CO₂ interactions, illustrate the relevance of monitoring the water quality at future industrial CO₂ storage sites.     

Multi-Element Reference Samples of Black Shale

This paper describes the creation of two new reference samples of black shale, SCHS-1 and SLg-1, collected in the Bodaibo Region of Irkutsk District (Siberia, Russia). These samples have been developed as State Standard Samples to meet the requirements of Russia's Normative Documents, so that they can be recognised by the Russian State Certification Agency and included on the State Register of certified reference materials. SCHS-1 was produced from shale of the Khomolkho suite; SLg-1 represents the average composition of the Sukhoy Log gold deposit (within Au mineralisation). The homogeneity was assessed by applying the concept of potential inhomogeneity, and the inhomogeneity uncertainties (σ_inhom) for the indicator elements that were estimated experimentally by silicate XRF analysis. It was found that a majority of the indicator elements was characterised by σ_inhom/Δ_p < 1/8 (Δ_p is the permissible certification uncertainty), and so σ_inhom was neglected for the indicator elements and all other elements that were included in their respective group. The value of σ_inhom calculated for Fe (σ_inhom,Fe/Δ_p > 1 /8) was taken into consideration for the noble metals and their associated elements (e.g., Ag, As, Au, etc.). More than fifty laboratories (thirty-seven from Russia) were involved in the interlaboratory experiment (ILE) to analyse these reference samples. The data obtained from the ILE allowed SCHS-1 and SLg-1 to be certified for, fifty-one and forty-three elements, respectively, as State Standard Samples.     

Focal mechanism of earthquakes in the Kamchatka-Commander region and heterogeneities of the active seismic zone

Summary Focal mechanisms of 74 shallow and 16 intermediate earthquakes (1964–1970) in the Kamchatka-Commander region are discussed. Regional stress systems and the nature of faulting are analyzed. Complex stress fields in the Kamchatka-Commander region and variable aftershock processes of strong Kamchatka earthquakes show the existence of heterogeneities of the active seismic zone in the Region.     

Hot Neutral H in the Heliosphere: Elastic H-H, H-p Collisions

Williams et al. (1997) have suggested that a population of hot hydrogen atoms is created in the heliosphere through elastic H-H collisions between energetic `solar' atoms (neutralized solar wind) and interstellar atoms. They used a BGK-like approximation (Bhatnagar et al., 1954) for the Boltzmann collision term and the collision cross sections suggested by Dalgarno (1960). We show that both assumptions result in a significant overestimation of the the H-H collision effect. On the basis of calculated momentum transfer cross-sections for elastic H-H collisions, we argue that elastic H-H and H-p collisions cannot produce hot H atoms in the heliosphere. This revised version was published online in July 2006 with corrections to the Cover Date.     

Faunal communities at sites of gas- and oil-bearing fluids in Lake Baikal

Macro- and meiofaunal communities were examined at four geomorphologically distinct sites with different gas- and oil-bearing fluid characteristics in the northern, central and southern basins of Lake Baikal. All sites had elevated concentrations of bicarbonate, nitrate, sulphate and chloride ions in pore fluids, with highest values at the Frolikha vent. Elevated levels of iron ions were found in pore waters of the St. Petersburg methane seep and the Gorevoy Utes oil seep. The chemical composition of pore waters at the Malenky mud volcano was similar to that reported in earlier work. Consistent with published data, the Frolikha vent (northern basin) and the St. Petersburg methane seep (central basin) were characterised by methane of mixed genesis (thermogenic + biogenic), whereas the methane source was mainly thermogenic at the Gorevoy Utes oil seep (central basin) and biogenic at the Malenky mud volcano (southern basin). In contrast to marine seep ecosystems, the macrofauna was dominated only by amphipods, giant planarians and oligochaetes, whereas bivalves were absent; the meiofauna was similar to its marine counterpart, being dominated by nematodes, cyclops, harpacticoids and ostracods. A statistically significant positive relationship was revealed between faunal abundance and the availability of bacterial mats on seep sediments. Moreover, ANOVA tests showed significant increases in both meiozoobenthic and macrozoobenthic densities at “hot spot” vent/seep sites relative to discharge-free reference sites. The isotopic composition of carbon and nitrogen at various trophic levels of these benthic vent/seep communities was found to differ markedly from that reported by earlier studies for the pelagic and other benthic food webs in Lake Baikal. As in marine seeps, the macrofauna had variable isotopic signatures. Light δ¹³C and δ¹⁵N values suggest the utilization of chemosynthetically fixed and/or methane-derived organic matter. By contrast, the heavy δ¹³C signatures of some mobile amphipods likely reflect consumption of photosynthetically derived carbon. These findings would at least partly explain why Lake Baikal is a notable outlier in global temperature–biodiversity patterns, exhibiting the highest biodiversity of any lake worldwide but at an extremely cold average temperature.     

Rotational evolution of exoplanets under the action of gravitational and magnetic perturbations

We investigate the evolution of the rotational axes of exoplanets under the action of gravitational and magnetic perturbations. The planet is assumed to be dynamically symmetrical and to be magnetised along its dynamical-symmetry axis. By qualitative methods of the bifurcation theory of multiparametric PDEs, we have derived a gallery of 69 phase portraits. The portraits illustrate evolutionary trajectories of the angular momentum of a planet for a variety of the initial conditions, for different values of the ratio between parameters describing gravitational and magnetic perturbations, and for different rates of the orbital evolution. We provide examples of the phase portraits, that reveal the differences in topology and the evolutionary track of in the vicinity of an equilibrium state. We determine the bifurcation properties, i.e., the way of reorganisation of phase trajectories in the vicinities of equilibria; and we point out the combinations of parameters’ values that permit ip-overs from a prograde to a retrograde spin mode.