The Third GABLS Intercomparison Case for Evaluation Studies of Boundary-Layer Models. Part B: Results and Process Understanding

We describe and analyze the results of the third global energy and water cycle experiment atmospheric boundary layer Study intercomparison and evaluation study for single-column models. Each of the nineteen participating models was operated with its own physics package, including land-surface, radiation and turbulent mixing schemes, for a full diurnal cycle selected from the Cabauw observatory archive. By carefully prescribing the temporal evolution of the forcings on the vertical column, the models could be evaluated against observations. We focus on the gross features of the stable boundary layer (SBL), such as the onset of evening momentum decoupling, the 2-m minimum temperature, the evolution of the inertial oscillation and the morning transition. New process diagrams are introduced to interpret the variety of model results and the relative importance of processes in the SBL; the diagrams include the results of a number of sensitivity runs performed with one of the models. The models are characterized in terms of thermal coupling to the soil, longwave radiation and turbulent mixing. It is shown that differences in longwave radiation schemes among the models have only a small effect on the simulations; however, there are significant variations in downward radiation due to different boundary-layer profiles of temperature and humidity. The differences in modelled thermal coupling to the land surface are large and explain most of the variations in 2-m air temperature and longwave incoming radiation among models. Models with strong turbulent mixing overestimate the boundary-layer height, underestimate the wind speed at 200 m, and give a relatively large downward sensible heat flux. The result is that 2-m air temperature is relatively insensitive to turbulent mixing intensity. Evening transition times spread 1.5 h around the observed time of transition, with later transitions for models with coarse resolution. Time of onset in the morning transition spreads 2 h around the observed transition time. With this case, the morning transition appeared to be difficult to study, no relation could be found between the studied processes, and the variation in the time of the morning transition among the models.     

Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes,Boundary Conditions and Spin-Up

We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about \(2\) K) and substantially underestimates it at night (about \(4\) K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( \(\approx \) 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

Evaluation of the Diurnal Cycle in the Atmospheric Boundary Layer Over Land as Represented by a Variety of Single-Column Models: The Second GABLS Experiment

We present the main results from the second model intercomparison within the GEWEX (Global Energy and Water cycle EXperiment) Atmospheric Boundary Layer Study (GABLS). The target is to examine the diurnal cycle over land in today??s numerical weather prediction and climate models for operational and research purposes. The set-up of the case is based on observations taken during the Cooperative Atmosphere-Surface Exchange Study-1999 (CASES-99), which was held in Kansas, USA in the early autumn with a strong diurnal cycle with no clouds present. The models are forced with a constant geostrophic wind, prescribed surface temperature and large-scale divergence. Results from 30 different model simulations and one large-eddy simulation (LES) are analyzed and compared with observations. Even though the surface temperature is prescribed, the models give variable near-surface air temperatures. This, in turn, gives rise to differences in low-level stability affecting the turbulence and the turbulent heat fluxes. The increase in modelled upward sensible heat flux during the morning transition is typically too weak and the growth of the convective boundary layer before noon is too slow. This is related to weak modelled near-surface winds during the morning hours. The agreement between the models, the LES and observations is the best during the late afternoon. From this intercomparison study, we find that modelling the diurnal cycle is still a big challenge. For the convective part of the diurnal cycle, some of the first-order schemes perform somewhat better while the turbulent kinetic energy (TKE) schemes tend to be slightly better during nighttime conditions. Finer vertical resolution tends to improve results to some extent, but is certainly not the solution to all the deficiencies identified.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part II: Continuous Formulations of Mixing Rates for Mass-Flux Schemes

The conditional sampling of coherent structures in large-eddy simulations of the convective boundary layer (Couvreux et al. Boundary-layer Meteorol 134:441–458, 2010) is used to propose and evaluate formulations of fractional entrainment and detrainment rates for mass-flux schemes. The proposed formulations are physically-based and continuous from the surface to the top of clouds. Entrainment is related to the updraft vertical velocity divergence, while detrainment depends on the thermal vertical velocity, on buoyancy and on the moisture contrast between the mean plume and its environment. The proposed formulations are first directly evaluated in simulations of shallow clouds. They are then tested in single-column simulations with the thermal plume model, a mass-flux representation of boundary-layer thermals.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part II: Mesoscale Modelling

The importance of soil moisture inputs and improved model physics in the prediction of the daytime boundary-layer structure during the Southern Great Plains Hydrology Experiment 1997 (SGP97) is investigated using the non-hydrostatic fifth-generation Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) Mesoscale Model MM5. This is Part II of a two-part study examining the relationship of surface heterogeneity to observed boundary-layer structure. Part I focuses on observations and utilizes a simple model while Part II uses observations and MM5 modelling. Soil moisture inputs tested include a lookup table based on soil type and season, output from an offline land-surface model (LSM) forced by atmospheric observations, and high-resolution ( 800 m) airborne microwave remotely sensed data. Model physics improvements are investigated by comparing an LSM directly coupled with the MM5 to a simpler force-restore method at the surface. The scale of land surface heterogeneities is compared to the scale of their effects on boundary-layer structure.The use of more detailed soil moisture fields allowed the MM5 to better represent the large-scale (hundreds of km) and small-scale (tens of km) horizontal gradients in surface-layer weather and, to a lesser degree, the atmospheric boundary-layer (ABL) height, which was evaluated against observations measured by differential absorption lidar (DIAL). The benefits of coupling an LSM to the MM5 were not readily evident in this summertime case, with the model having particular difficulty simulating the timing of maximum surface fluxes while underestimating the depth of the mixed layer.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part III: Derivation of a Statistical Scheme for Cumulus Clouds

We present a statistical cloud scheme based on the subgrid-scale distribution of the saturation deficit. When analyzed in large-eddy simulations (LES) of a typical cloudy convective boundary layer, this distribution is shown to be bimodal and reasonably well-fitted by a bi-Gaussian distribution. Thanks to a tracer-based conditional sampling of coherent structures of the convective boundary layer in LES, we demonstrate that one mode corresponds to plumes of buoyant air arising from the surface, and the second to their environment, both within the cloud and sub-cloud layers. According to this analysis, we propose a cloud scheme based on a bi-Gaussian distribution of the saturation deficit, which can be easily coupled with any mass-flux scheme that discriminates buoyant plumes from their environment. For that, the standard deviations of the two Gaussian modes are parametrized starting from the top-hat distribution of the subgrid-scale thermodynamic variables given by the mass-flux scheme. Single-column model simulations of continental and maritime case studies show that this approach allows us to capture the vertical and temporal variations of the cloud cover and liquid water.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple Prognostic Model

The atmospheric boundary-layer (ABL) depth was observed by airborne lidar and balloon soundings during the Southern Great Plains 1997 field study (SGP97). This paper is Part I of a two-part case study examining the relationship of surface heterogeneity to observed ABL structure. Part I focuses on observations. During two days (12–13 July 1997) following rain, midday convective ABL depth varied by as much as 1.5 km across 400 km, even with moderate winds. Variability in ABL depth was driven primarily by the spatial variation in surface buoyancy flux as measured from short towers and aircraft within the SGP97 domain. Strong correlation was found between time-integrated buoyancy flux and airborne remotely sensed surface soil moisture for the two case-study days, but only a weak correlation was found between surface energy fluxes and vegetation greenness as measured by satellite. A simple prognostic one-dimensional ABL model was applied to test to what extent the soil moisture spatial heterogeneity explained the variation in north–south ABL depth across the SGP97 domain. The model was able to better predict mean ABL depth and variations on horizontal scales of approximately 100 km using observed soil moisture instead of constant soil moisture. Subsidence, advection, convergence/divergence and spatial variability of temperature inversion strength also contributed to ABL depth variations. In Part II, assimilation of high-resolution soil moisture into a three-dimensional mesoscale model (MM5) is discussed and shown to improve predictions of ABL structure. These results have implications for ABL models and the influence of soil moisture on mesoscale meteorology     

The energy balance experiment EBEX-2000. Part II: Intercomparison of eddy-covariance sensors and post-field data processing methods

The eddy-covariance method is the primary way of measuring turbulent fluxes directly. Many investigators have found that these flux measurements often do not satisfy a fundamental criterion—closure of the surface energy balance. This study investigates to what extent the eddy-covariance measurement technology can be made responsible for this deficiency, in particular the effects of instrumentation or of the post-field data processing. Therefore, current eddy-covariance sensors and several post-field data processing methods were compared. The differences in methodology resulted in deviations of 10% for the sensible heat flux and of 15% for the latent heat flux for an averaging time of 30 min. These disparities were mostly due to different sensor separation corrections and a linear detrending of the data. The impact of different instrumentation on the resulting heat flux estimates was significantly higher. Large deviations from the reference system of up to 50% were found for some sensor combinations. However, very good measurement quality was found for a CSAT3 sonic together with a KH20 krypton hygrometer and also for a UW sonic together with a KH20. If these systems are well calibrated and maintained, an accuracy of better than 5% can be achieved for 30-min values of sensible and latent heat flux measurements. The results from the sonic anemometers Gill Solent-HS, ATI-K, Metek USA-1, and R.M. Young 81000 showed more or less larger deviations from the reference system. The LI-COR LI-7500 open-path H₂O/CO₂ gas analyser in the test was one of the first serial numbers of this sensor type and had technical problems regarding direct solar radiation sensitivity and signal delay. These problems are known by the manufacturer and improvements of the sensor have since been made. The National Center for Atmospheric Research is supported by the National Science Foundation.     

Effect of Surface Heterogeneity on the Boundary-Layer Height: A Case Study at a Semi-Arid Forest

We investigate the effects of an isolated meso-\(\gamma \)-scale surface heterogeneity for roughness and albedo on the atmospheric boundary-layer (ABL) height, with a case study at a semi-arid forest surrounded by sparse shrubland (forest area: \(28~\text{ km }^2\), forest length in the main wind direction: 7 km). Doppler lidar and ceilometer measurements at this semi-arid forest show an increase in the ABL height over the forest compared with the shrubland on four out of eight days. The differences in the ABL height between shrubland and forest are explained for all days with a model that assumes a linear growth of the internal boundary layer of the forest through the convective ABL upwind of the forest followed by a square-root growth into the stable free atmosphere. For the environmental conditions that existed during our measurements, the increase in ABL height due to large sensible heat fluxes from the forest (\(600~\text {W~m}^{-2}\) in summer) is subdued by stable stratification in the free atmosphere above the ABL, or reduced by high wind speeds in the mixed layer.     

Intercomparison and Evaluation of Precipitation Forecasts for MAP Seasons 1995 and 1996

Summary  The precipitation forecasts of four operational numerical weather prediction models over the Alpine region are evaluated and intercompared for two periods of interest to the Mesoscale Alpine Programme (MAP). The new analysis of Alpine rainfall of Frei and Scha?r (1998) is used to validate the models. It is found that the models have a tendency to overestimate the total precipitation and the frequency of intense rain events over high orography. The skill scores show good consistency between models, except for the ability to forecast light rain or heavy rain events. The partition between convective and stratiform rainfall is rather variable between the models. Received March 15, 1999/Revised July 12, 1999     

Wind-Tunnel Simulation of the Wake of a Large Wind Turbine in a Stable Boundary Layer. Part 1: The Boundary-Layer Simulation

Measurements have been made in both a neutral and a stable boundary layer as part of an investigation of the wakes of wind turbines in an offshore environment, in the EnFlo stratified flow wind tunnel. The working section is long enough for the flow to have become very nearly invariant with streamwise distance. In order to be systematic, the flow profile generators of Irwin-type spires and surface roughness were the same for both neutral and stable conditions. Achieving the required profiles by adjusting the flow generators, even for neutral flow, is a highly iterative art, and the present results indicate that it will be no less iterative for a stable flow (as well as there being more conditions to meet), so this was not attempted in the present investigation. The stable-case flow conformed in most respects to Monin–Obukhov similarity in the surface layer. A linear temperature profile was applied at the working section inlet, resulting in a near-linear profile in the developed flow above the boundary layer and ‘strong’ imposed stability, while the condition at the surface was ‘weak’. Aerodynamic roughness length (mean velocity) was not affected by stability even though the roughness Reynolds number ${<}1$ , while the thermal roughness length was much smaller, as is to be expected. The neutral case was Reynolds-number independent, and by inference, the stable case was also Reynolds-number independent.     

A Case Study of the Near-Neutral Coastal Internal Boundary-Layer Growth: Aircraft Measurements Compared with Different Model Estimates

A simplified version of a model of the internalboundary layer evolution is compared with aircraftmeasurements across the southern part of the east coast ofSweden (Blekinge). The two flight occasions used hererepresent relatively high wind velocities and negligiblesurface heat flux. Both the vertical potential temperaturegradient and the surface temperaturedifference between land and sea are small. Modelpredictions of the internal boundary-layer (IBL)height agree well with data, for both short and largerdistances from the coastline. The contributions to theIBL growth from the mechanical turbulence and theZilitinkevich correction are discussed, as are thesensitivities of the growth to a number of physicalparameters. In near neutralconditions, the vertical potential temperaturegradient is found to be an important parameter for thegrowth of the IBL.A comparison with several simple model estimates forthe height of the IBL is presented. The majority ofthe models overestimate the IBL height. Some of themodels give reasonable results close to the coast, butgenerally the simple models perform less well atlarger distances.     

中国土壤湿度的分布与变化 I. 多种资料间的比较

土壤湿度是影响气候的重要因子之一, 但观测资料的欠缺制约着该领域研究工作的开展。本文汇总了目前国际上运用较为广泛的四套土壤湿度资料: ERA40和NCEP/NCAR再分析资料、全球土壤湿度计划资料(GSWP2)、以及NCAR最近完成的利用观测资料强迫“通用陆面模式”CLM所产生的土壤湿度资料。在此基础上, 利用中国区域观测的19年 (1981～1999年) 的土壤湿度和13年 (1992～2004年) 的土壤相对湿度资料, 对四套资料在中国区域的可靠性进行了分析和比较, 主要结论如下: 四套资料基本揭示出土壤湿度的空间分布, GSWP2揭示了四套资料最多的共性, 即东北、 华南湿, 华北、 西北干, 土壤湿度基本由西北向东北和东南呈梯度增加的特征; GSWP2较好地描述了土壤湿度的季节循环; ERA40土壤湿度的年际变化与观测相关最好; 观测资料和四套资料都表明前期降水会增加土壤湿度, 但土壤湿度异常对后期降水的影响则不显著; 气温与土壤湿度的关系较复杂, 不同的区域有不同的特征。     

Numerical Prediction of the Boundary-Layer Flow Over the Bolund Hill: Assessment of Turbulence Models and Advection Schemes

Boundary-Layer Meteorology - The WindStation software package is applied to simulate the wind field over the Bolund hill. The standard, ReNormalization Group (RNG), realizable, and limited-length...     

海洋环流模式的发展和应用Ⅰ.全球海洋环流模式

概述近10年来中国科学院大气物理研究所大气科学和地球流体动力学数值模拟国家重点实验室全球海洋环流模式的发展及其在全球海气耦合模式的发展和气候模拟方面的应用.重点是:一个30层、0.5°×0.5°的准全球海洋环流模式LICOM的建立及其模拟的热带太平洋海洋环流和印度尼西亚贯穿流;以20层海洋模式为海洋分量建立的全球海洋-大气-陆面系统耦合模式GOALS在气候变化模拟方面的应用,和以海洋模式L30T63为海洋分量建立的灵活的耦合环流模式FGCM-0在热带太平洋-印度洋海气相互作用及古海洋-古气候模拟方面的应用     

A Regional Ensemble Forecast System for Stratiform Precipitation Events in Northern China. Part I: A Case Study

A single-model,short-range,ensemble forecasting system (Institute of Atmospheric Physics,Regional Ensemble Forecast System,IAP REFS) with 15-km grid spacing,configured with multiple initial conditions,multiple lateral boundary conditions,and multiple physics parameterizations with 11 ensemble members,was developed using the Weather and Research Forecasting Model Advanced Research modeling system for prediction of stratiform precipitation events in northern China.This is the first part of a broader research project to develop a novel cloud-seeding operational system in a probabilistic framework.The ensemble perturbations were extracted from selected members of the National Center for Environmental Prediction Global Ensemble Forecasting System (NCEP GEFS) forecasts,and an inflation factor of two was applied to compensate for the lack of spread in the GEFS forecasts over the research region.Experiments on an actual stratiform precipitation case that occurred on 5-7 June 2009 in northern China were conducted to validate the ensemble system.The IAP REFS system had reasonably good performance in predicting the observed stratiform precipitation system.The perturbation inflation enlarged the ensemble spread and alleviated the underdispersion caused by parent forecasts.Centering the extracted perturbations on higher-resolution NCEP Global Forecast System forecasts resulted in less ensemble mean root-mean-square error and better accuracy in probabilistic quantitative precipitation forecasts (PQPF).However,the perturbation inflation and recentering had less effect on near-surface-level variables compared to the mid-level variables,and its influence on PQPF resolution was limited as well.     

海洋环流模式的发展和应用Ⅱ.近海和区域环流模式

回顾了近20年来中国科学院大气物理研究所大气科学和地球流体动力学数值模拟国家重点实验室近海和区域海洋环流模式发展和应用的历程.介绍了正压区域海流模式在中国近海海流数值模拟研究当中的应用以及一个24层、最高水平分辨率为0.25°×0.25°的三重嵌套的中国近海环流模式及其应用.在嵌套模式基础上,基于变分原理,建立了一个中国近海环流的资料同化分析系统.初步结果表明,资料同化提高了模式计算的真实性,所揭示的黑潮的强度和流幅、夏季黄海冷水团的位置等都更接近观测结果.     

Subfilter-Scale Fluxes over a Surface Roughness Transition. Part II: A priori Study of Large-Eddy Simulation Models

The ability of subfilter-scale (SFS) models to reproduce the statistical properties of SFS stresses and energy transfers over heterogeneous surface roughness is key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. In this study, several SFS models are evaluated a priori using experimental data acquired downwind of a rough-to-smooth transition in a wind tunnel. The SFS models studied include the eddy-viscosity, similarity, non-linear and a mixed model consisting of a combination of the eddy-viscosity and non-linear models. The dynamic eddy-viscosity model is also evaluated. The experimental data consist of vertical and horizontal planes of high-spatial-resolution velocity fields measured using particle image velocimetry. These velocity fields are spatially filtered and used to calculate SFS stresses and SFS transfer rates of resolved kinetic energy. Coefficients for each SFS model are calculated by matching the measured and modelled SFS energy transfer rates. For the eddy-viscosity model, the Smagorinsky coefficient is also evaluated using a dynamic procedure. The model coefficients are found to be scale dependent when the filter scales are larger than the vertical measurement height and fall into the production subrange of the turbulence where the flow scales are anisotropic. Near the surface, the Smagorinsky coefficient is also found to decrease with distance downwind from the transition, in response to the increase in mean shear as the flow adjusts to the smooth surface. In a priori tests, the ability of each model to reproduce statistical properties of the SFS stress is assessed. While the eddy-viscosity model has low spatial correlation with the measured stress, it predicts mean stresses with the same accuracy as the other models. However, the deficiency of the eddy-viscosity model is apparent in the underestimation of the standard deviation of the SFS stresses and the inability to predict transfers of kinetic energy from the subfilter scales to the resolved scales. Overall, the mixed model is found to have the best performance.