Age estimates of coastal terraces in the Andaman and Nicobar Islands and their tectonic implications     

Kusala Rajendran  C.P. Rajendran  Anil Earnest  G.V. Ravi Prasad  K. Dutta  D.K. Ray  R. Anu 《Tectonophysics》2008,455(1-4):53-60

The great Indian Ocean earthquake of December 26, 2004 caused significant vertical changes in its rupture zone. About 800 km of the rupture is along the Andaman and Nicobar Islands, which forms the outer arc ridge of the subduction zone. Coseismic deformation along the exposed land could be observed as uplift/subsidence. Here we analyze the morphological features along the coast of the Andaman and Nicobar Islands, in an effort to reconstruct the past tectonics, taking cues from the coseismic effects. We obtained radiocarbon dates from coastal terraces of the island belt and used them to compute uplift rates, which vary from 1.33 mm yr^− 1 in the Little Andaman to 2.80 mm yr^− 1 in South Andaman and 2.45 mm yr^− 1 in the North Andaman. Our radiocarbon dates converge on  600 yr and  1000 yr old coastal uplifts, which we attribute to the level changes due to two major previous subduction earthquakes in the region.  相似文献   

Integration and magnitude homogenization of the Egyptian earthquake catalogue   总被引：2，自引：0，他引：2  

H. M. Hussein  K. M. Abou Elenean  I. A. Marzouk  A. Peresan  I. M. Korrat  E. Abu El-Nader  G. F. Panza  M. N. El-Gabry 《Natural Hazards》2008,47(3):525-546

The aim of the present work is to compile and update a catalogue of the instrumentally recorded earthquakes in Egypt, with uniform and homogeneous source parameters as required for the analysis of seismicity and seismic hazard assessment. This in turn requires a detailed analysis and comparison of the properties of different available sources, including the distribution of events with time, the magnitude completeness, and the scaling relations between different kinds of magnitude reported by different agencies. The observational data cover the time interval 1900–2004 and an area between 22°–33.5° N and 25°–36° E. The linear regressions between various magnitude types have been evaluated for different magnitude ranges. Using the best linear relationship determined for each available pair of magnitudes, as well as those identified between the magnitudes and the seismic moment, we convert the different magnitude types into moment magnitudes M _W, through a multi-step conversion process. Analysis of the catalogue completeness, based on the M _W thus estimated, allows us to identify two different time intervals with homogeneous properties. The first one (1900–1984) appears to be complete for M _W ≥ 4.5, while the second one (1985–2004) can be considered complete for magnitudes M _W ≥ 3.  相似文献   

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

Accurate local upscaling with variable compact multipoint transmissibility calculations     

James V. Lambers  Margot G. Gerritsen  Bradley T. Mallison 《Computational Geosciences》2008,12(3):399-416

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.  相似文献   

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data     

Bill X. Hu  Changming He 《Mathematical Geosciences》2008,40(8):845-859

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.  相似文献   

井下记录的PP和SS震相     

乌统昱  赵惠君 《地震地磁观测与研究》1992,13(2):20-26

本文在系统分析静海台井下地震记录图的基础上,发现并解决了具有一定实际和理论意义的两个问题:其一是证明了最初误认为是干扰的“双脉冲”图形实际上均是微地震。从而指出研究静海周围地区的地震活动性应考虑微震活动问题。其二确认静海周围地区微震波形的4个主要震相是P、PP、S和SS,并对地面反射波PP和SS能够观测到的条件作了初步分析,其结果进一步说明井下地震记录较之地面记录更为“逼真”,复杂和丰富。  相似文献   

Seismic potential of Southern Italy   总被引：1，自引：2，他引：1  

S. Jenny  S. Goes  D. Giardini  H.-G. Kahle 《Tectonophysics》2006,415(1-4):81-101

To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.  相似文献   

Earthquake generation in Cyprus revealed by the evolving stress field     

Eleftheria E. Papadimitriou  Vassilios G. Karakostas   《Tectonophysics》2006,423(1-4):61

Strong earthquake occurrence (M ≥ 6.0) onshore and offshore the Cyprus Island constitutes significant seismic hazard because they occur close to populated areas. Seismicity is weak south of the Island along the Cyprean Arc and strong events are aligned along the Paphos transform fault and Larnaka thrust fault zone that were already known and the Lemessos thrust fault zone that defined in the present study. By combining the past history of strong (M ≥ 6.0) events and the long-term tectonic loading on these major fault zones, the evolution of the stress field from 1896 until the present is derived. Although uncertainties exist in the location, magnitude and fault geometries of the early earthquakes included in our stress evolutionary model, the resulting stress field provides an explanation of later earthquake triggering. It was evidenced that the locations of all the strong events were preceded by a static stress change that encouraged failure. The current state of the evolved stress field may provide evidence for the future seismic hazard. Areas of positive static stress changes were identified in the southwestern offshore area that can be considered as possible sites of future seismic activity.  相似文献   

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?     

M. Ruiz  O. Gasp  J. Gallart  J. Díaz  J.A. Pulgar  J. García-Sansegundo  C. Lpez-Fernndez  J.M. Gonzlez-Cortina 《Tectonophysics》2006,424(3-4):223

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.  相似文献   

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area     

P. Tullen  P. Turberg  A. Parriaux 《Engineering Geology》2006,87(3-4):195-204

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.  相似文献