有限元集中质量法在非饱和土壤水流中的应用

利用Ｇａｌｅｒｋｉｎ有限元集中质量法,讨论了非饱和土壤水流动的数值计算模型,用以模拟均质土壤,地下水埋藏很深且在不同的初始和边界条件下的水分运动。通过将第一类边界条件齐次化,而对于第二类非线性边界条件采用变分的办法将其直接转化为已知边界通量的计算,从而容易且较好地处理了带通量约束边界条件。对于入渗或蒸发问题,用通常的离散方法容易产生振荡非物理解。我们采用有限元集中质量法,将质量矩阵的非对角元集中到主对角元,避免了这种非物理解的产生。利用该方案统一地模拟如下情况：稳定入渗和蒸发、蒸发一段时间后达到风干率,以及波动入渗蒸发时土壤水势的变化状况。模拟表明,在已知有关土壤参数的情况下,该模型能较好地模拟非饱和流水势及含水量分布。     

陆面过程模式中地下水位的参数化及初步应用

田间研究表明地表水和地下水有重要的相互作用,它与土壤含水量密切相关.土壤含水量不仅在陆气相互作用系统水和能量平衡中,而且在干旱、洪水预报、水资源管理、生态系统研究中起十分重要的作用.因此,研究地表水和地下水的相互作用,建立陆面模式中地下水位的动态表示,对于气候与水资源研究具有重要意义.将地下水位的动态表示问题归结为饱和与非饱和流问题,发展了其数值计算方案,建立了地下水位的动态表示,并与陆面过程模型耦合,建立了陆气相互作用中地下水位的动态表示,并进行了数值模拟研究.     

非饱和土壤水流方程的极值原理

讨论了非饱和土壤水流问题及其数值计算方案的极值原理,在理论上证明了连续问题及其数值格式均满足极值原理,即在一定的条件下,预报变量的最大值和最小值只能在边界上达到,从而描述了特定的物理性质.这些问题在讨论入渗、蒸发问题的研究中是重要的.     

土壤湿度初值对边界层物理量预报影响的分析

使用全球中期数值预报模式T213L31,对不同的土壤湿度初值对边界层物理量预报的影响进行敏感试验分析。试验表明土壤湿度初值的准确性非常重要,其变化导致模式计算的感热通量、潜热通量、地表温度和2m温度发生较大变化。而当使用与观测值相近的土壤湿度初值时得到的边界层物理量(如2m温度)更加准确。通过使用固定值与6h预报的背景场加权平均作为土壤湿度的初值,改进原有采用固定值的方法,经过检验模式系统,2m温度的预报误差明显减小。     

脉冲降雨—蒸发对土壤湿度影响的动力机制分析

土壤湿度直接作用于植被的呼吸、蒸腾和各种化学反应, 是生态气象学的重要的环境因素。降雨和蒸发作为土壤湿度的两个重要扰动因子, 对其进行定性的理论分析及定量的数值分析具有重要意义, 尤其在半干旱和干旱地区。蒸发对土壤湿度具有持续的耗散作用, 可以视为连续的, 而降雨作为湿度补充作用具有脉冲效应。为阐释降雨和蒸发对土壤湿度的影响, 该文通过建立脉冲降雨—蒸发对土壤湿度响应模型, 并对其进行定性分析, 得到此微分方程系统的周期解的充分条件。在此基础上, 通过Forcal对方程参数拟合, 并以Maple作为数值分析工具, 对此系统进行数值模拟。      

干旱及半干旱区土壤水热传输模式研究

本文发展了一个干旱地区土壤中水分及热量传输模式。首先从土壤孔隙内水汽运动研究入手,阐明了在土壤深层孔隙中水汽压与液面饱和值之间处在平衡态,而土壤表层则处在非平衡态,因而必要在计算表层蒸发的公式中引入必要的土壤表层阻力,才可使模式简单实用。在此基础上建立了考虑液态水及汽态水运动耦合的多层模式,并用HEIFE沙漠站的资料对模式进行了验证,结果表明此模式较好地再现土壤内及地气界面上的水热交换过程,并且也表明干旱地区土壤中水蒸汽输送对水分平衡及蒸发的计算是重要的。这种模式很易推广到气候研究的干旱土壤下垫面的模式中去。     

集合Kalman滤波在土壤湿度同化中的应用

基于非饱和土壤水模型和集合卡尔曼滤波 (Ensemble Kalman Filter, 简称EnKF) 并结合陆面水文模型——可变下渗能力模型 (Variable Infiltration Capacity, 简称VIC模型) 发展了一个土壤湿度同化方案。利用1998年6~8月淮河流域能量和水循环试验 (HUBEX) 项目外场观测试验区——史灌河流域梅山站土壤湿度逐日观测资料及1986~1993年合肥和南阳两站点的土壤湿度旬观测资料进行同化试验, 结果表明该同化方案能完整估计土壤湿度廓线, 同化的土壤湿度与观测资料基本吻合, 反映了土壤湿度的日、 旬、 月、 季变化, 同化方案是合理的。与基于扩展卡尔曼滤波 (Extended Kalman Filter, 简称EKF) 的土壤湿度同化方案的结果比较, 基于EnKF的土壤湿度同化方案易于实现, 且通过选择恰当的集合样本数其同化效果总体上略优于EKF同化方案, 但前者同化时需要花费较多的计算时间。     

Estimating the Soil Moisture Profile by Assimilating Near-Surface Observations with the Ensemble Kalman Filter （EnKF）

The paper investigates the ability to retrieve the true soil moisture profile by assimilating near-surface soil moisture into a soil moisture model with an ensemble Kalman filter （EnKF） assimilation scheme, including the effect of ensemble size, update interval and nonlinearities in the profile retrieval, the required time for full retrieval of the soil moisture profiles, and the possible influence of the depth of the soil moisture observation. These questions are addressed by a desktop study using synthetic data. The ＂true＂ soil moisture profiles are generated from the soil moisture model under the boundary condition of 0.5 cm d^-1 evaporation. To test the assimilation schemes, the model is initialized with a poor initial guess of the soil moisture profile, and different ensemble sizes are tested showing that an ensemble of 40 members is enough to represent the covariance of the model forecasts. Also compared are the results with those from the direct insertion assimilation scheme, showing that the EnKF is superior to the direct insertion assimilation scheme, for hourly observations, with retrieval of the soil moisture profile being achieved in 16 h as compared to 12 days or more. For daily observations, the true soil moisture profile is achieved in about 15 days with the EnKF, but it is impossible to approximate the true moisture within 18 days by using direct insertion. It is also found that observation depth does not have a significant effect on profile retrieval time for the EnKF. The nonlinearities have some negative influence on the optimal estimates of soil moisture profile but not very seriously.     

土壤-大气系统中蒸发过程的数值模拟

本文利用一维土壤-大气耦合数值模式,对土面蒸发过程进行了模拟研究。模式较好反映了蒸发过程三阶段的趋势特征,并试验了不同土壤性质、太阳辐射强度和地转风速对蒸发过程影响。模式还揭示了蒸发过程中,土面热量平衡各分量相互转换过程。     

2009/2010年冬季云南干旱的进一步研究——前期土壤湿度影响的数值模拟

采用新一代中尺度数值模式WRFv3.2版本,模拟研究了前期(秋季)土壤湿度异常对云南冬季降水的影响。数值模拟试验结果和一系列分析清楚表明,前期(秋季)土壤湿度的异常偏低,会导致云南地区冬季(12月1日~2月28日)降水的显著减少;前期土壤湿度减少一半,可以使云南冬季的降水量平均减少30%以上,小部分区域减少达50%以上,影响十分明显。大气环流及其主要参量模拟结果的对比分析清楚表明,持续的西偏北气流和干气团的控制以及云南地区大气散度场和垂直运动场等的异常是导致降水量减少的直接原因。对降水过程的分析表明,前期土壤湿度减少对降水过程的频次和发生时间的影响较小,但对各次过程的降水强度影响明显。这是前期土壤湿度减少所导致的包括区域性蒸发量和热通量等大气物理过程的改变决定的。本研究数值模拟结果与关于区域性土壤湿度异常影响机理的已有结论基本一致。     

内蒙古典型草原区近40年气候变化及其对土壤水分的影响

分析气候变化对草原区土壤水分的影响对了解草原退化原因、恢复草原生态环境有重要的指导意义。根据近40年气象资料和近20年的土壤水分观测资料,利用线性趋势等数理统计方法,分析了内蒙古典型草原区气候变化趋势和对土壤水分变化的影响,得出内蒙古典型草原区近40年气候变化趋势与全球气候变化规律相似;影响土壤湿度的气象因子主要是降水和蒸发,温度通过影响蒸发而间接影响土壤湿度,蒸降差是分析气候变化对土壤水分影响的直观指标。气候变暖导致蒸发加剧,在降水增加不明显的条件下,加速了土壤干旱化程度。     

Effect of surface layer thickness on simultaneous transport of heat and water in a bare soil and its implications for land surface schemes

Abstract

A physically‐based multi‐layer numerical model is developed to determine the coupled transport of heat and water in the soil and in the soil‐atmosphere boundary layer. Using inputs of standard weather data and initial soil conditions the model is capable of predicting the surface energy balance components as well as water content and temperature profiles in the soil. It is used to predict these variables for a bare silt loam soil under two tillage treatments, viz. culti‐packed and left loose after disc‐harrowing, and the predicted results are compared with measurements. Very good agreement between the model predictions and measured evaporation and heat fluxes and soil water and temperatures for a ten‐day period shows that the model is capable of simulating the coupled transport of soil heat and soil water and their transfer across the soil surface‐atmosphere interface adequately.

Model predictions were compared with those of CLASS (Canadian Land Surface Scheme). It is shown that CLASS, version 2.6, provides good estimates of evaporation and hence the latent heat flux density, Q_E, under wetter soil conditions, but overestimates Q_E at moderately wet soil conditions and underestimates it under dry soil conditions. Under dry to moderately wet soil conditions the calculation of evaporation from bare soil is very sensitive to the thickness of the top layer particularly as the thickness approaches 10 cm.     

Interaction between soil hydrology and boundary-layer development

A two-layer model of soil hydrology and thermodynamics is combined with a one-dimensional model of the planetary boundary layer to study various interactions between evolution of the boundary layer and soil moisture transport. Boundary-layer moistening through surface evaporation reduces the potential and actual surface evaporation as well as the boundary-layer growth. With more advanced stages of soil drying, the restricted surface evaporation allows greater sensible heat flux which enhances boundary-layer growth and entrainment drying.Special individual cases are studied where the wind speed is strong, solar radiation is reduced, transpiration is important, the soil is thin, or the soil is covered with organic debris.     

A method for determining thermophysical properties of organic material in aqueous solutions: Succinic acid

A method for determining evaporation rates and thermodynamic properties of aqueous solution droplets is introduced. The method combines evaporation rate measurements using modified TDMA technique with data evaluation using an accurate evaporation model. The first set of data has been collected and evaluated for succinic acid aqueous solution droplets.Evaporation rates of succinic acid solution droplets have been measured using a TDMA system at controlled relative humidity (65%) and temperature (298 K). A temperature-dependent expression for the saturation vapour pressure of pure liquid phase succinic acid at atmospheric temperatures has been derived by analysing the evaporation rate data with a numerical model. The obtained saturation vapour pressure of liquid phase succinic acid is ln(p) = 118.41 − 16204.8/T − 12.452ln(T). The vapour pressure is in unit of Pascal and the temperature in Kelvin. A linear expression for the enthalpy of vaporization for liquid state succinic acid is also presented.According to the results presented in the following, a literature expression for the vapour pressure of liquid phase succinic acid defined for temperatures higher than 461 K [Yaws, C.L., 2003. Yaws' Handbook of Thermodynamic and Physical Properties of Chemical Compounds, Knovel] can be extrapolated to atmospheric temperatures with very good accuracy. The results also suggest that at 298 K the mass accommodation coefficient of succinic acid is unity or very close to unity.     

中国地区土壤湿度记忆性及其与降水特征变化的关系

本文首先利用中国气象局国家气象信息中心提供的中国732个站点观测的土壤体积含水量,评估了CLM4.5（Community Land Model version 4.5）在CFSR（Climate Forecast System Reanalysis）近地面大气数据驱动下模拟的逐月土壤湿度（记为CLM4.5-CFSR）,然后基于CLM4.5-CFSR比较了皮尔逊相关法和自相关法计算得到的1980～2009年中国地区土壤湿度记忆性的区域及季节分布特征,量化了土壤湿度的记忆能力,研究了降水频率、降水强度和近地表气温分别对土壤湿度记忆性的影响。结果表明:CLM4.5-CFSR能较好地反映出大部分地区月时间尺度上土壤湿度的变化特征。两种方法描述的土壤湿度记忆性的空间分布特征相似,但季节特征不同。不同深度土壤湿度的记忆时长相差不大,在0.85～2.2个月不等,其中内蒙古东北部较大,新疆西南部较小。春季,较湿的土壤记忆性也较强。当降水频率较低时,其对蒸发速率较大的地区土壤湿度的记忆性影响很小,当降水强度较大时,它会迅速补充土壤散失的水分,破坏初始时刻土壤的干湿状态,引起其记忆性减弱。近地表气温变化主要通过影响土壤的蒸发过程减弱土壤湿度的记忆性。未来可利用气候模式开展数值敏感性试验对本文得到的结论进行机理研究,为进一步提高季节和季节内尺度的降水预报提供依据。     

Experiments using new initial soil moisture conditions and soil map in the Eta model over La Plata Basin

An effort towards a more accurate representation of soil moisture and its impact on the modeling of weather systems is presented. Sensitivity tests of precipitation to soil type and soil moisture changes are carried out using the atmospheric Eta model for the numerical simulation of the development of a mesoscale convective system over northern Argentina. Modified initial soil moisture conditions were obtained from a hydrological balance model developed and running operationally at INPE. A new soil map was elaborated using the available soil profile information from Brazil, Paraguay, Uruguay, and Argentina and depicts 18 different soil types. Results indicate that more accurate initial soil moisture conditions and incorporating a new soil map with hydraulic parameters, more representative of South American soils, improve daily total precipitation forecasts both in quantitative and spatial representations.     

Impacts of Soil Moisture on the Numerical Simulation of a Post-Landfall Storm

Surface heat and moisture fluxes are important to the evolution of a tropical storm after its landfall. Soil moisture is one of the essential components that influence surface heating and moisture fluxes. In this study, the impact of soil moisture on a pre-landfall numerical simulation of Tropical Storm Bill(2015), which had a much longer lifespan over land, is investigated by using the research version of the NCEP Hurricane Weather Research and Forecasting(HWRF) model. It is found that increased soil moisture with SLAB scheme before storm's landfall tends to produce a weaker storm after landfall and has negative impacts on storm track simulation. Further diagnoses with different land surface schemes and sensitivity experiments indicate that the increase in soil moisture inside the storm corresponds to a strengthened vertical mixing within the storm boundary layer, which is conducive to the decay of storm and has negative impacts on storm evolution. In addition, surface diabatic heating effects over the storm environment are also found to be an important positive contribution to the storm evolution over land, but their impacts are not so substantial as boundary layer vertical mixing inside the storm. The overall results highlight the importance and uncertainty of soil moisture in numerical model simulations of landfalling hurricanes and their further evolution over land.     

Validation and sensitivity of the global hydrologic budget in stand-alone simulations with the ISBA land-surface scheme

 Global soil moisture data of high quality and resolution are not available by direct observation, but are useful as boundary and initial conditions in comprehensive climate models. In the framework of the GSWP (Global Soil Wetness Project), the ISBA land-surface scheme of Météo-France has been forced with meteorological observations and analyses in order to study the feasibility of producing a global soil wetness climatology at a 1°×1° horizontal resolution. A control experiment has been performed from January 1987 to December 1988, using the ISLSCP Initiative I boundary conditions. The annual mean, the standard deviation and the normalised annual harmonic of the hydrologic fields have been computed from the 1987 monthly results. The global maps which are presented summarise the surface hydrologic budget and its annual cycle. The soil wetness index and snow cover distributions have been compared respectively to the results of the ECMWF reanalysis and to satellite and in situ observations. The simulated runoff has been validated against a river flow climatology, suggesting a possible underestimation over some large river basins. Besides the control run, other simulations have been performed in order to study the sensitivity of the hydrologic budget to changes in the surface parameters, the precipitation forcing and the runoff scheme. Such modifications have a significant impact on the partition of total precipitation into evaporation and runoff. The sensitivity of the results suggests that soil moisture remains one of the most difficult climatological parameters to model and that any computed soil wetness climatology must be considered with great caution. Received: 3 January 1997 / Accepted: 19 August 1987     

三类陆面模式模拟土壤湿度廓线的对比研究

针对不同陆面模式对土壤湿度方程求解方法以及对土壤分层结构的差异, 本文选取了三类陆面模式（CLM, CABLE和ECMWF陆面模式）作比较研究。为了避免不同模式参数化方案引起的上边界误差, 上边界采用固定蒸发率、 入渗率和固定表层土壤湿度三类边界条件。土壤分层采用101层（细网格）和11层（粗网格）两种, 并考虑土壤性质沿深度变化。结果表明： 当土壤性质均匀时, 求解的差别主要在第三类边界条件下CLM 求出的水分入渗速度比其它两种快; 改用粗网格后由于土壤深层厚度加大无法与细网格得出的土壤湿度廓线相重合。当土壤性质非均匀时, 模拟结果间差别加大, 只有ECMWF模式模拟的土壤湿度廓线是严格连续的。对于模式和上边界的不同组合, 粗、 细网格模拟结果间均方根偏差不一致。一般而言, CABLE模拟的偏差除第一类条件较小外, 其它都是最大的。第二三类边界条件引起的偏差较大, 第一类最小。上述结果提示我们, 在比较不同陆面模式以及用观测资料来检验模拟结果时应充分考虑土壤分层及土壤性质非均匀性的可能影响。     

Optimization of Upstream Profiles in Modelled Flow Over Complex Terrain

This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.