Abstract A two-dimensional version of the Pennsylvania State University mesoscale model has been applied to Winter Monsoon Experiment data in order to simulate the diurnally occurring convection observed over the South China Sea. The domain includes a representation of part of Borneo as well as the sea so that the model can simulate the initiation of convection. Also included in the model are parameterizations of mesoscale ice phase and moisture processes and longwave and shortwave radiation with a diurnal cycle. This allows use of the model to test the relative importance of various heating mechanisms to the stratiform cloud deck, which typically occupies several hundred kilometers of the domain. Frank and Cohen's cumulus parameterization scheme is employed to represent vital unresolved vertical transports in the convective area. The major conclusions are: Ice phase processes are important in determining the level of maximum large-scale heating and vertical motion because there is a strong anvil componen...

Numerical study of convection observed during the Winter Monsoon Experiment using a mesoscale two-dimensional model [presentation]

The Twentieth Century Reanalysis (20CR) project is an international effort to produce a comprehensive global atmospheric circulation dataset spanning the twentieth century, assimilating only surface pressure reports and using observed monthly sea-surface temperature and sea-ice distributions as boundary conditions. It is chiefly motivated by a need to provide an observational dataset with quantified uncertainties for validations of climate model simulations of the twentieth century on all time-scales, with emphasis on the statistics of daily weather. It uses an Ensemble Kalman Filter data assimilation method with background ‘first guess’ fields supplied by an ensemble of forecasts from a global numerical weather prediction model. This directly yields a global analysis every 6 hours as the most likely state of the atmosphere, and also an uncertainty estimate of that analysis.

The 20CR dataset provides the first estimates of global tropospheric variability, and of the dataset's time-varying quality, from 1871 to the present at 6-hourly temporal and 2° spatial resolutions. Intercomparisons with independent radiosonde data indicate that the reanalyses are generally of high quality. The quality in the extratropical Northern Hemisphere throughout the century is similar to that of current three-day operational NWP forecasts. Intercomparisons over the second half-century of these surface-based reanalyses with other reanalyses that also make use of upper-air and satellite data are equally encouraging.

It is anticipated that the 20CR dataset will be a valuable resource to the climate research community for both model validations and diagnostic studies. Some surprising results are already evident. For instance, the long-term trends of indices representing the North Atlantic Oscillation, the tropical Pacific Walker Circulation, and the Pacific–North American pattern are weak or non-existent over the full period of record. The long-term trends of zonally averaged precipitation minus evaporation also differ in character from those in climate model simulations of the twentieth century. Copyright © 2011 Royal Meteorological Society and Crown Copyright.

/pdf/the-twentieth-century-reanalysis-project-5fx9pcs9x0.pdf

The Twentieth Century Reanalysis Project

The purpose of this paper is to provide a comprehensive presentation and interpretation of the Ensemble Kalman Filter (EnKF) and its numerical implementation. The EnKF has a large user group, and numerous publications have discussed applications and theoretical aspects of it. This paper reviews the important results from these studies and also presents new ideas and alternative interpretations which further explain the success of the EnKF. In addition to providing the theoretical framework needed for using the EnKF, there is also a focus on the algorithmic formulation and optimal numerical implementation. A program listing is given for some of the key subroutines. The paper also touches upon specific issues such as the use of nonlinear measurements, in situ profiles of temperature and salinity, and data which are available with high frequency in time. An ensemble based optimal interpolation (EnOI) scheme is presented as a cost-effective approach which may serve as an alternative to the EnKF in some applications. A fairly extensive discussion is devoted to the use of time correlated model errors and the estimation of model bias.

The Ensemble Kalman Filter: Theoretical formulation and practical implementation

High-resolution information on climatic conditions is essential to many applications in environmental and ecological sciences. Here we present the CHELSA (Climatologies at high resolution for the earth’s land surface areas) data of downscaled model output temperature and precipitation estimates of the ERA-Interim climatic reanalysis to a high resolution of 30 arc sec. The temperature algorithm is based on statistical downscaling of atmospheric temperatures. The precipitation algorithm incorporates orographic predictors including wind fields, valley exposition, and boundary layer height, with a subsequent bias correction. The resulting data consist of a monthly temperature and precipitation climatology for the years 1979–2013. We compare the data derived from the CHELSA algorithm with other standard gridded products and station data from the Global Historical Climate Network. We compare the performance of the new climatologies in species distribution modelling and show that we can increase the accuracy of species range predictions. We further show that CHELSA climatological data has a similar accuracy as other products for temperature, but that its predictions of precipitation patterns are better. Machine-accessible metadata file describing the reported data (ISA-Tab format)

/pdf/climatologies-at-high-resolution-for-the-earth-s-land-23ipmy16re.pdf

Climatologies at high resolution for the earth’s land surface areas

/pdf/efficient-data-assimilation-for-spatiotemporal-chaos-a-local-zlp3sn8d6w.pdf

Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter

Accurate depiction of meteorological conditions, especially within the planetary boundary layer (PBL), is important for air pollution modeling, and PBL parameterization schemes play a critical role in simulating the boundary layer. This study examines the sensitivity of the performance of the Weather Research and Forecast (WRF) model to the use of three different PBL schemes [Mellor–Yamada–Janjic (MYJ), Yonsei University (YSU), and the asymmetric convective model, version 2 (ACM2)]. Comparison of surface and boundary layer observations with 92 sets of daily, 36-h high-resolution WRF simulations with different schemes over Texas in July–September 2005 shows that the simulations with the YSU and ACM2 schemes give much less bias than with the MYJ scheme. Simulations with the MYJ scheme, the only local closure scheme of the three, produced the coldest and moistest biases in the PBL. The differences among the schemes are found to be due predominantly to differences in vertical mixing strength and entr...

/pdf/evaluation-of-three-planetary-boundary-layer-schemes-in-the-4q46nj8em4.pdf

Evaluation of Three Planetary Boundary Layer Schemes in the WRF Model

The ensemble Kalman filter (EnKF) uses an ensemble of short-range forecasts to estimate the flow-dependent background error covariances required in data assimilation. The feasibility of the EnKF for convective-scale data assimilation has been previously demonstrated in perfect-model experiments using simulated observations of radial velocity from a supercell storm. The present study further explores the potential and behavior of the EnKF at convective scales by considering more realistic initial analyses and variations in the availability and quality of the radar observations. Assimilation of simulated radial-velocity observations every 5 min where there is significant reflectivity using 20 ensemble members proves to be successful in most realistic observational scenarios for simulated supercell thunderstorms, although the same degree of success may not be readily expected with real observations and an imperfect model, at least with the present EnKF implementation. Even though the filter converges toward the truth simulation faster from a better initial estimate, an experiment with the initial estimate of the supercell displaced by 10 km still yields an accurate estimate of the storm for both observed and unobserved variables within 40 min. Similarly, radial-velocity observations below 2 km are certainly beneficial to capturing the storm (especially the detailed cold pool structure), but in their absence the assimilation scheme can still achieve a comparably accurate estimate of the state of the storm given a slightly longer assimilation period. An experiment with radar observations only above 4 km fails to assimilate the storm properly, but, with the addition of a hypothetical surface mesonet taking wind and temperature observations, the EnKF can again provide a good estimate of the storm. The supercell can also be successfully assimilated in the case of radar observations only below 4 km (such as those from the ground-based mobile radars). More frequent observations can help the storm assimilation initially, but the benefit diminishes after half an hour. Results presented here indicate that the vertical resolution and the uncertainty of observations, for the typical range of most of the observational radars, would have little impact on the overall performance of the EnKF in assimilating the storm.

/pdf/impacts-of-initial-estimate-and-observation-availability-on-3zg6vika7f.pdf

Impacts of Initial Estimate and Observation Availability on Convective-Scale Data Assimilation with an Ensemble Kalman Filter

Assimilation of Doppler radar data into cloud models is an important obstacle to routine numerical weather prediction for convective-scale motions; the difficulty lies in initializing fields of wind, temperature, moisture, and condensate given only observations of radial velocity and reflectivity from the radar. This paper investigates the potential of the ensemble Kalman filter (EnKF), which estimates the covariances between observed variables and the state through an ensemble of forecasts, to assimilate radar observations at convective scales. In the basic experiment, simulated observations are extracted from a reference simulation of a splitting supercell and assimilated using the EnKF and the same numerical model that produced the reference simulation. The EnKF produces accurate analyses, including the unobserved variables, after roughly 30 min (or six scans) of radial velocity observations. Additional experiments, in which forecasts are made from the ensemble-mean analysis, reveal that forecast errors grow significantly in this simple system, so that the ability of the EnKF to track the reference solution is not simply because of stable system dynamics. It is also found that the covariances between radial velocity and temperature, moisture, and condensate are important to the quality of the analyses, as is the initialization chosen for the ensemble members prior to assimilating the first observations. These results are promising, especially given the ease of implementing the EnKF. A number of important issues remain, however, including the initialization of the ensemble prior to the first observation, the treatment of uncertainty in the environmental sounding, the role of error in the forecast model (particularly the microphysical parameterizations), and the treatment of lateral boundary conditions.

/pdf/assimilation-of-simulated-doppler-radar-observations-with-an-2z6t0m2945.pdf

Assimilation of Simulated Doppler Radar Observations with an Ensemble Kalman Filter

This paper reviews the development of the ensemble Kalman filter (EnKF) for atmospheric data assimilation. Particular attention is devoted to recent advances and current challenges. The distinguishing properties of three well-established variations of the EnKF algorithm are first discussed. Given the limited size of the ensemble and the unavoidable existence of errors whose origin is unknown (i.e., system error), various approaches to localizing the impact of observations and to accounting for these errors have been proposed. However, challenges remain; for example, with regard to localization of multiscale phenomena (both in time and space). For the EnKF in general, but higher-resolution applications in particular, it is desirable to use a short assimilation window. This motivates a focus on approaches for maintaining balance during the EnKF update. Also discussed are limited-area EnKF systems, in particular with regard to the assimilation of radar data and applications to tracking severe storms ...

Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation

For several decades, jets and fronts have been known from observations to be significant sources of internal gravity waves in the atmosphere. Motivations to investigate these waves have included their impact on tropospheric convection, their contribution to local mixing and turbulence in the upper troposphere, their vertical propagation into the middle atmosphere, and the forcing of its global circulation. While many different studies have consistently highlighted jet exit regions as a favored locus for intense gravity waves, the mechanisms responsible for their emission had long remained elusive: one reason is the complexity of the environment in which the waves appear; another reason is that the waves constitute small deviations from the balanced dynamics of the flow generating them; i.e., they arise beyond our fundamental understanding of jets and fronts based on approximations that filter out gravity waves. Over the past two decades, the pressing need for improving parameterizations of nonorographic gravity waves in climate models that include a stratosphere has stimulated renewed investigations. The purpose of this review is to present current knowledge and understanding on gravity waves near jets and fronts from observations, theory, and modeling, and to discuss challenges for progress in coming years.

/pdf/internal-gravity-waves-from-atmospheric-jets-and-fronts-22oa2c4omk.pdf

Fuqing Zhang

Papers

Evaluation of Three Planetary Boundary Layer Schemes in the WRF Model

Impacts of Initial Estimate and Observation Availability on Convective-Scale Data Assimilation with an Ensemble Kalman Filter

Assimilation of Simulated Doppler Radar Observations with an Ensemble Kalman Filter

Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation

Internal gravity waves from atmospheric jets and fronts