地表非均匀性对区域平均水分通量参数化的影响

次网格尺度地表非均匀性对于网格区平均通量具有重要影响。若将网格区视为均一地表 ,并不能真实描述地 气通量交换过程 ,且可造成很大误差。文中从理论上证明 ,区域平均水分通量的变化率可分解为两部分 :第一部分为区域水分通量的算术平均变化率 ;第二部分为非均匀性所引起的水分通量变化率扰动 ,它与区域内土壤水分空间分布的变差系数有关。数值试验表明 ,地表土壤水分的水平空间变差系数集中反映了区域内土壤水分分布的非均匀程度 ,不同土壤对同样的非均匀程度其敏感性是不同的。变差系数愈大 ,非均匀性愈强 ,在相同的土壤水分平均值下 ,不同土壤类型对地表非均匀程度的敏感性并不相同。例如沙土和粘土受非均匀性的影响就可相差数十倍。     

A modified ISBA surface scheme for modeling the hydrology of Athabasca River Basin with GCM-scale data

A soil–vegetation–atmosphere transfer model (SVAT), interactions between the soil–biosphere–atmosphere (ISBA) of Météo France, is modified and applied to the Athabasca River Basin (ARB) to model its water and energy fluxes. Two meteorological datasets are used: the archived forecasts from the Meteorological Survey of Canada’s Global Environmental Multiscale Model (GEM) and the European Centre for Mid-range Weather Forecasts global re-analysis (ERA-40), representing spatial scales typical of a weather forecasting model and a global circulation model (GCM), respectively. The original treatment of soil moisture and rainfall in ISBA (OISBA) is modified to statistically account for sub-grid heterogeneity of soil moisture and rainfall to produce new, highly non-linear formulations for surface and sub-surface runoff (MISBA). These new formulations can be readily applied to most existing SVATs. Stand alone mode simulations using the GEM data demonstrate that MISBA significantly improves streamflow predictions despite requiring two fewer parameters than OISBA. Simulations using the ERA-40 data show that it is possible to reproduce the annual variation in monthly, mean annual, and annual minimum flows at GCM scales without using downscaling techniques. Finally, simulations using a simple downscaling scheme show that the better performance of higher resolution datasets can be primarily attributed to improved representation of local variation of land cover, topography, and climate.     

Models of sub-grid variability in numerical simulations of solute transport in heterogeneous porous formations: three-dimensional flow and effect of pore-scale dispersion

Modeling transport of contaminants in the earths subsurface relies on numerical solutions over grids with blocks larger than Darcys scale. The hydraulic conductivity is homogenized over the grid blocks and the plumes spreading is reduced as a consequence of the wiped-out variability. To compensate for this loss Rubin et al. (1999) proposed to augment mixing by block-effective dispersion coefficients, and Rubin et al. (2003) showed, by means of two dimensional simulations, how this concept can be applied in practice. In this paper, we present new solutions of the block-effective dispersion tensor for an axisymmetric exponential covariance model. In addition, we discuss the influence of pore-scale dispersion in both two- and three-dimensional applications.     

海冰非均匀条件下大气环流模式和海洋环流模式耦合的实现

通过采用一种较好地考虑了海气通量交换对下垫面特征非线性依赖的通量计算方案,利用较成熟的逐日通量距平耦合方案实现了次网格尺度海冰非均匀条件下一个全球大气环流模式和一个全球海洋环流模式在高纬地区的耦合,并完成了一个50年的长期积分。模拟结果表明,通量计算方案中水道的作用是合理的,耦合模式能较好地模拟出北半球高纬海冰的主要地理分布特征,但夏季北冰洋内部靠近欧亚大陆部分边缘海区海冰密集度偏大。     

Large-eddy simulations of convection-driven dynamos using a dynamic scale-similarity model

Dynamo simulations require sub-grid scale (SGS) models for the momentum and heat flux, the Lorentz force, and the magnetic induction. Previous large eddy simulations (LES) using the scale similarity model have represented many aspects of the SGS motion. However, discrepancies are observed due to interchanging the order of filtering operation and spatial differentiation. In this study, we implement a correction term for this commutation error specifically for the scale-similarity model. Furthermore, we implement a dynamic scheme to evaluate time-dependent coefficients for the SGS models. We perform dynamo simulations in a rotating plane layer with different spatial resolutions, and compare results for the time dependence of the large-scale magnetic field. Simulations are performed at two different Rayleigh numbers, using constant values for the other dimensionless numbers (Ekman, Prandtl, and magnetic Prandtl numbers). Both cases show that the dynamic LES can accurately represent the large-scale magnetic field, whereas the dynamo failed in the direct simulations without the SGS terms at the same spatial resolutions. We conclude that the dynamic versions of the SGS and commutation error correction are essential for successful dynamos on coarser grids.     

Review of Development and Application of High Resolution Global Climate System Model

Climate models have been used as an important tool to quantitatively study climate variability and to predict or project climate change in the future. One of the most important pathways for development and improvement of climate system model is to increase the spatial resolution and improve the corresponding physical parameterization schemes, which is very important for understanding climate change and improving climate prediction skill. Based on a brief introduction of the importance of developing high-resolution global climate system model, a review of recent progresses in the development and application of high-resolution models was summarized. The paper also introduced the current situation and problems for the development and evaluation of high-resolution models and focused on the key scientific and technical bottlenecks which restrict the development of high-resolution models, including the development of dynamic framework of the high-resolution ocean and atmospheric models and massive high performance parallel computing, the improvement of the sub-grid physical parameterization scheme, and mesoscale air-sea interaction. Meanwhile, the scientific objects and experiments design of the international high resolution climate model intercomparison project (HiResMIP) of the coupled model intercomparison project phase 6 (CMIP6) was introduced. Finally, we prospect the future developments and evaluations of high-resolution climate models in China was proposed.     

Optimisation of an idealised ocean model,stochastic parameterisation of sub-grid eddies

An optimisation scheme is developed to accurately represent the sub-grid scale forcing of a high dimensional chaotic ocean system. Using a simple parameterisation scheme, the velocity components of a 30 km resolution shallow water ocean model are optimised to have the same climatological mean and variance as that of a less viscous 7.5 km resolution model. The 5 day lag-covariance is also optimised, leading to a more accurate estimate of the high resolution response to forcing using the low resolution model.The system considered is an idealised barotropic double gyre that is chaotic at both resolutions. Using the optimisation scheme, we find and apply the constant in time, but spatially varying, forcing term that is equal to the time integrated forcing of the sub-grid scale eddies. A linear stochastic term, independent of the large-scale flow, with no spatial correlation but a spatially varying amplitude and time scale is used to represent the transient eddies. The climatological mean, variance and 5 day lag-covariance of the velocity from a single high resolution integration is used to provide an optimisation target. No other high resolution statistics are required. Additional programming effort, for example to build a tangent linear or adjoint model, is not required either.The focus of this paper is on the optimisation scheme and the accuracy of the optimised flow. However the forcing can provide insights in the design of deterministic and stochastic parameterisations. In the present study, we found that the stochastic parameterisation correcting the model variance is associated with the spatial pattern of eddy-decorrelation timescales rather than the spatial pattern of the amplitude of the variance. The method can be applied in future investigations into the physical processes that govern barotropic turbulence and it can perhaps be applied to help understand and correct biases in the mean and variance of a more realistic coarse or eddy-permitting ocean model. The method is complementary to current parameterisations and can be applied at the same time without modification.     

非均匀地表陆面过程参数化研究

地表固有的非均匀性影响近地层大气的垂直结构,甚至改变局地天气条件,亦使得大气数值模式大尺度网格面积水热通量的计算对其具有较强的敏感性。为了提高气候模式性能,非均匀陆面过程参数化已是当前大气边界层和陆面过程模式研究的热点和难点问题之一。本文在调研国内外大量文献的基础上,综述了近年来非均匀地表陆面过程参数化的研究现状,分析和比较了不同的参数化方法的优缺点以及数值模式模拟结果对它们的响应,提出了目前尚待继续探讨和解决的几个关键性问题。     

Mosaic方法在非均匀下垫面上的适用性研究

在非均匀下垫面情况下, Mosaic方法是目前国际上广泛运用于模式中计算地表通量的方法.大量的研究表明, 下垫面的非均匀分布会引发局地环流, 非均匀分布的空间尺度较大时, 所引起的环流甚至可以达到海陆风的强度.这种环流的存在直接影响到次网格地表通量的计算.次网格地表非均匀分布, 尤其是大尺度模式中的次网格非均匀分布, 必将影响地表通量的计算.本文针对次网格地表非均匀问题, 设计了高分辨率的Mosaic试验和非均匀试验, 开展了不同背景风情况下的一系列数值试验, 以探讨这种影响的程度.结果表明, 在土壤湿度空间分布不均匀的情况下, 运用Mosaic方法计算得到的地表潜热通量偏小, 背景风较小的时候偏差较大, 背景风增强时偏差减小.     

青藏高原地形重力波拖曳的初步分析及数值模拟研究

针对目前对青藏高原大地形激发的重力波拖曳相关问题还不十分清楚,在GRAPES_Meso模式中引入次网格地形重力波拖曳参数化方案,通过数值试验初步研究了青藏高原地区次网格地形重力波拖曳的一些相关参数,结果指出:(1)沿30°N地形重力波拖曳的垂直分布显示,阻塞拖曳主要存在于模式的低层(第1—5层),重力波拖曳主要存在于模式的第5—10层;从水平分布看,模式第3层以阻塞拖曳为主,主要位于青藏高原边缘地区,阻塞拖曳大值区沿喜马拉雅山脉走向和青藏高原东坡;模式第5层以重力波拖曳为主,主要位于青藏高原东部地区和云贵高原的北部边缘。(2)弗劳德数和气流绕流高度分析表明,在青藏高原喜马拉雅山脉一带和高原东部边缘地区,气流爬坡能力强,同时在这一地区绕流高度最高;弗劳德数越大的地区绕流高度距离地表越高。(3)采用次网格地形重力波拖曳参数化方案后,对于低层和高层地形重力波破碎的发生有更准确的描述,地形重力波是向上垂直传播的。(4)个例和批量试验检验结果表明,采用次网格地形重力波拖曳参数化方案对于风场和降水模拟有正效果,提高了模式预报的准确率。