Geochemistry constraints of Mesozoic–Cenozoic calc-alkaline magmatism in the South Shetland arc, Antarctica

Geochemical data from basalts, basaltic andesites, and andesites of the Mesozoic–Cenozoic (143–44 Ma) from Livingston, Greenwich, Robert, King George, and Ardley Islands of the South Shetland archipelago, Antarctica, are presented. The rocks have variable SiO₂ of approximately 46–61 wt%, Al₂O₃ of 15–26 wt%, and total alkali (K₂O+Na₂O) of 2–6 wt%. Most samples have low Mg#, Cr, and Ni, which indicates that they have undergone significant fractional crystallization from mantle-derived melts. The presence of olivine cumulatic in the samples from Livingston and Robert Islands explains some high MgO, Ni, and Cr values, whereas low Rb, Zr, and Nb values could be related to undifferentiated magmas. N-MORB-normalized trace element patterns show that South Shetland Islands volcanic rocks have a geochemical pattern similar to that found for other island arcs, with enrichment in LILE relative to HFSE and in LREE relative to HREE. The geochemistry pattern and presence of calcic plagioclase, orthopyroxene, Mg-olivine, and titanomagnetite phenocrysts suggest a source related to the subduction process. The geochemical data also suggest magma evolution from the tholeiitic to the calc-alkaline series; some samples show a transitional pattern. Samples from the South Shetland archipelago show moderate LREE/HREE ratios relative to N-MORB and OIB, depletion in Nb relative to Yb, and high Th/Yb ratios. These patterns probably reflect magma derived from a lithospheric mantle source previously modified by fluids and sediments from a subduction zone.

Resumo

Dados geoquímicos de basaltos, andesitos basálticos e andesitos mesozóicos–cenozóicos (143–44 Ma) das ilhas Livingston, Greenwich, Robert, King George e Ardley do Arquipélago Shetland do Sul, Antártica são discutidas neste artigo. As rochas tem conteúdos de SiO₂ variando de 46 a 61%, Al₂O₃ de 15 a 26% e álcalis (K₂O+Na₂O) de 2 a 6%. A maior parte das amostras tem conteúdos baixos de Mg#, Cr e Ni, indicando que sofreram significante cristalização fracionada de fusões derivadas do manto. A presença de fases cumuláticas nas amostras das ilhas Livingston e Robert explicaria os elevados valores de MgO, Ni, Cr, enquanto que baixos valores de Rb, Zr e Nb observados nas amostras destas ilhas poderiam estar relacionados a magmas não diferenciados. Os padrões de elementos-traço normalizados pelo N-MORB mostram que as rochas vulcânicas das Ilhas Shetland do Sul têm padrão geoquímico similar àqueles encontrados em outros arcos de ilhas com enriquecimento em LILE em relação aos HFSE e em ETRL em relação aos ETRP. O padrão geoquímico e a ocorrência de fenocristais de plagioclásio cálcico, ortopiroxênio, olivina magnesiana e titanomagnetita sugerem origem relacionada a processos de subducção. Dados geoquímicos obtidos para as amostras do arquipélago Shetland do Sul sugerem um magma evoluindo de toleítico para cálcico-alcalino, observando-se em algumas amostras um padrão transicional. As amostras do arquipélago Shetland do Sul mostram em relação ao N-MORB e OIB, moderadas razões ETRL/ETRP, empobrecimento em Nb relativo a Yb e elevada razão Th/Yb Estes padrões refletem, provavelmente, magma derivado de uma fonte mantélica litosférica, que foi modificada por fluídos e sedimentos da zona de subducção.     

Assessing the Vibrational Frequencies of the Holweide Hospital in the City of Cologne (Germany) by Means of Ambient Seismic Noise Analysis and FE modelling

Ambient seismic noise measurements were conducted inside the Holweide Hospital (Cologne) for checking whether its frequencies of vibration fall into the range where soil amplification is expected. If this is the case, damage might increase in case of an earthquake due to an amplified structural response of the building. Two different techniques were used: the ratio between the horizontal and vertical components of the spectra recorded at stations located inside the building and the ratio between the corresponding components of the spectra recorded simultaneously inside the building and at a reference station placed outside. While the former method might be preferred because of less equipment involved, the latter has the advantage of producing more stable results and deleting automatically the influence of the sedimentary cover, which might obscure some eigenfrequencies of vibration of the building. An independently performed finite-element analysis of the hospital showed a good correlation between measured and calculated eigenmodes.     

Ab initio study of the high-pressure behavior of CaSiO3 perovskite

Using density functional simulations, within the generalized gradient approximation and projector-augmented wave method, we study structures and energetics of CaSiO₃ perovskite in the pressure range of the Earths lower mantle (0–150 GPa). At zero Kelvin temperature the cubic CaSiO₃ perovskite structure is unstable in the whole pressure range, at low pressures the orthorhombic (Pnam) structure is preferred. At 14.2 GPa there is a phase transition to the tetragonal (I4/mcm) phase. The CaIrO₃-type structure is not stable for CaSiO₃. Our results also rule out the possibility of decomposition into oxides.

地质体三维可视化建模的技术方法研究

地质体三维可视化模型具有形象、直观、动态、准确、信息丰富等特点,在指导地质找矿和矿山开发等方面具有很大的优势及实用价值。文章对地质体三维可视化建模的整个技术方法流程作了总结性阐述,并分析其技术关键,同时还进一步探讨了模型的分类及其应用。     

Probing the Obscuring Medium Around Active Nuclei Using Masers: The Case of 3C 403

We report the first detection of a water megamaser in a radio-loud galaxy, 3C 403, and present a follow-up study using the VLA. 3C 403 has been observed as a part of a small sample of FR II galaxies with evidence of nuclear obscuration. The isotropic luminosity of the maser is 1200 L. With a recessional velocity of cz 17680 km s^–1 it is the most distant water maser so far reported. The line arises from the densest (> 10⁸ cm^–3) interstellar gas component ever observed in a radio-loud galaxy. Two spectral features are identified, likely bracketing the systemic velocity of the galaxy. Our interferometric data clearly indicate that these arise from a location within 0.1 (110 pc) from the active galactic nucleus. We conclude that the maser spots are most likely associated with the tangentially seen parts of a nuclear accretion disk, while an association with dense warm gas interacting with the radio jets cannot yet be ruled out entirely.     

Geometry and growth of sill complexes: insights using 3D seismic from the North Rockall Trough

Doleritic sill complexes, which are an important component of volcanic continental margins, can be imaged using 3D seismic reflection data. This allows unprecedented access to the complete 3D geometry of the bodies and an opportunity to test classic sill emplacement models. The doleritic sills associated with basaltic volcanism in the North Rockall Trough occur in two forms. Radially symmetrical sill complexes consist of a saucer-like inner sill at the base with an arcuate inclined sheet connecting it to a gently inclined, commonly ragged, outer rim. Bilaterally symmetrical sill complexes are sourced by magma diverted from a magma conduit feeding an overlying volcano. With an elongate, concave upwards, trough-like geometry bilaterally symmetrical sills climb away from the magma source from which they originate. Both sill complex types can appear as isolated bodies but commonly occur in close proximity and consequently merge, producing hybrid sill complexes. Radial sill complexes consist of a series of radiating primary flow units. With dimensions up to 3 km, each primary flow unit rises from the inner saucer and is fed by primary magma tube. Primary flow units contain secondary flow units with dimensions up to 2 km, each being fed by a secondary magma tube branching from the primary magma tube. Secondary flow units in turn are composed of 100-m scale tertiary flow units. A similar branching hierarchy of flow units can also be seen in bilaterally symmetrical sill complexes, with their internal architecture resembling an enlarged version of a primary flow unit from a radial sill complex. This branching flow pattern, as well as the interaction between flow units of varying orders, provides new insights into the origin of the structures commonly seen within sill complexes and the hybrid sill bodies produced by their merger. The data demonstrate that each radially symmetrical sill complex is independently fed from a source located beneath the centre of the inner saucer, grows by climbing from the centre outwards and that peripheral dyking from the upper surface is a common feature. These features suggest a laccolith emplacement style involving peripheral fracturing and dyking during inner saucer growth and thickening. The branching hierarchy of flow units within bilaterally symmetrical sill complexes is broadly similar to that of primary flow units within a radially symmetrical sill complex, suggesting that the general features of the laccolith emplacement model also apply.Editorial responsibility: J. Stix     

Marquardt算法磁测三维反演方法策略

在磁法三维反演中,反演对象数目往往巨大,受限于当今的计算机技术与计算方法,在超大规模方程组的求解时遇到极大的困难。所以,在解决实际问题时能力非常有限,对于Marquardt算法磁测三维反演问题尤为突出。为此,在这里提出了两个较为有效的方法策略,即将模型组体化和将组体分类,并在理论应用中得到了较好的效果。     

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.     

Building and Testing an <Emphasis Type="Italic">a priori</Emphasis> Geophysical Model for Western Eurasia and North Africa

We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q _P and Q _S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle _n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.