地壳对海洋潮汐的响应

应用三维动态有限元方法研究了中国北部地区的地壳对邻近的渤海与黄海海平面变化的响应。虽然此应力场过于微弱不足以引发地震，但发现应力集中的位置及应力场变化较大的位置恰好与某些现代地震的震中一致。这一结果表明研究地壳对广泛分布的载荷的响应对研究区域地震构造是有帮助的。     

冬半年青藏高原低涡发生发展的天气学分析

本文针对冬半年高原低涡几个不同源地的低涡生成背景,着重分析500hpa环流场,指出:地形对低涡生成是非常重要的条件,但一定范围内有利的环境流场也是低涡生成不可缺少的外部条件。并给出了不同涡源区有利于低涡发生,发展的500hpa环流形势。     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Automated analysis of eclipsing binary light curves – II. Statistical analysis of OGLE LMC eclipsing binaries

In the first paper of this series, we presented EBAS – Eclipsing Binary Automated Solver, a new fully automated algorithm to analyse the light curves of eclipsing binaries, based on the ebop code. Here, we apply the new algorithm to the whole sample of 2580 binaries found in the Optical Gravitational Lensing Experiment (OGLE) Large Magellanic Cloud (LMC) photometric survey and derive the orbital elements for 1931 systems. To obtain the statistical properties of the short-period binaries of the LMC, we construct a well-defined subsample of 938 eclipsing binaries with main-sequence B-type primaries. Correcting for observational selection effects, we derive the distributions of the fractional radii of the two components and their sum, the brightness ratios and the periods of the short-period binaries. Somewhat surprisingly, the results are consistent with a flat distribution in log P between 2 and 10 d. We also estimate the total number of binaries in the LMC with the same characteristics, and not only the eclipsing binaries, to be about 5000. This figure leads us to suggest that  (0.7 ± 0.4)  per cent of the main-sequence B-type stars in the LMC are found in binaries with periods shorter than 10 d. This frequency is substantially smaller than the fraction of binaries found by small Galactic radial-velocity surveys of B stars. On the other hand, the binary frequency found by Hubble Space Telescope ( HST ) photometric searches within the late main-sequence stars of 47 Tuc is only slightly higher and still consistent with the frequency we deduced for the B stars in the LMC.     

Comet nuclei: Morphology and implied processes of surface modification

Surface morphology and related issues for nuclei of three comets: Halley, Borrelly and Wild 2, are considered in the paper. Joint consideration of publications and results of our analysis of the comets’ images led to conclusions, partly new, partly repeating conclusions published by other researchers. It was found that typical for all three nuclei is the presence of rather flat areas: floors of craters and other depressions, mesas and terraces. This implies that flattening surfaces or planation is a process typical for the comet nuclei. Planation seems to work through the sublimation-driven slope collapse and retreat. This requires effective sublimation so this process should work only when a comet is close to the Sun and if on the nucleus there are starting slopes, steep and high enough to support the “long-distance” avalanching of the collapsing material. If the surface had no starting slopes, then instead of planation, the formation of pitted-and-hilly surfaces should occur. An example of this could be the mottled terrain of the Borelly nucleus. Both ways of the sublimational evolution on the nucleus surface should lead to accumulation of cometary regolith. The thickness of the degassed regolith is not known, but it is obvious that in surface depressions, including the flat-floor ones, it should be larger compared with nondepression areas. This may have implications for the in situ study of comets by the Deep Impact and Rosetta missions.Our morphological analysis puts constraints on the applicability of the popular “rubble-pile comet nucleus” hypothesis (Weissman, 1986. Are cometery nuclei primordial rubble piles? Nature 320, 242-244.). We believe that the rubble pile hypothesis can be applicable to the blocky Halley nucleus. The Borelly and Wild 2 nuclei also could be rubble piles. But in these cases the “rubbles” have to be either smaller than 30-50 m (a requirement to keep lineament geometry close to ideal), or larger than 1-2 km (a requirement to form the rather extended smooth, flat surfaces of mesa tops and crater floors). Another option is that the Borelly and Wild 2 nuclei are not rubble piles.In relation to surface morphology we suggest that three end-member types of the comet nuclei may exist: (1) impact cratered “pristine” bodies, (2) non-cratered fragments of catastrophic disruption, and (3) highly Sun-ablated bodies. In this threefold classification, the Wild 2 nucleus is partially ablated primarily cratered body. Borrelly is significantly ablated and could be either primarily cratered or not-cratered fragment. Halley is certainly partially ablated but with the available images it is difficult to say if remnants of impact craters do exist on it.Recently published observations and early results of analysis of the Tempel 1 nucleus images taken by Deep Impact mission are in agreement with our conclusions on the processes responsible for the Halley, Borrelly and Wild 2 nuclei morphologies. In particular, we have now more grounds to suggest that decrease in crater numbers and increase of the role of smooth flat surfaces in the sequence Wild 2?Tempel 1?Borelli reflects a progress in the sublimational degradation of the nucleus surface during comet passages close to the Sun.     

Single star scidar: atmospheric parameters profiling using the simulated annealing algorithm

In this paper, we present a new method to estimate, for each turbulent layer labelled i , the horizontal wind speed   v ( h _i )  , the standard deviation of the horizontal wind speed fluctuations  σ _v ( h_i )  and the integrated value of   C ² _n   over the thickness  Δ h_i   of the turbulent layer   C ² _n ( h_i )Δ h_i   , where   h_i   is the altitude of the turbulent layer. These parameters are extracted from single star scintillation spatiotemporal cross-correlation functions of atmospheric speckles obtained within the generalized mode. This method is based on the simulated annealing algorithm to find the optimal solution required to solve the problem. Astrophysics parameters for adaptive optics are also calculated using   C ² _n ( h_i )  and   v ( h_i )  values. The results of other techniques support this new method.