首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 560 毫秒
1.
During recent gravity surveys in Patagonia, we were unable to acquire data on the North Patagonian Icecap, which is completely inaccessible for land surveying. The Icecap is clearly of tectonic significance, but, because of its tectonic history, it is uplifted and rugged, and consequently inaccessible. Therefore the distribution of geophysical observations which can be acquired in a rugged field area is dependent on the tectonic history. This suggests that a tectonic history which is reconstructed from geophysical data may be systematically biased. We use models of local and flexural isostasy to estimate the gravity errors associated with interpolation across inaccessible topography. The gravity error is largest for Pratt isostasy, where the mass deficit which supports the topography is at relatively shallow depths. The gravity error is least for flexural isostasy, because in this case the inaccessible topography is supported regionally by a mass deficit which extends beyond the inaccessible region. An analytical flexural interpolation scheme is proposed for interpolation across data gaps associated with inaccessible topography. Flexural theory and Gauss’s theorem are used to predict the gravity anomaly due to the mass excess of the inaccessible topography. We apply this scheme to the North Patagonian Icecap, to demonstrate that flexural interpolation predicts a relative gravity low at the site of the Icecap, which would not be predicted by purely geometrical interpolation schemes.  相似文献   

2.
气象台站20 cm蒸发皿观测资料自然正交分解显示,1980～2000年中国区域气温显著增加期间,长江中游至河套、东北等区域地表年蒸发潜力呈增加趋势;相反在长江以南、东部和西南等地区年蒸发潜力呈下降趋势.辐射观测资料分析结果表明,自20世纪70年代中国区域太阳入射能整体呈下降趋势,因此对于蒸发潜力增加的地区,太阳辐射产生...  相似文献   

3.
Estimating transpiration of the trees in agroforestry system is important in water management of the site. Sap flow of intercropped fast‐growing young poplar trees and microclimate factors in semiarid northeastern China was measured in two growing seasons (2008 and 2009). Sapwood growth and water storage of wood and leaf increment during the growing season were involved in the calculation of sap flow. The results showed that diurnal variation of sap flow followed to that of short wave solar radiation. Sap flows both in 10 min mean and daily gross values mainly depended on solar radiation and vapor pressure deficit, and the relations well fit hyperbolic function. The regression coefficients of monthly window data indicated that the seasonal variation of sap flow capacities decreased gradually from June to September. Moderate soil water stress of upper soil layer (0–50 cm) did not constrain the sap flow because the trees could use the water at deeper soil layer. The daily sap flow per tree ranged 0.8 to 18.1 and 3.7 to 23.8 kg d?1 tree?1, with averages of 8.7 and 14.3 kg d?1 tree?1 in 2008 and 2009 respectively. An empirical model was established to estimate the sap flow of the poplar trees by solar radiation, vapor pressure deficit, leaf area index and Julian days. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
The numerical dispersion and computational cost are high for conventional Taylor series expansion staggered-grid finite-difference forward modeling owing to the high frequency of the wavelets and the large grid intervals. In this study, the cosine-modulated binomial window function (CMBWF)-based staggered-grid finite-difference method is proposed. Two new parameters, the modulated time and modulated range are used in the new window function and by adjusting these two parameters we obtain different characteristics of the main and side lobes of the amplitude response. Numerical dispersion analysis and elastic wavefield forward modeling suggests that the CMBWF method is more precise and less computationally costly than the conventional Taylor series expansion staggered-grid finitedifference method.  相似文献   

5.
The field deployment of a heated distributed temperature sensor (DTS) for over three years has revealed two obstacles to estimating soil moisture (θ) that may hamper subsurface DTS applications as well as use of other subsurface thermal probes. The first observed obstacle was a hysteretic response of the DTS sensor. The relationship between θ and the temperature response (?T) within the cable was not only dependent on θ of the soil, but also on the previous wetting and drying cycles leading to that state. The second observed obstacle was soil structure healing. Soil structure healing causes the relationship between ?T and θ to evolve through time; this calibration curve becomes flatter, or less sensitive, as the surrounding soil makes better contact with the cable. Effects of the hysteretic response of the instrument and soil structure healing are largely the result of small gaps between the cable and soil. These small gaps can be approximated by a contact resistance between the cable and soil. In this article we characterize the occurrence of hysteretic and soil structure healing effects from field data and parameterize contact resistance by simulating heat transfer using a numerical modelling approach Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.  相似文献   

7.
Multiscale analysis of vegetation surface fluxes: from seconds to years   总被引：2，自引：0，他引：2  
The variability in land surface heat (H), water vapor (LE), and CO2 (or net ecosystem exchange, NEE) fluxes was investigated at scales ranging from fractions of seconds to years using eddy-covariance flux measurements above a pine forest. Because these fluxes significantly vary at all these time scales and because large gaps in the record are unavoidable in such experiments, standard Fourier expansion methods for computing the spectral and cospectral statistical properties were not possible. Instead, orthonormal wavelet transformations are proposed and used. The are ideal at resolving process variability with respect to both scale and time and are able to isolate and remove the effects of missing data (or gaps) from spectral and cospectral calculations. Using the spectra, we demonstrated unique aspects in three appropriate ranges of time scales: turbulent time scales (fractions of seconds to minutes), meteorological time scales (hour to weeks), and seasonal to interannual time scales corresponding to climate and vegetation dynamics. We have shown that: (1) existing turbulence theories describe the short time scales well, (2) coupled physiological and transport models (e.g. CANVEG) reproduce the wavelet spectral characteristics of all three land surface fluxes for meteorological time scales, and (3) seasonal dynamics in vegetation physiology and structure inject strong correlations between land surface fluxes and forcing variables at monthly to seasonal time scales. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the flux variability is dominated by a narrow band or bands of modes, and these modes “resonate” with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE.  相似文献   

8.
Various numerical methods have been used in the literature to simulate single and multiphase flow in fractured media. A promising approach is the use of the discrete-fracture model where the fracture entities in the permeable media are described explicitly in the computational grid. In this work, we present a critical review of the main conventional methods for multiphase flow in fractured media including the finite difference (FD), finite volume (FV), and finite element (FE) methods, that are coupled with the discrete-fracture model. All the conventional methods have inherent limitations in accuracy and applications. The FD method, for example, is restricted to horizontal and vertical fractures. The accuracy of the vertex-centered FV method depends on the size of the matrix gridcells next to the fractures; for an acceptable accuracy the matrix gridcells next to the fractures should be small. The FE method cannot describe properly the saturation discontinuity at the matrix–fracture interface. In this work, we introduce a new approach that is free from the limitations of the conventional methods. Our proposed approach is applicable in 2D and 3D unstructured griddings with low mesh orientation effect; it captures the saturation discontinuity from the contrast in capillary pressure between the rock matrix and fractures. The matrix–fracture and fracture–fracture fluxes are calculated based on powerful features of the mixed finite element (MFE) method which provides, in addition to the gridcell pressures, the pressures at the gridcell interfaces and can readily model the pressure discontinuities at impermeable faults in a simple way. To reduce the numerical dispersion, we use the discontinuous Galerkin (DG) method to approximate the saturation equation. We take advantage of a hybrid time scheme to alleviate the restrictions on the size of the time step in the fracture network. Several numerical examples in 2D and 3D demonstrate the robustness of the proposed model. Results show the significance of capillary pressure and orders of magnitude increase in computational speed compared to previous works.  相似文献   

9.
Subslab soil gas sampling and analysis is a common line of evidence for assessing human health risks associated with subsurface vapor intrusion to indoor air for volatile organic compounds; however, conventional subslab sampling methods have generated data that show substantial spatial and temporal variability, which often makes the interpretation difficult. A new method of monitoring has been developed and tested that is based on a concept of integrating samples over a large volume of soil gas extracted from beneath the floor slab of a building to provide a spatially averaged subslab concentration. Regular field screening is also conducted to assess the trend of concentration as a function of the volume removed to provide insight into the spatial distribution of vapors at progressive distances away from the point of extraction. This approach minimizes the risk of failing to identify the areas of elevated soil vapor concentrations that may exist between discrete sample locations, and can provide information covering large buildings with fewer holes drilled through the floor. The new method also involves monitoring the extraction flow rate and transient vacuum response for mathematical analysis to help interpret the vapor concentration data and to support an optimal design for any subslab venting system that may be needed.  相似文献   

10.
Increases in atmospheric CO2 concentration not only affects climate variables such as precipitation and air temperature, but also affects intrinsic ecosystem physiological properties such as bulk stomatal conductance and intercellular CO2 concentration. De-convolving these two effects remains uncertain in biosphere–atmosphere water and carbon cycling. Using a simplified analytical net ecosystem CO2 exchange (NEE) model, tested with recently collected flux measurements in a humid grassland ecosystem in Ireland, we assess how much projected climate shifts affect net canopy photosynthesis (A) without physiological adjustments and contrast those findings with published field data on physiological adjustments for several grassland ecosystems. Our analysis suggests that the intrinsic grassland ecosystem physiological adjustment of A is about 45 times more important than the resulting climatic forcing shifts from the IS92a scenario (and a doubling of atmospheric CO2 concentration). Also, our analysis shows that increase in precipitation results in concomitant decrease in the two climate variables—net radiation and vapor pressure deficit, and these decreases have opposite (and almost canceling) effects on A. Implications to afforestation policy and future experimental efforts to quantify the carbon sink from humid grassland ecosystems are also discussed.  相似文献   

11.
Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

12.
Very early times in the order of 2–3 μs from the end of the turn‐off ramp for time‐domain electromagnetic systems are crucial for obtaining a detailed resolution of the near‐surface geology in the depth interval 0–20 m. For transient electromagnetic systems working in the off time, an electric current is abruptly turned off in a large transmitter loop causing a secondary electromagnetic field to be generated by the eddy currents induced in the ground. Often, however, there is still a residual primary field generated by remaining slowly decaying currents in the transmitter loop. The decay disturbs or biases the earth response data at the very early times. These biased data must be culled, or some specific processing must be applied in order to compensate or remove the residual primary field. As the bias response can be attributed to decaying currents with its time constantly controlled by the geometry of the transmitter loop, we denote it the ‘Coil Response’. The modelling of a helicopter‐borne time‐domain system by an equivalent electronic circuit shows that the time decay of the coil response remains identical whatever the position of the receiver loop, which is confirmed by field measurements. The modelling also shows that the coil response has a theoretical zero location and positioning the receiver coil at the zero location eliminates the coil response completely. However, spatial variations of the coil response around the zero location are not insignificant and even a few cm deformation of the carrier frame will introduce a small coil response. Here we present an approach for subtracting the coil response from the data by measuring it at high altitudes and then including an extra shift factor into the inversion scheme. The scheme is successfully applied to data from the SkyTEM system and enables the use of very early time gates, as early as 2–3 μs from the end of the ramp, or 5–6 μs from the beginning of the ramp. Applied to a large‐scale airborne electromagnetic survey, the coil response compensation provides airborne electromagnetic methods with a hitherto unseen good resolution of shallow geological layers in the depth interval 0–20 m. This is proved by comparing results from the airborne electromagnetic survey to more than 100 km of Electrical Resistivity Tomography measured with 5 m electrode spacing.  相似文献   

13.
Understanding the principal causes and possible solutions for groundwater depletion in India is important for its water security, especially as it relates to agriculture. A study was conducted in an agricultural watershed in Andhra Pradesh, India to assess the impacts on groundwater of current and alternative agricultural management. Hydrological simulations were used as follows: (1) to evaluate the recharge benefits of water‐harvesting tillage through a modified Soil and Water Assessment Tool (SWAT) model and (2) to predict the groundwater response to changing extent and irrigation management of rice growing areas. The Green–Ampt infiltration routine was modified in SWAT was modified to represent water‐harvesting tillage using maximum depression storage parameter. Water‐harvesting tillage in rainfed croplands was shown to increase basin‐scale groundwater recharge by 3% and decrease run‐off by 43% compared with existing conventional tillage. The groundwater balance (recharge minus irrigation withdrawals), negative 11 mm/year under existing management changed to positive (18–45 mm/year) when rice growing areas or irrigation depths were reduced. Groundwater balance was sensitive to changes in rice cropland management, meaning even small changes in rice cropland management had large impacts on groundwater availability. The modified SWAT was capable of representing tillage management of varying maximum depression storage, and tillage for water‐harvesting was shown to be a potentially important strategy for producers to enhance infiltration and groundwater recharge, especially in semi‐arid regions where rainfall may be becoming increasingly variable. This enhanced SWAT could be used to evaluate the landscape‐scale impacts of alternative tillage management in other regions that are working to develop strategies for reducing groundwater depletion. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
基于局部斜率属性的VSP波场分离研究   总被引：6，自引：2，他引：4       下载免费PDF全文
基于垂直地震剖面(VSP)资料中上、下行波视速度的差异,利用地震剖面同相轴局部斜率属性参数,提出了一种分离上、下行波场的新方法.首先利用Fourier变换初步分离上、下行波场,然后利用平面波分解滤波器(Plane Wave Destruction (PWD) Filter)技术估计初始分离波场的同相轴局部斜率属性参数,在此基础上对VSP原始资料波场分离.该方法是一种时间域最小平方优化分离波场的方法,不存在其他滤波方法阈值滤波器边界的影响,减少了因镶边问题带来的假象.模拟和实际资料处理结果表明,该方法与传统方法相比,分离出的上、下行波噪声假象少,振幅保持好,更好地消除了上、下行波的相互影响.  相似文献   

15.
The rough‐sea reflection‐response varies (1) along the streamer (2) from shot to shot and (3) with time along the seismic trace. The resulting error in seismic data can be important for time‐lapse imaging. One potential way of reducing the rough‐sea receiver error is to use conventional statistical deconvolution, but special care is needed in the choice of the design and application windows. The well‐known deconvolution problem associated with the non‐whiteness of the reflection series is exacerbated by the requirement of an unusually short design window – a requirement that is imposed by the non‐stationary nature of the rough‐sea receiver wavelet. For a synthetic rough‐sea data set, with a white 1D reflection series, the design window needs to be about 1000 ms long, with an application window about 400 ms long, centred within the design window. Although such a short design window allows the deconvolution operator to follow the time‐variation of the rough‐sea wavelet, it is likely to be too short to prevent the non‐whiteness of the geology from corrupting the operator when it is used on real data. If finely spatial‐sampled traces are available from the streamer, the design window can be extended to neighbouring traces, making use of the spatial correlations of the rough‐sea wavelet. For this ‘wave‐following’ approach to be fruitful, the wind (and hence the dominant wave direction) needs to be roughly along the line of the streamer.  相似文献   

16.
To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.  相似文献   

17.
The results of comprehensive field testing of on‐site vapor‐phase‐based groundwater monitoring methods are presented to demonstrate their utility as a robust and cost‐effective approach for rapidly obtaining volatile organic compounds (VOCs) concentration data from a monitoring well. These methods—which rely on sensitive, commercially available field equipment to analyze vapor in equilibrium with groundwater—proved easy to implement and can be tailored to site‐specific needs, including multilevel sampling. During field testing, low‐flow groundwater concentrations could be reasonably estimated using submerged passive vapor diffusion samplers or field equilibration of collected groundwater (R2 = 0.85 to 0.96). These two methods are not as reliant on in‐well mixing to overcome vertical stratification within wells as simpler headspace methods. The importance of well and aquifer‐specific factors on concentration data (and therefore method selection) is highlighted, including the effect of changing in‐well patterns due to seasonal temperature gradients. Results indicated that vertical stratification was relatively limited within the set of wells included in these studies, resulting in similar performance for short depth‐discrete passive vapor diffusion samplers (constructed from 40‐mL vials) and longer samplers (2.5 to 5 feet in length) designed to cover a larger portion of the screened interval. A year‐long, multi‐event evaluation demonstrated that vapor‐phase‐based monitoring methods are no more variable than conventional groundwater monitoring methods, with both types subject to similar spatial and temporal variability that can be difficult to reduce. Vapor sampling methods represent a promising approach for estimation of groundwater concentrations by reducing the cost liabilities associated with monitoring while providing a more sustainable approach.  相似文献   

18.
Generic indoor air:subslab soil gas attenuation factors (SSAFs) are important for rapid screening of potential vapor intrusion risks in buildings that overlie soil and groundwater contaminated with volatile chemicals. Insufficiently conservative SSAFs can allow high‐risk sites to be prematurely excluded from further investigation. Excessively conservative SSAFs can lead to costly, time‐consuming, and often inconclusive actions at an inordinate number of low‐risk sites. This paper reviews two of the most commonly used approaches to develop SSAFs: (1) comparison of paired, indoor air and subslab soil gas data in empirical databases and (2) comparison of estimated subslab vapor entry rates and indoor air exchange rates (IAERs). Potential error associated with databases includes interference from indoor and outdoor sources, reliance on data from basements, and seasonal variability. Heterogeneity in subsurface vapor plumes combined with uncertainty regarding vapor entry points calls into question the representativeness of limited subslab data and diminishes the technical defensibility of SSAFs extracted from databases. The use of reasonably conservative vapor entry rates and IAERs offers a more technically defensible approach for the development of generic SSAF values for screening. Consideration of seasonal variability in building leakage rates, air exchange rates, and interpolated vapor entry rates allows for the development of generic SSAFs at both local and regional scales. Limitations include applicability of the default IAERs and vapor entry rates to site‐specific vapor intrusion investigations and uncertainty regarding applicability of generic SSAFs to assess potential short‐term (e.g., intraday) variability of impacts to indoor air.  相似文献   

19.
Trend analysis of reference evapotranspiration (ET0), as a key factor in irrigation programming, has an important role in water resources management. Many parameters affect ET0 and their variations can change its values. In this paper, the effect of temporal variation of meteorological variables including wind speed, temperature, solar radiation and saturation vapor pressure deficit on temporal variations of ET0 was analyzed. Trend analysis of ET0 and its more effective meteorological parameters was accomplished in 30 synoptic stations which are located in Iran using Spearman’s Rho test. The multiple linear regressions were also used to determine the relationship between ET0 trend and the trend of its more effective parameters. Increasing and decreasing trends in ET0 were obtained at annual and seasonal scales. Many studied stations which had decreasing trend in the annual and seasonal periods have been located in the arid climates while all stations which have been located in humid and very-humid climates, had an increasing trend in annual and seasonal periods. The trend results in studied variables showed that annual and seasonal values of wind speed, temperature and saturation vapor pressure deficit decrease however the values of solar radiation increases in most studied stations. Multiple linear regressions results demonstrated that ET0 trend can be calculated by the trend of two more effective variables including wind speed and saturation vapor pressure deficit.  相似文献   

20.
快速有效的转换波共转换点叠加技术   总被引：19，自引：8，他引：19       下载免费PDF全文
转换波共转换点(Common Conversion Point简称CCP)叠加的关键在于CCP抽道集和非双曲线正常时差校正NMO(Normal MoveOut). 目前方法的精度限制了其在中-浅层或大炮检距情况下的应用. 我们对CCP叠加技术进行了系统研究,导出了新的CCP位置计算公式和非双曲线时距关系式,并给出了具体的CCP抽道集方法. 理论模型试验和实验资料处理表明本文方法精度高、简便易行,特别是对于中-浅层和大炮检距情况也能得到良好的叠加效果.  相似文献   

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

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