Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.

     

Gravity and magnetic data inversion for 3D topography of the Moho discontinuity in the northern Red Sea area, Egypt

The main goal of our study is to investigate 3D topography of the Moho boundary for the area of the northern Red Sea including Gulf of Suez and Gulf of Aqaba. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. To separate sources of gravity and magnetic field in depth, a method is suggested, based on upward and downward continuation. Both new methods are applied to isolate the contribution of the Moho interface to the total field and to find its 3D topography. At the first stage, we separate near-surface and deeper sources. According to the obtained field of shallow sources a model of the horizontal layer above the depth of 7 km is suggested, which includes a density interface between light sediments and crystalline basement. Its depressions and uplifts correspond to known geological structures. At the next stage, we isolate the effect of very deep sources (below 100 km) and sources outside the area of investigation. After subtracting this field from the total effect of deeper sources, we obtain the contribution of the Moho interface. We make inversion separately for the area of rifts (Red Sea, Gulf of Suez and Gulf of Aqaba) and for the rest of the area. In the rift area we look for the upper boundary of low-density, heated anomalous upper mantle. In the rest of the area the field is satisfied by means of topography for the interface between lower crust and normal upper mantle. Both algorithms are applied also to the magnetic field. The magnetic model of the Moho boundary is in agreement with the gravitational one.     

接收函数确定Moho面速度和密度跃变的方法研究及应用

通过对单层模型反射和透射系数的推导,提出了利用接收函数一次转换波和多次波确定Moho面速度和密度跃变的速度-密度跃变（δβ-δρ）扫描叠加方法.利用反射率法计算了不同模型的远震理论地震图,按照与处理实际观测波形一致的方法和流程计算了理论接收函数;根据不同模型数值试验结果,深入分析了界面速度和密度跃变对接收函数震相幅度的影响.利用（δβ-δρ）扫描叠加方法,对理论接收函数进行了数值试验,结果证明了该方法的可行性.最后将该方法应用于位于青藏高原东北缘的高台（GTA）台和兰州（LZH）台,确定了两个台站下方Moho面的速度跃变分别约为（19±1）%和（20±1）%,密度跃变最小值为（4±2）%和（6±2）%.

     

一种新的拟合地震动衰减关系的分步回归法

提出了拟合地震动衰减关系的一种新的分步回归法.该方法首先确定近场饱和项,再获得距离衰减项,最后确定震级项.作为算例,采用NGA项目美国西部地区的强震观测资料,对基岩场地自由地表加速度反应谱平台高度值衰减关系进行了回归.结果表明,采用此方法对实际资料的拟合不存在系统偏差,结果稳定,能反映强震地震动在近场的距离饱和与震级饱和特性.对远场数据考虑截尾效应,能使回归结果正确反映地震动的远场衰减特征.     

Fuzzy chance-constrained linear fractional programming approach for optimal water allocation

A fuzzy chance-constrained linear fractional programming method was developed for agricultural water resources management under multiple uncertainties. This approach improved upon the previous programming methods, and could reflect the ratio objective function and multiple uncertainties expressed as probability distributions, fuzzy sets, and their combinations. The proposed approach is applied to an agricultural water resources management system where many crops are considered under different precipitation years. Through the scenarios analyses, the multiple alternatives are presented. The solutions show that it is applicable to practical problems to address the crop water allocation under the precipitation variation and sustainable development with ratio objective function of the benefit and the irrigation amount. It also provides bases for identifying desired agriculture water resources management plans with reasonable benefit and irrigation schedules under crops.     

联合多代卫星测高和多源重力数据的局部大地水准面精化方法

本文研究了基于泊松小波径向基函数融合多代卫星测高及多源重力数据精化大地水准面模型的方法.分别以沿轨垂线偏差和大地水准面高高差作为卫星测高观测量,研究了使用不同类型测高数据对于大地水准面建模精度的影响.针对全球潮汐模型在浅水区域及部分开阔海域精度较低的问题,引入局部潮汐模型研究了不同潮汐模型对于大地水准面的影响.数值分析表明:相比于使用沿轨垂线偏差作为测高观测量,基于沿轨大地水准面高高差解算得到的大地水准面模型的精度更高,特别是在海域区域,其精度提高了2.3 cm.由于使用沿轨大地水准面高高差作为测高观测量削弱了潮汐模型长波误差的影响,采用不同潮汐模型对大地水准面解算的影响较小.总体而言,船载重力及测高观测数据在海洋重力场的确定中呈现互补性关系,联合两类重力场观测量可以提高局部重力场的建模精度.     

A Multiple Criteria Decision Modelling approach to selection of estimation techniques for fitting extreme floods

The problem of fitting a probability distribution, here log-Pearson Type III distribution, to extreme floods is considered from the point of view of two numerical and three non-numerical criteria. The six techniques of fitting considered include classical techniques (maximum likelihood, moments of logarithms of flows) and new methods such as mixed moments and the generalized method of moments developed by two of the co-authors. The latter method consists of fitting the distribution using moments of different order, in particular the SAM method (Sundry Averages Method) uses the moments of order 0 (geometric mean), 1 (arithmetic mean), –1 (harmonic mean) and leads to a smaller variance of the parameters. The criteria used to select the method of parameter estimation are:

Comparison of spectral and spatial methods for a Moho recovery from gravity and vertical gravity-gradient data

In global studies investigating the Earth’s lithospheric structure, the spectral expressions for the gravimetric forward and inverse modeling of the global gravitational and crustal structure models are preferably used, because of their numerical efficiency. In regional studies, the applied numerical schemes typically utilize the expressions in spatial form. Since the gravity-gradient observations have a more localized support than the gravity measurements, the gravity-gradient data (such as products from the Gravity field and steady-state Ocean Circulation Explorer - GOCE - gravity-gradiometry satellite mission) could preferably be used in regional studies, because of reducing significantly the spatial data-coverage required for a regional inversion or interpretation. In this study, we investigate this aspect in context of a regional Moho recovery. In particular, we compare the numerical performance of solving the Vening Meinesz-Moritz’s (VMM) inverse problem of isostasy in spectral and spatial domains from the gravity and (vertical) gravity-gradient data. We demonstrate that the VMM spectral solutions from the gravity and gravity-gradient data are (almost) the same, while the VMM spatial solutions differ from the corresponding spectral solutions, especially when using the gravity-gradient data. The validation of the VMM solutions, however, reveals that the VMM spatial solution from the gravity-gradient data has a slightly better agreement with seismic models. A more detailed numerical analysis shows that the VMM spatial solution formulated for the gravity gradient is very sensitive to horizontal spatial variations of the vertical gravity gradient, especially in vicinity of the computation point. Consequently, this solution provides better results in regions with a relatively well-known crustal structure, while suppressing errors caused by crustal model uncertainties from distant zones. Based on these findings we argue that the gravity-gradient data are more suitable than the gravity data for a regional Moho recovery.     

A Crustal Model for the Eastern Alps Region and a New Moho Map in Southeastern Europe

In the last two decades, south-central Europe and the Eastern Alps have been widely explored by many seismic refraction experiments (e.g., CELEBRATION 2000, ALP 2002, SUDETES 2003). Although quite detailed images are available along linear profiles, a comprehensive, three-dimensional crustal model of the region is still missing. This limitation makes this region a weak spot in continental-wide comprehensive representations of crustal structure. To improve on this situation, we select and collect 37 published active-source seismic lines in this region. After geo-referencing each line, we sample them along vertical profiles—every 50?km or less along the line—and derive P-wave velocities in a stack of homogeneous layers (separated by discontinuities: depth of crystalline basement, top of lower crust, and Moho). We finally merge the information using geostatistical methods, and infer S-wave velocity and density using empirical scaling relations. We present here the resulting crustal model for a region encompassing the Eastern Alps, Dinarides, Pannonian basin, Western Carpathians and Bohemian Massif, covering the region within $45^{\circ}\text{--}51^{\circ}\hbox{N}$ and $11^{\circ} \text{--} 22^{\circ}\hbox{E}$ with a resolution of $0.2^{\circ} \times 0.2^{\circ}.$ We are also able to extend and update the map of Moho depth in a wider region within $35^{\circ}\text{--}51^{\circ}\hbox{N}$ and $12^{\circ}\text{--}45^{\circ}\hbox{E},$ gathering Moho values from the collected seismic lines, other published dataset and using the European plate reference EPcrust as a background. All the digitized profiles and the resulting model are available online.     

High spatial resolution mapping of surface velocities and depths for shallow overland flow

Point measurements of flow rate, depth or velocity are not sufficient to validate overland flow models, particularly when the interaction of the water with the soil surface creates a complex flow geometry. In this study, we present the coupling of two techniques obtaining spatial data of flow depths and surface velocity measurements for water depths as low as 1 mm. Overland flow experiments were performed in the laboratory at various flow rates and slopes on two surfaces. The first surface was 120 cm by 120 cm showing three undulations of sinusoidal shape with an amplitude of 1 cm and a wavelength of 20 cm, while the second was a 60 cm by 60 cm moulded reproduction of a seedbed with aggregates up to 2 cm in size. Large scale particle image velocimetry (LSPIV) was used for velocity measurements with a sub‐centimetre spatial resolution. An instantaneous‐profile laser scanner was used to map flow depths with a sub‐millimetre spatial resolution. A sensitivity analysis of the image processing of the LSPIV showed good robustness of the method. Comparison with measurements performed with hot film anemometer and salt velocity gauge showed that LSPIV surface velocities were representative of the flow. Water depths measured with the laser scanner were also in good agreement with single‐point measurements performed with a dial indicator. Spatially‐distributed flow rates could be computed by combining both presented techniques with a mean relative error less than 20%. Copyright © 2012 John Wiley & Sons, Ltd.     

IOD-related optimal initial errors and optimal precursors for IOD predictions from reanalysis data

This study explored the spatial patterns of winter predictability barrier (WPB)-related optimal initial errors and optimal precursors for positive Indian Ocean dipole (IOD) events, and the associated physical mechanisms for their developments were analyzed using the Simple Ocean Data Assimilation dataset. Without consideration of the effects of model errors on “predictions,” it was assumed that different “predictions” are caused by different initial conditions. The two types of WPB-related optimal initial errors are almost opposite for the start months of July (–1) and July (0), although they both present a west-east dipole pattern in the tropical Indian Ocean, with the maximum errors located at the thermocline depth. Bjerknes feedback and ocean waves play important roles in the growth of prediction errors. These two physical mechanisms compete during July–December and ocean waves dominate during January–June. The spatial patterns of optimal precursors and the physical mechanisms for their developments are similar to those of WPB-related optimal initial errors. It is worth noting that large values of WPB-related optimal initial errors and optimal precursors are concentrated within a few locations, which probably represent the sensitive areas of targeted observations for positive IOD events. The great similarities between WPB-related optimal initial errors and optimal precursors suggest that were intensive observations performed over these areas, this would not only reduce initial errors and thus, prediction errors, but it would also permit the detection of the signal of IOD events in advance, greatly improving the forecast skill of positive IOD events.     

Focal depths for moderate-sized aftershocks of the Wenchuan M_S8.0 earthquake and their implications

Sliding-window cross-correlation method is firstly adopted to identify sPn phase, and to constrain focal depth from regional seismograms, by measuring the time separation between sPn and Pn phases. We present the focal depths of the 17 moderate-sized aftershocks (MS≥5.0) of the Wenchuan MS8.0 earth-quake, using the data recorded by the regional seismic broadband networks of Shaanxi, Qinghai, Gansu, Yunnan and Sichuan. Our results show focal depths of aftershocks range from 8 to 20 km, and tend to cluster at two average depths, separate at 32.5°N, i.e., 11 km to the south and 17 km to the north, indicating that these aftershocks are origin of upper-to-middle crust. Combined with other results, we suggest that the Longmenshan fault is not a through-going crustal fault and the Pingwu-Qingchuan fault may be not the northward extension of the Longmenshan thrust fault.     

Improved tilt-depth method for fast estimation of top and bottom depths of magnetic bodies

The tilt-depth method can be used to make fast estimation of the top depth of magnetic bodies. However, it is unable to estimate bottom depths and its every inversion point only has a single solution. In order to resolve such weaknesses, this paper presents an improved tilt-depth method based on the magnetic anomaly expression of vertical contact with a finite depth extent, which can simultaneously estimate top and bottom depths of magnetic bodies. In addition, multiple characteristic points are selected on the tilt angle map for joint computation to improve reliability of inversion solutions. Two- and threedimensional model tests show that this improved tilt-depth method is effective in inverting buried depths of top and bottom bodies, and has a higher inversion precision for top depths than the conventional method. The improved method is then used to process aeromagnetic data over the Changling Fault Depression in the Songliao Basin, and inversion results of top depths are found to be more accurate for actual top depths of volcanic rocks in two nearby drilled wells than those using the conventional tilt-depth method.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

用GPS数据反演海原断裂带断层滑动速率和闭锁深度

以"中国地壳运动观测网络"区域站在海原断裂带附近的所有观测数据及跨断裂GPS剖面观测数据作为约束,用Smith 3-D体力模型反演了海原断裂带断层滑动速率和断层闭锁深度.从西到东断裂共分为5段,采用遗传算法拟合GPS水平运动速度场,拟合残差均方根为1.1 mm/a.反演结果为:毛毛山断裂左旋走滑运动速率为3.5 mm/a,闭锁深度为22.0km;老虎山断裂左旋走滑速率为6.5 mm/a,闭锁深度为10.3 km;海原断裂带西段、中段和东段的滑动速率依次为4.5 mm/a、5.6 mm/a和5.5 mm/a,闭锁深度依次为8.4 km、3.6km和4.3 km.表明毛毛山断裂左旋走滑运动速率小,闭锁深度大,有利于应变能的积累,使得该断裂及附近地区存在发生强震的背景.     

Comparison of Different Methods for a Moho Modeling Under Oceans and Marginal Seas: A Case Study for the Indian Ocean

Surveys in Geophysics - Since marine seismic studies are relatively sparse and unevenly distributed, detailed tomographic images of the Moho geometry under large parts of the world’s oceans...     

Probabilistic seismic hazard model for Cairo,Egypt: estimates and uncertainties

A new seismic hazard model for Cairo, the capital city of Egypt is developed herein based on comprehensive consideration of uncertainties in various components of the probabilistic seismic hazard analysis. The proposed seismic hazard model is developed from an updated catalogue of historical and instrumental seismicity, geodetic strain rates derived from GPS-based velocity-field of the crust, and the geologic slip rates of active faults. The seismic source model consists of area sources and active faults characterised to forecast the seismic productivity in the region. Ground motion prediction models are selected to describe the expected ground motion at the sites of interest. The model accounts for inherent epistemic uncertainties of statistical earthquake recurrence; maximum magnitude; ground motion prediction models, and their propagation toward the obtained results. The proposed model is applied to a site-specific hazard analysis for Kottamiya, Rehab City and Zahraa-Madinat-Nasr (hereinafter referred to as Zahraa) to the East of Cairo (Egypt). The site-specific analysis accounts for the site response, through the parameterization of the sites in terms of average 30-m shear-wave velocity (Vs30). The present seismic hazard model can be considered as a reference model for earthquake risk mitigation and proper resilience planning.     

Rapid in situ assessment of luminescence-bleaching depths for deriving burial and exposure chronologies of rock surfaces

Recent developments in luminescence dating offer new ways to date exposure and burial durations of rocks. The new rock surface dating methods ideally require high-resolution data, faster sample preparation and measurement times, and field screening methods to select samples with appropriate luminescence characteristics and bleaching histories. Presented here is a demonstration of an EMCCD (electron multiplying charge coupled device) based system capable of imaging high-resolution infrared stimulated luminescence (IRSL) and infrared-photoluminescence (IRPL) from rock samples. The IRPL can be detected at both 880 nm and 955 nm. Using this instrument, the entire luminescence-depth profile can be reconstructed by imaging a single surface cut perpendicular to the exposed rock face. We demonstrate the possibility of reconstructing luminescence-depth profiles suitable for rock surface dating from large (cm-scale) rock samples, without using a regeneration dose for normalisation of the natural luminescence signals. Based on the different bleaching characteristics of the IRSL and IRPL emissions at 880 nm and 955 nm, we show that it is possible to gain reliable estimates of bleaching depths from measurement of as few as two images of the IRPL signal (one for each emission), or from measurement of the IRSL decay curve. We thus by-pass laborious sample preparations and the need for a gamma source to estimate the bleaching depth, thereby extending the 2D luminescence-depth profile imaging technique to other laboratories that lack access to ionising radiation facilities. This study also makes a significant progression towards development of a field instrument for in situ relative exposure dating, and sample screening for rock surface burial dating.