A technique for modeling the Earth’s gravity field on the basis of satellite accelerations   总被引：3，自引：2，他引：1  

P. Ditmar  A. A. van Eck van der Sluijs 《Journal of Geodesy》2004,78(1-2):12-33

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.  相似文献   

Occurrence and transport of Irgarol 1051 and its major metabolite in coastal waters from South Florida   总被引：1，自引：0，他引：1  

Gardinali PR  Plasencia MD  Maxey C 《Marine pollution bulletin》2004,49(11-12):1072-1083

Irgarol 1051, a boosting antifouling agent often used to supplement copper based paints was found in surface waters from South Florida at stations collected from the Miami River, Biscayne Bay and selected areas of the Florida Keys. Concentrations of the herbicide ranged from below the method detection limit (1 ng/L) to as high as 182 ng/L in a canal system in Key Largo. The herbicide was present at 93% of the stations and often found in conjunction with its descyclopropyl metabolite (M1) previously reported to be the major degradation product of Irgarol under natural environmental conditions. The 90th percentile concentration calculated for all South Florida samples was 57.6 ng/L. Based on available data on the toxicity of Irgarol to algae and coral, only two stations (approximately 3%) ranked above the LC50 of 136 ng/L reported for the marine algae Naviculla pelliculosa and above the 100 ng/L level reported to reversibly inhibit photosynthesis of intact corals. However, a basic dissipation model for Irgarol using the Key Largo Harbor station as a point source indicated that concentrations of the herbicide decreased rapidly and concentrations below the MDL are observed within 2000 m of the source. No major coral based benthic habitats are documented for all the stations surveyed at distances that Irgarol may pose a substantial risk. However, other types of submerged vegetation like seagrasses are common around the marinas and the effects of Irgarol to such endpoints should be investigated further.  相似文献   

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

A.?JoshiEmail author  S.?Midorikawa 《Journal of Seismology》2004,8(4):467-484

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.  相似文献   

Mean and eddy motion in the Skagerrak/northern North Sea: insight from a numerical model     

Lars Petter?R?edEmail author  Ingerid?Fossum 《Ocean Dynamics》2004,54(2):197-220

Insight regarding the mean and eddy motion in the Skagerrak/northern North Sea area is gained through an analysis of model-simulated currents, hydrography, kinetic energy and relative vorticity for the 2 years 2000 and 2001. In this a -coordinate ocean model is used. Since the tidal currents are generally strong in the area, care is exercised to distinguish the mesoscale (eddy) motion from higher-frequency motion such as tides, before computing the mean and eddy kinetic energy. The model-simulated response is first compared with available knowledge of the circulation in the area, and when available, also with sea-surface temperature obtained from satellite imagery. It is concluded that the model appears to faithfully reproduce most of what is known, in particularly the upper mixed layer circulation. An analysis of the mean and eddy kinetic energy reveals that many of the mesoscale structures found in the area are recurrent. This is particularly true for the structures off the southern tip of Norway. Also in general, areas of strong mean and eddy kinetic energy are co-located. The exception is the area off the southern tip of Norway, where the eddy kinetic energy is much larger than its mean counterpart. An analysis of the relative vorticity reveals that the variability found is due to the occurrence of recurrent anticyclonic eddies. It is hypothesized that these eddies are generated due to an offshore veering of the Norwegian coastal current (NCC) as it reaches the eastern end of the Norwegian Trench plateau. Here it becomes a free jet, which is then vulnerable to either barotropic instability caused by the horizontal shear in the jet-like structure of the NCC at this point, or a baroclinic (frontal) instability. The latter may come into play when the NCC veers offshore and its relatively fresh water meets the inflowing saline water of Atlantic origin, a frontogenesis that may become strong enough for cyclogenesis to take place. Due to the depth-independent nature of the model-generated eddies, the barotropic instability is the most likely candidate. It remains to resolve the reason for the offshore veering of the NCC. The most likely candidate mechanisms are vortex squeezing or simply that the coastline curvature is large enough for the NCC to separate from the coast in a hydraulic sense.Responsible Editor: Phil Dyke  相似文献   

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

Alexander?YakubchukEmail author  Anatoly?Nikishin 《Mineralium Deposita》2004,39(2):125-142

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  相似文献   

Bulk-rock Major and Trace Element Compositions of Abyssal Peridotites: Implications for Mantle Melting, Melt Extraction and Post-melting Processes Beneath Mid-Ocean Ridges   总被引：26，自引：2，他引：24  

NIU  YAOLING 《Journal of Petrology》2004,45(12):2423-2458

This paper presents the first comprehensive major and traceelement data for 130 abyssal peridotite samples from the Pacificand Indian ocean ridge–transform systems. The data revealimportant features about the petrogenesis of these rocks, mantlemelting and melt extraction processes beneath ocean ridges,and elemental behaviours. Although abyssal peridotites are serpentinized,and have also experienced seafloor weathering, magmatic signaturesremain well preserved in the bulk-rock compositions. The betterinverse correlation of MgO with progressively heavier rare earthelements (REE) reflects varying amounts of melt depletion. Thismelt depletion may result from recent sub-ridge mantle melting,but could also be inherited from previous melt extraction eventsfrom the fertile mantle source. Light REE (LREE) in bulk-rocksamples are more enriched, not more depleted, than in the constituentclinopyroxenes (cpx) of the same sample suites. If the cpx LREErecord sub-ridge mantle melting processes, then the bulk-rockLREE must reflect post-melting refertilization. The significantcorrelations of LREE (e.g. La, Ce, Pr, Nd) with immobile highfield strength elements (HFSE, e.g. Nb and Zr) suggest thatenrichments of both LREE and HFSE resulted from a common magmaticprocess. The refertilization takes place in the ‘cold’thermal boundary layer (TBL) beneath ridges through which theascending melts migrate and interact with the advanced residues.The refertilization apparently did not affect the cpx relicsanalyzed for trace elements. This observation suggests grain-boundaryporous melt migration in the TBL. The ascending melts may notbe thermally ‘reactive’, and thus may have affectedonly cpx rims, which, together with precipitated olivine, entrappedmelt, and the rest of the rock, were subsequently serpentinized.Very large variations in bulk-rock Zr/Hf and Nb/Ta ratios areobserved, which are unexpected. The correlation between thetwo ratios is consistent with observations on basalts that D_Zr/D_Hf< 1 and D_Nb/D_Ta < 1. Given the identical charges (5⁺ forNb and Ta; 4⁺ for Zr and Hf) and essentially the same ionicradii (R_Nb/R_Ta = 1·000 and R_Zr/R_Hf = 1·006–1·026),yet a factor of 2 mass differences (M_Zr/M_Hf = 0·511 andM_Nb/M_Ta = 0·513), it is hypothesized that mass-dependentD values, or diffusion or mass-transfer rates may be importantin causing elemental fractionations during porous melt migrationin the TBL. It is also possible that some ‘exotic’phases with highly fractionated Zr/Hf and Nb/Ta ratios may existin these rocks, thus having ‘nugget’ effects onthe bulk-rock analyses. All these hypotheses need testing byconstraining the storage and distribution of all the incompatibletrace elements in mantle peridotite. As serpentine containsup to 13 wt % H₂O, and is stable up to 7 GPa before it is transformedto dense hydrous magnesium silicate phases that are stable atpressures of 5–50 GPa, it is possible that the serpentinizedperidotites may survive, at least partly, subduction-zone dehydration,and transport large amounts of H₂O (also Ba, Rb, Cs, K, U, Sr,Pb, etc. with elevated U/Pb ratios) into the deep mantle. Thelatter may contribute to the HIMU component in the source regionsof some oceanic basalts. KEY WORDS: abyssal peridotites; serpentinization; seafloor weathering; bulk-rock major and trace element compositions; mantle melting; melt extraction; melt–residue interaction; porous flows; Nb/Ta and Zr/Hf fractionations; HIMU mantle sources  相似文献   

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

Żelaźniewicz  Andrzej  Dörr  Wolfgang  Bylina  Paweł  Franke  Wolfgang  Haack  Udo  Heinisch  Helmut  Schastok  Janina  Grandmontagne  Klaus  Kulicki  Cyprian 《International Journal of Earth Sciences》2004,93(5):773-781

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.  相似文献   

Yaxcopoil-1 and the Chicxulub impact     

Wolfgang?StinnesbeckEmail author  Gerta?Keller  Thierry?Adatte  Markus?Harting  Doris?Stüben  Georg?Istrate  Utz?Kramar 《International Journal of Earth Sciences》2004,93(6):1042-1065

CSDP core Yaxcopoil-1 was drilled to a depth of 1,511 m within the Chicxulub crater. An organic-rich marly limestone near the base of the hole (1,495 to 1,452 m) was deposited in an open marine shelf environment during the latest Cenomanian (uppermost Rotalipora cushmani zone). The overlying sequence of limestones, dolomites and anhydrites (1,495 to 894 m) indicates deposition in various carbonate platform environments (e.g., sabkhas, lagoons). A 100-m-thick suevite breccia (894–794 m) identifies the Chicxulub impact event. Above the suevite breccia is a dolomitic limestone with planktic foraminiferal assemblages indicative of Plummerita hantkeninoides zone CF1, which spans the last 300 ky of the Maastrichtian. An erosional surface 50 cm above the breccia/dolomite contact marks the K/T boundary and a hiatus. Limestones above this contact contain the first Tertiary planktic foraminifera indicative of an upper P. eugubina zone P1a(2) age. Another hiatus 7 cm upsection separates zone P1a(2) and hemipelagic limestones of planktic foraminiferal Zone P1c. Planktic foraminiferal assemblages of Zone Plc to P3b age are present from a depth of 794.04 up to 775 m. The Cretaceous carbonate sequence appears to be autochthonous, with a stratigraphic sequence comparable to late Cretaceous sediments known from outside the Chicxulub crater in northern and southern Yucatan, including the late Cenomanian organic-rich marly limestone. There is no evidence that these sediments represent crater infill due to megablocks sliding into the crater, such as major disruption of sediments, chaotic changes in lithology, overturned or deep dipping megablocks, major mechanical fragmentation, shock or thermal alteration, or ductile deformation. Breccia units that are intercalated in the carbonate platform sequence are intraformational in origin (e.g., dissolution of evaporites) and dykes are rare. Major disturbances of strata by the impact therefore appear to have been confined to within less than 60 km from the proposed impact center. Yaxcopoil-1 may be located outside the collapsed transient crater cavity, either on the upper end of an elevated and tilted horst of the terrace zone, or even outside the annular crater cavity. The Chicxulub site thus records a large impact that predates the K/T boundary impact and mass extinction.  相似文献   

Palaeogeography and tectonic structure of allochthonous units in the German part of the Rheno-Hercynian Belt (Central European Variscides)     

H.?HuckriedeEmail author  K.?Wemmer  H.?Ahrendt 《International Journal of Earth Sciences》2004,93(3):414-431

New information on palaeogeography, orogenic evolution, tectonic structure, and boundaries of allochthonous units in the Rheno-Hercynian Belt is based on provenance analyses of clastic sediments and field studies. ⁴⁰K/⁴⁰Ar dating of detrital muscovites proved to be a particularly useful method because Cadomian, Caledonian and Early Variscan provenances of detrital material can be distinguished. Cadomian muscovite cooling ages are restricted to allochthonous units whereas Caledonian ages dominate within par-autochthonous and shortly displaced allochthonous units. The largest and uppermost preserved nappe, the Gießen-Harz Nappe, is derived from an oceanic flysch basin, which was not reached by Caledonian detritus. The other allochthonous units form a duplex-like structure sandwiched between the Gießen-Harz Nappe and par-autochthonous units at its base. The thick and heterogeneous roof- and floor-thrusts of this structure were previously often misinterpreted as olistostromes. The northern margin of allochthonous units is the steeply dipping Hörre-Gommern Zone. It consists of three sub-units derived from deep-water areas between the shelf at the southern margin of the Old Red Sandstone Continent and an oceanic basin to the south. The southeastern part of the duplex-structure (Harzgerode Zone) shows close affinities to Armorican terranes.  相似文献   

基于COM/DCOM的企业级GIS服务器研究     

熊伟  邰茜  蔡中祥 《测绘科学技术学报》2006,23(3):178-181

企业级GIS服务器可提供更高级、更完整的GIS服务,一般基于分布式对象技术构建.AO是ESRI完全基于COM工业标准开发的组件GIS,AGS建立在AO基础上,是企业级GIS服务器中的典范.深入研究COM组件技术、AGS的GIS服务器实现机制,对研发自主产权的企业级GIS服务器很有帮助.文中从应用角度,首先对COM/DCOM概念、本质及实现方式进行了研究;其次分析了AGS的体系结构,得出GIS服务器如何宿主和管理AO组件对象是AGS的核心所在,并对此进行了深入研究;最后提出基于COM/DCOM的PosaGIS服务器架构并进行了相关实践.  相似文献