共查询到20条相似文献，搜索用时 468 毫秒
1.
At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. To enable us to give our ‘measurements in the field’ a more quantitative interpretation than can be deduced from these alone, a series of experiments has been started under well-defined conditions. A cylindrical vessel containing water simulates wet soil. Mounted horizontally above the water surface is a pulsed triangular half-wave dipole which is used as a transmitting antenna (TA). It has a carrier-frequency of about 160 MHz and a pulse repetition-frequency of about 50 kHz. A movable receiving dipole (‘probe’) in the water measures the transverse, mutually orthogonal Eφ,- and Eθ-components of the pulses as a function of probe-position (r, θ, φ) and of the height h of the TA above the water surface. When these pulses are Fourier-transformed, the transverse electric fields Eφ and Eθ at 200 MHz are obtained. The resulting field patterns are compared with computational results on the basis of the theory of the continuous wave, infinitesimal electric dipole (‘point dipole’). It can be concluded that:
- 1 Far-field conditions have not fully developed at a depth of about 2.50 m, the largest value of the radius r at which field patterns were measured, although it represents a distance of about 15 wavelengths.
- 2 The attenuation constant of the tapwater used, as deduced from E-field measurements for θ= 0, 2.50 m < r < 2.75 m, is slightly less than the value measured using a network analyser and air line combination, in agreement with (1).
- 3 E φ field patterns calculated using the value of the conductivity σ corresponding to the former value of the attenuation constant agree reasonably well with the measured patterns for r≤ 2.50 m and for θ < 20° at all antenna heights considered. Calculated Eφ patterns do not agree so well with the measured patterns when h is close to zero. With increasing height the agreement inproves.
- 4 In accordance with the theory of the point-dipole, the angular distribution of the radiation patterns of the TA becomes wider as the frequency decreases.
- 5 The normalized underwater pulse-spectra shift to lower frequencies with increasing r. This can be explained since the attenuation constant of the water rises with rising frequency.
2.
Laboratory experiments using simulated tilting rill channels with different bed roughnesses have been conducted to determine hydraulic conditions for incipient motion of single rock fragments lying on rill beds. These threshold conditions can be described by a modified Shields' entrainment parameter (θm) and relative depth (Z), adjusted for rill conditions. The negative relation found between θm and Z also fits steep (stream) channel data collected by other investigators. Plots of the classical Shields' entrainment parameter (θ) versus both relative depth and relative particle size also showed clear negative trends. Special care has been taken to avoid false relations that can arise between compound parameters. This occurs when the latter are not independent from each other because of their definitions. 相似文献
3.
We analyze the properties of the scattering matrix associated with the incident and scattered electric fields used in GPR. The elements of the scattering matrix provide information produced by different polarizations of the incident wave field. Rotationally invariant quantities such as trace, determinant and Frobenius norm lead to images that combine the information contained in the four elements of the scattering matrix in a mathematically simple and sound manner. The invariant quantities remove the directional properties implicit in the dipolar field used in GPR allowing the application of standard processing techniques designed for scalar fields, such as those used in seismic data processing. We illustrate the non-directional properties of the invariants using a 3D simulation of the wavefield produced by a point scatterer. The estimation of the azimuth angle of elongated targets is also explored using rotation transformations that maximize alternatively the co-polarized or the cross-polarized responses. The angle estimation is essentially an unstable process, particularly if low amplitudes or noisy data are involved. We apply the Frobenius norm 6S6F as a criterion for selection of the best amplitudes to use for a more stable and significant angle estimation. The performance of our formulation was tested with synthetic data produced by a 3D model of an air-filled metal pipe buried in a homogeneous halfspace. The images resulting from the invariants show a clear diffraction hyperbola suitable for a scalar wavefield migration, while the azimuth of the pipe is neatly resolved for amplitudes selected with 6S6F ≥ 0.4. A field experiment conducted above an aqueduct pipe illustrates the proposed methods with real data. The images obtained from the invariants are better than those from the individual elements of the scattering matrix. The azimuth estimated using our formulation is in agreement with the probable orientation of the aqueduct. Finally, a field experiment above a buried air-filled barrel shows that combining the information in the way proposed in this work may lead to an improved image of the subsurface target, the cost to pay is the lost of directional information contained in the scattering matrix. In general, we claim that the methods proposed in this work can be useful to analyze the information acquired by multicomponent GPR surveys using standard scalar wavefield algorithms. 相似文献
4.
AbstractSoil water content (θ) and saturated hydraulic conductivity (Ks) vary in space. The objective of this study was to examine the effects of initial soil water content (θi) and Ks variability on runoff simulations using the LImburg Soil Erosion Model (LISEM) in a small watershed in the Chinese Loess Plateau, based on model parameters derived from intensive measurements. The results showed that the total discharge (TD) and peak discharge (PD) were underestimated when the variability of θi and Ks was partially considered or completely ignored compared with those when the variability was fully considered. Time to peak (TP) was less affected by the spatial variability compared to TD and PD. Except for TP in some cases, significant differences were found in all hydrological variables (TD, PD and TP) between the cases in which spatial variability of θi or Ks was fully considered and those in which spatial variability was partially considered or completely ignored. Furthermore, runoff simulations were affected more strongly by Ks variability than by θi variability. The degree of spatial variability influences on runoff simulations was related to the rainfall pattern and θi. Greater rainfall depth and instantaneous rainfall intensity corresponded to a smaller influence of the spatial variability. Stronger effects of the θi variability on runoff simulation were found in wetter soils, while stronger effects of the Ks variability were found in drier soils. For accurate runoff simulation, the θi variability can be completely ignored in cases of a 1-h duration storm with a return period greater than 10 years, while Ks variability should be fully considered even in the case of a 1-h duration storm with a return period of 20 years.
Editor D. Koutsoyiannis; Associate editor A. Fiori 相似文献
5.
The inductive response of a conducting horizontal cylinder embedded in a uniform earth is studied using numerical results obtained for an analytical solution for the problem of a conducting cylinder buried in a homogeneous earth for the case of a uniform inducing field. A check of the validity of the numerical results is made by a comparison with analogue model measurements for a number of cases. Numerical results for a range of cylinder radii (a = 1–10 km), depths of burial (d= 0–4 km), conductivity contrasts (σ2= 10?2-10 Sm?1), and source frequencies (f= 10?1-10?4 Hz) of interest in the interpretation of magnetotelluric field measurements are presented. The results indicate that for a uniform inducing field the conductivity and depth of burial of a horizontal cylindrical inhomogeneity are best determined through a measurement of the amplitudes Hy, Hz and Ex and the phases φy and Ψx. 相似文献
6.
Soil pipe flow tracer experiments: 2. Application of a streamflow transient storage zone model 
下载免费PDF全文
下载免费PDF全文 Soil pipes are important subsurface flow pathways in many soil erosion phenomena. However, limited research has been performed on quantifying and characterizing their flow and transport characteristics. The objectives of this research were to determine the applicability of a streamflow model with transient storage in deriving flow and transport characteristics of soil pipes. Tracer data from pulse inputs were collected in four different soil pipes after a fluorescein dye was injected in the upstream end of each soil pipe network in three branches (west, middle, and east) of a main catchment and a back catchment in Goodwin Creek Experimental Watershed in Mississippi. Multiple sampling stations were positioned along each soil pipe network. The transient storage zone model OTIS‐P was executed inversely to estimate transport parameters by soil pipe reach such as the soil pipe cross‐sectional area (A), soil storage zone cross‐sectional area (As), and exchange rate between the soil pipe and the soil storage zone (αs). Model convergence was achieved, and simulated breakthrough curves of the reaches were in good agreement with actual tracer data for eight of the nine reaches of the three branches of the Main Catchment and five of the seven reaches of the Back Catchment soil pipe. Simulation parameters for the soil pipe networks were similar to the range of values reported for flow and transport characteristics commonly observed in streams. Inversely, estimated soil pipe flow velocities were higher with increased tortuosity, which led to a smaller cross‐sectional areas predicted for the soil pipe flowpaths, while other parameters were not sensitive to tortuosity. In general, application of One‐Dimensional Transport with Inflow and Storage‐P to this unique soil pipe condition suggested larger transient storage (As and αs) compared with most stream systems. This was hypothesized to be because of relatively higher ratio of the wetted perimeter to flow area in the soil pipe, the hydraulic roughness of the soil pipe, potential retention in collapsed portions of the pipe, and interaction with smaller preferential flow systems. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
7.
根据场源不同可将地电场E分为大地电场ET和自然电场ESP。空间Sq电流系和固体潮一般被认为是大地电场ET的起源;自然电场ESP源于地下介质的物理化学作用,其局部性变化相对稳定。2013年南北地震带相继发生了芦山MS7.0地震（30.3°N,103.0°E）和岷县、漳县MS6.6地震（34.5°N,104.2°E）。两次地震周边约400 km范围内的地电场台站（成都、汉王）数据表明震前自然电场ESP出现明显的小幅度突跳或大幅度跃变,同时在时间上具有准同步性,其他地电场台站也有类似变化。基于大地电场岩体裂隙水（电荷）渗流（移动）模型,对两次地震周边自然电场ESP的变异机理进行探讨,尝试解析其物理过程。 相似文献
8.
Tsunami generated by submarine slumps and slides are investigated in the near-field, using simple source models, which consider the effects of source finiteness and directivity. Five simple two-dimensional kinematic models of submarine slumps and slides are described mathematically as combinations of spreading constant or slopping uplift functions. Tsunami waveforms for these models are computed using linearized shallow water theory for constant water depth and transform method of solution (Laplace in time and Fourier in space). Results for tsunami waveforms and tsunami peak amplitudes are presented for selected model parameters, for a time window of the order of the source duration.The results show that, at the time when the source process is completed, for slides that spread rapidly (cR/cT≥20, where cR is the velocity of predominant spreading), the displacement of the free water surface above the source resembles the displacement of the ocean floor. As the velocity of spreading approaches the long wavelength tsunami velocity
the tsunami waveform has progressively larger amplitude, and higher frequency content, in the direction of slide spreading. These large amplitudes are caused by wave focusing. For velocities of spreading smaller than the tsunami long wavelength velocity, the tsunami amplitudes in the direction of source propagation become small, but the high frequency (short) waves continue to be present. The large amplification for cR/cT1 is a near-field phenomenon, and at distances greater than several times the source dimension, the large amplitude and short wavelength pulse becomes dispersed.A comparison of peak tsunami amplitudes for five models plotted versus L/h (where L is characteristic length of the slide and h is the water depth) shows that for similar slide dimensions the peak tsunami amplitude is essentially model independent. 相似文献
9.
A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?2-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10?3 to 5.0× 10?5, for?α?=α 0cos?θ?sin?π?(r?ri ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?0. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos. 相似文献
10.
Conventional processes of extracting magnetotelluric signals from noisy records are reviewed: instrument noises and noises that are generated close to the detectors can be eliminated by the usual auto- and crosscorrelation processes. Identification of coherent noises, such as pulses due to field sources that are not uniform over at least 100 km in oil exploration or 1000 km in crustal studies, is much more tedious. The 5 components Hx, Hy, Hz, Ex, Ey, of the magnetotelluric field have been recorded in many areas in France at different periods of the year, (a) in non-uniform field sources in the vicinity of electric railways and of 50 cycle power lines, and (b) in areas of strong inhomogeneity at depth on the flanks of steep structures and near the sea shore. Means for detecting non-uniformity are reviewed. Measuring the vertical component of magnetic pulses is a good way of estimating field uniformity: if H vertical/H horizontal <10%, the uniform field assumption is valid, and the classical restitution formulas can be used; if H vertical/H horizontal > 10%, uniformity can not be assumed and there is some difficulty in deciding whether non-uniformity is due to the field source or to anisotropy or inhomogeneities at depth. Several ways to solve this difficulty are described. The reliability of calculation of actual resistivity at various depths is examined as a function of the precision of apparent resistivity measurements. 相似文献
11.
《国际泥沙研究》2019,34(6):564-576
Many parameters are involved in the deposition of a known mass of sand particles discharge from oblique pipes used in a range of engineering applications. Existing knowledge is not sufficient to accurately predict sediment mound dimensions and development from vertical and horizontal pipes. To better understand deposition patterns, laboratory experiments were done to investigate the development of subaqueous sand deposits from oblique pipes in stagnant water. Factors including nozzle diameter, release angle, release height, and sand mass were evaluated. It was found that nozzle size plays the most important role in shape formation and mound development. Release angle and release height were secondarily important. Five different shape patterns were observed and denoted as circular, ellipse, circular-ring, ellipse-ring, and pear-shaped. The scour hole diameter in the middle of the sediment mound was predicted using engineering assumptions and available semi-empirical correlations from the literature. It was found that the existing formulations can accurately predict the size of the scour hole for large release angles (θ = 60°). An aspect ratio of L/do was defined to combine the effect of nozzle size, do, and sand mass where L is the equivalent length of sand particles in a pipe. It was found that the mass balance equation can predict the deposition for L/do < 100. For higher aspect ratios, semi-empirical formulations were developed to predict the deposition length, width, height, and area with ±10% error. 相似文献
12.
Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θg), although no models have been proposed to reproduce the variation of WDPT with θg in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θg measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θg curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)?1. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
13.
Soil water content vertical profiles under natural conditions: matching of experiments and simulations by a conceptual model 总被引：1，自引：0，他引：1 下载免费PDF全文
下载免费PDF全文 The prediction of soil moisture content, θ, as a function of depth, z, and time, t, is of fundamental importance for applications in many hydrological processes. The main objective of this paper is to provide an approach to solve this problem at a local scale in soils with vegetation. The matching of soil moisture vertical profiles observed under natural conditions in grassy plots and their simulations by a conceptual model is presented. Experimental measurements were performed in a plot located in Central Italy, complete with hydrometeorological sensors specifically set up and equipped with a time domain reflectometry system providing the water content, θe(z, t). A conceptual model framework earlier proposed for two‐layered soil vertical profiles was modified and adopted for simulations. The changes concern the incorporation of evapotranspiration, the reduction of the original model for applications also to homogeneous soil vertical profiles, and a correction for the differences existing between assumed and observed initial moisture contents. In the model calibration, it was found that the effects of vegetation could be represented adequately by a fictitious soil vertical profile with a more permeable upper layer of saturated hydraulic conductivity, Ks, independent of time. Then, for the validation events, the model simulations in the stages of both infiltration and redistribution/evapotranspiration reproduced appropriately θe(z, t) with typical values of root mean square error in the range 0.0017–0.0657. Similar results were obtained by applying the modified two‐layered model for simulations of experimental data observed in three other plots located in Northern Italy and Germany. For all four vegetated sites, the two‐layer profile better matched the experimental data than the assumption of a homogeneous profile. Thus, the conceptual approach based on a two‐layered scheme for representing θ(z, t) in soils with vegetation appears to be appropriate for many hydrological applications. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
14.
M. Poddar 《Pure and Applied Geophysics》1981,119(5):1024-1036
FollowingDmitriev (1960) a rigorous theoretical solution for the problem of scattering by a perfectly conducting inclined half-plane buried in a uniform conductive half-space has been obtained for plane wave excitation. The resultant integral equation for the Laplace transform of scattering current in the half-plane is solved numerically by the method of successive approximation. The scattered fields at the surface of the half-space are found by integrating the half-space Green's function over the transform of the scattering current.The effects of depth of burial and inclination, of the half-plane on the scattered fields are studied in detail. An increase in the depth of burial leads to attenuation of the fields. Inclination introduces asymmetry in the field profiles beside affecting its magnitude. Depth of exploration is greater for quadrature component. An interpretation scheme based on a phasor diagram is presented for the VLF-EM method of exploration for rich vein deposits in a conductive terrain.List of symbols
x, y, z
Space co-ordinates
-
Half-space conductivity
- 0
Free-space permeability
-
Excitation frequency (angular)
-
T
Time
-
h
Depth of the half-plane
-
a
Inclination of the half-plane
-
E
x
x-Directed total electric field
-
E
x
p
x-Directed primary electric field
-
E
xo
p
x-Directed primary electric field atz=0 directly over the half-plane
-
H
y
y-Component of total magnetic field
-
H
y
p
y-Component of primary magnetic field
-
H
y0
p
y-Component of primary magnetic field atz=0 directly over the half-plane
-
H
z
z-Component of total magnetic field
-
H
z
p
z-Component of primary magnetic field
-
J
x
Surface density ofx-directed scattering current
-
G
Green's function
-
k
0,K
Wave numbers
-
u,u
0,u
1,u
2
Functions
-
Space co-ordinate
-
s
Variable in transform domain
-
Variable of integration
-
Normalized scattering current
-
Laplace transform of
-
N
Normalized
- , 0, 1, 2
Functions
-
t
Variable of integration
-
Skin depth
-
H
Total magnetic field
-
H
p
Primary magnetic field
-
H
0
p
Primary magnetic field atz=0 directly over the half-plane
-
M,Q,R,S,U,V
Functions
-
N
1,N
2
Functions 相似文献
15.
Field measurements were conducted to study the influence of aquatic vegetation on flow structures in floodplains under combined currents and wind-driven waves. Wave and turbulent velocities were decomposed from the time series of instantaneous velocity and analysed separately. In the present study, the wind waves were small, leading to the ratios of wave excursion (Ew) to stem spacing (S) for all cases tested here were less than 0.5. This caused the vertical distributions of time-averaged velocity (Uhoriz) and turbulent kinetic energy (TKE) impacted by vegetation similar with the vegetated flow structures under pure current conditions. For emergent vegetation, Uhoriz and TKE distributed uniformly through the entire water column or increased slightly from bed to water surface. Similar distributions were present in the lower part of submerged vegetation. In the upper part of submerged vegetation, Uhoriz and TKE increased rapidly toward water surface and TKE reached its maximum near the top of vegetation. The measured wave orbital velocity (Uw) fitted linear wave theory well through the entire water depth for both the emergent and submerged cases, so that with small Ew/S the wave velocity was not attenuated within vegetation and Uw within canopy can be predicted by the linear wave theory under combined currents and waves. However, wind-driven waves made the turbulence generated near the top of canopy penetrate a deeper depth into vegetation than predictions under pure current conditions. 相似文献
16.
The present study proposes a method for evaluating the effectiveness of road‐crossing drainage culverts in ephemeral streams. This approach is focused on estimating the culvert capacity in road–stream crossings and the probable runoff generated on the road from hydrological thresholds associated with hydromorphological criteria. In particular, discharges at bankfull and flood‐prone stages have been used in combination with 2.5 and 100‐year peak discharges. Different hydraulic variables have been considered for calculating the discharge through culverts under these conditions (e.g. tailwater and headwater depth, inlet control, pipe roughness, pipe cross‐area and slope, pipe outlet velocity, critical water depth, and flow rate over the road). Geomorphological factors such as bed stability, bed load transport, and channel roughness have also been considered because of their potential for obstructing the drains in this type of channel. In addition, a potential obstruction index (PIOBSTR) has been calculated, as a dependent parameter of the obstacle index (IOBST) and the potential build‐up of coarse sediments (PBCS). The study has been carried out on the Mediterranean coast in the region of Murcia (Spain), where there are numerous examples of road–stream crossings equipped with culverts in ephemeral channels that could cause highly dangerous situations for road traffic. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
17.
Factors affecting the measurement of the vertical hydraulic conductivity of a streambed sediment using standpipe tests 下载免费PDF全文
下载免费PDF全文 The vertical hydraulic conductivity (Kv) of a stream or lake sediment is often determined in the field using standpipe tests. Calculation of Kv is based on the assumption that the hydraulic head in the pipe is equal to that of the stream or lake stage. In this work, a modified equation for Kv is developed for the standpipe test which is applicable when this assumption is not valid. The equation involves not only the hydraulic head at different times but also the difference in the hydraulic head (a) between the groundwater level and river stage. The effects of certain factors on Kv, such as the ratio of the hydraulic head at different times (h1/h2), the difference a, and the initial water table height (h0), are also discussed. The results show that when h1/h2 is constant, the relative error (Er) in Kv increases with the ratio a/h2. Furthermore, if a/h2 < 0.05, then for any value of h1/h2, Er is less than 5% using the modified equation. Also, if a/h2 is large, hydraulic head readings with larger h1/h2 ratios must be used to avoid large Er values. The results of a field test also indicate that the error in Kv decreases as the value of h0 increases. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
18.
Max Sieber Sebastian Klar Michalis F. Vassiliou Ioannis Anastasopoulos 《地震工程与结构动力学》2020,49(14):1388-1405
Recognizing the beneficial effect of nonlinear soil–foundation response has led to a novel design concept, termed ‘rocking isolation’. The analysis and design of such rocking structures require nonlinear dynamic time history analyses. Analyzing the entire soil–foundation–structure system is computationally demanding, impeding the application of rocking isolation in practice. Therefore, there is an urgent need to develop efficient simplified analysis methods. This paper assesses the robustness of two simplified analysis methods, using (i) a nonlinear and (ii) a bilinear rocking stiffness combined with linear viscous damping. The robustness of the simplified methods is assessed by (i) one-to-one comparison with a benchmark finite element (FE) analysis using a selection of ground motions and (ii) statistical comparison of probability distributions of response quantities, which characterize the time history response of rocking systems. A bridge pier (assumed rigid) supported on a square foundation, lying on a stiff clay stratum, is used as an illustrative example. Nonlinear dynamic FE time history analysis serves as a benchmark. Both methods yield reasonably accurate predictions of the maximum rotation θmax. Their stochastic comparison with respect to the empirical cumulative distribution function of θmax reveals that the nonlinear and the bilinear methods are not biased. Thus, both can be used to estimate probabilities of exceeding a certain threshold value of θ. Developed in this paper, the bilinear method is much easier to calibrate than the nonlinear, offering similar performance. 相似文献
19.
Characteristics of coseismic water level changes at Tangshan well for the Wenchuan MS8.0 earthquake and its larger aftershocks 下载免费PDF全文
下载免费PDF全文 Baojun Yin Li Ma Huizhong Chen Jianping Huang Chaojun Zhang Wuxing Wang 《地震科学（英文版）》2009,22(2):149-157
Coseismic water level changes which may have been induced by the Wenchuan MS8.0 earthquake and its 15 larger aftershocks (MS≥?5.4) have been observed at Tangshan well. We analyze the correlation between coseismic parameters (maximum amplitude, duration, coseismic step and the time when the coseismic reach its maximum amplitude) and earthquake parameters (magnitude, well-epicenter distance and depth), and then compare the time when the coseismic oscillation reaches its maximum amplitude with the seismogram from Douhe seismic station which is about 16.3 km away from Tangshan well. The analysis indicates that magnitude is the main factor influencing the induced coseismic water level changes, and that the well-epicenter distance and depth have less influence. MS magnitude has the strongest correlation with the coseismic water level changes comparing to MW and ML magnitudes. There exists strong correlation between the maximum amplitude, step size and the oscillation duration. The water level oscillation and step are both caused by dynamic strain sourcing from seismic waves. Most of the times when the oscillations reach their maximum amplitudes are between S and Rayleigh waves. The coseismic water level changes are due to the co-effect of seismic waves and hydro-geological environments. 相似文献
20.
ABSTRACTSourcing subsurface evaporation (Ess) into groundwater (Eg) and unsaturated zone (Eu) components has received little scientific attention so far, despite its importance in water management and agriculture. We propose a novel sourcing framework, with its implementation in dedicated post-processing software called SOURCE (used along with the HYDRUS1D model), to study evaporation sourcing dynamics, define quantitatively “shallow” and “deep” water table conditions and test the applicability of water table fluctuation (WTF) and “bucket” methods for estimation of Eg and Eu separately.For the “shallow” and “deep” water table we propose Eg?>?0.95Ess and Eg = 0 criteria, respectively. Assessment of the WTF method allowed sourcing of very small fluxes otherwise neglected by standard hydrological methods. Sourcing with SOURCE software was more accurate than the standard “bucket” method mainly because of greater flexibility in spatio-temporal discretization. This study emphasized the dry condition relevance of groundwater evaporation which should be analysed by applying coupled flow of heat, vapour and liquid water.
Editor D. Koutsoyiannis; Associate editor S. Kanae 相似文献