The Tropospheric Humidity Trends of NCEP/NCAR Reanalysis Before the Satellite Era (original) (raw)

A Climatology of Midlatitude Continental Clouds from the ARM SGP Central Facility. Part II: Cloud Fraction and Surface Radiative Forcing

Journal of Climate, 2006

Data collected at the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Central Facility (SCF) are analyzed to determine the monthly and hourly variations of cloud fraction and radiative forcing between January 1997 and December 2002. Cloud fractions are estimated for total cloud cover and for single-layered low (0-3 km), middle (3-6 km), and high clouds (Ͼ6 km) using ARM SCF ground-based paired lidar-radar measurements. Shortwave (SW) and longwave (LW) fluxes are derived from up-and down-looking standard precision spectral pyranometers and precision infrared radiometer measurements with uncertainties of ϳ10 W m Ϫ2. The annual averages of total and single-layered low-, middle-, and high-cloud fractions are 0.49, 0.11, 0.03, and 0.17, respectively. Both totaland low-cloud amounts peak during January and February and reach a minimum during July and August; high clouds occur more frequently than other types of clouds with a peak in summer. The average annual downwelling surface SW fluxes for total and low clouds (151 and 138 W m Ϫ2 , respectively) are less than those under middle and high clouds (188 and 201 W m Ϫ2 , respectively), but the downwelling LW fluxes (349 and 356 W m Ϫ2) underneath total and low clouds are greater than those from middle and high clouds (337 and 333 W m Ϫ2). Low clouds produce the largest LW warming (55 W m Ϫ2) and SW cooling (Ϫ91 W m Ϫ2) effects with maximum and minimum absolute values in spring and summer, respectively. High clouds have the smallest LW warming (17 W m Ϫ2) and SW cooling (Ϫ37 W m Ϫ2) effects at the surface. All-sky SW cloud radiative forcing (CRF) decreases and LW CRF increases with increasing cloud fraction with mean slopes of Ϫ0.984 and 0.616 W m Ϫ2 % Ϫ1 , respectively. Over the entire diurnal cycle, clouds deplete the amount of surface insolation more than they add to the downwelling LW flux. The calculated CRFs do not appear to be significantly affected by uncertainties in data sampling and clear-sky screening. Traditionally, cloud radiative forcing includes not only the radiative impact of the hydrometeors, but also the changes in the environment. Taken together over the ARM SCF, changes in humidity and surface albedo between clear and cloudy conditions offset ϳ20% of the NET radiative forcing caused by the cloud hydrometeors alone. Variations in water vapor, on average, account for 10% and 83% of the SW and LW CRFs, respectively, in total cloud cover conditions. The error analysis further reveals that the cloud hydrometeors dominate the SW CRF, while water vapor changes are most important for LW flux changes in cloudy skies. Similar studies over other locales are encouraged where water and surface albedo changes from clear to cloudy conditions may be much different than observed over the ARM SCF.

The relationship between tropospheric wave forcing and tropical lower stratospheric water vapor

Atmospheric Chemistry and Physics, 2008

Using water vapor data from HALOE and SAGE II, an anti-correlation between planetary wave driving (here expressed by the mid-latitude eddy heat flux at 50 hPa added from both hemispheres) and tropical lower stratospheric (TLS) water vapor has been obtained. This appears to be a manifestation of the inter-annual variability of the Brewer-Dobson (BD) circulation strength (the driving of which is generally measured in terms of the mid-latitude eddy heat flux), and hence amount of water vapor entering the stratosphere. Some years such as 1991 and 1997 show, however, a clear departure from the anti-correlation which suggests that the water vapor changes in TLS can not be attributed solely to changes in extratropical planetary wave activity (and its effect on the BD circulation). After 2000 a sudden decrease in lower stratospheric water vapor has been reported in earlier studies based upon satellite data from HALOE, SAGE II and POAM III indicating that the lower stratosphere has become drier since then. This is consistent with a sudden rise in the combined mid-latitude eddy heat flux with nearly equal contribution from both hemispheres as shown here and with the increase in tropical upwelling and decrease in cold point temperatures found by . The low water vapor and enhanced planetary wave activity (in turn strength of the BD circulation) has persisted until the end of the satellite data records. From a multi-variate regression analysis applied to 27 years of NCEP and HadAT2 (radiosonde) temperatures (up to 2005) with contributions from solar cycle, stratospheric aerosols and QBO removed, the enhancement wave driving after 2000 is estimated to contribute up to 0.7 K cooling to the overall TLS temperature change during the period 2001-2005 when compared to the period 1996-2000. NCEP cold point temperature show an average decrease of nearly 0.4 K from changes in the wave driving, which is consistent with observed mean TLS water vapor changes of about −0.2 ppm after 2000.

Relationships Between Cloud Radiative Forcing and Local Meteorological Variables in Several Global Climate Models

Journal of Climate, 2006

The relations between local monthly mean shortwave cloud radiative forcing and aspects of the resolvedscale meteorological fields are investigated in hindcast simulations performed with 12 of the global coupled models included in the model intercomparison conducted as part of the preparation for Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). In particular, the connection of the cloud forcing over tropical and subtropical ocean areas with resolved midtropospheric vertical velocity and with lower-level relative humidity are investigated and compared among the models. The model results are also compared with observational determinations of the same relationships using satellite data for the cloud forcing and global reanalysis products for the vertical velocity and humidity fields. In the analysis the geographical variability in the long-term mean among all grid points and the interannual variability of the monthly mean at each grid point are considered separately. The shortwave cloud radiative feedback (SWCRF) plays a crucial role in determining the predicted response to large-scale climate forcing (such as from increased greenhouse gas concentrations), and it is thus important to test how the cloud representations in current climate models respond to unforced variability.

On the decadal increase in the tropical mean outgoing longwave radiation for the period 1984-2000

Atmospheric Chemistry and Physics, 2004

In the present paper, we have calculated the outgoing longwave radiation at the top of the atmosphere (OLR at TOA) using a deterministic radiation transfer model, cloud data from ISCCP-D, and atmospheric temperature and humidity data from NCEP/NCAR reanalysis, for the seventeen-year period 1984-2000. We constructed anomaly timeseries of the OLR at TOA, as well as of all of the key input climatological data, averaged in the tropical region between 20 • N and 20 • S. We compared the anomaly time-series of the model calculated OLR at TOA with that obtained from the ERBE S-10N (WFOV NF edition 2) non-scanner measurements. The model results display very similar seasonal and inter-annual variability as the ERBS data, and indicate a decadal increase of OLR at TOA of 1.9±0.2 Wm −2 /decade, which is lower than that displayed by the ERBS time-series (3.5±0.3 Wm −2). Analysis of the inter-annual and long-term variability of the various parameters determining the OLR at TOA, showed that the most important contribution to the observed trend comes from a decrease in high-level cloud cover over the period 1984-2000, followed by an apparent drying of the upper troposphere and a decrease in low-level cloudiness. Opposite but small trends are introduced by a decrease in low-level cloud top pressure, an apparent cooling of the lower stratosphere (at the 50 mbar level) and a small decadal increase in mid-level cloud cover.

Evaluation of upper tropospheric humidity forecasts from ECMWF using AIRS and CALIPSO data

Atmospheric Chemistry and Physics Discussions, 2008

An evaluation of the upper tropospheric humidity from the European Centre of Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) is presented. We first make an analysis of the spinup behaviour of ice supersaturation in weather forecasts. It shows that a spinup period of at least 12 h is necessary before using forecast humidity data from the upper troposphere. We compare the forecasted upper tropospheric humidity with coincident relative humidity fields retrieved from the Atmospheric InfraRed Sounder (AIRS) and with cloud vertical profiles from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). The analysis is made over one year, from October 2006 to September 2007, and we discuss how relative humidity and cloud features appear both in the IFS and in the observations. In a last part, we investigate the presence of ice supersaturation within low vertical resolution pressure layers by comparing the IFS outputs for high-resolution and low-resolution humidity profiles and by simulating the interpolation of humidity over radiosonde data. A new correction method is proposed and tested with these data.

Model-simulated humidity bias in the upper troposphere and its relation to the large-scale circulation

Journal of Geophysical Research, 2011

1] The depiction of water vapor in the upper troposphere in Geophysical Fluid Dynamics Laboratory (GFDL) climate model and ERA-40 reanalysis is evaluated through a model-to-radiance approach. Brightness temperatures of High-Resolution Infrared Radiation Sounder (HIRS) 6.7 mm channel, Special Sensor for Microwave Water Vapor Profiler (SSM/T-2) 183.31 ± 1 GHz channel, and Microwave Sounding Unit (MSU) 60 GHz channel simulated with data from the GFDL climate model and ERA-40 reanalysis show a distinct cold and moist bias in the upper troposphere compared to satellite observations, particularly over the subtropics. Temperature biases are a common feature in many climate models and complicate the interpretation of radiance-based comparisons with satellite data. We introduce a new method for evaluating the water vapor distribution which combines both HIRS 6.7 mm and SSM/T-2 183.31 ± 1 GHz channels and is much less sensitive to tropospheric temperature biases. Using this method, we show that GFDL climate model has a more humid upper troposphere over dry subtropical area than ERA-40 reanalysis. The geographical distribution of the humidity bias is found to exhibit a close association with differences in the 500 hPa vertical pressure velocity, suggesting that much of the bias in tropical upper tropospheric relative humidity can be attributed to errors in simulating the intensity of large-scale tropical circulation. Given the strong dependence of upper tropospheric water vapor on the large-scale circulation, these results suggest that long-term monitoring of upper tropospheric water vapor from satellites may also offer insight into variations in the large-scale atmospheric circulation. Citation: Chung, E.-S., B. J. Soden, B.-J. Sohn, and J. Schmetz (2011), Model-simulated humidity bias in the upper troposphere and its relation to the large-scale circulation,

The surface downward longwave radiation in the ECMWF forecast system (2002)

The surface downward longwave radiation, computed by the European Centre for Medium-Range Weather Forecasts (ECMWF) forecast system used for the ECMWF 40-yr reanalysis, is compared with surface radiation measurements for the April–May 1999 period, available as part of the Baseline Surface Radiation Network (BSRN), Surface Radiation (SURFRAD), and Atmospheric Radiation Measurement (ARM) programs. Emphasis is put on comparisons on a 1-h basis, as this allows discrepancies to be more easily linked to differences between model description and observations of temperature, humidity, and clouds. It is also possible to compare the model and observed temporal variability in the surface radiation fluxes. Comparisons are first carried out at locations for which the spectral model orography differs from the actual station height. Sensitivity of the model fluxes to various algorithms to correct for this discrepancy is explored. A simple interpolation–extrapolation scheme for pressure, temperature, and specific humidity allows the improvement of model calculations of the longwave surface fluxes in most cases. Intercomparisons of surface longwave radiation are presented for the various longwave radiation schemes operational since the 15-yr ECMWF Re-Analysis (ERA-15) was performed. The Rapid Radiation Transfer Model of Mlawer et al., now operational at ECMWF, is shown to correct for the major underestimation in clear-sky downward longwave radiation seen in ERA-15. Sensitivity calculations are also carried out to explore the role of the cloud optical properties, cloud effective particle size, and aerosols in the representation of the surface downward longwave radiation

Analysis of the decrease in the tropical mean outgoing shortwave radiation at the top of atmosphere for the period 1984-2000

Atmospheric Chemistry and Physics, 2005

A decadal-scale trend in the tropical radiative energy budget has been observed recently by satellites, which however is not reproduced by climate models. In the present study, we have computed the outgoing shortwave radiation (OSR) at the top of atmosphere (TOA) at 2.5 • longitudelatitude resolution and on a mean monthly basis for the 17year period 1984-2000, by using a deterministic solar radiative transfer model and cloud climatological data from the International Satellite Cloud Climatology Project (IS-CCP) D2 database. Anomaly time series for the mean monthly pixel-level OSR fluxes, as well as for the key physical parameters, were constructed. A significant decreasing trend in OSR anomalies, starting mainly from the late 1980s, was found in tropical and subtropical regions (30 • S-30 • N), indicating a decadal increase in solar planetary heating equal to 1.9±0.3 Wm −2 /decade, reproducing well the features recorded by satellite observations, in contrast to climate model results. This increase in solar planetary heating, however, is accompanied by a similar increase in planetary cooling, due to increased outgoing longwave radiation, so that there is no change in net radiation. The model computed OSR trend is in good agreement with the corresponding linear decadal decrease of 2.5±0.4 Wm −2 /decade in tropical mean OSR anomalies derived from ERBE S-10N nonscanner data (edition 2). An attempt was made to identify the physical processes responsible for the decreasing trend in tropical mean OSR. A detailed correlation analysis using pixel-level anomalies of model computed OSR flux and ISCCP cloud cover over the entire tropical and subtropical region (30 • S-30 • N), gave a correlation coefficient of 0.79, indicating that decreasing cloud cover is the main reason for the tropical OSR trend. According to the ISCCP-D2 data

Determination and significance of upper-tropospheric humidity

Atmospheric Chemistry and Physics Discussions, 2018

We present a novel retrieval for upper-tropospheric humidity (UTH) from HIRS channel 12 radiances that successfully bridges the wavelength change from 6.7 to 6.5 µm that occurred from HIRS 2 on NOAA 14 to HIRS 3 on NOAA 15. The jump in average brightness temperature (T 12) that this change caused (about −7 K) could be fixed with a statistical intercalibration method (Shi and Bates, 2011). Unfortunately, the retrieval of UTHi based on the intercalibrated data was not satisfying at the high tail of the distribution of UTHi. Attempts to construct a better intercalibration in the low T 12 range (equivalent to the high UTHi range) were either not successful (Gierens et al., 2018) or required additional statistically determined corrections to the measured brightness temperatures (Gierens and Eleftheratos, 2017). The new method presented here is based on the original one (Soden and Bretherton, 1993; Stephens et al., 1996; Jackson and Bates, 2001), but it extends linearisations in the formulation of water vapour saturation pressure and in the temperature-dependence of the Planck function to second order. To achieve the second-order formulation we derive the retrieval from the beginning, and we find that the most influential ingredient is the use of different optical constants for the two involved channel wavelengths (6.7 and 6.5 µm). The result of adapting the optical constant is an almost perfect match between UTH data measured by HIRS 2 on NOAA 14 and HIRS 3 on NOAA 15 on 1004 common days of operation. The method is applied to both UTH and UTHi, the upper-tropospheric humidity with respect to ice. For each case retrieval coefficients are derived. We present a number of test applications, e.g. on computed brightness temperatures based on high-resolution radiosonde profiles, on the brightness temperatures measured by the satellites on the mentioned 1004 common days of operation. Further we present time series of the occurrence frequency of high UTHi cases and we show the overall probability distribution of UTHi. The two latter applications expose clear indications of moistening of the upper troposphere over the last 35 years. Finally, we discuss the significance of UTH. We state that UTH algorithms cannot be judged for their correctness or incorrectness, since there is no true UTH. Instead, UTH algorithms should fulfill a number of usefulness-postulates, that we suggest and discuss. In the course of this discussion an alternative method to estimate the weighting function is presented.