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.     

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.     

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

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.     

Complex Site Effects in Thessaloniki (Greece): I. Soil Structure and Comparison of Observations with 1D Analysis

This paper presents a 2D model of the geological structure of Thessaloniki city and results of empirical and theoretical approaches for the evaluation of site response due to complex site effects. The construction of the 2D model is based on the available geophysical and geotechnical data in terms of the most important parameters needed to model site response. The well-known subsoil structure, despite the existence of some local uncertainties, gave the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from the analysis of one event in frequency and time domains and 1D numerical modelling. The study of site response shows the effect of the lateral variations on ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In this case history, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects.