Note of Reviewer on the Paper “On the Inductive Proof of Legendre Addition Theorem” by K. Maleček and Z. Nádeník

Studia Geophysica et Geodaetica -     

小尺度地磁场勒让德多项式建模方法

根据地面实测地磁场数据,提出一种基于勒让德多项式的对小尺度地磁场的模型建立方法,结合数据预处理、最小二乘拟合方法,采用高阶的模型描述典型小尺度区域的地磁场,绘制出地磁图.通过实验验证,模型拟合实验点的数据与真实数据的磁场总强度平均误差小于30 nT,并且反映出小尺度区域中磁场的精细结构,有效地抑制了边界效应,大大提高了地磁场建模能力.另外讨论了勒让德多项式建模方法的优缺点以及使用范围.这为地磁场在空间延拓、探矿、地下结构分析提供了重要的方法和理论基础.     

小尺度地磁场勒让德多项式建模方法

根据地面实测地磁场数据,提出一种基于勒让德多项式的对小尺度地磁场的模型建立方法,结合数据预处理、最小二乘拟合方法,采用高阶的模型描述典型小尺度区域的地磁场,绘制出地磁图.通过实验验证,模型拟合实验点的数据与真实数据的磁场总强度平均误差小于30 nT,并且反映出小尺度区域中磁场的精细结构,有效地抑制了边界效应,大大提高了地磁场建模能力.另外讨论了勒让德多项式建模方法的优缺点以及使用范围.这为地磁场在空间延拓、探矿、地下结构分析提供了重要的方法和理论基础.     

海洋局域地磁正常场勒让德多项式模型的建立

研究了利用勒让德多项式构建局部海域地磁正常场模型的方法,给出了一种趋势面地磁数据滤波方法和模型截止阶数的确定方法,研究和分析了参与建模的地磁总强度数据的代表性和密度对地磁正常场模型精度的影响,最后提出了一种分区建模的思想.试验表明分区建模相对整体建模具有较高的精度.     

不同规格化的引潮位展开及其转换

引入完全规格化的缔合勒让德函数,阐述国际上两类不同规格化的引潮位展开,使引潮位函数表达更加简明,大地系数定义更加规范. 由此得到由Doodson规格化与Hartmann & Wenzel规格化之间的转换系数.     

关于地磁场超低频脉动(ULF)的研究

本文是关于地磁场超低频脉动研究的综述,内容包括：（1）历史和内容提要;（2）脉动分类的探究;（3）关于“四多”（多观测仪器、多卫星和飞船、多地面和网络、多分析方法）;（4）其它行星ULF波的探求.     

地震空间分布函数的确定方法研究

选取了几个反央西南地区地震活动的时间和空间特征的因子，因子量化全部建立在因子与实际地震发生的统计关系上，使因子具有一定的概率预报意义和相互对比的基础，另外还对部分因子的预报效能了作评估和修正，采用３种方法分别求出了３套地震空间分布函数，并作了要相互比较。     

Meixner functions for derivation of the unit hydrograph

The Meixner functions are utilized to relate the effective rainfall, the direct runoff and the unit hydrograph through linkage equations. The linkage equations are then employed to derive the unit hydrograph for given rainfall-runoff data on a small agricultural watershed. These functions are tested with regard to their ability to reproduce and predict the direct runoff hydrograph. The Meixner functions are found to be an effective analytical tool for hydrograph synthesis. Further, they compare well with the least squares and linear programming methods of the unit hydrograph derivation.     

The BEAM-project: prediction and assessment of mixture toxicities in the aquatic environment

Freshwater and marine ecosystems are exposed to various multi-component mixtures of pollutants. Nevertheless, most ecotoxicological research and chemicals regulation focus on hazard and exposure assessment of individual substances only, the problem of chemical mixtures in the environment is ignored to a large extent. In contrast, the assessment of combination effects has a long tradition in pharmacology, where mixtures of chemicals are specifically designed to develop new products, e.g. human and veterinary drugs or agricultural and non-agricultural pesticides. In this area, two concepts are frequently used and are thought to describe fundamental relationships between single substance and mixture effects: Independent Action (Response Addition) and Concentration Addition. The question, to what extent these concepts may also be applied in an ecotoxicological and regulatory context may be considered a research topic of major importance, as the concepts would allow to make use of already existing single substance toxicity data for the predictive assessment of mixture toxicities. Two critical knowledge gaps are identified: (a) There is a lack of environmental realism, as a huge part of our current knowledge about the applicability of the concepts is restricted to artificial situations with respect to mixture composition or biological effect assessment. (b) The knowledge on what exactly is needed for using the concepts as tools for the predictive mixture toxicity assessment is insufficient. Both gaps seriously hamper the necessary, scientifically sound consideration of mixture toxicities in a regulatory context.In this paper, the two concepts will be briefly introduced, the necessity of considering the toxicities of chemical mixtures in the environment will be demonstrated and the applicability of Independent Action and Concentration Addition as tools for the prediction and assessment of mixture toxicities will be discussed. An overview of the specific aims and approaches of the BEAM project to fill in the identified knowledge gaps is given and first results are outlined.     

太湖富营养化问题及其综合控制对策

On the basis of nine-time current situation investigation for eutrophication of Taihu Lake during 1991-1995, this paper evaluated the trophic levels in the different periods and analyzed the development of the main nutrient content in the nearest 35 years. The results show that the trophic level of Taihu Lake is in the transition state from meso-eutrophic to eutrophic. The eutrophhic and hypereutrophic waters account for 10% or so. The limiting nutrient, P, rises most rapidly, which causes the ratio of N:P to decrease. The increase of P content is still one of the main factors giving rise to the eutrophication of Taihu. LakeSome proposals of comprehensive countermeasure for the eutrophication are put forward. They include the pollution source control of the basin, the littoral multiple management, optimal dispatch of water conservancy facilities, the engineering of helping Taihu Lake with diversion of the Changjiang River, and as well as the setting of the water quality protection and legal system.     

Theorem of turbulent intensity and macroscopic mechanism of the turbulence development

Turbulence is one of the most common nature phenomena in everyday experience, but that is not adequately understood yet. This article reviews the history and present state of development of the turbulence theory and indicates the necessity to probe into the turbulent features and mechanism with the different methods at different levels. Therefore this article proves a theorem of turbulent transpor- tation and a theorem of turbulent intensity by using the theory of the nonequilibrium thermodynamics, and that the Reynolds turbulence and the Rayleigh-Bénard turbulence are united in the theorems of the turbulent intensity and the turbulent transportation. The macroscopic cause of the development of fluid turbulence is a result from shearing effect of the velocity together with the temperature, which is also the macroscopic cause of the stretch and fold of trajectory in the phase space of turbulent field. And it is proved by the observed data of atmosphere that the phenomenological coefficient of turbulent in- tensity is not only a function of the velocity shear but also a function of temperature shear, viz the sta- bility of temperature stratification, in the atmosphere. Accordingly, authenticity of the theorem, which is proved by the theory of nonequilibrium thermodynamics, of turbulent intensity is testified by the facts of observational experiment.     

Short-term precursor of the geomagnetic horizontal field transfer functions

IntroductionThe history of studying the changes in deep geoelectric conductivity of an earthquake preparation using the parameters of geomagnetic short-period variations has been about 30 years. At first, Japanese calculated the real transfer functions as well as the a and b coefficients of Parkinson vector by taking the measurement of graph and using empirical formula. There are some examples, such as the 1923 Kanto M=7.9 earthquake (Yanagihara, 1972), the 1966 Tashkent M=5.5 earthquake, t…     

Modelling considerations in probabilistic performance‐based seismic evaluation: case study of the I‐880 viaduct

Limitations associated with deterministic methods to quantify demands and develop rational acceptance criteria have led to the emergence of probabilistic procedures in performance‐based seismic engineering. The Pacific Earthquake Engineering Research performance‐based methodology is one such approach. In this paper, the impact of certain modelling decisions made at different stages of the evaluation process on the performance assessment of a typical multi‐bent viaduct is examined. Modelling, in the context of this paper, covers hazard modelling, structural modelling and loss modelling. The specific application considered in this study is a section of an existing viaduct in California: the I‐880 interstate highway. Several simulation models of the viaduct are developed, a series of nonlinear time‐history analyses are carried out to predict demands, measures of damage are evaluated and the probability of closure of the viaduct is estimated using the specified hazard for the site. It is concluded that the methodology offers several advantages over existing deterministic performance‐based procedures. Results of the investigation indicate that the assessment methodology is particularly sensitive to the reliability of decisions made by bridge inspectors following a seismic event, and to the dispersion in the demand estimation, which in turn is influenced by several factors including soil–structure interaction effects and ground motion scaling procedures. Copyright © 2005 John Wiley & Sons, Ltd.     

A conceptual model of the relationship between channel characteristics and discharge

Regression equations with the form y = f(q), where y is a channel characteristic and q is a discharge of a particular frequency are widely used in fluvial geomorphology. They imply that channel geometry is related to a dominant discharge. For this to be so, the river must be ‘in regime’ i.e. its characteristics must fluctuate about a mean condition which is effectively stationary at the time scale considered. The static regression-dominant discharge approach does not fully describe the relationship between channel characteristics and discharge and an alternative model is proposed with the form y(t) = ∞ h(u)- q(t-u)du. This model states that channel form is a product of the whole series of discharges experienced by the channel. The dominant discharge-static regression approach provides a working approximation of the relationship between channel characteristics and discharge only when temporal variations in q are small or when y is insensitive to variations in q.     

唐贞观晋州地震考

对唐朝3次灾害最严重的地震作了研究,认为949年晋州(今山西临汾)地震死亡人数应为5000余人。从震亡人数、余震序列,对比734年天水、749年河套地震,认为这次地震烈度可达9度,震级可达7级。     

Rupture process of large earthquakes in the northern Mexico subduction zone

The Cocos plate subducts beneath North America at the Mexico trench. The northernmost segment of this trench, between the Orozco and Rivera fracture zones, has ruptured in a sequence of five large earthquakes from 1973 to 1985; the Jan. 30, 1973 Colima event (M _s 7.5) at the northern end of the segment near Rivera fracture zone; the Mar. 14, 1979 Petatlan event (M _s 7.6) at the southern end of the segment on the Orozco fracture zone; the Oct. 25, 1981 Playa Azul event (M _s 7.3) in the middle of the Michoacan gap; the Sept. 19, 1985 Michoacan mainshock (M _s 8.1); and the Sept. 21, 1985 Michoacan aftershock (M _s 7.6) that reruptured part of the Petatlan zone. Body wave inversion for the rupture process of these earthquakes finds the best: earthquake depth; focal mechanism; overall source time function; and seismic moment, for each earthquake. In addition, we have determined spatial concentrations of seismic moment release for the Colima earthquake, and the Michoacan mainshock and aftershock. These spatial concentrations of slip are interpreted as asperities; and the resultant asperity distribution for Mexico is compared to other subduction zones. The body wave inversion technique also determines theMoment Tensor Rate Functions; but there is no evidence for statistically significant changes in the moment tensor during rupture for any of the five earthquakes. An appendix describes theMoment Tensor Rate Functions methodology in detail.The systematic bias between global and regional determinations of epicentral locations in Mexico must be resolved to enable plotting of asperities with aftershocks and geographic features. We have spatially shifted all of our results to regional determinations of epicenters. The best point source depths for the five earthquakes are all above 30 km, consistent with the idea that the down-dip edge of the seismogenic plate interface in Mexico is shallow compared to other subduction zones. Consideration of uncertainties in the focal mechanisms allows us to state that all five earthquakes occurred on fault planes with the same strike (N65°W to N70°W) and dip (15±3°), except for the smaller Playa Azul event at the down-dip edge which has a steeper dip angle of 20 to 25°. However, the Petatlan earthquake does prefer a fault plane that is rotated to a more east-west orientation—one explanation may be that this earthquake is located near the crest of the subducting Orozco fracture zone. The slip vectors of all five earthquakes are similar and generally consistent with the NUVEL-predicted Cocos-North America convergence direction of N33°E for this segment. The most important deviation is the more northerly slip direction for the Petatlan earthquake. Also, the slip vectors from the Harvard CMT solutions for large and small events in this segment prefer an overall convergence direction of about N20°E to N25°E.All five earthquakes share a common feature in the rupture process: each earthquake has a small initial precursory arrival followed by a large pulse of moment release with a distinct onset. The delay time varies from 4 s for the Playa Azul event to 8 s for the Colima event. While there is some evidence of spatial concentration of moment release for each event, our overall asperity distribution for the northern Mexico segment consists of one clear asperity, in the epicentral region of the 1973 Colima earthquake, and then a scattering of diffuse and overlapping regions of high moment release for the remainder of the segment. This character is directly displayed in the overlapping of rupture zones between the 1979 Petatlan event and the 1985 Michoacan aftershock. This character of the asperity distribution is in contrast to the widely spaced distinct asperities in the northern Japan-Kuriles Islands subduction zone, but is somewhat similar to the asperity distributions found in the central Peru and Santa Cruz Islands subduction zones. Subduction of the Orozco fracture zone may strongly affect the seismogenic character as the overlapping rupture zones are located on the crest of the subducted fracture zone. There is also a distinct change in the physiography of the upper plate that coincides with the subducting fracture zone, and the Guerrero seismic gap to the south of the Petatlan earthquake is in the wake of the Orozco fracture zone. At the northern end, the Rivera fracture zone in the subducting plate and the Colima graben in the upper plate coincide with the northernmost extent of the Colima rupture zone.     

利用P波接收函数研究青藏高原东南缘地幔转换带结构

研究青藏高原东南缘的深部结构对于理解印度-欧亚板块的碰撞机理和青藏高原的形成演化具有重要的科学意义.本研究对布设在研究区域内566个固定和流动地震台站的波形资料进行了处理,获得77853条高质量P波接收函数,应用接收函数共转换点（CCP）叠加技术获得了研究区域下方精细的地幔转换带间断面起伏形态及转换带厚度变化图像.结果表明：研究区域南北方向上具有两个明显的转换带增厚异常区,南侧异常区位于滇中次级块体与印支块体下方,可能是新特提斯洋板片与上部印度板块间断离并部分滞留在转换带底部的结果;北侧川西地区异常增厚可能与上方岩石圈拆沉并降至转换带有关;腾冲火山起源可能是板块俯冲过程中发生断离造成软流圈物质部分熔融,湿热物质上涌所致.

     

Transitions in the style of mantle convection at high Rayleigh numbers

The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 10⁴ and 10⁸), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈10⁵. A second transition occurs at higher Rayleigh numbers, Ra≈5×10⁶, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRa^β) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.     

六盘山盆地热历史的裂变径迹证据

研究盆地的热历史将为确定生烃过程和探勘目标提供重要制约因素。磷灰石裂变径迹研究表明，六盘山盆地白垩系地层经历两次埋深加热事件。第一次在白垩经末之前达到最高古地温，第二次在晚新生代约8MaB.P.之前达到最高古地温.第一次最高古地温要高于第二次最高古地温.晚新生代六盘山盆地古地温梯度为16℃/km.从白垩纪到新生代,六盘山盆地可能发生古地温梯度降低事件.三叠纪、中侏罗统烃源岩达到或超过生油高峰温度.白垩纪乃家河组和马东山组虽然进入生油窗温度范围,但未达到生油峰温度.按照古地温资料推断,三叠纪、中侏罗统烃源岩应为六盘山盆地主要生烃源岩,晚白垩纪之前应为六盘山盆地烃源岩的主要生烃阶段.     

Structure of the crust and uppermost mantle from broadband seismic array in the western Dabie Mountains,east central China

A broadband seismic array of 7 stations was set up in the western Dabie Mountains (31°20′-31°50′N, 114°30′-115°E). Teleseismic events from May 2001 to November 2001 were collected and analyzed by radial receiver function to determine the S-wave velocity structure of the crust and uppermost mantle. The crustal thickness is 32-38 km beneath the array. The crust-mantle boundary appears as a gently north-dipping velocity discontinuity, but turns to be a velocity gradient beneath a station near the Qiliping shea...