Stability and Grid Dispersion of the P-SV 4th-Order Staggered-Grid Finite-Difference Schemes     

Moczo  Peter  Kristek  Jozef  Bystrický  Erik 《Studia Geophysica et Geodaetica》2000,44(3):381-402

Stability and grid dispersion in the P-SV 4th-order in space, 2nd-order in time, displacement-stress staggered-grid finite-difference scheme is investigated in the case of a homogeneous unbounded medium. All results, however, also apply to the velocity-stress and displacement- velocity-stress finite-difference schemes. Independent stability conditions for the P and S waves are obtained by exact separation of equations for the two types of waves. Since the S-wave group velocity can differ from the actual velocity as much as 5% for the sampling ratio 1/5, commonly used in numerical modelling, the sampling of the minimum S wavelength by 6 grid spacings (with the velocity difference not larger than 2.5%) is recommended. Grid dispersion is strongest for a wave propagating in a direction of a coordinate axis and weakest for a wave propagating along a plane diagonal. Grid dispersion in the 4 ^th -order scheme for the sampling ratios s = 1/5 and s = 1/6 is smaller than grid dispersion in the 2 ^nd -order scheme for s = 1/10 and s = 1/12, respectively.  相似文献   

On the uncertainty of the transmission function of the optically thick AGB dust shells     

Miroslav Kocifaj  Jozef Klačka  Thomas Posch 《Astrophysics and Space Science》2008,317(1-2):31-38

The uncertainty of the spectral transmission function due to the nonsphericity of cosmic dust particles is analysed for optically thick C-rich and O-rich dust shells. The transmission function directly prescribes the intensity of radiation that passes through a dust shell. It is shown that nonspherical particles affect the stellar spectra in a different way than a system of equally sized spherical particles. Discrepancies in the stellar spectra for both morphological models (spheres and irregular targets) grow with optical thickness of the dust shell and the spectral behaviour of the optical thickness depends on the particle size distribution function. In particular, two most frequently used distrubutions, a power law and Dirac’s delta function, are considered in the presented numerical runs. Light transmission through C-rich dust shells is strongly influenced by absorption which dominates in case of carbonaceous particles. Irregularly shaped carbonaceous particles reduce the intensity of stellar spectra more efficiently than spherical particles of the same composition (the difference is about 10–30%). MgO particles which may be present in O-rich dust shells are almost pure scatterers, thus the ratio of transmission functions for irregularly shaped particles and spherical ones show specific resonant features (especially at wavelengths below 0.4 μm).  相似文献   

Motion of a Meteoroid Released from an Asteroid     

Peter Vereš  Jozef Klačka  Ladislav Kómar  Juraj Tóth 《Earth, Moon, and Planets》2008,102(1-4):47-51

Evidence of asteroid surface features as regolith grains and larger boulders implies resurfacing possibility due to external forces such as gravitational tidal force during close planet encounters. Motion of a meteoroid released from an asteroid in the gravitational fields of the asteroid and the Earth is modeled. We are interested mainly in a distance between the meteoroid and the asteroid as a function of the time. Applications to Itokawa and some close approaching NEAs are presented.  相似文献   

Efficient Methods to Simulate Planar Free Surface in the 3D 4th-Order Staggered-Grid Finite-Difference Schemes   总被引：1，自引：0，他引：1  

Kristek  Jozef  Moczo  Peter  Archuleta  Ralph J. 《Studia Geophysica et Geodaetica》2002,46(2):355-381

We numerically tested accuracy of two formulations of Levander's (1988) stress-imaging technique for simulating a planar free surface in the 4^th-order staggered-grid finite-difference schemes. We have found that both formulations (one with normal stress-tensor components at the surface, the other with shear stress-tensor components at the surface) require at least 10 grid spacings per minimum wavelength ( _min÷h = 10) if Rayleigh waves are to be propagated without significant grid dispersion in the range of epicentral distances up to 15 _dom ^S.Because interior 4^th-order staggered-grid schemes usually do not require more than 6 grid spacings per minimum wavelength, in the considered range of epicentral distances, it was desirable to find alternative techniques to simulate a planar free surface, which would not require denser spatial sampling than _min÷h = 6. Therefore, we have developed and tested new techniques: 1. Combination of the stress imaging (with the shear stress-tensor components at the surface) with Rodrigues' (1993) vertically refined grid near the free surface. 2. Application of the adjusted finite-difference approximations to the z-derivatives at the grid points at and below the surface that uses no virtual values above the surface and no stress imaging. The normal stress-tensor components are at the surface in one formulation, while the shear stress-tensor components are at the surface in the other formulation.The three developed formulations give for the spatial sampling _min÷h = 6 results very close to those obtained by the discrete-wavenumber method. Because, however, the technique with the vertically refined grid near the free surface requires 3 times smaller time step (due to the refined grid), the technique with adjusted finite-difference approximations is the most accurate and efficient technique from the examined formulations in the homogeneous halfspace.  相似文献   

Nongravitational forces and meteoroid streams     

Jozef Klačka 《Earth, Moon, and Planets》1994,65(2):191-196

The action of the solar electromagnetic radiation on the moving interplanetary dust particles in its more complete form than the special case known as the Poynting-Robertson effect is theoretically discussed in application to meteoroid stream of comet Encke.Normal and transversal components of the perturbing nongravitational force are used due to the action of the solar electromagnetic radiation. It is shown that the normal component of the force is negligible. However, transversal component is very important: it can probably completely explain all the observed meteoroid streams situated along the orbit of comet Encke (and, possibly, some asteroids) as the product of the comet Encke alone. Much shorter time is required for producing such a meteoroid stream than is a general conception.If the idea about the significance of the transversal component of the nongravitational force (may be, not produced by electromagnetic radiation) is correct, it may have important consequences for our understanding of ageing of comets, global evolution of the cometary (and, partially, asteroidal) system, and, of course, for a long-term evolution of small interplanetary particles.  相似文献   

Estimating baking temperatures in a Roman pottery kiln by rock magnetic properties: implications of thermochemical alteration on archaeointensity determinations     

Simo Spassov  Jozef Hus 《Geophysical Journal International》2006,167(2):592-604

  相似文献   

Rotational overthrusting of the northwestern Himalaya: further palaeomagnetic evidence from the Riasi thrust sheet, Jammu foothills, India     

Chris T. Klootwijk  Madan Lal Sharma  Jozef Gergan  S.K. Shah  B.K. Gupta 《Earth and Planetary Science Letters》1986,80(3-4):375-393

Thermal demagnetization results (316 samples) are presented for the Tertiary succession of the Riasi thrust sheet (Jammu foothills, northwestern Himalaya). Primary and secondary magnetization directions of Murree Group red beds (Miocene to Upper Eocene) sampled northeast of Jammu indicate, for this part of the Riasi thrust sheet, a clockwise rotation over about 45° with respect to the Indian shield since Late Eocene/Early Miocene time. This accords with clockwise rotations of similar magnitude observed in the Panjal Nappe and the Krol Belt, and is interpreted as representative for the northwestern Himalaya. Results from the western part of the Kalakot inlier, sampled northwest of Jammu, i.e. basal Murree claystone (Middle Eocene) and carbonate from the Subathu Group (lower Middle to Lower Eocene), indicate an aberrant 20–25° counterclockwise rotation which is of local importance only. Available observations on rotation of Himalayan thrust sheets with respect to the Indian shield, indicate that the Himalayan Arc has formed through oroclinal bending. This supports Powell and Conaghan's and Veevers et al.'s model of Greater India with large-scale intracontinental underthrusting along the Main Central Thrust beneath the Tibetan Plateau. Minimal magnitudes of underthrusting of 550 km in the Krol Belt and 650 km in the Thakkhola region are concluded. Palaeolatitude observations (herein and in [1[) agree with absolute positioning of the Indian plate based on India-Africa relative movement data fixed to a hotspot frame in the Atlantic Ocean, and with palaeolatitude observations from DSDP cores on the Indian plate. Collision-related secondary magnetic components observed both to the north and to the south of the Indus-Tsangpo Suture zone show palaeolatitudes between the equator and 7°N. Comparison of both datasets indicates that initial contact between Greater India and south-central Asia had been established in the Hindu Kush—Karakorum region by about 60 Ma ago whereas eastwards progressive suturing had advanced to the Lhasa Block segment of the Indus-Tsangpo Suture zone before 50 Ma ago.  相似文献   

Composition of waters and sediments in the Antwerp harbor     

Jozef R. Senten  Roger H. Charlier 《Environmental Geology》1983,5(3):159-165

Industrial expansion in the harbor areas of Antwerp has altered environmental conditions to a great extent. This study examines the relation between the water economy of the port of Antwerp docks and the composition of dock waters, sediments, and interstitial waters, with emphasis on trace metals and chlorides. Some conclusions are reached as to the possible use of the dredged deposits in agriculture.  相似文献   

Squeezing information from regional climate change projections - results from a synthesis of CMIP5 results for south-east Queensland, Australia     

Ian Smith  Jozef Syktus  Clive McAlpine  Kenneth Wong 《Climatic change》2013,121(4):609-619

We present a synthesis of CMIP5 model results for projected rainfall changes for a single region (south-east Queensland, Australia) and note that, as was evident in CMIP3 results, the multi-model mean projected changes for the late 21st century are not statistically significant for any season nor annually. Taking account of the number of statistically significant changes to mean rainfall, we find some evidence favouring a decrease in both spring and annual rainfall, but this is not compelling. In almost all cases the most frequent result is for no significant change. However, if we consider the number of results where there is a statistically significant change in the distributions of rainfall amounts, there appears to be slightly more information available for risk assessment studies. These numbers suggest an increase in the frequency of both wet and dry events during summer and spring, and a shift towards more frequent dry events during winter. There is no evidence for any significant changes to the distributions for either autumn or annually. The findings suggest that, in one respect, multi-model rainfall projections may contain more information than is evident from syntheses which focus on changes to the means and that, for some regions where changes in the frequency of wet and dry seasons/years have known impacts, the model projections may be more valuable than previously thought.  相似文献   

Key structural parameters affecting earthquake ground motion in 2D and 3D sedimentary structures     

Peter Moczo  Jozef Kristek  Pierre-Yves Bard  Svetlana Stripajová  Fabrice Hollender  Zuzana Chovanová  Miriam Kristeková  Deborah Sicilia 《Bulletin of Earthquake Engineering》2018,16(6):2421-2450

Alluvial valleys generate strong effects on earthquake ground motion (EGM). These effects are rarely accounted for even in site-specific studies because of (a) the cost of the required geophysical surveys to constrain the site model, (b) lack of data for empirical prediction, and (c) poor knowledge of the key controlling parameters. We performed 3D, 2D and 1D simulations for six typical sedimentary valleys of various width and depth, and for a variety of modifications of these 6 “nominal models” to investigate sensitivity of EGM characteristics to impedance contrast, attenuation, velocity gradient and geometry. We calculated amplification factors, and 2D/1D and 3D/2D aggravation factors for 10 EGM characteristics, using a representative set of recorded accelerograms to account for input motion variability. The largest values of the amplification and aggravation factors are found for the Arias intensity and cumulative absolute velocity, the lowest for the root-mean-square acceleration. The aggravation factors are largest for the vertical component. For each model, at least one EGM characteristic exhibits a significant 2D/1D aggravation factor, while all EGM characteristics exhibit significant 2D/1D aggravation factor on the vertical component. For all investigated sites, there is always an area in the valley for which 1D estimates are not sufficient. 2D estimates are insufficient at several sites. The key structural parameters are the shape ratio and overall geometry of the sediment-bedrock interface, impedance contrast at the sediment-bedrock interface, and attenuation in sediments. The amplification factors may largely exceed the values that are usually considered in GMPEs between soft soils and rock sites.  相似文献