The position of the lowest levels in the boundary layer of atmospheric circulation models

Abstract

The error associated with the position of the lowest wind level in atmospheric boundary‐layer modelling is studied in connection with vertical resolutions typical of parametrization schemes for atmospheric circulation models. The test case is the neutrally stratified steady‐state boundary layer. Finite‐difference and finite‐element schemes of two types are used: in one case the lowest wind level is centred with respect to the surface and the lowest internal level where the shear stress is calculated, in the other case the lowest wind level is set very close to the surface (5 m). It is found that schemes of the latter type underestimate the friction force at the lowest level and consequently overestimate the wind and the surface stress. This error is largest at low resolution since it is due to the uncentring of the lowest wind level with respect to the stress levels. The error in schemes of the former type is different, and is associated with the determination of the surface stress from a wind at a height that may exceed the extent of the surface layer. For all neutral cases, this error can be made small by adding a corrective term to the traditional logarithmic formulation. This paper shows that considerably more error is generated by uncentred differencing than by deepening the surface layer.     

Testing asynchronous coupling on simple “ocean-atmosphere” dynamic systems

 Simulations of long-term climatic variations require the use of coupled oceanic and atmospheric models. Given their different scales in space and time and present computing limitations, asynchronous coupling schemes are often used. Comparing the evolutions of simple dynamic systems which we couple alternatively with synchronous and asynchronous schemes, we show that the behaviour can be modified, not only in minute details, but drastically, if asynchronous coupling is used with insufficient care. We study how it happens, then suggest methods to insure that the average behaviour of more realistic models is not modified with asynchronous schemes. Other branches of physics with coupled non-linear systems could make use of these remarks. Received: 19 February 1997 / Accepted: 17 June 1998     

CMIP6模式大气中南北半球际涛动的季节循环

利用1958—2014年47个CMIP6模式输出资料和NCEP/NCAR再分析资料,研究了模式大气中南北涛动(InterHemispheric Oscillation,IHO)的季节变化特征,且评估了CMIP6对IHO季节特征的模拟能力。结果表明:47个CMIP6模式都能模拟出IHO的季节演变特征,但模式间存在一定差异。通过比较,筛选出模拟IHO季节循环较好的16个模式,它们能成功模拟出半球大气质量的时间演变和空间结构。进一步分析表明,水汽对IHO季节变化有抵消作用且半球内部水汽质量变化可驱动越赤道质量流的产生;地表净短波辐射夏高冬低,其加热造成的水汽蒸发在水汽质量变化中起到重要作用;地表净长波辐射在春秋变化幅度较大,与大气质量逐月变化吻合。对比再分析资料表明,CMIP6模式模拟的半球大气质量的峰谷值变化有明显的月份偏差,且CMIP6模式模拟的地表气压异常值的偏差主要出现在北太平洋、欧亚大陆、南半球中纬度地区和两极极区,模拟的南北半球的蒸发和降水量、赤道风场、地表净长波和短波辐射通量等均存在明显的偏差。     

Prairie agroclimate boundary‐layer model: A simulation of the atmosphere/crop‐soil interface

Abstract

This paper describes a 1‐D agroclimatic model of the atmosphere/crop‐soil interface. Vertical profiles of wind, potential temperature and water vapour are constructed twice daily for the overnight‐low and maximum temperature times by combining 1200 and 0000 UTC upper‐air standard‐level grid‐point data with climatological observations. The vertical structure of the atmospheric boundary layer has a surface constant‐flux layer that is usually topped by a mixed layer by day but not at night. The crop‐soil boundary layer consists of a shallow top‐zone and a growing root‐zone. Vegetation cover and root depth depend upon crop type and phenological stage. Water‐balance accounting tracks the moisture contents of both the top‐ and root‐zones. Evapotranspiration or the vertical flux of water vapour in the atmospheric boundary layer is tied to the evolution of the crop‐soil boundary layer.

The model was calibrated using field data from the Regional Evaporation Study's primary site in an agricultural area of central Saskatchewan. The evolution of 1991's wheat‐soil boundary layer from the crop's heading to ripe stages was then successfully simulated at two additional sites in the same geographical area.     

Towards an equitable global climate change regime: compatibility with Article 2 of the Climate Change Convention and the link with sustainable development

Abstract

In this paper we argue that the discussion on how to get to an equitable global climate change regime requires a long-term context. Some key dimensions of this discussion are responsibility, capability and development needs. Each of these, separately or in combination, has been used in designing schemes for differentiation of commitments to limit or reduce greenhouse gas emissions. Many implementation problems of these proposals are often side-stepped. In particular, some proposals may be incompatible with Article 2 of the United Nations Framework Convention on Climate Change (UNFCCC), i.e. they are unlikely to keep the option open of long-term stabilisation of atmospheric greenhouse gas concentrations at relatively low levels. We present some evidence that shifting the emphasis from emission reduction to sustainable development needs can contribute significantly to relieving the threat of human-induced climate change.     

Atmospheric moisture budget estimates of regional evapotranspiration from RES‐91

Abstract

Diurnal changes in the local atmospheric moisture budget over the Canadian Prairies are computed using sequential radiosonde soundings from the 1991 Regional Evaporation Study (RES‐91). Previous attempts to estimate evapotranspiration with radiosonde data have used either similarity theory or a moisture budget, but have been confined to the boundary layer in either case. These studies, as well as semi‐empiric operational techniques which use surface‐based data, exclude the effects of moisture advection and energy exchanges between the boundary layer and the free atmosphere, assuming negligible effects on evapotranspiration. The moisture budget method adopted here includes horizontal advection explicitly, and treats vertical fluxes implicitly through a total tropospheric moisture budget.

Comparison of the evapotranspiration estimates with those of other techniques are positive only when results are averaged over several days to weeks. While the advection estimates are a major source of error for the “daily” estimates in this particular study, it is shown that neither advection nor moisture flux through the boundary layer can be ignored in estimating daily evapotranspiration, regardless of the technique used. The results also suggest that evapotranspiration is more variable on a daily basis than other techniques have indicated. With an improved synoptic database now available for advection estimates, the moisture budget technique may provide an excellent ground‐truth method for fine‐tuning techniques for remote sensing of evapotranspiration, and could lead to improved parametrization schemes for both NWP models and GCMs.     

Simulation of the Asian Summer Monsoon in Five European General Circulation Models

A comparison is made of the mean monsoon climatology in five different general circulation models (GCMs) which have been used by the participants of a project, funded by the European Union, entitled Studies of the Influence, Hydrology and Variability of the Asian summer monsoon (SHIVA). The models differ considerably, in horizontal and vertical resolution, numerical schemes and physical parametrizations, so that it is impossible to isolate the cause of differences in their monsoon simulations. Instead, the purpose of this comparison is to document and compare the representation of the mean monsoon in models which are being used to investigate the characteristics of the monsoon, its variability and its response to different boundary forcings. All of the models produce a reasonable representation of the monsoon circulation, although there are regional variations in the magnitude and pattern of the flow at both 850 hPa and 200 hPa. Considerable differences between the models are seen in the amount and distribution of precipitation. The models all reproduce the basic monsoon seasonal variation, although the timing of the onset and retreat, and the maxima in the winds and precipitation during the established phase, differ between them. There are corresponding differences in the evolution of the atmospheric structure between the pre-monsoon season and its established phase. It is hoped that this study will set in context the investigations of the monsoon system and its impacts carried out using these models, both during SHIVA and in the future.     

UV attenuation in the cloudy atmosphere

Ultraviolet (UV) energy absorption plays a very important role in the Earth–atmosphere system. Based on observational data for Beijing, we suggest that some atmospheric constituents utilize or transfer UV energy in chemical and photochemical (C&P) reactions, in addition to those which absorb UV energy directly. These constituents are primarily volatile organic compounds (VOCs) emitted from both vegetative and anthropogenic sources. The total UV energy loss in the cloudy atmosphere for Beijing in 1990 was 78.9 Wm⁻². This attenuation was caused by ozone (48.3 Wm⁻²), other compounds in the atmosphere (26.6 Wm⁻²) and a scattering factor (4.0 Wm⁻²). Our results for a cloudy atmosphere in the Beijing area show that the absorption due to these other compounds occurs largely through the mediation of water vapor. This fraction of energy loss has not been fully accounted for in previous models. Observations and previous models results suggest that 1) a cloudy atmosphere absorbs 25∼30 Wm⁻² more solar shortwave radiation than models predict; and 2) aerosols can significantly decrease the downward mean UV-visible radiation and the absorbed solar radiation at the surface by up to 28 and 23 Wm⁻², respectively. Thus, quantitative study of UV and visible absorption by atmospheric constituents involved in homogeneous and heterogeneous C&P reactions is important for atmospheric models.     

Sensitivity of transient eddies to climate change in the CCC general circulation model

Abstract

We use eddy life‐cycle simulations to evaluate the response of atmospheric transient eddies to a global warming caused by CO₂ doubling in the CCC general circulation model. In simulations using Northern Hemisphere winter conditions, transient waves attain larger kinetic energy and encompass a wider range of latitudes in the warmer climate. This behaviour contrasts with a previous investigation that used output from the NCAR and GFDL models. Our analysis indicates two primary factors for the difference between model responses: (1) a smaller change in the mid‐latitude temperature gradient in the CCC model, which allows (2) increased atmospheric water vapour in mid‐latitudes to catalyze a more rapidly evolving life‐cycle.     

The Response of Arctic Sea Ice to Global Change

The sea ice-covered polar oceans have received wider attention recently for two reasons. Firstly, the global conveyor belt circulation of the ocean is believed to be forced in the North and South Atlantic through deep water formation, which to a large degree is controlled by the variations of the sea ice margin and especially by the sea ice export to lower latitudes. Secondly, CO₂ response experiments with coupled climate models show an enhanced warming in polar regions for increased concentrations of atmospheric greenhouse gases. Whether this large response in high latitudes is due to real physical feedback processes or to unrealistic simplifications of the sea ice model component remains to be determined. Coupled climate models generally use thermodynamic sea ice models or sea ice models with oversimplified dynamics schemes. Realistic dynamic-thermodynamic sea ice models are presently implemented only at a few modeling centers. Sensitivity experiments with thermodynamic and dynamic-thermodynamic sea ice models show that the more sophisticated models are less sensitive to perturbations of the atmospheric and oceanic boundary conditions. Because of the importance of the role of sea ice in mediating between atmosphere and ocean an improved representation of sea ice in global climate models is required. This paper discusses present sea ice modeling as well as the sensitivity of the sea ice cover to changes in the atmospheric boundary conditions. These numerical experiments indicate that the sea ice follows a smooth response function: sea ice thickness and export change by 2% of the mean value per 1 Wm^-2 change of the radiative forcing.     

Multi-species chemical data assimilation with the Danish Eulerian hemispheric model: system description and verification

Satellite retrievals of atmospheric composition provide a wealth of data on a global scale. These complement results from atmospheric chemistry-transport models (CTMs), and can be combined using data assimilation. We present two assimilation schemes coupled to the Danish Eulerian Hemispheric Model (DEHM), a three-dimensional, off-line CTM with full photochemistry: a variant on the ensemble Kalman filter and the three-dimensional variational scheme. The aim of this paper is to describe the two schemes and present an initial assessment of their impacts on model skill. Retrievals of multiple atmospheric trace gases are assimilated, namely: NO₂ tropospheric column densities, CH₄ total column densities, and partial column concentrations of O₃, CO and CH₄; these data are retrieved from four satellite sensors. Data for each species are assimilated independently of one another, and other species are only adjusted indirectly via the model’s chemistry and dynamics. Assimilation results are compared with measurements from surface monitoring stations and other satellite retrievals, and preliminary validation results are presented.Reference simulations (without assimilation) grossly underestimate surface CO concentrations, and both assimilation schemes eliminate this large and systematic model bias. The assimilation improves the spatial correlation of modelled CO with surface observations, and improves the spatial correlation between forecasts and retrievals for CO, NO₂ and O₃. Results for CH₄ show a loss of skill due to a mismatch in model bias between two assimilated CH₄ data-sets. Finally, we discuss differences in methodology and results between this paper and a recent study on multi-species chemical data assimilation. Joint optimisation of initial conditions and emission rates offers a promising direction for improving modelled boundary-layer concentrations.     

A coupled atmosphere‐ocean model for transient climate change studies

Abstract

A new coupled atmosphere‐ocean model has been developed for climate predictions at decade to century scales. The atmospheric model is similar to that of Hansen et al. (1983) except that the atmospheric dynamic equations for mass and momentum are solved using Arakawa and Lamb's (1977) C grid scheme and the advection of potential enthalpy and water vapour uses the linear upstream scheme (Russell and Lerner, 1981). The new global ocean model conserves mass, allows for divergent flow, has a free surface and uses the linear upstream scheme for the advection of potential enthalpy and salt. Both models run at 4° × 5° resolution, with 9 vertical layers for the atmosphere and 13 layers for the ocean. Twelve straits are included, allowing for subgrid‐scale water flow. Runoff from land is routed into appropriate ocean basins. Atmospheric and oceanic surface fluxes are of opposite sign and are applied synchronously. Flux adjustments are not used. Except for partial strength alternating binomial filters (Shapiro, 1970), which are applied to the momentum components in the atmosphere and oceans, there is no explicit horizontal diffusion.

A 120‐year simulation of the coupled model starting from the oceanic initial conditions of Levitus (1982) is discussed. The model dynamics stabilize after several decades. The maximum northward ocean heat flux is 1.4 × 10¹⁵ W at 16°N. The model appears to maintain the vertical gradients characterizing the separation between the upper and deep ocean spheres. Inadequacies in the coupled model simulation lead to decreasing temperature and salinity in the high latitude North Atlantic and to a poor simulation of the northern North Atlantic thermohaline circulation. The mass transport of the Gulf Stream is about half of observed values, while the transports of the Kuroshio and Antarctic Circumpolar Currents are similar to observations. Additional deficiencies include a climate drift in the surface air temperature of 0.006°C year^‐1 due to a radiation imbalance of 7.4 Wm^‐2 at the top of the atmosphere and too warm temperatures in the eastern portions of tropical oceans. The coupled model should be useful for delineating modelling capabilities without the use of flux adjustments and should serve as a benchmark for future model improvements.     

Notes and correspondence

Abstract

This paper is an expansion of the lectures presented by the authors at the Canadian Middle Atmosphere Model (CMAM) summer school, 25–29 August 1997. It attempts to fulfill several goals as did the curriculum for the summer school. Firstly, it has a large tutorial component allowing newcomers a rapid introduction to the field, secondly it surveys the field of stratospheric chemistry, and thirdly it attempts to bring the reader to the forefront of some of the problem areas in stratospheric chemistry. As such, it includes much background information such as the magnitude of rate coefficients, photolysis rates and their calculation and time constants. The intended audience is individuals about to commence study or research in the field of atmospheric chemistry or established atmospheric scientists in other areas wishing to expand their knowledge of this topic.

The basic chemistry of the ozone layer is described with an emphasis on the catalytic cycles and the equilibria between “reservoir” and “active” species. To this end, a compendium of the important aspects of heterogeneous chemistry is also provided. The application of box (0‐dimensional) and more elaborate (1‐, 2‐ and 3‐dimensional) models to compare and validate with measurements is discussed. Also, factors affecting the present and future trends in ozone concentration, both at mid‐latitudes and in the polar region, are presented.     

Skew eddy fluxes as signatures of non‐linear tidal current interactions,with application to Georges Bank

Abstract

The eddy flux of a conservative scalar in a time‐dependent rotary velocity field may have a component that is normal to the scalar gradient. This component is the “skew flux”, which consists of the scalar transport by the Stokes velocity and a part that is always non‐divergent (and hence does not affect scalar evolution). Since tidal velocity fields usually have rotary features, tidal‐band eddy scalar fluxes may include a skew component that can be useful in indicating the occurrence of non‐linear current interactions.

The skew temperature flux associated with the semidiurnal tide in a continental shelf region is demonstrated using simple models, and moored current and temperature observations from Georges Bank. The observed fluxes on the Bank are largely directed along isobaths, with apparent contributions from the topographic rectification of the barotropic tidal current over the Bank's side and from the rotary tidal ellipses in a frontal region. Simple models indicate that the weaker cross‐isobath fluxes can arise through the influence of frictionally induced vertical structure on topographic tidal rectification, a baroclinic tidal current interaction, or the interaction of baroclinic and barotropic tidal currents. In some cases, the simple models show qualitative agreement with the observed fluxes and currents but, in general, more realistic models and better estimates of the background mean temperature field are required to obtain quantitative estimates of the relative importance of these interactions and other processes. Nevertheless, the observations and models suggest that non‐linear interactions involving both barotropic and baroclinic tidal currents are occurring on Georges Bank.     

Near-source particulate emissions and plume dynamics from agricultural field operations

Particulate matter emissions generated by agricultural field preparation and harvesting operations were measured remotely via aerosol lidar and sampled simultaneously with a variety of aerosol point samplers in order to quantify dust plume space and time dynamics and particulate mass and number concentrations. Data for two cotton operations (disking, harvesting) in a flood-irrigated field in New Mexico are presented. Dust plume dynamics varied with boundary layer meteorological conditions, especially atmospheric stability, with plume maximum height significantly lower under stable conditions. Plume tracking indicated little change in plume area with height under unstable conditions and plume movement depended on wind speed and direction. Particle mass distributions indicate approximately 50% of the measured PM₁₀ mass was PM_2.5, significantly higher than previously reported values, possibly due to the near-source nature of the samples collected here. Variability in plume movement matched the variability in short-term wind fluctuations and this variability helps explain why models that utilize long-term averages perform poorly when trying to capture plume dynamics for nonsteady sources such as tractor operations.     

Comparison of different boundary layer surface schemesusing single point micrometeorological field data

Summary  In the last decade, a vast number of land surface schemes has been designed for use in global climate models, atmospheric weather prediction, mesoscale numerical models, ecological models, and models of global changes. Since land surface schemes are designed for different purposes they have various levels of complexity in the treatment of bare soil processes, vegetation, and soil water movement. This paper is a contribution to a little group of papers dealing with intercomparison of differently designed and oriented land surface schemes. For that purpose we have chosen three schemes for classification: i) global climate models, BATS (Dickinson et al., 1986; Dickinson et al., 1992); ii) mesoscale and ecological models, LEAF (Lee, 1992) and iii) mesoscale models, LAPS (Mihailović, 1996; Mihailović and Kallos, 1997; Mihailović et al., 1999) according to the Shao et al. (1995) classification. These schemes were compared using surface fluxes and leaf temperature outputs obtained by time integrations of data sets derived from the micrometeorological measurements above a maize field at an experimental site in De Sinderhoeve (The Netherlands) for 18 August, 8 September, and 4 October 1988. Finally, comparison of the schemes was supported applying a simple statistical analysis on the surface flux outputs. Received November 23, 1999 Revised August 18, 2000     

The deep water stratification of ocean general circulation models

Abstract

A central problem in climate and ocean modelling is the accurate simulation of the climatological state of the oceanic density field. A constant vertical diffusivity for heat and salt is frequently employed in ocean general circulation models (OGCMs) and it is usually assigned a value designed to optimize the depth of the pycnocline. One undesired consequence of this choice is a poor representation of the deep water, which is usually insufficiently stratified. In contrast to the uniform diffusivity of many models, some observational studies suggest that the vertical diffusivity is not constant but increases with depth, possibly in inverse proportion to the local buoyancy frequency. Numerical experiments with an OGCM are presented that demonstrate that allowing the vertical diffusivity to increase below the pycnocline substantially increases the stratification of the abyssal water mass of these models without significantly affecting the pycnocline depth, and hence may lead to a better representation of the vertical density structure.     

Single-Column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, having a significant impact on medium-range weather forecasts and climate integrations. To pursue this further, a moderately stratified Arctic case is simulated by nineteen single-column turbulence schemes. Statistics from a large-eddy simulation intercomparison made for the same case by eleven different models are used as a guiding reference. The single-column parameterizations include research and operational schemes from major forecast and climate research centres. Results from first-order schemes, a large number of turbulence kinetic energy closures, and other models were used. There is a large spread in the results; in general, the operational schemes mix over a deeper layer than the research schemes, and the turbulence kinetic energy and other higher-order closures give results closer to the statistics obtained from the large-eddy simulations. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.     

ENSO in a hybrid coupled model. Part II: prediction with piggyback data assimilation

A hybrid coupled model (HCM) for the tropical Pacific ocean-atmosphere system is employed for ENSO prediction. The HCM consists of the Geophysical Fluid Dynamics Laboratory ocean general circulation model and an empirical atmospheric model. In hindcast experiments, a correlation skill competitive to other prediction models is obtained, so we use this system to examine the effects of several initialization schemes on ENSO prediction. Initialization with wind stress data and initialization with wind stress reconstructed from SST using the atmospheric model give comparable skill levels. In re-estimating the atmospheric model in order to prevent hindcast-period wind information from entering through empirical atmospheric model, we note some sensitivity to the estimation data set, but this is considered to have limited impact for ENSO prediction purposes. Examination of subsurface heat content anomalies in these cases and a case forced only by the difference between observed and reconstructed winds suggests that at the current level of prediction skill, the crucial wind components for initialization are those associated with the slow ENSO mode, rather than with atmospheric internal variability. A “piggyback” suboptimal data assimilation is tested in which the Climate Prediction Center data assimilation product from a related ocean model is used to correct the ocean initial thermal field. This yields improved skill, suggesting that not all ENSO prediction systems need to invest in costly data assimilation efforts, provided the prediction and assimilation models are sufficiently close. Received: 17 April 1998 / Accepted: 22 July 1999