Establishment and analysis of a High-Resolution Assimilation Dataset of the water-energy cycle in China

For better prediction and understanding of land-atmospheric interaction, in-situ observed meteorological data acquired from the China Meteorological Administration (CMA) were assimilated in the Weather Research and Forecasting (WRF) model and the monthly Green Vegetation Coverage (GVF) data, which was calculated using the Normalized Difference Vegetation Index (NDVI) of the Earth Observing System Moderate-Resolution Imaging Spectroradiometer (EOS-MODIS) and Digital Elevation Model (DEM) data of the Shuttle Radar Topography Mission (SRTM) system. Furthermore, the WRF model produced a High-Resolution Assimilation Dataset of the water-energy cycle in China (HRADC). This dataset has a horizontal resolution of 25 km for near surface meteorological data, such as air temperature, humidity, wind vectors and pressure (19 levels); soil temperature and moisture (four levels); surface temperature; downward/upward short/long radiation; 3-h latent heat flux; sensible heat flux; and ground heat flux. In this study, we 1) briefly introduce the cycling 3D-Var assimilation method and 2) compare results of meteorological elements, such as 2 m temperature and precipitation generated by the HRADC with the gridded observation data from CMA, and surface temperature and specific humidity with Global Land Data Assimilation System (GLDAS) output data from the National Aeronautics and Space Administration (NASA). We find that the simulated results of monthly 2 m temperature from HRADC is improved compared with the control simulation and has effectively reproduced the observed patterns. The simulated special distribution of ground surface temperature and specific humidity from HRADC are much closer to GLDAS outputs. The spatial distribution of root mean square errors (RMSE) and bias of 2 m temperature between observations and HRADC is reduced compared with the bias between observations and the control run. The monthly spatial distribution of surface temperature and specific humidity from HRADC is consistent with the GLDAS outputs over China. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations, and the simulated results could be used in further research on the long-term climatic effects and characteristics of the water-energy cycle over China.     

Diurnal cycle of surface energy fluxes in high mountain terrain: High-resolution fully coupled atmosphere-hydrology modelling and impact of lateral flow

Water and energy fluxes are inextricably interlinked within the interface of the land surface and the atmosphere. In the regional earth system models, the lower boundary parameterization of land surface neglects lateral hydrological processes, which may inadequately depict the surface water and energy fluxes variations, thus affecting the simulated atmospheric system through land-atmosphere feedbacks. Therefore, the main objective of this study is to evaluate the hydrologically enhanced regional climate modelling in order to represent the diurnal cycle of surface energy fluxes in high spatial and temporal resolution. In this study, the Weather Research and Forecasting model (WRF) and coupled WRF Hydrological modelling system (WRF-Hydro) are applied in a high alpine catchment in Northeastern Tibetan Plateau, the headwater area of the Heihe River. By evaluating and intercomparing model results by both models, the role of lateral flow processes on the surface energy fluxes dynamics is investigated. The model evaluations suggest that both WRF and coupled WRF-Hydro reasonably represent the diurnal variations of the near-surface meteorological fields, surface energy fluxes and hourly partitioning of available energy. By incorporating additional lateral flow processes, the coupled WRF-Hydro simulates higher surface soil moisture over the mountainous area, resulting in increased latent heat flux and decreased sensible heat flux of around 20–50 W/m² in their diurnal peak values during summertime, although the net radiation and ground heat fluxes remain almost unchanged. The simulation results show that the diurnal cycle of surface energy fluxes follows the local terrain and vegetation features. This highlights the importance of consideration of lateral flow processes over areas with heterogeneous terrain and land surfaces.     

Numerical simulations of land surface physical processes and land-atmosphere interactions over oasis-desert/Gobi region

A land-surface physical process model was coupled with a mesoscale atmospheric model. This cou- pled model was then used to simulate the interactions between land and the atmosphere, including surface temperature, net radiation, sensible heat flux and latent heat flux over a desert/Gobi with an oasis in northwestern semiarid regions in China. Comparisons between observations and simulations were made over the oasis and the desert/Gobi, respectively. Both cold island effect and wet island ef- fect, the so-called oasis effect, were observed and simulated. Lower temperature, higher specific hu- midity and weaker turbulent transfer were present over the oasis than the desert/Gobi. A subsidence occurred over the oasis, leading to a thermally-generated mesoscale circulation.     

南极中山站至Dome A考察断面湍流参数特征近地层湍流参数特征

利用南极中山站至Dome A考察断面上3个自动气象站2005～2007年的观测资料和2008年夏季在中山站附近冰盖获取的湍流观测资料,应用空气动力学方法和涡动相关法计算分析了中山站至Dome A断面上近地层各种湍流参数(感热通量,潜热通量,湍流温度、湿度和速度尺度,地表粗糙度,大气稳定度及动量输送系数)的季节变化、日变...     

CLM4-LISSS模型对太湖多时间尺度水热通量模拟性能的评估

湖泊模型为数值天气预报模型提供热量通量、水汽通量和动量通量等下边界条件,但是不同时间尺度上湖泊水热通量变化的控制因子不同,因此有必要对湖泊模型进行多时间尺度上的离线评估.本文利用2012-2016年太湖中尺度通量网避风港站的气象资料和辐射数据驱动CLM4-LISSS模型(Community Land Model version 4-Lake,Ice,Snow and Sediment Simulator),并与涡度相关观测(Eddy Covariance,EC)结果进行对比,以年平均潜热通量模拟结果最佳为目标调整了模式中的消光系数、粗糙度长度方案,研究了该模型从半小时到年尺度上对湖表温度和水热通量的模拟性能.结果表明:模型对湖表温度的模拟在各时间尺度上均比较理想,但是模拟的日较差较小;从半小时到年尺度上潜热通量的变化趋势都能被很好地模拟出来,但在季节尺度上,潜热通量的模拟出现了秋冬季偏高、春夏季偏低的情况,季节变化模拟不准确.湖表温度和潜热通量模拟偏差的原因可能是消光系数的参数化方案.相比之下,感热通量尽管年际变化趋势的模拟值与观测值一致,但是从半小时到年尺度均被高估.特别地,冷锋过境期间,模型能较好地模拟出潜热通量和感热通量的变化趋势,但对于高风速条件下的感热通量模拟效果不佳.本文的研究结果能为湖泊模式的应用与发展提供有用信息.     

Relative importance of advective heat transport and boundary layer decoupling in the melt dynamics of a patchy snow cover

For a patchy snow cover the advective heat transport and the near-surface boundary layer decoupling, which have an opposite effect on sensible heat transport onto the snow surface, are both expected to increase in magnitude. The main aim of this study is to investigate the effects of locally developing atmospheric stratification over a discontinuous snow cover which can result in a decoupling from the warm atmosphere. We are particularly interested in the effect of boundary layer decoupling on the net sensible heat flux into the snow. We therefore applied local eddy flux measurements over snow patches at three different heights above the snow surface. We identified wind velocity, turbulence intensity, fetch distance and topographical curvature as the main factors driving the boundary layer depth and the efficiency of advective heat transport to contribute to snow ablation. The atmospheric decoupling is thus shown to be a key mechanism in snow patch survival.     

Simulation and evaluation of 2-m temperature over Antarctica in polar regional climate model

The European Centre for Medium-Range Weather Forecasts Reanalysis ERA40,National Centers for Environmental Prediction(NCEP) 20th-century reanalysis,and three station observations along an Antarctic traverse from Zhongshan to Dome-A stations are used to assess 2-m temperature simulation skill of a regional climate model.This model(HIRHAM) is from the Alfred Wegener Institute for Polar and Marine Research in Germany.Results show:(1) The simulated multiyear averaged 2-m temperature field pattern is close to that of ERA40 and NCEP;(2) the cold bias relative to ERA40 over all of Antarctic regions is 1.8℃,and that to NCEP reaches 5.1℃;(3) bias of HIRHAM relative to ERA40 has seasonal variation,with a cold bias mainly in the summer,as much as 3.4℃.There is a small inland warm bias in autumn of 0.3℃.Further analysis reveals that the reason for the cold bias of 2-m temperature is that physical conditions of the near-surface boundary layer simulated by HIRHAM are different from observations:(1) During the summer,observations show that near-surface atmospheric stability conditions have both inversions and non-inversions,which is due to the existence of both positive and negative sensible heat fluxes,but HIRHAM almost always simulates a situation of inversion and negative sensible heat flux;(2) during autumn and winter,observed near-surface stability is almost always that of inversions,consistent with HIRHAM simulations.This partially explains the small bias during autumn and winter.     

鄱阳湖夏季水热通量特征及环境要素影响分析

气候变化加速了全球水文循环过程,然而,气候变化如何影响水体蒸发及其水热通量交换仍然不清楚.基于涡度相关系统观测鄱阳湖水体水热通量过程,在小时和日尺度分析了水热通量的变化规律及其主要影响因子.研究表明,潜热通量日变化波动剧烈,大部分为正值,变化范围在-50~580 W/m2之间.而感热通量数值较小,变化范围在-50~50 W/m2之间.8月份潜热通量和感热通量均呈波动下降趋势,均值分别为167.4和15.9 W/m2.8月份日平均潜热通量和感热通量之和大于净辐射,这是由于这一时段储存在水体中的热量释放并补充潜热通量和感热通量.小时尺度上潜热通量日变化在相位上与净辐射无显著相关性,而与风速显著相关.在日尺度变化趋势上,8月份日平均潜热通量仍主要受到风速和水温的影响,感热通量则主要受到风速和饱和水汽压差的影响.     

Spatial distribution of surface energy fluxes over the Loess Plateau in China and its relationship with climate and the environment

China's Loess Plateau is located at the edge of the Asian summer monsoon in a transition zone of climate and ecology. In the Loess Plateau, climate and environments change along with space, which has an obvious impact on the spatial distribution of surface energy fluxes. Because of scarce land-surface observation sites and short observation time in this area, previous studies have failed to fully understand the land-surface energy balance characteristics over the entire the Loess Plateau and their effect mechanisms. In this paper, we first test the simulation ability of the Community Land Model(CLM) model by comparing its simulated data with observed data. Based on the simulation data for the Loess Plateau over the past thirty years, we then analyze the spatial distribution of surface energy fluxes and compare the pattern differences between the area averages for the driest year and wettest year. Furthermore, we analyze the relationship between the spatial distribution of the components of the surface energy balance with longitude, latitude, altitude, precipitation and temperature. The main results are as follows: the spatial distribution of surface energy fluxes are significantly different, with the surface net radiation and sensible heat flux increasing from south to north and latent heat flux and soil heat flux decreasing from southeast to northwest. The sensible heat flux at the driest point is nearly twice as high as that at the wettest point, whereas the latent heat flux and soil heat flux at the driest point are half as much as that at the wettest point. The impact of variations of annual precipitation on the components of the surface energy balance is also obvious, and the maximum magnitude of the changes to the sensible heat flux and latent heat flux is nearly 30%. To a certain extent, geographical factors(including longitude, latitude, and altitude) and climate factors(including temperature and precipitation) affect the surface energy fluxes. However, the surface net radiation is more closely related to latitude and altitude, sensible heat flux is more closely related to the monsoon rainfall and latitude, and latent heat flux and soil heat flux are more closely related to the monsoon rainfall.     

Improving WRF-Hydro runoff simulations of heavy floods through the sea surface temperature fields with higher spatio-temporal resolution

This study investigates the impact of the spatio-temporal accuracy of four different sea surface temperature (SST) datasets on the accuracy of the Weather Research and Forecasting (WRF)-Hydro system to simulate hydrological response during two catastrophic flood events over the Eastern Black Sea (EBS) and the Mediterranean (MED) regions of Turkey. Three time-variant and high spatial resolution external SST products (GHRSST, Medspiration and NCEP-SST) and one coarse-resolution and time-invariant SST product (ERA5- and GFS-SST for EBS and MED regions, respectively) already embedded in the initial and the boundary conditions datasets of WRF model are used in deriving near-surface atmospheric variables through WRF. After the proper event-based calibration is performed to the WRF-Hydro system using hourly and daily streamflow data in both regions, uncoupled model simulations for independent SST events are conducted to assess the impact of SST-triggered precipitation on simulated extreme runoff. Some localized and temporal differences in the occurrence of the flood events with respect to observations depending on the SST representation are noticeable. SST products represented with higher cross-correlations (GHRSST and Medspiration) revealed significant improvement in flood hydrographs for both regions. The GHRSST dataset shows a substantial improvement in NSE (~70%), RMSE reduction up to 20%, and an increase in correlation from 0.3 to 0.8 with respect to the invariable SST (ERA5) in simulated runoffs over the EBS region. The use of both GHRSST and Medspiration SST data characterized with high spatio-temporal correlation resulted in runoff simulations exactly matching the observed runoff peak of 300 m³/s by reducing the overestimation seen in invariable SST (GFS) in the MED region. Improved precipitation simulation skills of the WRF model with the detailed SST representation show that the hydrographs of GHRSST and Medspiration simulations show better performance compared to the simulated hydrographs by observed precipitation.     

Determination of regional land surface heat fluxes over a heterogeneous landscape of the Jiddah area of Saudi Arabia by using Landsat‐7 ETM data

In this study, a parameterization method based on Landsat‐7 Enhanced Thematic Mapper (ETM) data and field observations is presented and tested for deriving the regional land surface variables, vegetation variables and land surface heat fluxes over a heterogeneous landscape. As a case study, the method and two Landsat‐7 ETM images are applied to the Jiddah area of Saudi Arabia. The regional distribution maps of surface reflectance, normalized difference vegetation index, modified soil adjusted vegetation index (MSAVI), vegetation coverage, leaf area index, surface temperature, net radiation flux, soil heat flux, sensible heat flux and latent heat flux have been determined over the Jiddah area. The derived results have been validated by using the ‘ground truth’. The results show that the more reasonable regional distributions of land surface variables (surface reflectance, surface temperature), vegetation variables (MSAVI and vegetation coverage), net radiation, soil heat flux and sensible heat flux can be obtained by using the method proposed in this study. Further improvement of the method is also discussed. Copyright © 2006 John Wiley & Sons, Ltd.     

The effects of small water bodies on the atmospheric heat and water budgets over the MacKenzie River Basin

The effects of small water bodies or lakes on the surface sensible and latent heat fluxes and the transport of heat and water vapour in the atmospheric boundary layer (ABL) over the Mackenzie River Basin (MRB) are studied from two cases, which occurred on 2 and 8 June 1999 during the warm season. The synoptic condition for the cases is representative of about 33% of the synoptic situation over the MRB. The two events are simulated using the Canadian mesoscale compressible community (MC2) model. A one‐way nesting grid approach is employed with the highest resolution of 100 m over a domain of 100 km². Experiments were conducted with (LAKE) and without (NOLAKE) the presence of small water bodies, whose size distribution is obtained through an inversion algorithm using information of their linear dimension determined from aircraft measurement of surface temperature during MAGS (the Mackenzie GEWEX (Global Energy and Water Cycle Experiment) Study) in 1999. The water bodies are assumed to be distributed randomly in space with a fractional area coverage of 10% over the MRB. The results show that, in the presence of lakes, the domain‐averaged surface sensible heat flux on 2 June 1999 (8 June 1999) decreases by 9·3% (6·6%). The surface latent heat flux is enhanced by 18·2% (81·5%). Low‐level temperature advection and the lake surface temperature affect the air–land/lake temperature contrast, which in turn controls the sensible heat flux. In the absence of lakes the surface wind speed impacts the latent heat flux, but in the presence of lakes the moisture availability and the atmospheric surface layer stability control the latent heat flux. The enhancement is smaller on 2 June 1999 as a result of a stable surface layer caused by the presence of colder lake temperatures. The domain‐averaged apparent heat source and moisture sink due to turbulent transports were also computed. The results show that, when lakes are present, heating and drying occur in the lowest 100 m from the surface. Above 100 m and within the ABL, there was apparent cooling. However, the apparent moistening profiles reveal that lakes tend to moisten the ABL through transfer of moisture from the lowest 50–100 m layer. Copyright © 2004 John Wiley & Sons, Ltd.     

海面动量、感热和水汽输送的动力学分析

基于大气边界层动量、感热和水汽通量的基本方程,定量地计算了波动海面的动量、感热和水汽通量。首先,应用Prandtl的混合长概念,推导出贴海面大气层中风速、位温和比湿的涡动交换率及其贴海面层厚度,并且证明了波面上位温或比湿贴海面层厚度与速度贴海面层厚度的比值,和平面上的相应比值完全相等。随后,利用空气动力学理论讨论了贴海面动量、感热和水汽输送的参数化问题。最后,对现有五种理论模式进行了比较说明。     

YSU边界层参数化方案中不同地形订正方法对地面风速及温度模拟的影响

为了修正中尺度气象模式WRF(Weather Research and Forecasting model)对低层风速模拟的系统性误差,有学者在新版本WRF模式的YSU(Yonsei University)边界层参数化方案中加入了两个地形订正方法:Jiménez方法和UW(University of Washington method)方法.本文利用这两个地形订正方法,选取了两个时间段,对北京地区的地面气象要素以及气象要素垂直廓线进行了个例模拟研究,模拟结果和观测数据的比对表明在北京地区:是否采用地形订正,对地面温度的模拟几乎没有影响;采用地形订正后,模式对地面风速的模拟有明显的改进,两种方法对风速模拟的差别主要体现在山/丘陵地区;Jiménez方法在山/丘陵地区的模拟风速明显偏大,而采用UW方法进行订正后,模拟的风速减小,更接近观测值;两种方法在山谷地区对风速均有一定的过度订正.通过分析气象要素的垂直廓线发现,不同地形订正方法主要影响的是2000m以下的低层风速.总体而言,UW地形订正方法在北京地区更为适用,采用UW方法后,模拟得到的地面气象要素的各项统计参数基本达到了统计基准值.     

TOGA-COARE IOP期间的海气通量观测结果

在TOGA－COAREIOP期间用涡度相关法对海气热通量进行了船载直接观测．对垂直风速、温度和湿度湍流脉动观测数据的谱分析显示它们在高频区基本满足"－2／3次方律"．对船体简谐震荡影响的讨论从理论上证明该影响在热通量计算中可被忽略．根据以此方法得到的通量求出了中性层结条件下感热和潜热的整体输送系数分别为2．25×10^－3和1．26×10^－3．对海气边界层特性的分析表明该海域的近海层主要呈中性或弱不稳定层结．海气通量的变化与背景环流形势密切相关,潜热通量主要受海面风场强度的影响,而感热通量变化除了风场的影响外,层结变化也是一个重要因素．用整体输送法计算TOGA－COAREIOP期间以及TOGA期间8个航次的通量结果而得到的Bowen比约为0．1,显示潜热通量是暖池大气的主要热源．     

Assimilation of Doppler Weather Radar Data in WRF Model for Simulation of Tropical Cyclone Aila

For the accurate and effective forecasting of a cyclone, it is critical to have accurate initial structure of the cyclone in numerical models. In this study, Kolkata Doppler weather radar (DWR) data were assimilated for the numerical simulation of a land-falling Tropical Cyclone Aila (2009) in the Bay of Bengal. To study the impact of radar data on very short-range forecasting of a cyclone's path, intensity and precipitation, both reflectivity and radial velocity were assimilated into the weather research and forecasting (WRF) model through the ARPS data assimilation system (ADAS) and cloud analysis procedure. Numerical experiment results indicated that radar data assimilation significantly improved the simulated structure of Cyclone Aila. Strong influences on hydrometeor structures of the initial vortex and precipitation pattern were observed when radar reflectivity data was assimilated, but a relatively small impact was observed on the wind fields at all height levels. The assimilation of radar wind data significantly improved the prediction of divergence/convergence conditions over the cyclone's inner-core area, as well as its wind field in the low-to-middle troposphere (600–900 hPa), but relatively less impact was observed on analyzed moisture field. Maximum surface wind speed produced from DWR–Vr and DWR–ZVr data assimilation experiments were very close to real-time values. The impact of radar data, after final analysis, on minimum sea level pressure was relatively less because the ADAS system does not adjust for pressure due to the lack of pressure observations, and from not using a 3DVAR balance condition that includes pressure. The greatest impact of radar data on forecasting was realized when both reflectivity and wind data (DWR–ZVr and DWR–ZVr00 experiment) were assimilated. It is concluded that after final analysis, the center of the cyclone was relocated very close to the observed position, and simulated cyclone maintained its intensity for a longer duration. Using this analysis, different stages of the cyclone are better captured, and cyclone structure, intensification, direction of movement, speed and location are significantly improved when both radar reflectivity and wind data are assimilated. As compared to other experiments, the maximum reduction in track error was noticed in the DWR–ZVr and DWR–ZVr00 experiments, and the predicted track in these experiments was very close to the observed track. In the DWR–ZVr and DWR–ZVr00 experiments, rainfall pattern and amount of rainfall forecasts were remarkably improved and were similar to the observation over West Bengal, Orissa and Jharkhand; however, the rainfall over Meghalaya and Bangladesh was missed in all the experiments. The influence of radar data reduces beyond a 12-h forecast, due to the dominance of the flow from large-scale, global forecast system models. This study also demonstrates successful coupling of the data assimilation package ADAS with the WRF model for Indian DWR data.     

2002年南海季风爆发前后西沙海区海-气通量交换及其变化

利用2002年4月24日至6月20日在西沙海区进行的第三次南海海-气通量观测试验资料,采用涡相关法和TOGA COARE25b版本通量计算方案,计算了西南季风爆发前后海洋-大气间的通量交换,讨论了辐射、动量、感热通量、潜热通量、海洋热量净收支的时间变化特征及其与气象要素变化的关系.结果表明：西南季风爆发前后,太阳短波辐射、海面净辐射、潜热通量和海洋热量净收支变化特别强烈;通量变化受不同环境要素的影响：感热通量与海-气温差呈正相关关系,与气温呈明显的负相关关系.潜热通量与风速、海-气温差及海面水温均有正相关关系,其中与风速的关系最密切.动量通量（τ）主要随风速变化,它与风速（V）的关系可以表示为τ=000185V2-000559V+001248.     

降水条件下的典型干旱区陆面特征模拟验证

利用“我国西北干旱区陆－气相互作用观测试验"在敦煌双墩子戈壁站取得的观测资料及最近的一些研究成果对陆面模式中反照率、相似性函数及地表粗糙度(动量粗糙度，标量粗糙度)3个方面的参数化方案进行了改进，然后对一次典型降水过程的陆面特征及近地层的风、温、湿进行了模拟.结果表明：改进的模式能对降水条件下的干旱区陆面特征进行较好的模拟，其中对辐射、地表温度的模拟相当好，而对能量的模拟虽然还有要改进的地方，但总体结果令人满意；另外改进后的模式对近地层的温度和湿度的模拟也有明显的改善.     

On the performance of surface renewal analysis to estimate sensible heat flux over two growing rice fields under the influence of regional advection

High-frequency temperature data were recorded at one height and they were used in Surface Renewal (SR) analysis to estimate sensible heat flux during the full growing season of two rice fields located north–northeast of Colusa, CA (in the Sacramento Valley). One of the fields was seeded into a flooded paddy and the other was drill seeded before flooding. To minimize fetch requirements, the measurement height was selected to be close to the maximum expected canopy height. The roughness sub-layer depth was estimated to discriminate if the temperature data came from the inertial or roughness sub-layer. The equation to estimate the roughness sub-layer depth was derived by combining simple mixing-length theory, mixing-layer analogy, equations to account for stable atmospheric surface layer conditions, and semi-empirical canopy–architecture relationships. The potential for SR analysis as a method that operates in the full surface boundary layer was tested using data collected over growing vegetation at a site influenced by regional advection of sensible heat flux. The inputs used to estimate the sensible heat fluxes included air temperature sampled at 10 Hz, the mean and variance of the horizontal wind speed, the canopy height, and the plant area index for a given intermediate height of the canopy. Regardless of the stability conditions and measurement height above the canopy, sensible heat flux estimates using SR analysis gave results that were similar to those measured with the eddy covariance method. Under unstable cases, it was shown that the performance was sensitive to estimation of the roughness sub-layer depth. However, an expression was provided to select the crucial scale required for its estimation.     

巢湖对冬季陆面辐射和热量过程的影响

从物理气候学的观点出发,分析了湖泊对水体和陆地辐射平衡和热量平衡各分量的影响。以巢湖地区冬季观测资料为例,揭示了在晴稳天气湖陆风环流对陆面显热输送的影响,以及在冷平流天气,湖泊对上下风方陆面显热输送的影响。