全球地震纬度分布特征和北纬35°峰值异常浅析

本文取1900年─1980年全球Ms≥6．0地震7936次分析后发现，地震发生在南北回归线（23°5N—23°5S）之间占49．69％，发生在南北极因（63°5N—63°5S）之间占99．17％，发生在南北极圈外只占0．83％，即很少发生地震；其中8级地震98次，在35°N线为高峰。又分析了公元648年─1979年全球M≥8级地震257次主要亦集中在35°N线附近。尤其是历史上破坏性极大的地震在35°N线附近者多。这些现象可能和地球自转有一定关系。     

Optimal weighting in the finite difference solution of the convection-dispersion equation

Oscillation and numerical dispersion limit the reliability of numerical solutions of the convection-dispersion equation when finite difference methods are used. To eliminate oscillation and reduce the numerical dispersion, an optimal upstream weighting with finite differences is proposed. The optimal values of upstream weighting coefficients numerically obtained are a function of the mesh Peclet number used. The accuracy of the proposed numerical method is tested against two classical problems for which analytical solutions exist. The comparison of the numerical results obtained with different numerical schemes and those obtained by the analytical solutions demonstrates the possibility of a real gain in precision using the proposed optimal weighting method. This gain in precision is verified by interpreting a tracer experiment performed in a laboratory column.     

Conductivity of NaCl solution at 0.4–5.0 GPa and 25–500 °C

NaCI-H₂O is the most fundamental ternary system in geology. Until now, the measurements of electrical conductivity of NaCl solutions are still little at high pressures (> O.5 GPa) We measured the conductivity of 0.01 m NaCl solution at 0.4–5.0 GPa and 25-500°C. The results are consistent with that of Quist and Marshall (1968) at 0.4 GPa. The conductivity of NaCl solution increases with increasing temperature. The results also show that the conductivity of NaCl solution changes little with increasing pressure below 1.5 GPa and changes rapidly with increasing pressure above 1.5 GPa. The rapid increase of the conductivity of NaCl solution may play an important role in many geological processes (such as the genesis of ore deposits under hydrothermal condition) and other fields. Project supported by the National Natural Science Foundation of China.     

Forest operations,tree species composition and decline in rainfall explain runoff changes in the Nacimiento experimental catchments,south central Chile

Few long-term studies have explored how intensively managed short rotation forest plantations interact with climate variability. We examine how prolonged severe drought and forest operations affect runoff in 11 experimental catchments on private corporate forest land near Nacimiento in south central Chile over the period 2008–2019. The catchments (7.7–414 ha) contain forest plantations of exotic fast-growing species (Pinus radiata, Eucalyptus spp.) at various stages of growth in a Mediterranean climate (mean long-term annual rainfall = 1381 mm). Since 2010, a drought, unprecedented in recent history, has reduced rainfall at Nacimiento by 20%, relative to the long-term mean. Pre-drought runoff ratios were <0.2 under 8-year-old Eucalyptus; >0.4 under 21-year-old Radiata pine and >0.8 where herbicide treatments had controlled vegetation for 2 years in 38% of the catchment area. Early in the study period, clearcutting of Radiata pine (85%–95% of catchment area) increased streamflow by 150 mm as compared with the year before harvest, while clearcutting and partial cuts of Eucalyptus did not increase streamflow. During 2008–2019, the combination of emerging drought and forestry treatments (replanting with Eucalyptus after clearcutting of Radiata pine and Eucalyptus) reduced streamflow by 400–500 mm, and regeneration of previously herbicide-treated vegetation combined with growth of Eucalyptus plantations reduced streamflow by 1125 mm (87% of mean annual precipitation 2010–2019). These results from one of the most comprehensive forest catchment studies in the world on private industrial forest land indicate that multiple decades of forest management have reduced deep soil moisture reservoirs. This effect has been exacerbated by drought and conversion from Radiata pine to Eucalyptus, apparently largely eliminating subsurface supply to streamflow. The findings reveal tradeoffs between wood production and water supply, provide lessons for adapting forest management to the projected future drier climate in Chile, and underscore the need for continued experimental work in managed forest plantations.     

福建仙游地震序列的震源机制解

1,209 earthquakes occurred in Xianyou, Fujian from August 4, 2010 to October 4, 2013. The largest earthquake was ML5.0 on September 4, 2013. In order to study the Xianyou earthquake sequence and understand the causative structure and stress field of Xianyou, the focal mechanism solutions of six earthquakes （ ML 〉 3. 5 ） in the Xianyou earthquake sequence are calculated using the broadband digital data of the Fujian Seismic Network with the seismic moment tensor inverse method. The results show that the focal faults of the six earthquakes are similar, which are all strike-slip faults striking to the northwest with high dip angles. The direction of the principal compressive stress axes is near SN, which is different from the stress field of Fujian region. The Xianyou earthquake sequence may have been induced by the stress adjustment after the impoundment of Jinzhong reservoir.     

Sea‐level rise impact on fresh groundwater lenses in two‐layer small islands

The fresh groundwater lenses (FGLs) of small islands can be highly vulnerable to climate change impacts, including sea‐level rise (SLR). Many real cases of atoll or sandy islands involve two‐layer hydrogeological conceptualizations. In this paper, the influential factors that affect FGLs in two‐layer small islands subject to SLR are investigated. An analytical solution describing FGLs in circular islands, composed of two geological layers, is developed for the simplified case of steady‐state and sharp‐interface conditions. An application of the developed model is demonstrated to estimate the FGL thickness of some real‐world islands by comparison with existing FGL thickness data. Furthermore, numerical modelling is applied to extend the analysis to consider dispersion effects and to confirm comparable results for both cases. Sensitivity analyses are used to assess the importance of land‐surface inundation caused by SLR, relative to other parameters (i.e. thickness of aquifer layers, hydraulic conductivity, recharge rate and land‐surface slope) that influence the FGL. Dimensionless parameters are used to generalize the findings. The results demonstrate that land‐surface inundation has a considerable impact on a FGL influenced by SLR, as expected, although the FGL volume is more sensitive to recharge, aquifer thickness and hydraulic conductivity than SLR impacts, considering typical parameter ranges. The methodology presented in this study provides water resource managers with a rapid‐assessment tool for evaluating the likely impacts of SLR and accompanying LSI on FGLs.     

Reconstruction of Conceptual Prediction Model for the Three Rainfall Patterns in the summer of eastern China under global warming

With the influence of global warming,the global climate has undergone significant inter-decadal variation since the late 1970s.Although El Nio-Southern Oscillation(ENSO)has been the strongest signal for predicting global climate inter-annual variability,its relation with the summer rainfall in China has significantly changed,and its indicative function on the summer rainfall in China has weakened.This has led to a significant decrease in the accuracy rate of early conceptual prediction models for the Three Rainfall Patterns in the summer of eastern China.On the basis of the difference analysis of atmospheric circulation system configuration in summer,as well as the interaction of ocean and atmospheric in previous winter between two phases,i.e.before and after the significant global warming(1951 to 1978 and 1979 to 2012,respectively),we concluded that(1)Under different inter-decadal backgrounds,the atmospheric circulations that impacted the Three Rainfall Patterns in the summer of eastern China showed consistency,but in the latter phase of the global warming,the Western Pacific Subtropical High(WPSH)was on the strong side,the position of which was in the south,and the blocking high in the Eurasia mid-high latitudes was active,while the polar vortex extended to the south,and meridional circulation intensified.This circulation background may have been conducive to the increase of the circulation frequency of Patterns II and III,and the decrease of the circulation frequency of Pattern I,thus leading to more Patterns II and III and fewer Pattern I in the summer rainfall of eastern China.(2)In the former phase,the corresponding previous winter SST fields of different rainfall patterns showed visible differences.The impact of ENSO on North Pacific Oscillation(NPO)was great,and the identification ability of which on Patterns I and II of summer rainfall was effective.In the latter phase,this identification ability decreased,while the impact of ENSO on the Pacific/North American(PNA)teleconnection pattern increased,and the identification ability of the PNA on Patterns II and III also increased.Based on the new inter-decadal climate background,this study reconstructs the conceptual prediction model for the Three Rainfall Patterns in summer of eastern China by using the previous winter PNA and the Eurasian(EU)teleconnection indexes.The fitting effect was satisfying,though it is necessary to be further tested.     

Elastic and anelastic anomalies in (Ca,Sr)TiO3 perovskites: Analogue behaviour for silicate perovskites

Values of bulk modulus (K), shear modulus (G) and mechanical quality factor (Q) have been determined for polycrystalline samples across the CaTiO₃ (CST0)–SrTiO₃ (CST100) solid solution by resonant ultrasound spectroscopy. Because of similarities with low frequency elastic and anelastic anomalies due to twin wall motion reported in previous studies, a working hypothesis is developed in which dissipation processes are interpreted in terms of twin wall displacements. At high temperatures in CST50 the stability field of the I4/mcm structure is marked by the disappearance of all resonance peaks (superattenuation). This is attributed to anelastic domain wall sliding. At room temperature the I4/mcm phase of CST70 and CST80 has values of G which are lower than those of cubic or orthorhombic phases, and a concomitant drop in Q is interpreted as implying that the domain wall pinning process reported elsewhere to occur below 400–450 K is only partial. A similar drop in G and Q was found in CST95 below the transition at 238 K. The I4/mcm ↔ Pbcm transition in CST70 at 230 K is marked by an abrupt increase in Q, suggesting that mobile twins in crystals with the I4/mcm structure become effectively immobile in antiferroelectric crystals with the Pbcm structure. The I4/mcm ↔ Pnma transition in CST50 is marked by a similarly abrupt increase in Q, consistent with twin walls becoming effectively immobile also in crystals with the Pnma structure. A fall in Q below 800 K in CST0, however, could imply that a degree of twin wall mobility might develop in Pnma crystals if the tetragonal spontaneous strain departs significantly from zero. The remarkable attenuation behaviour of crystals with the I4/mcm structure at the relatively low stress conditions which apply during resonances of a parallelepiped with edge dimensions of 2–4 mm, is consistent with the view that a characteristic signature for tetragonal CaSiO₃ in the Earth's lower mantle should be a marked attenuation of seismic waves.     

Comparison of two different solutions to Molodensky’s G1 term

The main purpose of this paper is comparison of two different approaches of solution to the Simple Molodensky’s Problem, the Molodensky’s Approach and the Analytical Continuation Approach, based on numerical computation. Although these approaches have been described theoretically by several authors, e.g. Molodensky et al. (1960), Heiskanen and Moritz (1967), Vaníček (1974), Moritz (1980) and Holota (1991, unpublished results) and theoretical proof of equivalence was given by Heiskanen and Moritz (1967), Moritz (1971), Ecker (1971) and Pellinen (1972, unpublished results), only very few practical experiences about the differences between particular solutions and computational efficiency exist. In this paper we compare the above two mentioned approaches in terms of the G₁-effect on quasigeoid. Both quasigeoid solutions are tested by the independent GPS/levelling approach and are also compared with the previous quasigeoid model of Slovakia where the G₁-term has been approximated using the classical terrain correction. The effect of the G₂-term is also numerically estimated, revealing that it might be significant for the precise quasigeoid determination. Some practical comments and recommendations are given at the end of the paper.     

On the existence of atmospheric attractors

In this paper, we consider the long-time dynamics for the primitive equations of large-scale dry at- mosphere. First, by energy methods, we obtain the existence and uniqueness of global strong solu- tions of the problem. Second, by studying the long-time behavior of strong solutions, we construct a global attractor which captures all the trajectories.