A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

Boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice

During the field experiment ARKTIS 1993 ten cases of boundary-layer modification in wintertime cold-air outbreaks from the Arctic sea ice in the Spitsbergen region were observed by aircraft over a distance ranging from about 50 km over the ice to about 300 km over the water. The modification depends decisively on the initial conditions over the ice, the boundary conditions at the bottom and top of the boundary layer and on the conditions of the large-scale flow. The modification of the bulk boundary-layer characteristics in relation to these conditions is presented.Besides the air-sea temperature contrast, the most important role for the boundary-layer modification is played by the stability on top of the boundary layer and by the divergence of the large-scale flow. According to the high variability of these conditions the observed boundary-layer modifications were very variable ranging from 100 to 300 m thick boundary layers with air temperatures between -32 and -22 °C over the ice to thicknesses between 900 and 2200 m and air temperatures between -15 and -5 °C after 300 km fetch over the open water. In most cases the large-scale flow was anticyclonic and divergent over the ice and changed to cyclonic and convergent over the water and an ice-sea breeze was superimposed on it.The sensible and latent heat fluxes are the dominant terms in the surface energy budget over the open water and ranged between 200 and 700 W m^-2 whereas the net longwave radiation is the dominating term over the ice with the heat fluxes only about 10 W m^-2.     

Wind-Tunnel Study Of Atmospheric Stable Boundary Layers Over A Rough Surface

Wind-tunnel simulations of theatmospheric stable boundary layer (SBL) developedover a rough surface were conducted by using athermally stratified wind tunnel at the Research Institutefor Applied Mechanics (RIAM), Kyushu University. Thepresent experiment is a continuation of the workcarried out in a wind tunnel at Colorado StateUniversity (CSU), where the SBL flows were developed over asmooth surface. Stably stratified flows were createdby heating the wind-tunnel airflow to a temperature ofabout 40–50°and by cooling the test-section floor toa temperature of about 10°. To simulate therough surface, a chain roughness was placed over thetest-section floor. We have investigated the buoyancyeffect on the turbulent boundary layer developed overthis rough surface for a wide range of stability,particularly focusing on the turbulence structure andtransport process in the very stable boundary layer.The present experimental results broadly confirm theresults obtained in the CSU experiment with the smoothsurface, and emphasizes the following features: thevertical profiles of turbulence statistics exhibitdifferent behaviour in two distinct stability regimes with weak and strong stability,corresponding to the difference in the verticalprofiles of the local Richardson number. The tworegimes are separated by the critical Richardsonnumber. The magnitudes in turbulence intensities andturbulent fluxes for the weak stability regime aremuch greater than those of the CSU experiments becauseof the greater surface roughness. For the very stableboundary layer, the turbulent fluxes of momentum andheat tend to vanish and wave-like motions due to theKelvin–Helmholtz instability and the rolling up andbreaking of those waves can be observed. Furthermore,the appearance of internal gravity waves is suggestedfrom cross-spectrum analyses.     

General circulation model simulation of mild nuclear winter effects

The climatic effects of an elevated uniform global layer of purely absorbing smoke of absorption optical depth 0.2 have been simulated using a version of the 9-level spectral model of McAvaney et al. (1978). The model was run at rhomboidal wave number 21 with convective adjustment, prognostic precipitation and soil hydrology, but fixed zonally averaged climatological cloud and fixed sea surface temperature, for constant January and July conditions with and without smoke absorption. Results show a reduction in convective rainfall in the tropics and monsoonal regions of the order of 50%, with diurnal average soil surface coolings of several degrees C except in those locations where the reduction in soil moisture is sufficient to effectively stop evaporation at the surface. In that case, small increases in temperature may occur. Results over Australia are consistent with the zonal mean picture. Run in a diurnal cycle mode, the model shows that daily maximum temperatures are more strongly affected, with soil surface coolings of the order of 2°–3° C in summer (with some local warmings) and 4°–6° C in winter. Overninght minimum temperatures cool by only 1°–2° C in both summer and winter. Possible effects of a lowering of sea surface temperature, variations in cloud cover, neglect of scattering by smoke, and infrared absorption and emission by the smoke are discussed.     

Improvement of a soil-atmosphere-transfer model for the simulation of bare soil surface energy balances in semiarid areas

A soil-atmosphere-transfer model (SATM) was evaluated using observational data from the Tongyu Cropland Station and Audubon Research Ranch in semiarid areas, where the land cover was nearly bare soil during the simulation period. Simulations by the SATM at both sites were conducted using the new and original surface thermal roughness length parameterization schemes, respectively. Comparisons of simulations and observations have demonstrated that using the new surface thermal roughness length scheme in this model made sound improvements in the simulation of soil surface temperatures, sensible heat fluxes and net radiation fluxes in the daytime at both sites, compared to the original scheme, because the new scheme produced a larger aerodynamic resistance for turbulent heat transfer in the daytime. With respect to latent heat fluxes, the improvement was not as obvious as that attained for soil surface temperature since the soil water content in the surface layer in a semiarid area is a more important factor than surface soil temperature in controlling evaporation rate. Accordingly, it can be concluded that the new surface thermal roughness length parameterization scheme could improve the ability of the SATM to simulate bare soil surface energy budget with latent heat flux component being innegligible in semiarid areas.     

Air mass modification over Lake Michigan

The modification of a relatively cold air mass over the warm water of Lake Michigan is studied by using a two-dimensional nonlinear mesoscale model. Considerable amounts of heat and water vapor are supplied from the water surface to the lower atmosphere by turbulent eddies. A convective mixed layer develops and grows toward the downwind region with stratocumulus clouds over the lake.The model simulates the warming and moistening of the mixed layer, the development of a boundary layer, the divergence and convergence of wind near the coastlines, and the turbulent fluxes.The model warming of the mixed layer across the lake was about 2.2 °K and the moistening of the mixed layer was about 0.8 g kg^–1, which are comparable to 2.7 °K and 0.8 g kg^–1 observed by Lenschow (1973). The convective boundary layer, which includes the cloud layer, subcloud layer, and superadiabatic layer near the water surface, is well simulated. The tilt of the inversion which coincides with the cloud top is also well reproduced. When a prescribed cooling rate is applied at the cloud top, stronger turbulence and a deeper cloud layer are generated. Without the cooling, the cloud is shallow and the shape of the cloud base is determined by surface conditions. The rise of the inversion is due to upward vertical motion, and deepening of the convective layer in the downwind region.     

Estimation of bulk transfer coefficient for latent heat flux (Ce)

The bulk transfer coefficient for latent heat flux (Ce) has been estimated over the Arabian Sea from the moisture budget during the pre-monsoon season of 1988.The computations have been made over two regions (A: 0–8 ° N; 60–68 ° E; B: 0–4 ° N; 56–60 ° E) with the upper computational boundary fixed at the 300 mb level. The precipitation amount (P) was negligible for region A while the observed values of P have been used for region B. The Ce estimates have been compared with those obtained with other available schemes (Kondo, 1975: Bunker, 1976). which are based on wind speed and atmospheric stability within the surface layer. Our value of Ce is higher in region A and lower in region B than the other estimates.     

Fluxes of gases and particles above a deciduous forest in wintertime

Eddy-correlation measurements of the vertical fluxes of ozone, carbon dioxide, fine particles with diameter near 0.1 m, and particulate sulfur, as well as of momentum, heat and water vapor, have been taken above a tall leafless deciduous forest in wintertime. During the experimental period of one week, ozone deposition velocities varied from about 0.1 cm s^–1 at night to more than 0.4 cm s^-1 during the daytime, with the largest variations associated primarily with changes in solar irradiation. Most of the ozone removal took place in the upper canopy. Carbon dioxide fluxes were directed upward due to respiration and exhibited a strong dependence on air temperature and solar heating. The fluxes were approximately zero at air temperatures less than 5 °C and approached 0.8 mg m^–2 s^–1 when temperatures exceeded 15 °C during the daytime. Fine-particle deposition rates were large at times, with deposition velocities near 0.8 cm s^–1 when turbulence levels were high, but fluxes directed upward were found above the canopy when the surface beneath was covered with snow. Diffusional processes seemed to dominate fine-particle transfer across quasilaminar layers and subsequent deposition to the upper canopy. Deposition velocities for particulate sulfur were highly variable and averaged to a value small in magnitude as compared to similar measurements taken previously over a pine forest in summer.     

The spray contribution to net evaporation from the sea: A review of recent progress

The part that sea spray plays in the air-sea transfer of heat and moisture has been a controversial question for the last two decades. With general circulation models (GCMs) suggesting that perturbations in the Earth's surface heat budget of only a few W m^–2 can initiate major climatic variations, it is crucial that we identify and quantify all the terms in that heat budget. Thus, here we review recent work on how sea spray contributes to the sea surface heat and moisture budgets. In the presence of spray, the near-surface atmosphere is characterized by a droplet evaporation layer (DEL) with a height that scales with the significant-wave amplitude. The majority of spray transfer processes occur within this layer. As a result, the DEL is cooler and more moist than the atmospheric surface layer would be under identical conditions but without the spray. Also, because the spray in the DEL provides elevated sources and sinks for heat and moisture, the vertical heat fluxes are no longer constant with height. We use Eulerian and Lagrangian models and a simple analytical model to study the processes important in spray droplet dispersion and evaporation within this DEL. These models all point to the conclusion that, in high winds (above about 15 m/s), sea spray begins to contribute significantly to the air-sea fluxes of heat and moisture. For example, we estimate that, in a 20-m/s wind, with an air temperature of 20°C, a sea surface temperature of 22°C, and a relative humidity of 80%, the latent and sensible heat fluxes resulting from the spray alone will have magnitudes of order 150 and 15 W/m², respectively, in the DEL. Finally, we speculate on what fraction of these fluxes rise out of the DEL and, thus, become available to the entire marine boundary layer.     

Modelling Regional Scale Surface Energy Exchanges And Cbl Growth In A Heterogeneous, Urban-Rural Landscape

Over the last decade, simple models of theconvective boundary layer (CBL) have beensuggested as an approach to inferring regionallyaveraged land-air exchanges of heat, water and tracegases, because the properties of the CBL respond toan average of the underlying small-scaleheterogeneity. This paper explores the use of anintegral CBL method to infer regionally averagedfluxes in a landscape that has at least three majorsources of heterogeneity – irrigated andnon-irrigated rural land use and a large urban area(Sacramento region, California).The first part of the paper assesses the validity ofthe simple slab model of the CBL – this isintegrated forwards in time using local-scalemeasured heat and water vapour fluxes, to predictmixed-layer depth, temperature and humidity. Of thefour different CBL growth schemes used, the Tennekesand Driedonks model is found to give the bestperformance. Evaluation of the model performancewith different weightings of heat and water vapourfluxes based on the land use characteristics in theregion suggest that the source area for theboundary-layer sonde measurements is larger thanphysically-based estimates would suggest.Finally, measured time series of potentialtemperature are used to infer regionally averagedsensible heat fluxes using an integral CBL (ICBL)method. These ICBL fluxes are compared with thosemeasured at the local scale over the three land usetypes that comprise the region of interest. They arefound to be closest to the heat fluxes calculated byappropriately weighting the measured heat fluxes inthe source area calculated for the ICBL. We concludethat the integral CBL budget method providesadequate estimates of regionally-averaged surfaceheat fluxes in a landscape that is characterised bysurface types with distinctly different surfaceenergy budgets.     

Summer boundary-layer height at the plateau site of Dome’C,antarctica

Measurements of the mean and turbulent structure of the planetary boundary layer using a sodar and a sonic anemometer, and radiative measurements using a radiometer, were carried out in the summer of 1999–2000 at the Antarctic plateau station of Dome C during a two-month period. At Dome C strong ground-based inversions dominate for most of the year. However, in spite of the low surface temperatures (between −50 and −20 °C), and the surface always covered by snow and ice, a regular daytime boundary-layer evolution, similar to that observed at mid-latitudes, was observed during summertime. The mixed-layer height generally reaches 200–300 m at 1300–1400 LST in high summer (late December, early January); late in the summer (end of January to February), as the solar elevation decreases, it reduces to 100–200 m. A comparison between the mixed-layer height estimated from sodar measurements and that calculated using a mixed-layer growth model shows a rather satisfactory agreement if we assign a value of 0.01–0.02 m s⁻¹ to the subsidence velocity at the top of the mixed layer, and a value of 0.003–0.004 K m⁻¹ to the potential temperature gradient above the mixed layer.     

Permafrost Thaw Accelerates in Boreal Peatlands During Late-20th Century Climate Warming

Permafrost covers 25% of the land surface in the northern hemisphere, where mean annual ground temperature is less than 0°C. A 1.4–5.8 °C warming by 2100 will likely change the sign of mean annual air and ground temperatures over much of the zones of sporadic and discontinuous permafrost in the northern hemisphere, causing widespread permafrost thaw. In this study, I examined rates of discontinuous permafrost thaw in the boreal peatlands of northern Manitoba, Canada, using a combination of tree-ring analyses to document thaw rates from 1941–1991 and direct measurements of permanent benchmarks established in 1995 and resurveyed in 2002. I used instrumented records of mean annual and seasonal air temperatures, mean winter snow depth, and duration of continuous snow pack from climate stations across northern Manitoba to analyze temporal and spatial trends in these variables and their potential impacts on thaw. Permafrost thaw in central Canadian peatlands has accelerated significantly since 1950, concurrent with a significant, late-20th-century average climate warming of +1.32 °C in this region. There were strong seasonal differences in warming in northern Manitoba, with highest rates of warming during winter (+1.39 °C to +1.66 °C) and spring (+0.56 °C to +0.78 °C) at southern climate stations where permafrost thaw was most rapid. Projecting current warming trends to year 2100, I show that trends for north-central Canada are in good agreement with general circulation models, which suggest a 4–8 °C warming at high latitudes. This magnitude of warming will begin to eliminate most of the present range of sporadic and discontinuous permafrost in central Canada by 2100.     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

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.     

A new method to estimate ice age temperatures

On glacial time scales, the waxing and waning of the Eurasian and North American ice sheets depend largely on variations in atmospheric temperature. As global sea level is primarily determined by the volume of these ice sheets, there is a direct (yet complex) relation between global sea level and the northern hemispheric (NH) temperature. This relation is essentially represented by a model of the NH ice sheets. We use a thermomechanical ice-sheet–ice shelf–bedrock model in conjunction with an inverse method to deduce a time series of NH temperature (from 120 kyr BP until present) that is consistent with the observed global sea level record. The advantage of this method is that it provides the annual mean surface air temperature averaged over the NH continents north of 40°N. The results reveal that ice age temperatures were 4–10°C lower than today, which agrees with other temperature reconstructions. However, reconstructed temperatures are comparitively low during the early stages of the glacial, a feature that is consistent with the rapid growth of the ice sheets. The sensitivity of the results for uncertainties in precipitation rate, in observed sea level and in some other model parameters is examined to quantify the error in reconstructed temperature. During the glacial period (120–15 kyr BP), surface air temperatures in the NH (north of 40°N) were 7.2±1.5°C lower than todays (interglacial) temperatures.     

A study of mean boundary-layer structures over the Arabian Sea and the Bay of Bengal during active and break monsoon periods

Observations from research ships which took part in the Indo-Soviet Monsoon Experiment of 1977 (MONSOON 77) and the International Monsoon Experiments (MONEX 79) over the central Arabian Sea and the north central Bay of Bengal were analyzed to study the mean wind and temperature structure of the monsoon boundary layer during active and break conditions. Mean profiles of wind speed and direction along with virtual potential temperature obtained by averaging data from several research ships during 1977 and 1979 indicate that onset conditions were associated with substantial increases in wind speed over the Arabian Sea and a shift to strong southwest flow. Monsoon onset was also characterized by near-neutral to slightly unstable temperature profiles in the lowest kilometer. Break conditions in 1977 in which the monsoon trough moved northward and substantial (5 mb) pressure rises were noted over the Arabian Sea show wind speeds typically decreasing from approximately 18 m s^–1 during active conditions to roughly 8 m s ^–1. Temperature profiles during break conditions are similar to those observed in pre-monsoon conditions in that the boundary layer is observed to be generally much more stable up to 900 mb. Above 900 mb, profiles of virtual potential temperature show little variation.Analysis of latent and sensible heat fluxes during June 1977 calculated by the bulk aerodynamic method indicates values of latent heat flux during active conditions to be roughly two to three times larger than those during break conditions. Sensible heat flux shows an increase from approximately 20 to 80 W m ^–1 during the onset of the monsoon. Surface fluxes of water vapor indicate the importance of water vapor transport over the ship observation region in the central Arabian Sea during active conditions. Onset of the monsoon over the Arabian Sea is accompanied by an increase in the surface moisture flux by a factor of about two. Time histories of precipitable water show decreases of approximately 15% from active to break periods.     

Local climates simulated by two generations of Canadian GCM land surface schemes

Abstract

The performance of two Canadian land surface schemes of widely differing complexity is compared and contrasted in a pair of year‐long simulations using the GCM developed at Atmospheric Environment Service, Canada. The old land surface model incorporates the force‐restore method for soil temperatures and the bucket approximation for soil moisture; the new model, CLASS (Canadian Land Surface Scheme) features three soil layers, an explicitly modelled snow layer, a thermally separate vegetation canopy, and physically‐based calculations of heat and moisture transfers between all of the land surface components and the atmosphere.

It was reported in previous papers that compared with observations, the old scheme tends to generate a climate which is characterized by anomalously high precipitation rates and cold temperatures over land. In this paper, by reference to field measurements and to the energy fluxes and temperatures generated by the two models at local scales, the hypotheses earlier postulated as to the underlying reasons for this are validated. The main factor contributing to the climate anomalies observed with the old scheme is found to be its generation of excessive evaporation rates; this is caused by the fact that the evaporation rate is never directly energy‐limited, the fact that the scaling of the evaporation rale with decreasing soil moisture content underestimates the effect of vegetation stomatal resistance, and the fact that the evaporation rate over bare soil depends not on the surface soil moisture, but on the moisture content of whole modelled soil column. The cold surface temperatures are additionally attributed to systematic errors incurred by the forward‐stepping temperature scheme, and to the fact that soils subjected to subzero temperature forcing in the winter are modelling as freezing completely. Finally, the inability of the old scheme to simulate partially frozen soils means that it proves unable to handle either shallow frost penetration at temperature latitudes, or the development of an active layer in permafrost.     

Temperature stratification of the atmospheric boundary layer over the Deccan Plateau,India, during the summer monsoon

The vertical and horizontal temperature structure of the atmospheric boundary layer (ABL) were studied using aircraft observations made in the lowest 2.4 km above ground level during the summer monsoon.The vertical temperature structure of the ABL in the region may be classified into the following four categories.Category The ABL consisted of two layers of thickness 700–900 m separated by a thin transition layer. The lapse rates in the former two layers were dry adiabatic.Category The lowest layer of the ABL of thickness 400–600 m was adiabatically stratified and the overlying layer was stable with gradients of potential temperature 4–5°C km^–1. The stable layer contained a thin adiabatic stratified layer of 200–300 m thickness at a height of 1.5 km.Category The lowest 200–400 m layer of the ABL was adiabatically stratified and the overlying layer was stable with potential temperature gradients of 5–6 °C km¹.Category The ABL was mainly stable with potential temperature gradients of 6 °C km^–1 or greater. Occasionally thin layers with adiabatic stratification were found embedded in the ABL.The temperature distribution of the horizontal temperature at 900 m was mainly normal. The high-frequency portion of the spectra lying between 0.05 and 0.16 Hz (corresponding to wave length 1 km to 300 m) oscillated around the –\2/3 power law line. The spectral curve showed a significant peak at 0.011 Hz having a wave-length of 5 km.Department of Geoscience, North Carolina State University, Raleigh, NC, 27650, U.S.A.     

Observations of the convective plume of a lake under cold-air advective conditions

Moderating effects of Lake Apopka, Florida on downwind surface temperatures were evaluated under cold-air advective conditions. Point temperature measurements north and south of the lake and data obtained from a thermal scanner flown at 1.6 km indicate that surface temperatures directly downwind may be higher than surrounding surface temperatures by as much as 5 °C under conditions of moderate winds (~4 m s^–1). No substantial temperature effects were observed with surface wind speed less than 1 m s^–1. Fluxes of sensible and latent heat from Lake Apopka were calculated from measurements of lake temperature, net radiation, relative humidity and air temperature above the lake. Bulk transfer coefficients and the Bowen ratio were calculated and found to be in agreement with reported data for non-advective conditions.IFAS Journal Series No. 1006.     

Thermodynamic structure of the Atmospheric Boundary Layer over the Arabian sea as revealed by MONSOON-77 data

The thermodynamic structure of the Atmospheric Boundary Layer (ABL) over the Arabian sea region has been studied with the help of 135 aerological observations obtained during MONSOON-77 in the region (10–14° N, 64–68° E) by USSR research vessels. Low-level inversions were observed over the western Arabian sea region (west of 66° E) in association with suppressed convection. The different sublayers of the ABL, viz. the mixed layer, the cloud layer and the inversion/isothermal/stable layer were identified. The low-level stability analysis indicated that in the region east of 66° E, conditions were favourable for deep convection. The thermodynamic transformation of the boundary layer after precipitation was documented.