Validation of ECMWF and NCEP–NCAR Reanalysis Data in Antarctica

The European Center for Medium-Range Weather Forecast (ECMWF) Re-Analysis (ERA-40) and the National Centers for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) ECMWF (ERA-40) and NCEP-NCAR reanalysis data were compared with Antarctic station observations, including surface-layer and upper-layer atmospheric observations, on intraseasonal and interannual timescales. At the interannual timescale, atmospheric pressure at different height levels in the ERA-40 data are in better agre...     

Dynamical downscaling of ERA-40 in complex terrain using the WRF regional climate model

Results from a first-time employment of the WRF regional climate model to climatological simulations in Europe are presented. The ERA-40 reanalysis (resolution 1°) has been downscaled to a horizontal resolution of 30 and 10?km for the period of 1961?C1990. This model setup includes the whole North Atlantic in the 30?km domain and spectral nudging is used to keep the large scales consistent with the driving ERA-40 reanalysis. The model results are compared against an extensive observational network of surface variables in complex terrain in Norway. The comparison shows that the WRF model is able to add significant detail to the representation of precipitation and 2-m temperature of the ERA-40 reanalysis. Especially the geographical distribution, wet day frequency and extreme values of precipitation are highly improved due to the better representation of the orography. Refining the resolution from 30 to 10?km further increases the skill of the model, especially in case of precipitation. Our results indicate that the use of 10-km resolution is advantageous for producing regional future climate projections. Use of a large domain and spectral nudging seems to be useful in reproducing the extreme precipitation events due to the better resolved synoptic scale features over the North Atlantic, and also helps to reduce the large regional temperature biases over Norway. This study presents a high-resolution, high-quality climatological data set useful for reference climate impact studies.     

Vertical structure of recent arctic warming from observed data and reanalysis products

Spatiotemporal patterns of recent (1979–2008) air temperature trends are evaluated using three reanalysis datasets and radiosonde data. Our analysis demonstrates large discrepancies between the reanalysis datasets, possibly due to differences in the data assimilation procedures as well as sparseness and inhomogeneity of high-latitude observations. We test the robustness of arctic tropospheric warming based on the ERA-40 dataset. ERA-40 Arctic atmosphere temperatures tend to be closer to the observed ones in terms of root mean square error compared to other reanalysis products used in the article. However, changes in the ERA-40 data assimilation procedure produce unphysical jumps in atmospheric temperatures, which may be the likely reason for the elevated tropospheric warming trend in 1979–2002. NCEP/NCAR Reanalysis data show that the near-surface upward temperature trend over the same period is greater than the tropospheric trend, which is consistent with direct radiosonde observations and inconsistent with ERA-40 results. A change of sign in the winter temperature trend from negative to positive in the late 1980s is documented in the upper troposphere/lower stratosphere with a maximum over the Canadian Arctic, based on radiosonde data. This change from cooling to warming tendency is associated with weakening of the stratospheric polar vortex and shift of its center toward the Siberian coast and possibly can be explained by the changes in the dynamics of the Arctic Oscillation. This temporal pattern is consistent with multi-decadal variations of key arctic climate parameters like, for example, surface air temperature and oceanic freshwater content. Elucidating the mechanisms behind these changes will be critical to understanding the complex nature of high-latitude variability and its impact on global climate change.     

CAM3 bias over the Arctic region during northern winter studied with a linear stationary model

This study builds upon two prior papers, which examine Arctic region bias of CAM3 (NCAR Community Atmosphere Model version 3) simulations during winter. CAM3 output is compared with ECMWF (European Centre for Medium-Range Weather Forecasts) 40?year reanalysis (ERA-40) data. Our prior papers considered the temperature and the vorticity equation terms and demonstrated that diabatic, transient, and linear terms dominate nonlinear bias terms over most areas of interest. Accordingly, this paper uses a linearized form of the model??s dynamical core equations to study aspects of the forcing that lead to the CAM3 biases. We treat the model??s long term winter bias as a solution to a linear stationary wave model (LSWM). Key features of the bias in the vorticity, temperature, and ln of surface pressure (=q) fields are shown at medium resolution. The important features found at medium resolution are captured at the much lower LSWM resolution. The Arctic q bias has two key features: excess q over the Barents Sea and a missing Beaufort High (negative maximum q bias) to the north of Alaska and eastern Siberia. The forcing fields are calculated by the LSWM. Horizontal advection tends to create multi-polar combinations of negative and positive extrema in the forcing. The positive and negative areas of forcing approximately match corresponding areas in the bias. There is a broad relation between cold bias with elevated q bias, as expected from classical theory. Forcing in related quantities: near surface vorticity and surface pressure combine to produce the sea level pressure bias.     

WRF模式对青藏高原那曲地区大气边界层模拟适用性研究

采用WRF（Weather Research and Forecasting）模式4种边界层参数化方案对青藏高原那曲地区边界层特征进行了数值模拟,并利用"第三次青藏高原大气科学试验"在青藏高原那曲地区5个站点的观测资料对模拟结果进行验证,分析不同参数化方案在那曲地区的适用性。研究表明,YSU、MYJ、ACM2和BouLac方案对2 m气温和地表温度的模拟偏低。BouLac方案模拟的地表温度偏差较小。通过对能量平衡各分量的对比分析发现,温度模拟偏低可能是向下长波辐射模拟偏低以及感热通量和潜热通量交换过强导致的。对于边界层风、位温和相对湿度垂直结构的模拟,局地方案的模拟效果均优于非局地方案。BouLac方案对那曲地区近地层温度、边界层内位温和相对湿度的垂直分布模拟效果较好。      

Evaluation of an ensemble of Arctic regional climate models: spatiotemporal fields during the SHEBA year

Simulations of eight different regional climate models (RCMs) have been performed for the period September 1997–September 1998, which coincides with the Surface Heat Budget of the Arctic Ocean (SHEBA) project period. Each of the models employed approximately the same domain covering the western Arctic, the same horizontal resolution of 50 km, and the same boundary forcing. The models differ in their vertical resolution as well as in the treatments of dynamics and physical parameterizations. Both the common features and differences of the simulated spatiotemporal patterns of geopotential, temperature, cloud cover, and long-/shortwave downward radiation between the individual model simulations are investigated. With this work, we quantify the scatter among the models and therefore the magnitude of disagreement and unreliability of current Arctic RCM simulations. Even with the relatively constrained experimental design we notice a considerable scatter among the different RCMs. We found the largest across-model scatter in the 2 m temperature over land, in the surface radiation fluxes, and in the cloud cover which implies a reduced confidence level for these variables. An erratum to this article can be found at     

Storm track sensitivity to sea surface temperature resolution in a regional atmosphere model

A high resolution regional atmosphere model is used to investigate the sensitivity of the North Atlantic storm track to the spatial and temporal resolution of the sea surface temperature (SST) data used as a lower boundary condition. The model is run over an unusually large domain covering all of the North Atlantic and Europe, and is shown to produce a very good simulation of the observed storm track structure. The model is forced at the lateral boundaries with 15–20 years of data from the ERA-40 reanalysis, and at the lower boundary by SST data of differing resolution. The impacts of increasing spatial and temporal resolution are assessed separately, and in both cases increasing the resolution leads to subtle, but significant changes in the storm track. In some, but not all cases these changes act to reduce the small storm track biases seen in the model when it is forced with low-resolution SSTs. In addition there are several clear mesoscale responses to increased spatial SST resolution, with surface heat fluxes and convective precipitation increasing by 10–20% along the Gulf Stream SST gradient.     

U.S. Climate Sensitivity Simulated with the NCEP Regional Spectral Model

10-year continuous U.S. climate simulations were conducted with the Regional Spectral Model (RSM) using boundary conditions from the National Centers for Environmental Prediction/Dept. of Energy reanalyses and the global PCM (Parallel Climate Model) simulations for present day (1986–1996) andfuture (2040–2050) CO₂ concentrations (about a 36% increasedCO₂). In order to examine the influence of physical parameterization differences as well as grid-resolution, fine resolution RSM simulations (50 km) were compared to coarse resolution (180 and 250 km) RSM simulations, which had resolutions comparable to the T62 reanalysis and PCM simulations. During the winter, the fine resolution RSM simulations provided more realistic detail over the western mountains. During the summer, large differences between the RSM and driving PCM simulations were found. Our results with presentCO₂ suggest that most of the differences between the regionalclimate model simulations and the climate simulations driven by the global model used to drive the regional climate model were not due to the finer resolution of the regional climate model but to the different treatment of the physical processes in the two models, especially when the subgrid scale physics was important, like during summer. Compared to the coarse resolution RSM simulation results, on the other hand, the fine resolution RSM simulations did show improved simulation skills especially when a good boundary condition such as the reanalysis was used to drive the RSM. Under increased CO₂, the driving PCM and downscaled RSM simulations exhibitedwarming over all vertical layers and all regions. Both the RSM and PCM had increased precipitation during the winter, but during the summer, the PCM simulation had an overall precipitation increase mainly due to increased subgrid scale convective activity, whereas the RSM simulations exhibited precipitation decreases and the resulting RSM soil moisture became dryer, especially in the U.S. Southwest. Most of differences in the simulated climate change signals were produced by the distinct model physics rather than by differences in grid resolution.     

Uncertainties in simulating regional climate of Southern Africa: sensitivity to physical parameterizations using WRF

This study aims at quantifying seasonal biases of regional climate model outputs during southern African summer, against a dense in situ measurement network (daily rain-gauge and surface air temperature records, and 12?h UTC radiosondes), and uncertainties associated with some physical parameterizations. Using the non-hydrostatic Advanced Research Weather Forecast (WRF) laterally forced by ERA40 reanalysis, twenty-seven experiments configured with three schemes of cumulus (CU), planetary boundary layer (PBL) and microphysics (MP), are performed at 35?km horizontal resolution during the core of a summer rainy season (December 1993 to February 1994 season) representative of the South African rainfall climatology. WRF simulates accurately seasonal large-scale rainfall patterns, as well as seasonal gradients of South African rainfall and 2-m temperature, and seasonal vertical profiles of the air temperature and humidity. However seasonal biases fluctuate strongly from an experiment to another, denoting considerable uncertainties generated by the physical package. Rainfall amounts are the most sensitive parameter to the tested schemes. Their geography, intensity, and intraseasonal characteristics are predominantly sensitive to CU schemes, and much less to PBL and MP schemes. Some CU-PBL combinations produce additive effects, which can dramatically either reduce or increase biases. Satisfactory configurations are found for South African climate, which would not have been possible without testing numerous physical parameterizations.     

Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

Numerical Simulation Experiment of Land Surface Physical Processes and Local Climate Effect in Forest Underlying Surface

Based on the basic principles of atmospheric boundary layer and plant canopy micrometeorology, a forest underlying surface land surface physical process model and a two-dimensional atmospheric boundary layer numerical model are developed and numerical simulation experiments of biosphere and physiological processes of vegetation and soil volumetric water content have been done on land surface processes with local climate effect. The numerical simulation results are in good agreement with realistic observations, which can be used to obtain reasonable simulations for diurnal variations of canopy temperature, air temperature in canopy, ground surface temperature, and temporal and spatial distributions of potential temperature and vertical wind velocity as well as relative humidity and turbulence exchange coefficient over non-homogeneous underlying surfaces. It indicates that the model developed can be used to study the interaction between land surface process and atmospheric boundary layer over various underlying surfaces and can be extended to local climate studies. This work will settle a solid foundation for coupling climate models with the biosphere.     

Reduction of systematic biases in regional climate downscaling through ensemble forcing

Simulations of the East Asian summer monsoon for the period of 1979–2001 were carried out using the Weather Research and Forecast (WRF) model forced by three reanalysis datasets (NCEP-R2, ERA-40, and JRA-25). The experiments forced by different reanalysis data exhibited remarkable differences, primarily caused by uncertainties in the lateral boundary (LB) moisture fluxes over the Bay of Bengal and the Philippine Sea. The climatological mean water vapor convergence into the model domain computed from ERA-40 was about 24% higher than that from the NCEP-R2 reanalysis. We demonstrate that using the ensemble mean of NCEP-R2, ERA-40, and JRA-25 as LB forcing considerably reduced the biases in the model simulation. The use of ensemble forcing improved the performance in simulated mean circulation and precipitation, inter-annual variation in seasonal precipitation, and daily precipitation. The model simulated precipitation was superior to that in the reanalysis in both climatology and year-to-year variations, indicating the added value of dynamic downscaling. The results suggest that models having better performance under one set of LB forcing might worsen when another set of reanalysis data is used as LB forcing. Use of ensemble mean LB forcing for assessing regional climate model performance is recommended.     

Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.     

The impact of vertical resolution upon GCM simulations of marine stratocumulus

 The impact of increased vertical resolution in the Hadley Centre Climate Model upon the simulation of stratocumulus is investigated in experiments using single column (SCM) and general circulation (GCM) model configurations. A threefold enhancement of vertical resolution in the boundary layer leads to improvements in the vertical structure of the cloud-topped boundary layer produced by the SCM and GCM in both well-mixed and decoupled situations. However, single and decoupled mixed layers in the marine stratocumulus subsidence regions are still too shallow and, despite increasing, layer cloud amounts remain generally too low. Moreover, closer examination of GCM data and SCM timeseries reveals an underlying sensitivity to vertical resolution in model interactions between boundary layer and convection processes which appears unrealistic. Stratocumulus simulation is thus unlikely to improve significantly as a result of enhanced resolution alone and further work is being undertaken to improve the Hadley Centre model’s boundary layer scheme and, in particular, its interaction with the convection scheme. Nevertheless, this study shows that the full benefit of an improved boundary layer scheme will not be realized if the boundary layer structure is constrained by the rather poor lower troposphere resolution of the standard 19-level climate model. Future Hadley Centre model versions will seek to combine the added flexibility of a better resolved structure with improvements to the subgrid boundary layer parametrizations. Received: 14 April 1998 / Accepted: 5 November 1998     

High resolution re-analysis for the Baltic Sea region during 1965–2005 period

Regional reanalysis database BaltAn65+ comprising meteorological data for Baltic Sea region for the time period 1965–2005 is described. For data assimilation and hindcasts, the numerical weather prediction model HIRLAM 7.1.4 is applied, with 11 km horizontal and 60-layer vertical resolution. Reanalysis includes three-dimensional weather analysis data. Standard surface observations and meteorological soundings together with ship and buoy measurements from WMO observational network are used in analysis. Boundary fields are obtained from ECMWF ERA-40 global re-analysis. The BaltAn65+ can be considered as a regional refinement of ERA-40 for Baltic Sea region, providing the historical weather and climate data with enhanced spatial resolution, which is main motivation for creation of this novel reanalysis database.     

Mixed-phase clouds cause climate model biases in Arctic wintertime temperature inversions

Temperature inversions are a common feature of the Arctic wintertime boundary layer. They have important impacts on both radiative and turbulent heat fluxes and partly determine local climate-change feedbacks. Understanding the spread in inversion strength modelled by current global climate models is therefore an important step in better understanding Arctic climate and its present and future changes. Here, we show how the formation of Arctic air masses leads to the emergence of a cloudy and a clear state of the Arctic winter boundary layer. In the cloudy state, cloud liquid water is present, little to no surface radiative cooling occurs and inversions are elevated and relatively weak, whereas surface radiative cooling leads to strong surface-based temperature inversions in the clear state. Comparing model output to observations, we find that most climate models lack a realistic representation of the cloudy state. An idealised single-column model experiment of the formation of Arctic air reveals that this bias is linked to inadequate mixed-phase cloud microphysics, whereas turbulent and conductive heat fluxes control the strength of inversions within the clear state.     

Oceanic climatology in the coupled model FGOALS-g2: Improvements and biases

The present study examines simulated oceanic climatology in the Flexible Global Ocean-Atmosphere-Land System model, Grid-point Version 2 (FGOALS-g2) forced by historical external forcing data. The oceanic temperatures and circulations in FGOALS-g2 were found to be comparable to those observed, and substantially improved compared to those simulated by the previous version, FGOALS-g1.0. Compared with simulations by FGOALS-g1.0, the shallow mixed layer depths were better captured in the eastern Atlantic and Pacific Ocean in FGOALS-g2. In the high latitudes of the Northern Hemisphere, the cold biases of SST were about 1°C–5°C smaller in FGOALS-g2. The associated sea ice distributions and their seasonal cycles were more realistic in FGOALS-g2. The pattern of Atlantic Meridional Overturning Circulation (AMOC) was better simulated in FGOALS-g2, although its magnitude was larger than that found in observed data. The simulated Antarctic Circumpolar Current (ACC) transport was about 140 Sv through the Drake Passage, which is close to that observed. Moreover, Antarctic Intermediate Water (AAIW) was better captured in FGOALS-g2. However, large SST cold biases (>3°C) were still found to exist around major western boundary currents and in the Barents Sea, which can be explained by excessively strong oceanic cold advection and unresolved processes owing to the coarse resolution. In the Indo-Pacific warm pool, the cold biases were partly related to the excessive loss of heat from the ocean. Along the eastern coast in the Atlantic and Pacific Oceans, the warm biases were due to overestimation of shortwave radiation. In the Indian Ocean and Southern Ocean, the surface fresh biases were mainly due to the biases of precipitation. In the tropical Pacific Ocean, the surface fresh biases (>2 psu) were mainly caused by excessive precipitation and oceanic advection. In the Indo-Pacific Ocean, fresh biases were also found to dominate in the upper 1000 m, except in the northeastern Indian Ocean. There were warm and salty biases (3°C–4°C and 1–2 psu) from the surface to the bottom in the Labrador Sea, which might be due to large amounts of heat transport and excessive evaporation, respectively. For vertical structures, the maximal biases of temperature and salinity were found to be located at depths of >600 m in the Arctic Ocean, and their values exceeded 4°C and 2 psu, respectively.     

A Comparative Study and Evaluation of Mixing-Height Estimation Based on Sodar-RASS, Ceilometer Data and Numerical Model Simulations

The unmanned aerial system SUMO (Small Unmanned Meteorological Observer) has been used for the observation of the structure and behaviour of the atmospheric boundary layer above the Advent Valley, Svalbard during a two-week period in early spring 2009. Temperature, humidity and wind profiles measured by the SUMO system have been compared with measurements of a small tethered balloon system that was operated simultaneously. It is shown that both systems complement each other. Above 200?m, the SUMO system outperforms the tethered balloon in terms of flexibility and the ability to penetrate strong inversion layers of the Arctic boundary layer. Below that level, the tethered balloon system provides atmospheric profiles with higher accuracy, mainly due to its ability to operate at very low vertical velocities. For the observational period, a numerical mesoscale model has been run at high resolution and evaluated with SUMO profiles reaching up to a height of 1500?m above the ground. The sensitivity to the choice of atmospheric boundary-layer schemes and horizontal resolution has been investigated. A new scheme especially suited for stable conditions slightly improves the temperature forecast in stable conditions, although all schemes show a warm bias close to the surface and a cold bias above the atmospheric boundary layer. During one cold and cloudless night, the SUMO system could be operated nearly continuously (every 30?C45?minutes). This allowed for a detailed case study addressing the structure and behaviour of the air column within and above Advent Valley and its interaction with the local topography. The SUMO measurements in conjunction with a 10-m meteorological mast enabled the identification of a very stable nocturnal surface layer adjacent to the valley bottom, a stable air column in the valley and a strong inversion layer above the summit height. The results indicate the presence of inertial-gravity waves during the night, a feature not captured by the model.     

Modelling of Water Surface Temperature of Three Lakes on the Tibetan Plateau using a Physically Based Lake Model

ABSTRACT

This study was carried out to quantify the physical processes of lakes in the Tibetan Plateau using the Community Land Model, version 4.5 (CLM4.5), coupled with a physically based, 10-layer lake model developed by the National Center for Atmospheric Research. The CLM was forced with 10?km resolution reanalysis data to attempt to understand detailed lake processes and how these processes affect lake surface temperature. In this study, we simulated seasonal and interannual variations of lake surface temperature for Lake Qinghai, Zhaling Co, and Nam Co in the Tibetan Plateau and compared these simulations with observations. The results showed that the CLM4.5 lake model simulations reproduced the observed lake surface temperatures for Lake Qinghai and Zhaling Co well but reproduced those for Nam Co poorly. Through detailed analysis, we found that the simulated biases for Nam Co result largely from the unrealistic parameterization of eddy diffusivity. By expanding this parameter, the lake surface temperature simulations improved remarkably. In addition, erroneous lake ice cover simulations contributed to the simulated lake surface temperature bias in the cold seasons.