GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

Similarities between strike-slip faults at different scales and a simple age determining method for active faults

Abstract Several differently scaled strike‐slip faults were examined. The faults shared many geometric features, such as secondary fractures and linkage structures (damage zones). Differences in fault style were not related to specific scale ranges. However, it was recognized that differences in style may occur in different tectonic settings (e.g. dilational/contractional relays or wall/linkage/tip zones), different locations along the master fault or different fault evolution stages. Fractal dimensions were compared for two faults (Gozo and San Andreas), which supports the idea of self‐similarity. Fractal dimensions for traces of faults and fractures of damage zones were higher (D ～1.35) than for the main fault traces (D ～1.005) because of increased complexity due to secondary faults and fractures. Based on the statistical analysis of another fault evolution study, single event movements in earthquake faults typically have a maximum earthquake slip : rupture length ratio of approximately 10^?4, although this has only been established for large earthquake faults because of limited data. Most geological faults have a much higher maximum cumulative displacement : fault length ratio; that is, approximately 10^?2 to 10^?1 (e.g. Gozo, ～10^?2; San Andreas, ～10^?1). The final cumulative displacement on a fault is produced by accumulation of slip along ruptures. Hence, using the available information from earthquake faults, such as earthquake slip, recurrence interval, maximum cumulative displacement and fault length, the approximate age of active faults can be estimated. The lower limit of estimated active fault age is expressed with maximum cumulative displacement, earthquake slip and recurrence interval as T ? (d_max /u) · I(M).     

Investigation of upper crust anisotropy in Ghaen-Birjand region, east-central Iran

IntroductionWhenpropagatingthroughananisotropicmedium,ashearwavesplitsintotwo(quasi)shearwaveswithdifferentpropagationspeedsandpolarizedorthogonally.Owingtotherecentdevel-opmentofseismicobservationsystem,detectionofshearwavessplittingwithverysmalldelaytimesbetweenfasterandslowershearwavesbecameavailableandprovidedpowerfulapproachfordetectionofcrustalanisotropy.Crampin(1978)emphasizedtheroleofalignedmicrocracksasacauseofcrustalanisotropyandpointedoutthatforverticallyalignedmicrocracksthedirecti…     

红砂岩高陡岩质边坡稳定性分析及治理

充分利用边坡失稳前自身抗滑能力进行高陡边坡治理，结合工程实例，提出了主动性公路边坡治理的思路。分别利用极限平衡法和有限元数值法，综合考虑了两种方法的优缺点，对红砂岩高陡边坡进行了全面的稳定性分析，并针对性地提出了经济合理的治理方案。     

黄河流域分布式水文模型初步研究与进展

结合国家重点基础研究规划项目"黄河流域水资源演化规律与可再生性维持机理"的研究,综述了黄河流域分布式水文模型研究的主要进展,包括分布式水文模型的特点和基本结构、分布式水文模型中的数字高程模型(DEM)和水文学基础、地理信息系统(GIS)和遥感(RS)在分布式水文模拟中的应用等方面的内容,对未来分布式水文模型的研究提出几点认识。     

基于GIS/RS的流域水文过程分布式模拟——Ⅰ模型的原理与结构

面向流域水资源管理,提出了一个基于GIS/RS的流域分布式水文模型,模型主要包括单元水文模型与河网汇流模型两大部分。单元水文模型涉及到冠层截留、融雪、蒸散发、坡面流、非饱和土壤水运动和地下水出流等水文物理过程。产流计算考虑到地形坡度的影响采用基于地形指数的计算方法。汇流演算基于河网结构采用分段马斯京根方法。模型的大部分参数与输入信息可以利用GIS和RS技术获取,能够对气候变化和人类活动对下垫面的改变,做出快速的模拟与响应。     

基于拉格朗日插值法修正地形影响的分布式降水模型研究

降水量一直是进行水文分析计算的输入项，对其分布状况的模拟直接影响水文分析成果精度。降水中心位置及其中心降水量对暴雨分析尤为重要，而目前尚未见对降水中心位置方面的模型研究。在对目前国内外降水模型分析的基础上，根据天气系统降水如不受地形影响其降水量等值线在平面上的分布近似为一组同心椭圆这一原理，建立了一种能够模拟次降水过程的降水中心位置及其中心降水量的新型分布式降水数学模型，并对其进行地形影响因素修正。由于模型建立原理简单，易于实现对流域未设站研究点的实时降水量估计，同时由于模型能够指明降水中心位置及其中心降水量，因此在流域暴雨分析和洪水预报中具有实用价值。模型经实践检验具有较高精度。     

莺歌海盆地断裂活动与地质灾害分析

莺歌海盆地是南海西北部海域地质灾害类型较齐全的地区之一。根据野外调查资料,本文简要地分析了该区的活动断裂和地质灾害类型。区内活动断裂主要走向呈NW向,该方向的活动断裂主要有红河断裂带、中央拗陷东部断裂带和中央拗陷西部断裂带。红河断裂带在地震资料上显示为一倾向SW的正断层,表现为明显的重力梯级带和磁异常分界带,NE向断裂带分布于盆地中部,在重磁异常上表现为明显的异常带,在北部表现为被NW向的断裂切割。莺歌海盆地地质灾害类型主要有地震活动、泥底辟构造、滑坡与崩塌等。同时,还对该区地质灾害基本特征和分布规律进行了初步分析。     

Imaging the ramp–décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake

We use coseismic GPS data from the 1999 Chi-Chi, Taiwan earthquake to estimate the subsurface shape of the Chelungpu fault that ruptured during the earthquake. Studies prior to the earthquake suggest a ramp–décollement geometry for the Chelungpu fault, yet many finite source inversions using GPS and seismic data assume slip occurred on the down-dip extension of the Chelungpu ramp, rather than on a sub-horizontal décollement. We test whether slip occurred on the décollement or the down-dip extension of the ramp using well-established methods of inverting GPS data for geometry and slip on faults represented as elastic dislocations. We find that a significant portion of the coseismic slip did indeed occur on a sub-horizontal décollement located at 8 km depth. The slip on the décollement contributes 21% of the total modeled moment release. We estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth: homogeneous, layered, and layered with lateral material contrast across the fault. It is shown, however, that heterogeneity has little influence on our estimated fault geometry. We also investigate several competing interpretations of deformation within the E/W trending rupture zone at the northern end of the 1999 ground ruptures. We demonstrate that the GPS data require a 22- to 35-km-long lateral ramp at the northern end, contradicting other investigations that propose deformation is concentrated within 10 km of the Chelungpu fault. Lastly, we propose a simple tectonic model for the development of the lateral ramp.     

Neotectonic deformation in the transition zone between the Dead Sea Transform and the East Anatolian Fault Zone, Southern Turkey: a palaeomagnetic study of the Karasu Rift Volcanism

In southern Turkey ongoing differential impingement of Arabia into the weak Anatolian collisional collage resulting from subduction of the Neotethyan Ocean has produced one of the most complex crustal interactions along the Alpine–Himalayan Orogen. Several major transforms with disputed motions, including the northward extension of the Dead Sea Fault Zone (DSFZ), meet in this region. To evaluate neotectonic motion on the Amanos and East Hatay fault zones considered to be northward extensions of the DSFZ, the palaeomagnetism of volcanic fields in the Karasu Rift between these faults has been studied. Remanence carriers are low-Ti magnetites and all except 5 of 51 basalt lavas have normal polarity. Morphological, polarity and K–Ar evidence show that rift formation occurred largely during the Brunhes chron with volcanism concentrated at 0.66–0.35 Ma and a subsidiary episode at 0.25–0.05. Forty-four units of normal polarity yield a mean of D/I=8.8°/54.7° with inclination identical to the present-day field and declination rotated clockwise by 8.8±4.0°. Within the 15-km-wide Hassa sector of the Karasu Rift, the volcanic activity is concentrated between the Amanos and East Hatay faults, both with left lateral motions, which have rotated blocks bounded by NW–SE cross faults in a clockwise sense as the Arabian Block has moved northwestwards. An average lava age of 0.5 Ma yields a minimum cumulative slip rate on the system bounding faults of 0.46 cm/year according with the rate deduced from the Africa–Arabia Euler vector and reduced rates of slip on the southern extension of the DSFZ during Plio-Quaternary times. Estimates deduced from offsets of dated lavas flows and morphological features on the Amanos Fault Zone [Tectonophysics 344 (2002) 207] are lower (0.09–0.18 cm/year) probably because they are limited to surface fault breaks and do not embrace the seismogenic crust.Results of this study suggest that most strike slip on the DSFZ is taken up by the Amanos–East Hatay–Afrin fault array in southern Turkey. Comparable estimates of Quaternary slip rate are identified on other faults meeting at an unstable FFF junction (DSFZ, East Anatolian Fault Zone, Karatas Fault Zone). A deceleration in slip rate across the DSFZ and its northward continuation during Plio-Quaternary times correlates with reorganization of the tectonic regime during the last 1–3 Ma including tectonic escape within Anatolia, establishment of the North and East Anatolian Fault Zones bounding the Anatolian collage in mid–late Pliocene times, a contemporaneous transition from transpression to transtension and concentration of all basaltic magmatism in this region within the last 1 Ma.