首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 281 毫秒
1.
Abstract

Abstract Recently, substantial progress has been made in detection and observation of non-aqueous phase liquids (NAPLs) in the subsurface using different experimental techniques. However, there is still a lack of appropriate direct methods to measure the saturation of NAPL (θNAPL). This paper provides a guide for estimating θNAPL and water content (θ w ) in unsaturated and saturated sand based on direct measurements of soil dielectric constant (Ka ) and electrical conductivity (σ a ) using time domain reflectometry (TDR). The results show that the previously used dielectric mixing model fails to predict θNAPL in the case of a four-phase system. A new methodology is suggested and exemplified by showing that the measured Ka gives accurate estimation of θNAPL for a three-phase system while in a four-phase system, both θ a and Ka need to be measured. The results show that using the suggested methodology, accurate predictions of θ w (R 2 = 0.9998) and θNAPL lower than 0.20 m3 m-3 (average R 2 = 0.9756) are possible.  相似文献   

2.
Selecting three half orbits near the epicenter of Pu’er earthquake, we analyzed the Ne data recorded in their revis-ited orbits during a year before this earthquake, and extracted Ne precursors. The results show that: ① There are significant seasonal variations of ionospheric Ne in night time, which exhibit different shapes respectively in four seasons; ② There are three main shapes of Ne: single-peak, saddle-shaped and even-shaped, all of which may oc-cur in four seasons, but each season with its typical shape relatively; ③ Spatial images of Ne showed high values near the epicenter in 30 days before the earthquake, and there is a good correlation between anomaly and distribu-tion of earthquake in space and time, which reflects that these spatial anomalies were indeed concerned with the earthquake; ④There shows a certain similarity of the Ne curves among revisited orbits, which can provide back-ground information for distinguishing and identification of seismic anomaly.  相似文献   

3.
The objective of this paper is to simulate the progress of the soil water content distribution in the soil profile with a water table at the bottom of the soil profile during ponding irrigation. This simulation can be done by solving the two‐dimensional Richards's equation for the assimilation of the advancing water jet, which uses the conditions of the two exponential functional forms k = ks eαψ and θ = θr + (θs − θr) eαψ to represent the hydraulic conductivity and volumetric water content, with ψ the pressure as the third variable. We assume that the ground surface becomes ponded and saturated as soon as the water flux passes the dry ground surface. By the technique of transformation, the analytical solution of these two‐dimensional Richards' equations has enabled figures of volumetric water content distribution to be obtained in successive time periods after irrigation. For the example of loam soil, it can simulate the variation of volumetric water content during and after irrigation in the soil profile. The analytical solutions of this paper reflect the real situation simulated, and can be applied to verify those complicated solutions from other analytical models. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

4.
Solute transport experiments using a non-reactive tracer were conducted on short, undisturbed, saturated columns of a sandy loam soil. All columns, 20 cm in diameter and 20 cm long, were collected along a transect of 35 m. Most of the soil columns had pre-existing macropores. The columns were leached at a steady flow-rate under ponding conditions. The resulting breakthrough curves (BTCs) showed a large heterogeneity. Several of the BTCs displayed early breakthrough and long tailing. All the data were interpreted in terms of dimensional time moments, the classical convection-dispersion equation (CDE) and the mobile-immobile transport model (MIM). Experimental time moments were found to vary significantly among the different BTCs. Analysis of the time moments also revealed that the variance of the field-scale BTC was several times larger than the average of the local-scale variance. The pore water velocity v and dispersion coefficient D were obtained by fitting the CDE to the local-scale BTCs, resulting in an average dispersivity of 7·4 cm. Frequency distributions for the CDE parameters v and D were equally well described by a normal or log-normal probability density function (pdf). When a log-normal pdf for D is considered, the variance of the loge transformed D values (σln D2) was found to be 2·1. For the MIM model, two additional parameters were fitted: the fraction of mobile water, θm/θ, and the first-order mass transfer coefficient, α. The MIM was more successful in describing the data than the CDE transport model. For the MIM model, the average dispersivity was about 2 cm. The MIM parameters v, D and θm/θ were best described by a log-normal pdf rather than a normal pdf. Only the parameter α was better described by a normal pdf. Mobile water fractions, θm/θ ranged from 0·01 to 0·98, with a mean of 0·43 (based on a log-normal pdf). When the CDE and MIM were applied to the data, the fitted pore water velocities, v, compared favourably with the effective pore water velocities, veff, obtained from moment analysis.  相似文献   

5.
Intergrain contact density indices for granular mixes I: Framework   总被引:2,自引:0,他引:2  
Mechanical behavior such as stress-strain response, shear strength, resistance to liquefaction, modulus, and shear wave velocity of granular mixes containing coarse and fine grains is dependent on intergrain contact density of the soil. The global void ratio e is a poor index of contact density for such soils. The contact density depends on void ratio, fine grain content (Cv), size disparity between particles, and gradation among other factors. A simple analysis of a two-sized particle system with large size disparity is used to develop an understanding of the effects of Cv, e, and gradation of coarse and fine grained soils in the soil mix on intergrain contact density. An equivalent intergranular void ratio (ec)oq is introduced as a useful intergrain contact density for soils at fines content of less than a threshold value Crth. Beyond this value, an equivalent interfine void ratio (ef)eq is introduced as a primary intergrain contact density index. At higher values of Cv beyond a limiting value of fine grains content CVL, an interfine void ratio ef is introduced as the primary contact density index. Relevant equivalent relative density indices (Drc)eq and (Drf)eq are also presented. Experimental data show that these new indices correlate well with steady state strength, liquefaction resistance, and shear wave velocities of sands, silty sands, sandy silts, and gravelly sand mixes.  相似文献   

6.
Earthquakes result from continuous geodynamic processes. A topic of significant interest for the scientific community is to elaborate on the phenomena governing the faulting and fracturing of crustal rocks. Therefore, in this study, uniaxial compressive shear failure experiments were conducted on Fangshan marble rock samples with a prefabricated slot to simulate thrust faulting. The center of each marble plate (105 mm × 80 mm × 5 mm) was engraved with a 30-mm long double-sided nonpenetrating slot (depth: 2 mm, width: 0.5 mm). The deformation and destruction processes of the rock surface were recorded using a high-speed camera. The digital image correlation method was used to calculate the displacement and strain distribution and variation at different loading stages. The accumulative and incremental displacement fields u and v , strain field ex and ey, and shear strain exy were analyzed. When the loading level reached its ultimate value, the strain field was concentrated around the prefabricated slot. The concentration reached a maximum at the ends of the prefabricated slot. The magnitude of shear strain reached 0.1. This experiment contributes to our understanding of the dynamic process of active faulting.  相似文献   

7.
The water retention curve (θ(ψ)), which defines the relationship between soil volumetric water content (θ) and matric potential (ψ), is of paramount importance in characterizing the hydraulic behaviour of soils. However, few methods are so far available for estimating θ(ψ) in undisturbed soil samples. We present a new design of TDR‐pressure cell (TDR‐Cell) for estimating θ(ψ) in undisturbed soil samples. The TDR‐Cell consists of a 50‐mm‐long and 50‐mm internal diameter stainless steel cylinder (which constitutes the outer frame of a coaxial line) attached to a porous ceramic disc and closed at the ends with two aluminium lids. A 49‐mm‐long and 3‐mm‐diameter stainless steel rod, which runs longitudinally through the centre of the cylinder, constitutes the inner rod of a coaxial TDR probe. The TDR‐Cell was used to determine the θ(ψ) curves of a packed sand and seven undisturbed soil samples from three profiles of agricultural soils. These θ(ψ) curves were subsequently compared to those obtained from the corresponding 2‐mm sieved soils using the pressure plate method. Measurements of bulk electrical conductivity, σa, as a function of the water content, σa(θ), of the undisturbed soil samples were also performed. An excellent correlation (R2 = 0·988) was found between the θ values measured by TDR on the different undisturbed soils and the corresponding θ obtained from the soil gravimetric water content. A typical bimodal θ(ψ) function was found for most of the undisturbed soil samples. Comparison between the θ(ψ) curves measured with the TDR‐Cell and those obtained from the 2‐mm sieved soils showed that the pressure plate method overestimates θ at low ψ values. The σa(θ) relationship was well described by a simple power expression (R2 > 0·95), in which the power factor, defined as tortuosity, ranged between 1·18 and 3·75. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
We develop a parameterized theory of convection driven by radiogenic and tidal heating. The tidal heating depends on eccentricity e of a satellite’s orbit. Using parameterized theory we determine the intensity of convection as a function of e and satellite’s properties. The theory is used for 6 medium sized satellites of Saturn. We find that endogenic activity on Tethys and Dione is possible if e exceeds some critical values e cr . For Enceladus, e was probably close to the present value for billions of years. We cannot find constrains for e of Mimas and Iapetus. The theory successfully predicts the possibility of present endogenic activity in Dione and rules out such activity in Tethys. Both these facts were recently confirmed by Cassini mission.  相似文献   

9.
Results are presented of recent sensitivity calibration of 76 accelerographs (SMA-1) of the Los Angeles Strong Motion Array. These have pendulum-like transducers and optical recording system. One characteristic of their design is off-axis sensitivity, which is magnified by transducer misalignment. A new calibration procedure was applied, which considers off-axis sensitivity and measures the angles of misalignment (φ and ψ), as well as the incident angle of the light beam onto the film (θ0). These are required (1) for accurate estimation of sensitivity, and (2) for proper instrument correction of recorded accelerograms which considers also cross-axis sensitivity and misalignment. These effects are important near large acceleration peaks (approaching and exceeding 1g), e.g. like the ones recorded near the source of the 1994 Northridge earthquake (ML=6·4). This earthquake was recorded by 65 stations of the Los Angeles Strong Motion Array, at epicentral distances from 2 to 85 km. Histograms showing distribution of the misalignment angles, light beam incidence angle θ0 (for unloaded position) and the transducer sensitivities are presented. These indicate that the misalignment angles are typically 1–1·5°, but may also be 3–4°. Angle θ0 (usually neglected), is mostly between ±8°, but may reach ±12°. Assuming θ0=0 leads to systematically smaller values of the measured sensitivity (e.g. by ∼3% for θ0=8° and ∼4% for θ0=12°). Comparison of the newly measured sensitivities with those measured prior to installation (in 1979/1980), sold, shows that, in general, the new values are systematically smaller. The difference is typically within 5 per cent, but in some cases is as large as 10 per cent. Other principal sources of the observed differences and their mechanisms are discussed. Those include long-term changes in the transducers (e.g. change of stiffness, reflected in changes of the natural frequency) and differences in the calibration procedure (e.g. errors associated with manual reading film records with tilt test data, and with transducer and instrument housing misalignment). The presented results may be considered typical of similar strong motion arrays worldwide. © 1998 John Wiley & Sons, Ltd.  相似文献   

10.
In this paper a spatially distributed model of the hillslope sediment delivery processes, named the sediment delivery distributed (SEDD) model, is initially reviewed; the model takes into account the sediment delivery processes due to both the hillslope sediment transport and the effects of slope curvature. Then the rainfall and sediment yield events measured at the experimental SPA2 basin, in Sicily, are used both to calibrate the SEDD model and to verify the predictive capability of the distributed sediment delivery approach at event scale. For the SPA2 basin discretized into morphological units and stream tubes, the SEDD model is calibrated at event scale using the measurements carried out at the outlet of the experimental basin in the period December 2000–January 2001. The model calibration is used to determine a relationship useful for estimating the unique coefficient βe of the model by rainfall erosivity factor Re at event scale. To test the predictive capability of the βe = f(Re) relationship, 20 events measured in the period September 2002–December 2005 are used; the comparison between measured sediment yield values and calculated ones for all monitored events shows that the sediment delivery distributed approach has a good predictive ability at event scale. The analysis also shows that estimating βe by the relationship βe = f(Re) gives a better agreement between measured and calculated sediment yields than obtained with the median value βe,m of all 27 calculated βe values. Finally the analysis at annual scale, for the period 2000–2005, allows the estimation of the median value βa,m representative of the annual behaviour. This analysis shows that the sediment delivery distributed approach also has a good predictive ability at annual scale. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

11.
The eddy covariance technique and the cuvette method were used to investigate water use efficiency in an irrigated winter wheat (Triticum asetivum L.)/summer maize (Zea mays L.) rotation system in the North China Plain. The results show that ecosystem water use efficiency (WUEe) changed diurnally and seasonally. Daily maximal WUEe appeared in the morning. WUEe generally peaked in late April in wheat field and in late July/early August in maize field. From 2003 to 2006, seasonal mean WUEe was 6.7–7.4 mg CO2 g−1 H2O for wheat and 8.4–12.1 mg CO2 g−1 H2O for maize. WUEe was much lower than canopy water use efficiency (WUEc) under small leaf area index (LAI) but very close to WUEc under large LAI. With the increase in LAI, WUEe enlarged rapidly under low LAI but slowly when LAI was higher than one. WUEe was greater on the cloudy days than on the sunny days. Under the same solar radiation, WUEe was higher in the morning than in the afternoon. The ratio of internal to ambient CO2 partial pressure (Ci/Ca) decreased significantly with the increase in photosynthetically active radiation (PAR) when PAR was lower than the critical values (around 500 and 1000 μmol m−2 s−1 for wheat and maize, respectively). Beyond critical PAR, Ci/Ca was approximately constant at 0.69 for wheat and 0.42 for maize. Therefore, when LAI and solar radiation was large enough, WUEe has negative correlation with vapor pressure deficit in both of irrigated wheat and maize fields.  相似文献   

12.
Abstract

Magnetic instabilities play an important role in the understanding of the dynamics of the Earth's fluid core. In this paper we continue our study of the linear stability of an electrically conducting fluid in a rapidly rotating, rigid, electrically insulating spherical geometry in the presence of a toroidal basic state, comprising magnetic field BMB O(r, θ)1ø and flow UMU O(r, θ)1ø The magnetostrophic approximation is employed to numerically analyse the two classes of instability which are likely to be relevant to the Earth. These are the field gradient (or ideal) instability, which requires strong field gradients and strong fields, and the resistive instability, dependent on finite resistivity and the presence of a zero in the basic state B O(r,θ). Based on a local analysis and numerical results in a cylindrical geometry we have established the existence of the field gradient instability in a spherical geometry for very simple basic states in the first paper of this series. Here, we extend the calculations to more realistic basic states in order to obtain a comprehensive understanding of the field gradient mode. Having achieved this we turn our attention to the resistive instability. Its presence in a spherical model is confirmed by the numerical calculations for a variety of basic states. The purpose of these investigations is not just to find out which basic states can become unstable but also to provide a quantitative measure as to how strong the field must become before instability occurs. The strength of the magnetic field is measured by the Elsasser number; its critical value c describing the state of marginal stability. For the basic states which we have studied we find c 200–1000 for the field gradient mode, whereas for the resistive modes c 50–160. For the field gradient instability, c increases rapidly with the azimuthal wavenumber m whereas in the resistive case there is no such pronounced difference for modes corresponding to different values of m. The above values of c indicate that both types of instability, ideal and resistive, are of relevance to the parameter regime found inside the Earth. For the resistive mode, as is increased from c, we find a shortening lengthscale in the direction along the contour BO = 0. Such an effect was not observable in simpler (for example, cylindrical) models.  相似文献   

13.
Whether the presence of non-plastic silt in a granular mix soil impact its liquefaction potential and how to evaluate liquefaction resistance of sand containing different amounts of silt contents are both controversial issues. This paper presents the results of an experimental evaluation to address these issues. Two parameters, namely, equivalent intergranular void ratio (ec)eq and equivalent interfine void ratio (ef)eq, proposed in a companion paper (Thevanayagam, 2007) as indices of active grain contacts in a granular mix, are used to characterize liquefaction resistance of sands and silty sands. Results indicate that, at the same global void ratio (e), liquefaction resistance of silty sand decreases with an increase in fines content (Cv) up to a threshold value (Crth). This is due to a reduction in intergrain contact density between the coarse grains. Beyond Crth, with further addition of fines, the interfine contacts become significant while the inter-coarse grain contacts diminish and coarse grains become dispersed. At the same e, the liquefaction resistance increases and the soil becomes stronger with a fttrther increase in silt content. Beyond a limiting fines content (CrL), the liquefaction resistance is controlled by interfine contacts only. When Cr〈Crth, at the same (e)eq, the liquefaction resistance of silty sand is comparable to that of the host clean sand at a void ratio equal to (ec)eq. When CF〉CFth, at the same (ef)eq, the cyclic strength of a sandy silt is comparable to the host silt at a void ratio equal to (ef)eq.  相似文献   

14.
Abstract

We investigate the influence of differential rotation on magnetic instabilities for an electrically conducting fluid in the presence of a toroidal basic state of magnetic field B 0 = BMB0(r, θ)1 φ and flow U0 = UMU0 (r, θ)1φ, [(r, θ, φ) are spherical polar coordinates]. The fluid is confined in a rapidly rotating, electrically insulating, rigid spherical container. In the first instance the influence of differential rotation on established magnetic instabilities is studied. These can belong to either the ideal or the resistive class, both of which have been the subject of extensive research in parts I and II of this series. It was found there, that in the absence of differential rotation, ideal modes (driven by gradients of B 0) become unstable for Ac ? 200 whereas resistive instabilities (generated by magnetic reconnection processes near critical levels, i.e. zeros of B0) require Ac ? 50. Here, Λ is the Elsasser number, a measure of the magnetic field strength and Λc is its critical value at marginal stability. Both types of instability can be stabilised by adding differential rotation into the system. For the resistive modes the exact form of the differential rotation is not important whereas for the ideal modes only a rotation rate which increases outward from the rotation axis has a stabilising effect. We found that in all cases which we investigated Λc increased rapidly and the modes disappeared when Rm ≈ O(ΛC), where the magnetic Reynolds number Rm is a measure of the strength of differential rotation. The main emphasis, however, is on instabilities which are driven by unstable gradients of the differential rotation itself, i.e. an otherwise stable fluid system is destabilised by a suitable differential rotation once the magnetic Reynolds number exceeds a certain critical value (Rm )c. Earlier work in the cylindrical geometry has shown that the differential rotation can generate an instability if Rm ) ?O(Λ). Those results, obtained for a fixed value of Λ = 100 are extended in two ways: to a spherical geometry and to an analysis of the effect of the magnetic field strength Λ on these modes of instability. Calculations confirm that modes driven by unstable gradients of the differential rotation can exist in a sphere and they are in good agreement with the local analysis and the predictions inferred from the cylindrical geometry. For Λ = O(100), the critical value of the magnetic Reynolds number (Rm )c Λ 100, depending on the choice of flow U0 . Modes corresponding to azimuthal wavenumber m = 1 are the most unstable ones. Although the magnetic field B 0 is itself a stable one, the field strength plays an important role for this instability. For all modes investigated, both for cylindrical and spherical geometries, (Rm )c reaches a minimum value for 50 ≈ Λ ≈ 100. If Λ is increased, (Rm )c ∝ Λ, whereas a decrease of Λ leads to a rapid increase of (Rm )c, i.e. a stabilisation of the system. No instability was found for Λ ≈ 10 — 30. Optimum conditions for instability driven by unstable gradients of the differential rotation are therefore achieved for ≈ Λ 50 — 100, Rm ? 100. These values lead to the conclusion that the instabilities can play an important role in the dynamics of the Earth's core.  相似文献   

15.
Relationships between gravimetric soil moisture content (w) and matric potential (ϕ), and between volumetric soil moisture content (θv) and pressure head (h) were approximated for the unsaturated zone on Long Island, New York. Soil samples were collected from two sites using a hand auger. The soil moisture content was determined using the filter‐paper (wf) and gravimetric (w) methods, respectively. The wf was then used in an empirical equation to estimate ϕm. Each set of ϕm and w was combined with a straight‐line empirical model to obtain a wm) relationship. Soil ϕm was converted to h, and w to the volumetric moisture content θv, in order to produce a θv(h) curve. Graphical and statistical comparison showed that the resulting θv(h) curves fell within one order of magnitude of similar curves generated by a more sophisticated non‐linear model developed previously. The simplicity and low cost of the filter‐paper approach described in this study recommends it for preliminary studies of hydraulic properties in the unsaturated zone. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

16.
The elastic thickness of the lithosphere in the Pacific Ocean   总被引:1,自引:0,他引:1  
In this study, we present determinations of the effective elastic thicknessTe of the oceanic lithosphere along Pacific chains or archipelagoes.Te is determined by computing the deflection of a continuous elastic plate under the load of volcanoes, and constrained by geoid heights provided by SEASAT. In the South Central Pacific, estimates of 14 km for the Marquesas and 6 km or less for the Pitcairn-Mururoa-Gloucester chain are in good agreement with a previous work in this region (Cook-Austral and Society chains). Around the Line Islands chain, SEASAT data reveal that the bathymetry is poorly known, preventing fine analysis. Meanwhile,Te looks globally very low ( 6 km), except for three volcanoes but these results may be unreliable. The Easter chain features lowTe values ( 6 km), with no noticeable variation along the chain. Higher values are found for a Samoan island, Manuae (24 km), and along the Hawaiian-Emperor seamounts chain (from 32 km at the eastern end of the chain to 21.5 km for the Hawaiian volcanoes, and from 25.5 to 15 km for the Emperor seamounts). The large number ofTe estimates obtained in this study points out a noticeable difference between North and South Pacific results. Those from the North Pacific agree with the general trend (increase with the square root of age plate at loading time), while those from the South Central Pacific are much lower, according to their plate age. These lowTe results from the South Pacific are only partly explained by taking account of thermal perturbations using the rejuvenation model. Therefore, these results then point out a regional difference in oceanic lithosphere.  相似文献   

17.
The responses, re, given by several multicomponent combination rules used in seismic codes for determining peak responses to three ground motion components are evaluated for elastic systems and compared with the critical response rcr; this is defined as the largest response for all possible incident angles of the seismic components and obtained by means of the CQC3‐rule when a principal seismic component is vertical, or the GCQC3‐rule when it departs from the vertical direction. The combination rules examined are the SRSS‐, 30%‐, 40%‐ and IBC‐rules, considering different alternatives for the design horizontal spectrum. Assuming that a principal seismic component is along the vertical direction, the upper and lower bounds of the ratio re/rcr for each combination rule are determined as a function of the spectral intensity ratio of the horizontal seismic components and of the responses to one seismic component acting alternately along each structural axis. Underestimations and overestimations of the critical response are identified for each combination rule and each design spectrum. When a component departs from the vertical direction, the envelopes of the bounds of the ratio re/rcr for each combination rule are calculated, considering all possible values of the spectral intensity ratios. It is shown that the inclination of a principal component with respect to the vertical axis can significantly reduce the values of re/rcr with respect to the case when the component is vertical. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

18.
The effective relative dielectric constant ?e, r and the effective conductivity σe have each been determined as a function of frequency in the range 1–3000 MHz at volumetric water contents of up to approximately 0.74 for clays, 0.83 for a peat and 0.56 for a silt. At frequencies above about 25 MHz (depending on soil type), ?e, rincreases with water content for all samples. However, at lower frequencies, ?e, ronly increases with water content as long as the wet density also increases, which is the case for water contents up to a critical value lying between 0.35 and 0.48. At higher water contents, ?e, rand the wet density decrease with increasing water content. Consequently, curves of ?e, rversus frequency for two wet samples with different water contents, at least one of them higher than the critical value, are seen to cross at about 25 MHz. Below the critical value the curve of the sample with the lower water content is below the other curve at all freqencies applied. At a given frequency, σe has a maximum as a function of water content. This is tentatively explained by assuming that σe is the sum of pore water conductivity (increasing with water content until all salts in the soil are dissolved into the water and then decreasing) and surface water conductivity (increasing with wet density and therefore increasing with water content up to the critical value and then decreasing). At frequencies higher than 1000 MHz, ?e, rdepends only weakly on salinity (which is represented by the measured conductivity). It shows an increasing dependence if the frequency is decreased towards 1 MHz. The highest values of ?e, rand σe, measured in this work, occur for a sample of wet, nearly saturated silt originating from a location below sea-level near to the Dead Sea, Israel: ?e, rdecreases continuously from a value of about 104 at 3 MHz to about 102 at 200 MHz, while σe rises from about 4 S/m to 5 S/m at these respective frequencies. The dependence of the wavelength on the loss-tangent is strong and the wavelength is considerably smaller than it would be in a dielectric. This is the only sample for which σe increases with water content, even if the latter is above its critical value. Therefore it is assumed that the pore water conductivity is greater than the surface water conductivity if the volumetric water content is lower than 0.564, the maximum value applied. The measurements give evidence for the presence of a relaxation at about 3 MHz for all samples examined.  相似文献   

19.
Abstract

Time-domain reflectometry (TDR) is an electromagnetic technique for measurements of water and solute transport in soils. The relationship between the TDR-measured dielectric constant (Ka ) and bulk soil electrical conductivity ([sgrave]a) to water content (θW) and solute concentration is difficult to describe physically due to the complex dielectric response of wet soil. This has led to the development of mostly empirical calibration models. In the present study, artificial neural networks (ANNs) are utilized for calculations of θw and soil solution electrical conductivity ([sgrave]w) from TDR-measured Ka and [sgrave]a in sand. The ANN model performance is compared to other existing models. The results show that the ANN performs consistently better than all other models, suggesting the suitability of ANNs for accurate TDR calibrations.  相似文献   

20.
Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θF. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θf; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α2Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θf—θF) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α2Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R α, explicitly ε½ ? R 2 α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θf—θF small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D min| ≤ |D c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R 2 α ? ε½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号