Stability Analysis and Numerical Simulation of Differential Frost Heave     

Rorik A. Peterson 《Mathematical Geosciences》2008,40(3):277-298

Differential frost heave is often implicated in the formation of patterned ground in regions subject to recurrent freezing and thawing. A linear stability analysis (LSA) indicates that a continuum model of frost heave is linearly unstable under typical natural freezing conditions of silty-clay soils. A two-dimensional non-linear numerical analysis corroborates the frozen time LSA results, and also indicates the importance of non-linear and time-dependent terms that ultimately lead to a preferred mode, which the LSA fails to predict. Instability of the one-dimensional solution occurs at shallow freezing depths and near-zero surface loads when positive perturbations in the ice content at the freezing front lead to a concomitant increase in thermomolecular pressure and upward ice velocity. Differential frost heave can then occur because of the increased heat flux from the perturbed surfaces. A three-dimensional model using random initial surface perturbations indicates that regular surface patterns will evolve with a length scale in the order of 2–4 meters, which corresponds quite closely with naturally-occurring non-sorted patterned ground.  相似文献   

Hybrid seismic modeling based on discrete-wave number and finite-difference methods     

Jiří Zahradník  Peter Moczo 《Pure and Applied Geophysics》1996,148(1-2):21-38

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.  相似文献   

Accelerating ground warming in the Canadian Prairie provinces: Is it a result of global warming?     

Jacek A. Majorowicz 《Pure and Applied Geophysics》1996,147(1):1-24

The results of precision temperature logs made to depths of several hundred meters in some 80 wells in Western Canada, most of which are located in the Prairie Provinces, show evidence of warming at the ground surface in the 0.5 K to 3.5 K range (average=2.2±0.7 K, for 80 unevenly distributed sites). Modeling shows that this warming mostly pertains to this century and it has been most substantal in the last four decades if the ramp function of the linear increase of surface temperature is assumed. Using the step function model's increase of surface temperature (land clearing, forest fires, etc.) the calculated onset of warming would pertain mostly to the last two decades. Contour maps of ground temperatures currently and previously and a contour map of the ground warming magnitude dilineate a large regional character of the ground temperature change at the southern marigin of permafrost for the large area of the Prairie Provinces. In many cases however, the magnitude of ground warming is much larger than the magnitude of air warming. This is especially evident for the northern areas of Alberta in the boreal forest ecoprovince. The magnitude of ground warming is equal to the magnitude of surface air warming in southern Alberta in the grassland and aspen parkland ecoprovinces. The analysis of the temperature depth response to the surface warming from well data shows the integrated effect of surface air warming together with the increases in ground temperature due to natural terrain effects and other anthropogenical changes to the surface of the earth.  相似文献   

Energy balance of simple elastodynamic sources     

Václav Vavryčuk 《Pure and Applied Geophysics》1994,143(4):563-586

Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic  相似文献   

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality     

J. Zahradník  P. Moczo  F. Hron 《Natural Hazards》1994,10(1-2):149-170

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.  相似文献   

The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy     

P. D. Cole  J. E. Guest  A. M. Duncan 《Bulletin of Volcanology》1993,55(7):467-480

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.  相似文献   

The arrangement in mean elements space of the periodic orbits close to that of the Moon     

G. B. Valsecchi  E. Perozzi  A. E. Roy  B. A. Steves 《Celestial Mechanics and Dynamical Astronomy》1993,56(1-2):373-380

In a simplified model of the Earth-Moon-Sun system based on the restricted circular 3-dimensional 3-body problem, it is possible to find numerically a set of 8 periodic orbits whose time evolutions closely resemble that of the Moon's orbit. These orbits have a period of 223 synodic months (i.e. the period of the Saros cycle known for more than two millennia as a means of predicting eclipses), and are characterized by a secular rotation of the argument of perigee . Periodic orbits of longer durations exhibiting this last feature are very abundant in Earth-Moon-Sun dynamical models. Their arrangement in the space of the mean orbital elements- for various values of the lunar mean motion is presented.  相似文献   

Generalized Lyapunov exponents indicators in Hamiltonian dynamics: An application to a double star system     

E. Lohinger  C. Froeschlé  R. Dvorak 《Celestial Mechanics and Dynamical Astronomy》1993,56(1-2):315-322

The Lyapunov characteristic numbers (LCNs) which are defined as the mean value of the distribution of the local variations of the tangent vectors to the flow (=ln _k ⁱ ) (see Froeschlé, 1984) have been found to be sensitive indicators of stochasticity. So we computed the distribution of these local variations and determined the moments of higher order for the integrable and stochastic regions in a binary star system with =0.5.  相似文献   

Asteroid motion near the 3 : 1 resonance     

John D. Hadjidemetriou 《Celestial Mechanics and Dynamical Astronomy》1993,56(4):563-599

A mapping model is constructed to describe asteroid motion near the 3 : 1 mean motion resonance with Jupiter, in the plane. The topology of the phase space of this mapping coincides with that of the real system, which is considered to be the elliptic restricted three body problem with the Sun and Jupiter as primaries. This model is valid for all values of the eccentricity. This is achieved by the introduction of a correcting term to the averaged Hamiltonian which is valid for small values of the ecentricity.We start with a two dimensional mapping which represents the circular restricted three body problem. This provides the basic framework for the complete model, but cannot explain the generation of a gap in the distribution of the asteroids at this resonance. The next approximation is a four dimensional mapping, corresponding to the elliptic restricted problem. It is found that chaotic regions exist near the 3 : 1 resonance, due to the interaction between the two degrees of freedom, for initial conditions close to a critical curve of the circular model. As a consequence of the chaotic motion, the eccentricity of the asteroid jumps to high values and close encounters with Mars and even Earth may occur, thus generating a gap. It is found that the generation of chaos depends also on the phase (i.e. the angles andv) and as a consequence, there exist islands of ordered motion inside the sea of chaotic motion near the 3 : 1 resonance. Thus, the model of the elliptic restricted three body problem cannot explain completely the generation of a gap, although the density in the distribution of the asteroids will be much less than far from the resonance. Finally, we take into account the effect of the gravitational attraction of Saturn on Jupiter's orbit, and in particular the variation of the eccentricity and the argument of perihelion. This generates a mixing of the phases and as a consequence the whole phase space near the 3 : 1 resonance becomes chaotic. This chaotic zone is in good agreement with the observations.  相似文献   

地铁施工对济南白泉泉群流量的影响     

束龙仓  王小博  李虎  倪寒茜  李罡  余亚飞  王鑫  张曼琦 《吉林大学学报(地球科学版)》2021,51(1):192-200

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。  相似文献