Traction image method for irregular free surface boundaries in finite difference seismic wave simulation   总被引：16，自引：1，他引：16  

Wei Zhang  Xiaofei Chen 《Geophysical Journal International》2006,167(1):337-353

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.  相似文献   

Shear wave splitting changes associated with the 2001 volcanic eruption on Mt Etna     

Francesca Bianco  Luciano Scarf&#;  Edoardo Del Pezzo  Domenico Patanè 《Geophysical Journal International》2006,167(2):959-967

  相似文献   

H/V ratio: a tool for site effects evaluation. Results from 1-D noise simulations   总被引：9，自引：1，他引：9  

Sylvette Bonnefoy-Claudet  Cécile Cornou  Pierre-Yves Bard  Fabrice Cotton  Peter Moczo  Jozef Kristek  Donat Fäh 《Geophysical Journal International》2006,167(2):827-837

  相似文献   

On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces     

Thorsten W. Becker 《Geophysical Journal International》2006,167(2):943-957

  相似文献   

Three-dimensional structure across the Tintina strike-slip fault, northern Canadian Cordillera, from seismic refraction and reflection tomography     

Colin A. Zelt  Robert M. Ellis  Barry C. Zelt 《Geophysical Journal International》2006,167(3):1292-1308

  相似文献   

Late Holocene linkages between decade–century scale climate variability and productivity at Lake Tanganyika, Africa     

Andrew S. Cohen  Kiram E. Lezzar  Julia Cole  David Dettman  Geoffrey S. Ellis  Meagan Eagle Gonneea  Pierre-Denis Plisnier  Victor Langenberg  Maarten Blaauw  Derrick Zilifi 《Journal of Paleolimnology》2006,36(2):189-209

Microlaminated sediment cores from the Kalya slope region of Lake Tanganyika provide a near-annually resolved paleoclimate record between ∼∼2,840 and 1,420 cal. yr B.P. demonstrating strong linkages between climate variability and lacustrine productivity. Laminae couplets comprise dark, terrigenous-dominated half couplets, interpreted as low density underflows deposited from riverine sources during the rainy season, alternating with light, planktonic diatomaceous ooze, with little terrigenous component, interpreted as windy/dry season deposits. Laminated portions of the studied cores consist of conspicuous dark and light colored bundles of laminae couplets. Light and dark bundles alternate at decadal time scales. Within dark bundles, both light and dark half couplets are significantly thinner than within light bundles, implying slower sediment accumulation rates during both seasons over those intervals.Time series analyses of laminae thickness patterns demonstrate significant periodicities at interannual–centennial time scales. Longer time scale periodicities (multidecadal to centennial scale) of light and dark half couplet thicknesses are coherent and in some cases are similar to solar cycle periods on these time scales. Although laminae thickness cycles do not strongly covary with the actual Δ¹⁴C record for this same time period, two large Δ¹⁴C anomalies are associated with substantial decreases in both light and dark laminae thickness. In contrast to the multidecadal– centennial time scale, significant annual to decadal periodicities, which are broadly consistent with ENSO/PDO forcing and their impact on East African climate, are not coherent between light and dark half couplets. The coherency of light–dark couplets at decadal–centennial time scales, but not at shorter time scales, is consistent with a model of a long-term relationship between precipitation (recorded in wet season dark laminae thickness) and productivity (light laminae thickness), which is not manifest at shorter time scales. We hypothesize that this coupling results from long-term recharging of internal nutrient loading during wet periods (higher erosion of soil P) and reduced loading during drought intervals. The relationship is not expressed on short time scales during which the dominant control on productivity is wind-driven, dry season upwelling, which is uncorrelated with wet-season precipitation. Our record greatly extends the temporal record of this quasi-periodic behavior throughout the late Holocene and provides the first evidence linking decade- to century-scale episodes of enhanced productivity to enhanced precipitation levels and nutrient recharge in a productive tropical lake.  相似文献   

Verification of the motion fields over Southeast Asia by the Florida State University (FSU) superensemble forecast     

Mastura Mahmud  Robert S. Ross 《Singapore journal of tropical geography》2006,27(2):176-190

The Florida State University (FSU) multimodel superensemble forecast is evaluated against several other operational weather models for the Southeast Asia region. The superensemble technique has demonstrated its exceptional skills in forecasting precipitation, motion and mass fields compared to either individual global operational or ensemble mean forecasts. The motion field investigation for the season of 2001 reveals that the superensemble forecasts are closer to the observed data compared to the other global member operational models through its low systematic errors at the 850 hPa level. The FSU multimodel superensemble forecasts exhibit the lowest root mean square errors (RSMEs), the highest correlation against the best observed data and the lowest systematic errors compared to the other operational model members. These forecasts have the potential to provide better daily weather predictions over the Southeast Asia region, particularly during the early northeast monsoon that often causes heavy rainfall in the equatorial part of the Southeast Asia region.  相似文献   

Surface waves image the top of the Eifel plume   总被引：1，自引：0，他引：1  

Jan P. Mathar  Joachim R. R. Ritter  Wolfgang Friederich 《Geophysical Journal International》2006,164(2):377-382

  相似文献   

The detection of low magnitude seismic events using array-based waveform correlation   总被引：12，自引：0，他引：12  

Steven J. Gibbons  Frode Ringdal 《Geophysical Journal International》2006,165(1):149-166

  相似文献   

Computation of Large Anisotropic Seismic Heterogeneities (CLASH)     

Éric Beucler  Jean-Paul Montagner 《Geophysical Journal International》2006,165(2):447-468

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.  相似文献