Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

三维高阶深度偏移方程及其数值求解方法

从三维非均匀介质中的波动方程出发,利用拟微分算子理论,Ｐａｄｅ逼近方法及因式分解技巧,获得了非均匀介质的三维高阶深度偏移方程,相应地提出了逐次低阶方程方法、低阶方程组方法及分裂方法等３种求解方法．与二维情形不同,以上每一种方法在数值求解时均存在由测线坐标ｙ的出现而带来的困难．为了克服这一困难,我们提出了差分算子分解方法,避免了近年来人们竞相研究的ｘ－ｙ方向微分算子分裂带来的分裂误差,保持了应有的相容性,解决了这一令人烦恼的问题．     

地质雷达方法在上海市轨道交通明珠线工程中的应用

地质雷达方法在上海市轨道交通明珠线工程中的应用表明,此方法在探测地下障碍物中能依靠雷达图象识别其类型、材质、深度等,有利于快速、准确定位,有一定的优越性。     

A search for a problem free cosmology

We explore features of a nonminimally coupled theory of scalar fields with an effective potential that supports non topological soliton solutions. It is suggested that a problem free cosmology results.     

人为因素导致的地面沉降

天然气、石油、地下水、地热热水和卤水等流体运移出地下储存地层时,地层耗尽了流体而产生压缩变形,这些变形传递到地表表面就形成了人为的地面沉降。在本文中:(1)列举了世界上主要的地面沉降区域;(2)论述了引发地面沉降的力学机制,这些可以量测的地面沉降发生在地下含水层体系、天然气气田和石油油田之上;(3)描述了目前可以用于地面沉降测量和岩石变形原位测试的技术手段;(4)简要介绍了几个地面沉降预测的数值模型;(5)说明了用于控制地面沉降发展和减轻相关环境影响的几种防治措施。     

Statistics of extremes and estimation of extreme rainfall: II. Empirical investigation of long rainfall records / Statistiques de valeurs extrêmes et estimation de précipitations extrêmes: II. Recherche empirique sur de longues séries de précipitations

Abstract

Abstract In the first part of this study, theoretical analyses showed that the Gumbel distribution is quite unlikely to apply to hydrological extremes and that the extreme value distribution of type II (EV2) is a more consistent choice. Based on these theoretical analyses, an extensive empirical investigation is performed using a collection of 169 of the longest available rainfall records worldwide, each having 100–154 years of data. This verifies the theoretical results. In addition, it shows that the shape parameter of the EV2 distribution is constant for all examined geographical zones (Europe and North America), with value κ = 0.15. This simplifies the fitting and the general mathematical handling of the distribution, which become as simple as those of the Gumbel distribution.     

Dynamic Behaviour of a Caisson Quaywall Under Strong Earthquake

—Effective stress analysis is performed to evaluate the residual displacement of a caissonquaywall during 1994 Hokkaido-Toho-Oki Earthquake and 1993 Kushiro-Oki Earthquake.Theconstitutive model used in this study is a multiple shear mechanism type defined in strain space and cantake into account the effect of rotation of principal stress axis.The earthquake accelerations recorded atthe outcropping rock during the earthquake are used as input bedrock motion.The results of finite ele-ment analysis are in good agreement with the observed behaviour of the quaywall.The analysis also indi-cates that liquefaction and high excess pore water pressure have a significant effect on the deformation ofthe caisson.Soil improvement is speculated as the most reliable measures against liquefaction.The influ-ence of soil improvement and the reasonable improved area are discussed in the paper.     

Toward an improved global network for determination of tropospheric ozone climatology and trends

An examination of typical tropospheric ozone variability on daily, monthly, annual and interannual timescales and instrumental precision indicates that the current ozonesonde network is insufficient to detect a trend in tropospheric ozone of 1% per year at the 2 level even at stations with records a decade in length. From a trend prediction analysis we conclude that in order to detect a 1% per year trend in a decade or less it will be necessary to decrease the time between observations from its present value of 3–7 days to 1 day or less. The spatial distribution of the current ozonesonde stations is also inadequate for determining the global climatology of ozone. We present a quantitative theory taking into account photochemistry, surface deposition, and wind climatology to define the effectively sampled region for an observing station which, used in conjunction with the instrumental precision and the above prediction analysis, forms the basis for defining a suitable global network for determining regional and global ozone climatology and trends. At least a doubling of the present number of stations is necessary, and the oceans, most of Asia, Africa, and South America are areas where more stations are most needed. Differential absorption lidar ozone instruments have the potential for far more frequent measurements of ozone vertical profiles and hence potentially more accurate climatology and trend determinations than feasible with ozonesondes but may produce a (fair weather) biased data set above the cloud base. A strategy for cloudy regions in which either each station utilizes both lidars and sondes or each station is in fact a doublet comprised of a near-sea-level lidar and a proximal-mountain-top lidar could serve to minimize this bias.     

The Darwin Rise demise: the western Pacific guyot heights trace the trans-Pacific Mendocino fracture zone

The Darwin Rise has been proposed so many times and in so many forms and places that the time has come to make a more comprehensive examination of the region. Lying on the NW Pacific Plate between the Geisha Guyots, the Mid-Pacific Mountains, the equator, and the trenches, the region is roughly bounded by magnetic anomaly M20 (147 Ma). It was subjected to a massive outpouring of lava about 105 to 120 Ma, which created the guyots and seamounts in that region. Guyots are excellent tools for studying events of long ago because they eroded in the same lowstand in the Cretaceous and guyot relief, therefore, is a surrogate for paleo-sealevel. The relief is derived by subtracting the break depth of the summit plateau of a guyot from the regional depth. Guyot relief would necessarily be less in the center than to the periphery if the feature formed on a pre-existing rise, as has been postulated. The existence of a paleo-Darwin Rise would give concentric contours for the region in question. Of the sixty guyots used in this study, thirty-seven of these guyots were surveyed using SASS multibeam in the Marcus-Wake seamount group. Twenty-three guyots were surveyed using random track single-beam sonar surveys. An entirely different scenario is shown. Data revealed a major fracture passing through the area coevally or after the guyots formed. Because the depths to the summit are not the same now, vertical tectonics occurred after subaerial erosion. This means the fracture formed during and after the erosion (roughly 105 Ma) and influenced the normal sequence of events in guyot formation. Depending on how one deciphers trends through the Hess Rise morass, SASS bathymetry shows a continuation of the Surveyor/Mendocino fracture zone swarm inside the M20 region to the NE of these data. The fracture swarm continues to the western Pacific trench system. Based on this information, if the Darwin Rise ever existed, it had to have done so elsewhere.     

Trend detection in river flow series: 1. Annual maximum flow / Détection de tendance dans des séries de débit fluvial: 1. Débit maximum annuel

Abstract

Results of a study on change detection in hydrological time series of annual maximum river flow are presented. Out of more than a thousand long time series made available by the Global Runoff Data Centre (GRDC) in Koblenz, Germany, a worldwide data set consisting of 195 long series of daily mean flow records was selected, based on such criteria as length of series, currency, lack of gaps and missing values, adequate geographical distribution, and priority to smaller catchments. The analysis of annual maximum flows does not support the hypothesis of ubiquitous growth of high flows. Although 27 cases of strong, statistically significant increase were identified by the Mann-Kendall test, there are 31 decreases as well, and most (137) time series do not show any significant changes (at the 10% level). Caution is advised in interpreting these results as flooding is a complex phenomenon, caused by a number of factors that can be associated with local, regional, and hemispheric climatic processes. Moreover, river flow has strong natural variability and exhibits long-term persistence which can confound the results of trend and significance tests.