Analysis of A New Model for Unsaturated Flow in Porous Media Including Hysteresis and Dynamic Effects     

A.Y. Beliaev  R.J. Schotting 《Computational Geosciences》2002,5(4):345-368

A new model for unsaturated flow in porous media, including capillary hysteresis and dynamic capillary effects, is analyzed. Existence and uniqueness of solutions are established and qualitative and quantitative properties of (particular) solutions are analyzed. Some results of numerical computations are given. The model under consideration incorporates simple play-type hysteresis and a dynamic term (time-derivative with respect to water content) in the capillary relation. Given an initial water content distribution, the model determines which parts of the flow domain are in drainage and which parts are in imbibition. The governing equations can be recast into an elliptic problem for fluid pressure and an evolution equation for water content. Standard methods are used to obtain numerical results. A comparison is given between J.R. Philip's semi-explicit similarity solution for horizontal redistribution in an infinite one-dimensional domain and solutions of the new model.  相似文献   

ON BROADENING MECHANISMS OF THE Ca  XIX RESONANCE LINE IN SOLAR FLARES     

W.Q. Gan  Y.P. Li 《Solar physics》2002,205(1):117-121

There are two explanations for the broadening mechanisms of the Caxix resonance line in the rising phase of a solar flare. One is due to the velocity dispersion of the chromospheric evaporation along the loop; the other is due to the non-thermal turbulence resulting from the flare energy release. In this paper we distinguish between these two possibilities by studying the influence of the loop orientation on the broadening of the Caxix resonance line, based on the data observed with SXT and BCS on Yohkoh. The results seem to support that the broadening of the Caxix resonance line is caused by non-thermal turbulence, rather than by the velocity dispersion.  相似文献   

Electron–Whistler Interaction in Coronal Loops and Radiation Signatures     

Stepanov  A.V.  Tsap  Y.T. 《Solar physics》2002,211(1-2):135-154

Interaction of the 30–300 keV electrons with whistlers in solar coronal loops is studied using a quasi-linear approach. We show that the electron–whistler interaction may play a dominant role in the formation of fast electron spectra within the solar flare loops with the plasma temperature 10⁷ K and plasma density 10¹¹ cm^–3. It is found that Landau damping of whistlers provides weak and intermediate pitch-angle diffusion regimes of fast electrons in coronal loops. The level of whistler turbulence in the weak diffusion regime under flare conditions is estimated as 10^–7 of the energy density in the thermal particles. The `top – footpoint' relations between the hard X-ray flux densities and spectra are derived. The reason for a `broken' spectrum of the flare microwave emission is discussed.  相似文献   

The Generating Region of Bidirectional Electron Beams in the Corona   总被引：1，自引：0，他引：1  

Xie  R.X.  Fu  Q.J.  Wang  M.  Liu  Y.Y. 《Solar physics》2000,197(2):375-385

Metric and decimetric type III bursts and microwave spike emissions with negative and positive frequency drift rates which were observed with radio spectrometers at Yunnan and Beijing Observatories are presented. The frequencies and heights at which the bidirectional electron beams originated are estimated. Three events reveal a separatrix frequency (at 250, 1300, and 2900 MHz) between normal- and reverse-drifting radio bursts, indicating a compact acceleration source where electron beams are injected in both upward and downward directions. These cases may indicate that the changeover frequencies of bidirectional electron beams are within a large band from 250 to 2900 MHz and the frequency bands of separatrices are in very small (4 to 100 MHz) and different bands. These type III bursts appear to be a plasma emission phenomenon from a beam of electrons which seem to have widely separated acceleration regions from the high to the low corona. These cases suggest that current sheets that separate open and closed magnetic fluxes in the low corona, and oppositely directed open field lines in the high corona are possible sites for bidirectional electron acceleration. The regions of magnetic topology from closed to open magnetic field structures should be very large (from about 20000 to 107000 km above the photosphere).  相似文献   

The Control Of Coherent Eddies In Vegetation Canopies: Streamwise Structure Spacing, Canopy Shear Scale And Atmospheric Stability   总被引：9，自引：9，他引：0  

Y. Brunet  M. R. Irvine 《Boundary-Layer Meteorology》2000,94(1):139-163

An analogy has been established between a plane mixing layer and the atmospheric flow near the top of a vegetation canopy. It is based on a common feature, a strong inflection in the mean velocity profile, responsible for hydrodynamical instabilities that set the pattern for the coherent eddies and determine the turbulence length scales. In an earlier study, this analogy was tested using a small data set from thirteen experiments, all in near-neutral conditions. It provided a good prediction of the streamwise spacing _w of the dominant canopy eddies (evaluated from time series of vertical velocity) that appears to depend on a shear length scale L_s = U(h)/U'(h), where h is canopy height, U is mean velocity and U' the vertical gradient dU/dz. The present analysis utilizes an extensive data set of approximately 700 thirty-minute runs, from six experiments on two forest sites and a maize crop, with a large range of stability conditions. _w was estimated for each run using the wavelet transform as an objective, automated detection method. First, the variations of _w and L_s with atmospheric stability are discussed. Neutral and unstable values exhibit a large scatter whereas in stable conditions both variables decrease with increasing stability. It is subsequently found that _w is directly related to L_s, in a way close to the neutral prediction _w /h = 8.1L_s/h.The Strouhal number S_tr = L_s /_w is then shown to vary with atmospheric stability, weakly in unstable conditions, more significantly in stable conditions. Altogether these results suggest that, to some extent, the plane mixing-layer analogy can be extended to non-neutral conditions. It is argued that the primary effect of atmospheric stability, at least in stable conditions, is to modify the shear length scale L_s through changes in U(h) and U'(h), which in turn determines the streamwise spacing of the active, coherent motions.  相似文献   

The Gravity Field of the Vogtland and NW Bohemia: Presentation of a Project     

Hofmann  Y.  Jahr  T.  Jentzsch  G.  Bankwitz  P.  Bram  K. 《Studia Geophysica et Geodaetica》2000,44(4):608-610

The Vogtland and NW Bohemia are characterized geoscientifically by periodically occurrence of swarm earthquakes. The basic geophysical mechanism is not yet sufficiently clarified, just like detail questions to geology in especially the deeper underground. Complex geophysical investigations in the seismoactive region indicate geodynamic phenomena like mass redistribution or stress accumulation and release (Spiák et al., 1998). According to Grünthal (1989) a weakness zone is suggested in the region of the swarm earthquakes. This zone can be caused by fluid-tectonics (Kämpf et al., 1992), a mantle plume (pers. com. J. Svancara, 1999) and/or by the geometry of the geological structures (Neunhöfer & Güth, 1988). A three-dimensional gravimetric model can clear up the underground situation. By means of high-resolution gravimetry a three-dimensional model will be developped for the Vogtland and NW Bohemia region. In the first step a homogeneous Bouguer map of the Vogtland and NW Bohemia was created (fig. 1) containing gravity structures analysed by Ibrmajer & Suk (1989) and Blízkovsky et al. (1985). The used gravimetric data were made available by the Saxonian National Office for Environment and Geology, by the Czech Geological Survey, Prague and by the GGA Hannover. In the context with the interpretation of the deep-seismic profile MVE 90 a two-dimensional gravimetric modeling was carried out (Behr et al., 1994), too. Anomaly-producing source bodies apparently do not offer themselves in a two-dimensional model, because after Jung (1961) the length of a gravimetric source structure must be about four times larger than it's width. The technique of the three-dimensional gravimetric modeling by means of any polyhedrons was developed by Götze (1976, 1984). Gravimetry is a potential method and supplies an infinite number of solutions, so the model has to be developed close to other geoscientific results. The aim is to construct a high-resolution three-dimensional underground model, which includes the upper earth's crust and the deep-seated structures of the middle and lower crust, too. The determination of the mass distribution in the underground supplies contradicting or supporting facts for geodynamic views in the Vogtland and NW Bohemia for example of Bankwitz et al. (1993). The interpretation of the Bouguer map of the Vogtland and a three-dimensional gravimetric model ought to contribute a substantial, also geodynamic part to understand the origin and the emergence of the swarm earthquakes in this region.  相似文献   

Straightforward methods to detect non-linear response of the soil. Application to the recordings of the Kobe earthquake (Japan, 1995)     

C. Lacave-Lachet  P.Y. Bard  J.C. Gariel  K. Irikura 《Journal of Seismology》2000,4(2):161-173

Simple straightforward methods are applied to testtheir ability to detect the non-linear response of thesoil. Recordings of the main shock and aftershocks ofthe 1995, Hyogo-ken Nanbu (Kobe) earthquake are used.Non-linear effects are investigated using twodifferent techniques, on a collection of data for 12sites situated on different geological structures inthe Kobe and Osaka areas. The first method used is theso-called receiver functions technique (Langston,1979), which consists of computing the spectral ratiobetween horizontal and vertical components of motion.This ratio has been shown to reveal the fundamentalfrequency of a site (Lachet and Bard, 1994; Lachet etal., 1996; Theodulidis et al., 1995, 1996). For eachsite, recordings of the main shock and a set ofaftershocks are considered. The variation of thisspectral ratio for different values of the maximumacceleration recorded at a site is investigated. Bothvariations of the amplitude of the H/V ratio (due tonon-linear behavior, on the horizontal components inparticular) and of the frequency position of theamplified band-width are observed. The secondtechnique used in this study is related to thevariation of the high frequency content of therecordings during the main-shock and its aftershocks.The high frequency spectral decay of the motion,characterized by parameter, is assumed to berelated mainly to the near-surface attenuation. Itshould then increase with increasing peak velocity, incase of non-linearity. The value of kappa iscalculated for the 12 sites in the Kobe area, fordifferent types of soil conditions, and againdifferent values of peak ground acceleration.Variations of kappa are then related to non-linearbehavior of the soil during the Kobe earthquake.  相似文献   

Hydrogeochemical Precursors in Kamchatka (Russia) Related to the Strongest Earthquakes in 1988–1997     

Biagi  P. F.  Ermini  A.  Cozzi  E.  Khatkevich  Y. M.  Gordeev  E. I. 《Natural Hazards》2000,21(2-3):263-276

The Kamchatka peninsula, located in the far east of Russia, is a geologically active margin where the Pacific plate subducts beneath the North American and Eurasia plates. This area is characterised by frequent and strong seismic activity (magnitudes up to 8.5) and epicentres are generally distributed offshore along the eastern coast of the peninsula. For many years, hydrogeochemicals have been collected with a mean sampling frequency of three days in the form of the flow rate and the most common ions and gases in the groundwater of three deep wells in the southern area of the Kamchatka peninsula, where the capital city Petropavlovsk is located. Beginning in 1988, five earthquakes with M > 6.5 occurred in this area. These earthquakes were powerful enough to be considered as potential precursor sources in the sense that the stresses and strains building up before them might be expected to cause precursory activity. In order to reveal any possible precursors of these earthquakes, we analysed the hydrogeochemical data collected. We considered any signal having an amplitude three times the standard deviation to be an irregularity and we defined as an anomaly the existence of an irregularity occurring simultaneously in more than one parameter at each well. Then, on the basis of the worldwide past results and the time recurrence of the quoted earthquakes, we chose 158 days as the maximum temporal window between a possible anomaly and the subsequent earthquake. We identified some premonitory anomalies in hydrogeochemical parameters at different wells. On the basis of these results some earthquake forecasting criteria in southern Kamchatka may be tentatively formulated  相似文献   

An interval-parameter multi-stage stochastic programming model for water resources management under uncertainty   总被引：1，自引：0，他引：1  

Y.P. Li  G.H. Huang  S.L. Nie 《Advances in water resources》2006

In this study, an interval-parameter multi-stage stochastic linear programming (IMSLP) method has been developed for water resources decision making under uncertainty. The IMSLP is a hybrid methodology of inexact optimization and multi-stage stochastic programming. It has three major advantages in comparison to the other optimization techniques. Firstly, it extends upon the existing multi-stage stochastic programming method by allowing uncertainties expressed as probability density functions and discrete intervals to be effectively incorporated within the optimization framework. Secondly, penalties are exercised with recourse against any infeasibility, which permits in-depth analyses of various policy scenarios that are associated with different levels of economic consequences when the promised water-allocation targets are violated. Thirdly, it cannot only handle uncertainties through constructing a set of scenarios that is representative for the universe of possible outcomes, but also reflect dynamic features of the system conditions through transactions at discrete points in time over the planning horizon. The developed IMSLP method is applied to a hypothetical case study of water resources management. The results are helpful for water resources managers in not only making decisions of water allocation but also gaining insight into the tradeoffs between environmental and economic objectives.  相似文献   

Effects of variable attribute weights on landform classification     

Y. X. Deng  J. P. Wilson  J. Sheng 《地球表面变化过程与地形》2006,31(11):1452-1462

This paper focuses on the attribute weight issue and advocates use of modi?able attribute weights in terrain‐based environmental analysis and classi?cation. A question was asked: ‘How much will the result of a terrain‐based environmental analysis be affected if the weights of used terrain attributes are changed?’ The literature on landform classi?cation and the fuzzy k‐means method was reviewed in particular to help clarify the background and importance of this weight assignment issue. As an example, the effects of modifying attribute weights were evaluated for fuzzy k‐means landform classi?cation in a case study area. A total of 102 classi?cations were compared with each other and with a soil map, and comparison methods were speci?cally designed to evaluate the differences between these classi?cations. The results show that fuzzy k‐means landform classi?cation is sensitive to weight adjustments of adopted terrain attributes. The sensitivity is particularly high when the attribute weights started to be tuned away from the standard (i.e. uniform) weight of one. Better matching between landform classi?cation and a soil map may be produced when attribute weights are tuned. In all, we advocate the widespread adoption of an exploratory attitude in assigning attribute weights for environmental analysis and classi?cation. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献