首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and physical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length z0m and excess resistance κB-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB-1 at the two stations. Results indicate that z0m has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of z0m in the land surface model and that obtained from this study. κB-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of z0m and κB-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.  相似文献   

2.
Based on the spherical earth dislocation theory and a fault slip model of the Tohoku-Oki MW9.0 earthquake, the co-seismic Coulomb failure stress changes (ΔCFS) on the northern Tanlu fault zone at depths of 0–40 km are calculated. By comparing two sets of results from the spherical earth dislocation theory and the semi-infinite space one, the effect of earth curvature on the calculation results is analyzed quantitatively. First, we systematically summarize previous researches related to the northern Tanlu fault zone, divide the fault zone as detailed as possible, give the geometric parameters of each segment, and establish a segmented structural model of the northern Tanlu fault zone. Second, we calculate the Coulomb stress changes on the northern Tanlu fault zone by using the spherical earth dislocation theory. The result shows the Coulomb stress changes are no more than 0.003 MPa, which proves the great earthquake did not significantly change the stress state of the fault zone. Finally, we quantitatively analyze the disparities between the results of semi-infinite space dislocation theory and the spherical earth one. The average disparity between them is about 7.7% on the northern Tanlu fault zone and is 16.8% on the Fangzheng graben, the maximum disparity on this graben reaches up to 25.5%. It indicates that the effect of earth curvature can not be ignored. So it’s necessary to use the spherical earth dislocation theory instead of the semi-infinite space one to study the Coulomb stress change in the far field.  相似文献   

3.
The stability of cohesive sediments from Venice lagoon has been measured in situ using the benthic flume Sea Carousel. Twenty four stations were occupied during summertime, and a sub-set of 13 stations was re-occupied during the following winter. Erosion thresholds and first-order erosion rates were estimated and showed a distinct difference between inter-tidal and sub-tidal stations. The higher values for inter-tidal stations are the result of exposure that influences consolidation, density, and organic adhesion. The thresholds for each state of sediment motion are well established. However, the rate of erosion once the erosion threshold has been exceeded has been poorly treated. This is because normally a time-series of sediment concentration (C) and bed shear stress (τ0(t)) is used to define threshold stress or cohesion (τcrit,z) and erosion rate (E). Whilst solution of the onset of erosion, τcrit,0, is often reported, the evaluation of the erosion threshold variation through the process of erosion (eroded depth) is usually omitted or not estimated. This usually leads to assumptions on the strength profile of the bed which invariably has no credibility within the topmost mm of the bed where most erosion takes place. It is possible to extract this information from a time-series through the addition of a step in data processing. This paper describes how this is done, and the impact of this on the accuracy of estimates of the excess stress (τ0(t)–τcrit,z) on E.  相似文献   

4.
The hydrodynamics of a system where there is a coupled flow above and below a sediment–water interface (SWI) are not completely understood. We numerically simulate mean two-dimensional, unidirectional, steady, viscous flow in these systems using a sequentially coupled formulation. Simulations were conducted to determine fundamental relationships between bedform geometry, Reynolds number for the water-column flow (Re), interfacial exchange zone depth (dz) in the sediments, and flux through the SWI (qint); the latter two parameters play a significant role in biogeochemical and aquatic-life processes across the SWI. dz and Re are functionally related through an asymptotic growth-curve model while qint and Re follow a power function. These relationships are dynamically explained by the manner in which pressure gradients along the SWI develop due to current–bedform interactions at different Res and by Darcy’s Law. We found that the coupling between water column and exchange zone flow is controlled by the behavior of the water-column eddy. The eddy detaches at or near the point of minimum pressure along the interface, and reattaches near the point of maximum pressure. These two critical points determine the pressure gradient along the bed surface that controls the exchange zone flow field. Moreover, the reattachment point corresponds to flow divides within the sediments. Lastly, pore-water velocities drop with depth below the SWI, and are larger below the bedform crests than below the troughs.  相似文献   

5.
A very fast and efficient approach for gravity data inversion based on the regularized conjugate gradient method has been developed. This approach simultaneously inverts for the depth (z), and the amplitude coefficient (A) of a buried anomalous body from the gravity data measured along a profile. The developed algorithm fits the observed data by a class of some geometrically simple anomalous bodies, including the semi-infinite vertical cylinder, infinitely long horizontal cylinder, and sphere models using the logarithms of the model parameters [log(z) and log(|A|)] rather than the parameters themselves in its iterative minimization scheme. The presented numerical experiments have shown that the original (non-logarithmed) minimization scheme, which uses the parameters themselves (z and |A|) instead of their logarithms, encountered a variety of convergence problems. The aforementioned transformation of the objective functional subjected to minimization into the space of logarithms of z and |A| overcomes these convergence problems. The reliability and the applicability of the developed algorithm have been demonstrated on several synthetic data sets with and without noise. It is then successfully and carefully applied to seven real data examples with bodies buried in different complex geologic settings and at various depths inside the earth. The method is shown to be highly applicable for mineral exploration, and for both shallow and deep earth imaging, and is of particular value in cases where the observed gravity data is due to an isolated body embedded in the subsurface.  相似文献   

6.
  相似文献   

7.
A theoretical model of grain size variation of domain transitions in titanomagnetite (x = 0.6) as a function of oxidation (z) is presented. The superparamagnetic (SP) to single-domain (SD) transition ds, the SD to two-domain (TD) transition d0, the TD to three-domain (3D) transition and the pseudo-single domain (PSD) to multi-domain (MD) transition are calculated as a function of z. It is shown that all the transition grain sizes increase with z, except for the PSD-MD transition for z > 0.6. The calculations predict that ds increases from 0.044 to 0.197 μm, d0 increases from 0.54 to 13 μm, the TD-3D transition increases from 1.6 to 49 μm as z varies from 0 to 0.8. The PSD-MD transition increases from 42 μm at z = 0 to 150 μm at z = 0.6, whereas between z = 0.6 to z = 0.8, the PSD-MD transition decreases to 49 μm. Qualitatively, the model explains some of the trends in magnetic properties of submarine basalts with low-temperature oxidation. Quantitatively, the model does give reasonable estimates of the PSD-MD boundary and d0, which are close to the experimental values for x = 0.6 and z = 0. Furthermore, the model predicts that psarks or two-domain grains could be the major contributors to the remanence of oxidized submarine pillow basalts.  相似文献   

8.
In a previous work the authors have developed a model, providing Kp as a function of the interplanetary magnetic field Bz component. They introduced a modified Bz function (denoted as Bzm), exhibiting a delayed reaction to Bz changes. The modified function Bzm was defined by using the analogy with a damping RC-circuit output voltage. The delaying reaction of Bzm to Bz was characterized by two time constants, one for rising and one for decreasing parts of Bz. The cross-correlation between Kp and Bzm has increased to 0.7, compared with −0.4 between Kp and Bz. In this paper, new dependences of Kp on solar wind velocity and dynamic pressure are included in the model to improve its accuracy. These solar wind parameters are found to correlate best with Kp. The hourly interpolated values are also added to the 3-h Kp values to increase the statistics. The new Kp data set is denoted as Kp1. The mean dependence of Kp on Bzm and dynamic pressure are approximated with parabolas, while the dependence on the velocity is linear. The constants in the model expression are obtained by using ACE data (1998–2000). The overall model error is estimated at 0.63 units Kp. The improvement over the previous simpler dependence in terms of the model error is about 30%.  相似文献   

9.
  相似文献   

10.
For assessing earthquake hazard of metro cities, knowledge of soil amplification, thickness and properties of sedimentary layer are essential. In order to map the soil thickness using microtremor survey method, in Bangalore city, it is required to calibrate the relation between fundamental resonance frequency of the soil layer and its thickness for the region. For this purpose microtremor survey was carried out at 34 locations in the city where borehole log was available. The resonance frequency of the soil is evaluated from the microtremor recordings using the H/V ratio technique. A nonlinear regression relation between the thickness of sedimentary layer h (m), from the borehole logs, and the resonance frequency fr (Hz), was derived as h=(58.3±8.8)fr−(0.95)±0.1. Using the model of shear wave velocity increasing with depth at these locations, the derived average shear wave velocity and the corresponding soil thickness were used, to get an empirical relation between VS (m/s) and depth z(m), as Vs=(174±28)(1+z)0.16±0.07. This relation also compares reasonably with the fit obtained between simulated VS and depth from borehole logs for Bangalore city. The calibrated relations can be used at locations in Bangalore city where borehole logs are not available, for finding the thicknesses and shear wave velocities of the local soil layers at the survey locations.  相似文献   

11.
In this work the surface temperature anomaly (dTG) and sunspot number (Rz) time series in the period 1880–2000 are studied with wavelet multi-resolution analysis. We found a very low correlation of 0.11 between dTG and Rz in the 11-yr-solar cycle band. A higher correlation of 0.66 is found in the ～22-yr-band with zero lag correlation coefficient between dTG and Rz. Furthermore, the long-term trend is markedly different between dTG and Rz. This might occurs because of the long-term warming on the last century, which is attributed mainly to anthropogenic effects.  相似文献   

12.
In this paper, a two-dimensional integrated numerical model is developed to examine the influences of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of a pipeline buried in a porous seabed. In the wave model, the RANS (Reynolds Averaged Navier–Stokes) equation is used to govern the wave motion. In the seabed model, the residual soil response in the vicinity of an embedded pipeline is considered with the 2-D elasto-plastic solution, where the phase-resolved shear stress is used as a source for the build-up of residual pore pressure. Classical Biot׳s consolidation equation is used for linking the solid-pore fluid interaction. The validation of the proposed integrated numerical model is conducted by the comparisons with the previous experimental data. Numerical examples show that the pore pressures can accumulate to a large value, thus resulting in a larger area of liquefaction potential in the given anisotropic soil compared to that with isotropic solution. The influences of anisotropic parameters on the wave-induced residual soil response in the vicinity of pipeline are significant. A high rate of pore pressure accumulation and dissipation is observed and the liquefaction potential develops faster as the anisotropic parameters increase. Finally, a simplified approximation based on a detailed parametric investigations is proposed for the evaluation of maximum liquefaction depth (zL) in engineering application.  相似文献   

13.
Applications of the entrainment process to layers at the boundary, which meet the self similarity requirements of the logarithmic profile, have been studied. By accepting that turbulence has dominating scales related in scale length to the height above the surface, a layer structure is postulated wherein exchange is rapid enough to keep the layers internally uniform. The diffusion rate is then controlled by entrainment between layers. It has been shown that theoretical relationships derived on the basis of using a single layer of this type give quantitatively correct factors relating the turbulence, wind and shear stress for very rough surface conditions. For less rough surfaces, the surface boundary layer can be divided into several layers interacting by entrainment across each interface. This analysis leads to the following quantitatively correct formula compared to published measurements. 1 $$\begin{gathered} \frac{{\sigma _w }}{{u^* }} = \left( {\frac{2}{{9Aa}}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \left( {1 - 3^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} \frac{a}{k}\frac{{d_n }}{z}\frac{{\sigma _w }}{{u^* }}\frac{z}{L}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ = 1.28(1 - 0.945({{\sigma _w } \mathord{\left/ {\vphantom {{\sigma _w } {u^* }}} \right. \kern-\nulldelimiterspace} {u^* }})({z \mathord{\left/ {\vphantom {z L}} \right. \kern-\nulldelimiterspace} L})^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ \end{gathered} $$ where \(u^* = \left( {{\tau \mathord{\left/ {\vphantom {\tau \rho }} \right. \kern-0em} \rho }} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} \) , σ w is the standard deviation of the vertical velocity,z is the height andL is the Obukhov scale lenght. The constantsa, A, k andd n are the entrainment constant, the turbulence decay constant, Von Karman's constant, and the layer depth derived from the theory. Of these,a andA, are universal constants and not empirically determined for the boundary layer. Thus the turbulence needed for the plume model of convection, which resides above these layers and reaches to the inversion, is determined by the shear stress and the heat flux in the surface layers. This model applies to convection in cool air over a warm sea. The whole field is now determined except for the temperature of the air relative to the water, and the wind, which need a further parameter describing sea surface roughness. As a first stop to describing a surface where roughness elements of widely varying sizes are combined this paper shows how the surface roughness parameter,z 0, can be calculated for an ideal case of a random distribution of vertical cylinders of the same height. To treat a water surface, with various sized waves, such an approach modified to treat the surface by the superposition of various sized roughness elements, is likely to be helpful. Such a theory is particularly desirable when such a surface is changing, as the ocean does when the wind varies. The formula, 2 $$\frac{{0.118}}{{a_s C_D }}< z_0< \frac{{0.463}}{{a_s C_D (u^* )}}$$ is the result derived here. It applies to cylinders of radius,r, and number,m, per unit boundary area, wherea s =2rm, is the area of the roughness elements, per unit area perpendicular to the wind, per unit distance downwind. The drag coefficient of the cylinders isC D . The smaller value ofz o is for large Reynolds numbers where the larger scale turbulence at the surface dominates, and the drag coefficient is about constant. Here the flow between the cylinders is intermittent. When the Reynolds number is small enough then the intermittent nature of the turbulence is reduced and this results in the average velocity at each level determining the drag. In this second case the larger limit forz 0 is more appropriate.  相似文献   

14.
Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling twodimensional wave fields and therefore cannot accurately simulate three-dimensional (3-D) full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity–stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.  相似文献   

15.
The factors influencing soil erosion may vary with scale. It remains unclear whether the spatial variation in soil erosion resistance is controlled by regional variables (e.g. precipitation, temperature, and vegetation zone) or by local specific variables (e.g. soil properties, root traits, land use, and farming operations) when the study area enlarges from a hillslope or catchment to the regional scale. This study was performed to quantify the spatial variations in soil erosion resistance to flowing water under three typical land uses along a regional transect on the Loess Plateau and to identify whether regional or local specific variables are responsible for these changes. The results indicated that the measured soil detachment capacities (Dc) of cropland exhibited an irregular trend along the regional transect. The Dc of grassland increased with mean annual precipitation, except for two sites (Yijun and Erdos). The measured Dc of woodland displayed an inverted ‘U’ shape. The changes in rill erodibility (Kr) of three land uses were similar to Dc, whereas no distinguishable trend was found for critical shear stress (τc). No significant correlation was detected between Dc, Kr and τc, and the regional variables. The spatial variation in soil erosion resistance could be explained reasonably by changes in soil properties, root traits, land use, and farming operations, rather than regional variables. The adjustment coefficient of Kr for grassland and woodland could be well simulated by soil cohesion and root mass density (R2 = 0.70, P < 0.01), and the adjustment coefficient of critical shear stress could be estimated with aggregate stability (R2 = 0.57, P < 0.01). The results are helpful for quantifying the spatial variation in soil detachment processes by overland flow and to develop process‐based erosion model at a regional scale. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

16.
The use of the shear wave velocity data as a field index for evaluating the liquefaction potential of sands is receiving increased attention because both shear wave velocity and liquefaction resistance are similarly influenced by many of the same factors such as void ratio, state of stress, stress history and geologic age. In this paper, the potential of support vector machine (SVM) based classification approach has been used to assess the liquefaction potential from actual shear wave velocity data. In this approach, an approximate implementation of a structural risk minimization (SRM) induction principle is done, which aims at minimizing a bound on the generalization error of a model rather than minimizing only the mean square error over the data set. Here SVM has been used as a classification tool to predict liquefaction potential of a soil based on shear wave velocity. The dataset consists the information of soil characteristics such as effective vertical stress (σ′v0), soil type, shear wave velocity (Vs) and earthquake parameters such as peak horizontal acceleration (amax) and earthquake magnitude (M). Out of the available 186 datasets, 130 are considered for training and remaining 56 are used for testing the model. The study indicated that SVM can successfully model the complex relationship between seismic parameters, soil parameters and the liquefaction potential. In the model based on soil characteristics, the input parameters used are σ′v0, soil type, Vs, amax and M. In the other model based on shear wave velocity alone uses Vs, amax and M as input parameters. In this paper, it has been demonstrated that Vs alone can be used to predict the liquefaction potential of a soil using a support vector machine model.  相似文献   

17.
The negative and positive fronts of the IMF B z component arrived at intervals of 3 h during a strong magnetic storm of May 15, 2005. The occurrence of Pc5 pulsations at these three characteristic instants has been considered based on the WIND satellite magnetic data. Pulsations originated not only during sudden compression SC of the magnetosphere but also during the B z sign reversal from positive to negative. The IMF B z sign reversal from negative to positive did not affect the development of pulsations. It is assumed that Pc5 pulsations observed after the negative IMF B z front are related to the development of surface waves at the magnetopause as a result of impulsive reconnection of field lines.  相似文献   

18.
针对工程中大量存在的平面应变问题,依据平面应变条件和广义胡克定律,基于SMP、Lade-Duncan、AC-SMP和广义Mises强度准则,推导出考虑中主应力及泊松比影响的无黏性土主、被动土压力计算公式,并将其扩展至黏性土,讨论基于各强度准则土压力计算公式的适用范围。结果表明：考虑中主应力对土强度的贡献后,基于各强度准则所计算的主动土压力均小于朗肯主动土压力,被动土压力均大于朗肯被动土压力;主动土压力Pa随着泊松比的增大而减小,被动土压力PP随着泊松比的增大而增大,且泊松比越大,与实测数据更为接近;基于同一强度准则下得到的主、被动土压力适用的内摩擦角范围随着泊松比的增大而增大;基于各强度准则的土压力计算公式均能较好的描述挡土结构上土压力的大小,其中广义Mises强度准则计算结果与实际工程更为吻合,研究成果可为挡土结构上土压力的计算提供一定理论参考。  相似文献   

19.
The results of a laboratory electromagnetic analogue model study, which employs a horizontal inducing field over a simple model of the British Isles region, delineate the location and frequency dependence of the major coast effect induction anomalies of the Scotland region. Contours of amplitudes, amplitude ratios, and in-phase and quadrature parts of the model field measurements are presented. The model vertical magnetic fields for two orthogonal source field polarizations and field station values for two hypothetical events for corresponding polarizations are compared.While major discrepancies occur between model and field Hz amplitudes, the Hz gradients across Scotland, which can be attributed to the coast-effect, are comparable in value, although sometimes reversed in sign. Superimposed on this coast-effect, the field data indicate the existence of current concentrations associated with the Great Glen and the Southern Uplands faults and possibly also of currents within the Scottish mainland near the east and west coasts.  相似文献   

20.
兰州地区红层泥岩形成于干旱、半干旱环境,其物理力学特性与其他地区的红层泥岩不同。以G6京藏高速兰海养护维修工程的红层泥岩路基填料为研究对象,通过大量的常规三轴试验,研究含水率对其力学特性及邓肯-张模型参数的影响。结果表明:在不同含水率下兰州地区红层泥岩的应力-应变曲线关系符合双曲线模型。含水率小于最优含水率时,破坏形式为剪切破坏;含水率大于最优含水率,破坏形式为鼓状破坏;随着含水率的增加,黏聚力和内摩擦角逐渐减小,内摩擦角减小的幅度更大,究其原因是含水率对黏聚力和内摩擦角的影响机理不同而导致。采用数学模型进行拟合发现:黏聚力、内摩擦角与含水率分别呈二次抛物线与对数曲线关系。不同含水率下邓肯-张模型参数破坏比R_f在区间(0.869,0.984)内波动,K值随含水率增大而减小呈负线性关系,n值随含水率的增大而增大呈线性关系。  相似文献   

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

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