A technique for modeling the Earth’s gravity field on the basis of satellite accelerations

A technique is proposed for Earths gravity field modeling on the basis of satellite accelerations that are derived from precise orbit data. The functional model rests on Newtons second law. The computational procedure is based on the pre-conditioned conjugate-gradient (PCCG) method. The data are treated as weighted average accelerations rather than as point-wise ones. As a result, a simple three-point numerical differentiation scheme can be used to derive them. Noise in the orbit-derived accelerations is strongly dependent on frequency. Therefore, the key element of the proposed technique is frequency-dependent data weighting. Fast convergence of the PCCG procedure is ensured by a block-diagonal pre-conditioner (approximation of the normal matrix), which is derived under the so-called Colombo assumptions. Both uninterrupted data sets and data with gaps can be handled. The developed technique is compared with other approaches: (1) the energy balance approach (based on the energy conservation law) and (2) the traditional approach (based on the integration of variational equations). Theoretical considerations, supported by a numerical study, show that the proposed technique is more accurate than the energy balance approach and leads to approximately the same results as the traditional one. The former finding is explained by the fact that the energy balance approach is only sensitive to the along-track force component. Information about the cross-track and the radial component of the gravitational potential gradient is lost because the corresponding force components do no work and do not contribute to the energy balance. Furthermore, it is shown that the proposed technique is much (possibly, orders of magnitude) faster than the traditional one because it does not require the computation of the normal matrix. Hints are given on how the proposed technique can be adapted to the explicit assembling of the normal matrix if the latter is needed for the computation of the model covariance matrix.Acknowledgments. Professor R. Klees is thanked for support of the project and for numerous fruitful discussions. The authors are also thankful to Dr. J. Kusche for useful remarks and to Dr. E. Schrama, his solid background in satellite geodesy proved to be very helpful. A large number of valuable comments were made by Dr. S.-C. Han, Dr. P. Schwintzer, and an anonymous reviewer; their contribution is greatly acknowledged. The satellite orbits used in the numerical study were kindly provided by Dr. P. Visser (Aerospace Department, Delft University of Technology). Access to the SGI Origin 3800 computer was provided by Stichting Nationale Computerfaciliteiten (NCF), grant SG-027.     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.     

南海北部一百万年以来的表层古生产力变化：来自ODP1144站的蛋白石记录

通过南海北部ODP 1144站蛋白石含量测定及其堆积速率的计算,并结合氧同位素记录等相关资料,获得南海北部1050ka以来高分辨率的表层古生产力变化与冰期旋回和东亚季风的关系。约900ka以来,蛋白石含量及其堆积速率较900ka以前明显增加,反映了“中更新世革命”事件之后,全球气候变冷,并导致表层生产力的提高。由于第四纪冰期旋回中的冬、夏季风的加强,加上1144站特殊的地理位置,使该站在冰期时表层生产力增加,间冰期时表层生产力降低。浮游有孔虫氧同位素记录与蛋白石含量及其堆积速率的时间序列频谱分析结果显示,三者均出现了相对应的偏心率周期、斜率周期和岁差周期,说明该站表层生产力的变化主要受地球轨道周期的驱动。     

Amoeboid olivine aggregates and related objects in carbonaceous chondrites: records of nebular and asteroid processes

Amoeboid olivine aggregates (AOAs) are the most common type of refractory inclusions in CM, CR, CH, CV, CO, and ungrouped carbonaceous chondrites Acfer 094 and Adelaide; only one AOA was found in the CB_b chondrite Hammadah al Hamra 237 and none were observed in the CB_a chondrites Bencubbin, Gujba, and Weatherford. In primitive (unaltered and unmetamorphosed) carbonaceous chondrites, AOAs consist of forsterite (Fa_<2), Fe, Ni-metal (5-12 wt% Ni), and Ca, Al-rich inclusions (CAIs) composed of Al-diopside, spinel, anorthite, and very rare melilite. Melilite is typically replaced by a fine-grained mixture of spinel, Al-diopside, and ±anorthite; spinel is replaced by anorthite. About 10% of AOAs contain low-Ca pyroxene replacing forsterite. Forsterite and spinel are always ¹⁶O-rich (δ^17,18O∼−40‰ to −50‰), whereas melilite, anorthite, and diopside could be either similarly ¹⁶O-rich or ¹⁶O-depleted to varying degrees; the latter is common in AOAs from altered and metamorphosed carbonaceous chondrites such as some CVs and COs. Low-Ca pyroxene is either ¹⁶O-rich (δ^17,18O∼−40‰) or ¹⁶O-poor (δ^17,18O∼0‰). Most AOAs in CV chondrites have unfractionated (∼2-10×CI) rare-earth element patterns. AOAs have similar textures, mineralogy and oxygen isotopic compositions to those of forsterite-rich accretionary rims surrounding different types of CAIs (compact and fluffy Type A, Type B, and fine-grained, spinel-rich) in CV and CR chondrites. AOAs in primitive carbonaceous chondrites show no evidence for alteration and thermal metamorphism. Secondary minerals in AOAs from CR, CM, and CO, and CV chondrites are similar to those in chondrules, CAIs, and matrices of their host meteorites and include phyllosilicates, magnetite, carbonates, nepheline, sodalite, grossular, wollastonite, hedenbergite, andradite, and ferrous olivine.Our observations and a thermodynamic analysis suggest that AOAs and forsterite-rich accretionary rims formed in ¹⁶O-rich gaseous reservoirs, probably in the CAI-forming region(s), as aggregates of solar nebular condensates originally composed of forsterite, Fe, Ni-metal, and CAIs. Some of the CAIs were melted prior to aggregation into AOAs and experienced formation of Wark-Lovering rims. Before and possibly after the aggregation, melilite and spinel in CAIs reacted with SiO and Mg of the solar nebula gas enriched in ¹⁶O to form Al-diopside and anorthite. Forsterite in some AOAs reacted with ¹⁶O-enriched SiO gas to form low-Ca pyroxene. Some other AOAs were either reheated in ¹⁶O-poor gaseous reservoirs or coated by ¹⁶O-depleted pyroxene-rich dust and melted to varying degrees, possibly during chondrule formation. The most extensively melted AOAs experienced oxygen isotope exchange with ¹⁶O-poor nebular gas and may have been transformed into magnesian (Type I) chondrules. Secondary mineralization and at least some of the oxygen isotope exchange in AOAs from altered and metamorphosed chondrites must have resulted from alteration in the presence of aqueous solutions after aggregation and lithification of the chondrite parent asteroids.     

Noril’sk–Talnakh Cu–Ni–PGE deposits: a revised tectonic model

The Norilsk mining district is located at the northwest margin of the Tunguska basin, in the centre of the 3,000×4,000 km Siberian continental flood basalt (CFB) province. This CFB province was formed at the Permo-Triassic boundary from a superplume that ascended into the geometric centre of the Laurasian continent, which was surrounded by subducting slabs of oceanic crust. We suggest that these slabs could have reached the core–mantle boundary, and they may have controlled the geometric focus of the superplume. The resulting voluminous magma intruded and erupted in continental rifts and related extensive flood basalt events over a 2–4 Ma period. Cu–Ni–PGE sulfide mineralization is found in olivine-bearing differentiated mafic intrusions beneath the flood basalts at the northwestern margin of the Siberian craton and also in the Taimyr Peninsula, some 300 km east of a triple junction of continental rifts, now buried beneath the Mesozoic–Cenozoic sedimentary basin of western Siberia. The Norilsk-I and Talnakh-Oktyabrsky deposits occur in the Norilsk–Kharaelakh trough of the Tunguska CFB basin. The Cu–Ni–PGE-bearing mineralized intrusions are 2–3 km-wide and 20 km-long differentiated chonoliths. Previous studies suggested that parts of the magma remained in intermediate-level crustal chambers where sulfide saturation and accumulation took place before emplacement. The 5–7-km-thick Neoproterozoic to Palaeozoic country rocks, containing sedimentary Cu mineralization and evaporites, may have contributed additional metal and sulfur to this magma. Classic tectonomagmatic models for these deposits proposed that subvertical crustal faults, such as the northeast-trending Norilsk–Kharaelakh fault, were major trough-parallel conduits providing access for magmas to the final chambers. However, geological maps of the Norilsk region show that the Norilsk–Kharaelakh fault offsets the mineralization, which was deformed into folds and offset by related reverse faults, indicating compressional deformation after mineralization in the Late Triassic to Early Jurassic. In addition, most of the intrusions are sills, not dykes as should be expected if the vertical faults were major conduits. A revised tectonic model for the Norilsk region takes into account the fold structure and sill morphology of the dominant intrusions, indicating a lateral rather than vertical emplacement direction for the magma into final chambers. Taking into account the fold structure of the country rocks, the present distribution of the differentiated intrusions hosting the Norilsk-I and Talnakh–Oktyabrsky deposits may represent the remnants of a single, >60 km long, deformed and eroded palm-shaped cluster of mineralized intrusions, which are perceived as separate intrusions at the present erosional level. The original direction of sill emplacement may have been controlled by a northeast-trending paleo-rise, which we suggest is present at the southeastern border of the Norilsk–Kharaelakh trough based on analysis of the unconformity at the base of the CFB. The mineralized intrusions extend along this rise, which we interpret as a structure that formed above the extensionally tilted block in the metamorphic basement. Geophysical data indicate the presence of an intermediate magma chamber that could be linked with the Talnakh intrusion. In turn, this T-shaped flat chamber may link with the Yenisei–Khatanga rift along the northwest-trending Pyasina transform fault, which may have served as the principal magma conduit to the intermediate chamber. It then produced the differentiated mineralized intrusions that melted through the evaporites with in situ precipitation of massive, disseminated, and copper sulfide ore. The Norilsk–Kharaelakh crustal fault may not relate to mineralization and possibly formed in response to late Mesozoic spreading in the Arctic Ocean.Editorial handling: P. Lightfoot     

The granite hosted gold deposit of Moulin de Chéni (Saint-Yrieix district,Massif Central,France): petrographic,structural, fluid inclusion and oxygen isotope constraints

The Moulin de Chéni orogenic gold deposit is the only granite-hosted deposit of the Saint-Yrieix district, French Massif Central. It occurs in 338±1.5 Ma-old peraluminous leucogranites and is characterized by intense microfracturing and bleaching of the granite in relation to pervasive sulfide crystallization. Formation of quartz veins and gold deposition occurred in two successive stages: an early mesozonal stage of quartz-sulfide (Fe-As-S) deposition, usually devoid of gold and a late epizonal stage of base metal and gold deposition. Both stages postdate peak metamorphism and granite intrusion. The genesis of the deposit is the result of four successive fluid events: (1) Percolation of aqueous-carbonic metamorphic fluids under an assumed lithostatic regime of 400–450 °C, at a maximum depth of 13 km; (2) Formation of the main quartz lodes with coeval K-alteration and introduction of As and S from aqueous-carbonic fluids percolating along regional faults. Arsenopyrite and pyrite deposition was linked to the alteration of Fe-silicates into K-feldspar and phengite at near-constant iron content in the bulk granite. Temperature was similar to that of the preceding stage, but pressure decreased to 100–50 MPa, suggesting rapid uplift of the basement up to 7.5 km depth; (3) The resulting extensional tectonic leads to the deposition of gold, boulangerite, galena and sphalerite in brecciated arsenopyrite and pyrite from aqueous fluids during a mixing process. Temperature and salinity decrease from 280 to 140 °C and 8.1 wt% eq. NaCl to 1.6 wt% eq. NaCl, respectively; (4) Sealing of the late fault system by barren comb quartz which precipitated from dilute meteoric aqueous fluids (1.6 wt% eq. NaCl to 0.9 wt% eq. NaCl) under hydrostatic conditions at 200–150 °C.Editorial handling: B. Lehmann     

The eastern continuation of the Cadomian orogen: U–Pb zircon evidence from Saxo-Thuringian granitoids in south-western Poland and the northern Czech Republic

We report U–Pb single zircon ages from three pre-Variscan granitoids in the NE part of the Bohemian Massif. The Platerówka granodiorite from the Lausitz-Izera Unit, the Polish Sudetes, has been dated at 533±9 Ma. The Bitouchov granite form the SW part of the South Krkonoe Unit, the Czech Sudetes, gave an age of 540+11/–10 Ma, and the Wdroe granodiorite in the Fore-Sudetic Block yielded 548±9 Ma. All these latest Vendian/Early Cambrian granitoids represent the post-tectonic expression of a late Proterozoic Cadomian orogenic cycle and demonstrate the eastward extent of the Cadomian basement into the Variscan orogen. Granodiorites of similar age have so far been reported from Brittany and especially from the Saxo-Thuringian Terrane to the NE and SW of the Elbe Fault Zone. We conclude that the Saxo-Thuringian Terrane extends across the Elbe and Sudetic Marginal Fault Zones into the Fore-Sudetic Block.     

亳州盐渍土的分布及其工程地质特征

盐渍土是气候、地貌、地质、工程地质、水文地质等因素综合影响的工程地质灾害,本文在大量调查、分析与研究的基础上,采用了槽探、水样(水化学分析)、土样(土壤含盐量分析、土壤毛细水试验)等综合方法,分析了亳州盐渍土灾害的危害性及对农业生态系统、公路建设的影响,同时,考虑到人为和环境地质因素对盐渍土的影响,圈定了盐渍土的潜在威胁区,从而,预测了盐渍土的发展、演变趋势,因地制宜地提出了整治盐渍土的基本对策与措施。     

Rhinocerotid tooth enamel 18O/16O variability between 23 and 12 Ma in southwestern France

The relationship between the oxygen isotope ratio of mammal tooth enamel and that of drinking water was used to reconstruct changes in the Miocene oxygen isotope ratio of rainfall (meteoric water ${\delta }^{18}$O_MW). These, in turn, are related to climatic parameters (temperature, precipitation and evaporation rate). ${\delta }^{18}$O values of rhinocerotid teeth from the Aquitaine Basin (southwestern France) suggest a significant climatic change between 17 and 12 Ma, characterized by cooling together with precipitation increase, in agreement with other terrestrial and oceanic records. To cite this article: I. Bentaleb et al., C. R. Geoscience 338 (2006).     

The Trade and Industries of Western South America

Abstract

This research was based on an instructional module developed and used to investigate whether GIS can be used to enhance fourth grade students' geographic knowledge and map skills. Another goal was to identify challenges the teachers and the students face using GIS.

Findings from the study suggest that using GIS in the classroom helps students improve their geographic knowledge and maps skills. In addition, this study found that the students' geographic content knowledge tended to improve gradually and the knowledge transferred to the subsequent lessons, while certain map skills tended to develop in a more specific manner.