scholarly journals The Sensitivity of Rainfall Estimation to Error Assumptions in a Bayesian Passive Microwave Retrieval Algorithm

2015 ◽  
Vol 54 (2) ◽  
pp. 408-422 ◽  
Author(s):  
Gregory S. Elsaesser ◽  
Christian D. Kummerow

AbstractThe Goddard profiling algorithm (GPROF) uses Bayesian probability theory to retrieve rainfall over the global oceans. A critical component of GPROF and most Bayes theorem–based retrieval frameworks is the specification of uncertainty in the observations being utilized to retrieve the parameter of interest. In the case of GPROF, for any sensor, uncertainties in microwave brightness temperatures (Tbs) arise from radiative transfer model errors, satellite sensor noise and/or degradation, and nonlinear, scene-dependent Tb offsets added during sensor intercalibration procedures. All mentioned sources impact sensors in a varying fashion, in part because of sensor-dependent fields of view. It is found that small errors in assumed Tb uncertainty (ranging from 0.57 K at 10 GHz to 2.29 K at 85 GHz) lead to a 3.6% change in the retrieved global-average oceanic rainfall rate, and 10%–20% (20%–40%) shifts in the pixel-level (monthly) frequency distributions for given rainfall bins. A mathematical expression describing the sensitivity of retrieved rainfall to uncertainty is developed here. The strong global sensitivity is linked to rainfall variance scaling systematically as Tb varies. For ocean scenes, the same emission-dominated rainfall–Tb physics used in passive microwave rainfall retrieval is also responsible for the substantial underestimation (overestimation) of global rainfall if uncertainty is overestimated (underestimated). Proper uncertainties are required to quantify variability in surface rainfall, assess long-term trends, and provide robust rainfall benchmarks for general circulation model evaluations. The implications for assessing global and regional biases in active versus passive microwave rainfall products, and for achieving rainfall product agreement among a constellation of orbiting microwave radiometers [employed in the Global Precipitation Measurement (GPM) mission], are also discussed.

Sensors ◽  
2020 ◽  
Vol 20 (4) ◽  
pp. 1225
Author(s):  
Lanka Karthikeyan ◽  
Ming Pan ◽  
Dasika Nagesh Kumar ◽  
Eric F. Wood

Passive microwave sensors use a radiative transfer model (RTM) to retrieve soil moisture (SM) using brightness temperatures (TB) at low microwave frequencies. Vegetation optical depth (VOD) is a key input to the RTM. Retrieval algorithms can analytically invert the RTM using dual-polarized TB measurements to retrieve the VOD and SM concurrently. Algorithms in this regard typically use the τ-ω types of models, which consist of two third-order polynomial equations and, thus, can have multiple solutions. Through this work, we find that uncertainty occurs due to the structural indeterminacy that is inherent in all τ-ω types of models in passive microwave SM retrieval algorithms. In the process, a new analytical solution for concurrent VOD and SM retrieval is presented, along with two widely used existing analytical solutions. All three solutions are applied to a fixed framework of RTM to retrieve VOD and SM on a global scale, using X-band Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) TB data. Results indicate that, with structural uncertainty, there ensues a noticeable impact on the VOD and SM retrievals. In an era where the sensitivity of retrieval algorithms is still being researched, we believe the structural indeterminacy of RTM identified here would contribute to uncertainty in the soil moisture retrievals.


2012 ◽  
Vol 12 (1) ◽  
pp. 479-523 ◽  
Author(s):  
B. M. Monge-Sanz ◽  
M. P. Chipperfield ◽  
A. Untch ◽  
J.-J. Morcrette ◽  
A. Rap ◽  
...  

Abstract. A new linear parameterisation for stratospheric methane (CoMeCAT) has been developed and tested. The scheme is derived from a 3-D full chemistry transport model (CTM) and tested within the same chemistry model itself, as well as in an independent general circulation model (GCM). The new CH4/H2O scheme is suitable for any global model and here is shown to provide realistic profiles in the 3-D TOMCAT/SLIMCAT CTM and in the ECMWF (European Centre for Medium-Range Weather Forecasts) GCM. Simulation results from the new stratospheric scheme are in good agreement with the full-chemistry CTM CH4 field and with observations from the Halogen Occultation Experiment (HALOE). The CH4 scheme has also been used to derive a source for stratospheric water. Stratospheric water increments obtained in this way within the CTM produce vertical and latitudinal H2O variation in fair agreement with satellite observations. Stratospheric H2O distributions in the ECMWF GCM present realistic overall features although concentrations are lower than in the CTM run (up to 0.5 ppmv lower above 10 hPa). The potential of the new CoMeCAT scheme for evaluating long-term transport within the ECMWF model is exploited to assess the impacts of nudging the free running GCM to ERA-40 and ERA-Interim reanalyses. In this case, the nudged GCM shows similar transport patterns to the CTM forced by the corresponding reanalysis data, ERA-Interim producing better results than ERA-40. The impact that the new methane description has in the GCM radiation scheme is also explored. Compared to the default CH4 climatology used by the ECMWF model, CoMeCAT produces up to 2 K cooling in the tropical lower stratosphere. The effect of using the CoMeCAT scheme for radiative forcing (RF) calculations has been investigated using the off-line Edwards-Slingo (E-S) radiative transfer model. Compared to the use of a tropospheric global 3-D CH4 value, the CoMeCAT distributions produce an overall decrease in the annual mean net RF, with the largest decrease found over the Southern Hemisphere high latitudes. The effect of the new CH4 stratospheric distribution on these RF calculations is of up to 30 mW m−2, i.e. the same order of magnitude, and opposite sign, as the inclusion of aircraft contrails formation in the radiative model.


2006 ◽  
Vol 19 (16) ◽  
pp. 3792-3809 ◽  
Author(s):  
Martin Wild ◽  
Erich Roeckner

Abstract Radiative fluxes in the ECHAM5 general circulation model (GCM) are evaluated using both surface and satellite-based observations. The fluxes at the top of the atmosphere (TOA) are generally in good agreement with the satellite data. Larger deviations in simulated cloud forcing are found especially at lower latitudes where the shortwave component within the intertropical convergence zone is overestimated during boreal summer and underestimated in the marine stratocumulus regimes, especially during boreal winter. At the surface the biases in the radiative fluxes are significantly smaller than in earlier versions of the same model and in other GCMs. The shortwave clear-sky fluxes are shown to be in good agreement with newly derived observational estimates. Compared to the preceding model version, ECHAM4, the spurious absorption of solar radiation in the cloudy atmosphere disappears due to the higher resolution in the near-infrared bands of the shortwave radiation code. This reduces the biases with respect to collocated surface and TOA observations. It is illustrated that remaining biases in atmospheric absorption may be related to the crude aerosol climatology, which does not account for high loadings of absorbing aerosol such as from biomass burning, whereas the biases disappear in areas and seasons where aerosol effects are less important. In the longwave, the introduction of the Rapid Radiative Transfer Model (RRTM) radiation code leads to an increase in the longwave downward flux at the surface at high latitudes, thereby reducing biases typically found in GCMs. The considerable skill in the simulation of the fluxes at the earth’s surface underlines the suitability of ECHAM5 as an atmospheric component of an integrated earth system model.


2020 ◽  
Author(s):  
Zadvornykh Ilya ◽  
Gribanov Konstantin ◽  
Zakharov Vyacheslav ◽  
Denisova Nina ◽  
Imasu Ryoichi ◽  
...  

<p><span>Global monitoring of isotopic composition of water vapor in the atmosphere provides information regarding atmospheric hydgrological cycle in the Earth’s climate system. The observed HDO to H</span><sub><span>2</span></sub><span>O ratio (δD) in atmospheric water vapor gives information on origin and history of air masses in the atmosphere.</span></p><p><span>In this study we present a method, software tool and some results on δD ratio retrieval from spectra, measured simultaneously in the thermal (TIR) and short wave infrared (SWIR) spectral ranges. The TANSO-FTS high spectral resolution spectrometer on board GOSAT-2 satellite is unique to perform simultaneous measurements in TIR and SWIR spectral bands. A method of simultaneous using of these both bands can improve vertical resolution of δD retrieved profile.</span></p><p><span>We applied conventional optimal estimation method to solve inverse problem. The output data of atmospheric general circulation model ECHAM5-wiso were used as a statistical ensemble of HDO and H2O a priori profiles. The averaging kernels, a posteriori covariance matrices and degrees of freedom are calculated. The retrieval algorithm is implemented using original FIRE-ARMS software and VLIDORT radiation transfer model.</span></p><p><span>This study is supported by the Russian Science Foundation grant No. 18-11-00024.</span></p>


2013 ◽  
Vol 13 (18) ◽  
pp. 9641-9660 ◽  
Author(s):  
B. M. Monge-Sanz ◽  
M. P. Chipperfield ◽  
A. Untch ◽  
J.-J. Morcrette ◽  
A. Rap ◽  
...  

Abstract. This study evaluates effects and applications of a new linear parameterisation for stratospheric methane and water vapour. The new scheme (CoMeCAT) is derived from a 3-D full-chemistry-transport model (CTM). It is suitable for any global model, and is shown here to produce realistic profiles in the TOMCAT/SLIMCAT 3-D CTM and the ECMWF (European Centre for Medium-Range Weather Forecasts) general circulation model (GCM). Results from the new scheme are in good agreement with the full-chemistry CTM CH4 field and with observations from the Halogen Occultation Experiment (HALOE). The scheme is also used to derive stratospheric water increments, which in the CTM produce vertical and latitudinal H2O variations in fair agreement with satellite observations. Stratospheric H2O distributions in the ECMWF GCM show realistic overall features, although concentrations are smaller than in the CTM run (up to 0.5 ppmv smaller above 10 hPa). The potential of the new CoMeCAT tracer for evaluating stratospheric transport is exploited to assess the impacts of nudging the free-running GCM to ERA-40 and ERA-Interim reanalyses. The nudged GCM shows similar transport patterns to the offline CTM forced by the corresponding reanalysis data. The new scheme also impacts radiation and temperature in the model. Compared to the default CH4 climatology and H2O used by the ECMWF radiation scheme, the main effect on ECMWF temperatures when considering both CH4 and H2O from CoMeCAT is a decrease of up to 1.0 K over the tropical mid/low stratosphere. The effect of using the CoMeCAT scheme for radiative forcing (RF) calculations is investigated using the offline Edwards–Slingo radiative transfer model. Compared to the default model option of a tropospheric global 3-D CH4 value, the CoMeCAT distribution produces an overall change in the annual mean net RF of up to −30 mW m−2.


2020 ◽  
Vol 13 (1) ◽  
pp. 116
Author(s):  
Lucie Leonarski ◽  
Laurent C.-Labonnote ◽  
Mathieu Compiègne ◽  
Jérôme Vidot ◽  
Anthony J. Baran ◽  
...  

The present study aims to quantify the potential of hyperspectral thermal infrared sounders such as the Infrared Atmospheric Sounding Interferometer (IASI) and the future IASI next generation (IASI-NG) for retrieving the ice cloud layer altitude and thickness together with the ice water path. We employed the radiative transfer model Radiative Transfer for TOVS (RTTOV) to simulate cloudy radiances using parameterized ice cloud optical properties. The radiances have been computed from an ice cloud profile database coming from global operational short-range forecasts at the European Center for Medium-range Weather Forecasts (ECMWF) which encloses the normal conditions, typical variability, and extremes of the atmospheric properties over one year (Eresmaa and McNally (2014)). We performed an information content analysis based on Shannon’s formalism to determine the amount and spectral distribution of the information about ice cloud properties. Based on this analysis, a retrieval algorithm has been developed and tested on the profile database. We considered the signal-to-noise ratio of each specific instrument and the non-retrieved atmospheric and surface parameter errors. This study brings evidence that the observing system provides information on the ice water path (IWP) as well as on the layer altitude and thickness with a convergence rate up to 95% and expected errors that decrease with cloud opacity until the signal saturation is reached (satisfying retrievals are achieved for clouds whose IWP is between about 1 and 300 g/m2).


2007 ◽  
Vol 7 (20) ◽  
pp. 5391-5400 ◽  
Author(s):  
K. M. Nissen ◽  
K. Matthes ◽  
U. Langematz ◽  
B. Mayer

Abstract. We introduce the improved Freie Universität Berlin (FUB) high-resolution radiation scheme FUBRad and compare it to the 4-band standard ECHAM5 SW radiation scheme of Fouquart and Bonnel (FB). Both schemes are validated against the detailed radiative transfer model libRadtran. FUBRad produces realistic heating rate variations during the solar cycle. The SW heating rate response with the FB scheme is about 20 times smaller than with FUBRad and cannot produce the observed temperature signal. A reduction of the spectral resolution to 6 bands for solar irradiance and ozone absorption cross sections leads to a degradation (reduction) of the solar SW heating rate signal by about 20%. The simulated temperature response agrees qualitatively well with observations in the summer upper stratosphere and mesosphere where irradiance variations dominate the signal. Comparison of the total short-wave heating rates under solar minimum conditions shows good agreement between FUBRad, FB and libRadtran up to the middle mesosphere (60–70 km) indicating that both parameterizations are well suited for climate integrations that do not take solar variability into account. The FUBRad scheme has been implemented as a sub-submodel of the Modular Earth Submodel System (MESSy).


2018 ◽  
Vol 11 (8) ◽  
pp. 4707-4723 ◽  
Author(s):  
Norbert Glatthor ◽  
Thomas von Clarmann ◽  
Gabriele P. Stiller ◽  
Michael Kiefer ◽  
Alexandra Laeng ◽  
...  

Abstract. Discrepancies in ozone retrievals in MIPAS channels A (685–970 cm−1) and AB (1020–1170 cm−1) have been a long-standing problem in MIPAS data analysis, amounting to an interchannel bias (AB–A) of up to 8 % between ozone volume mixing ratios in the altitude range 30–40 km. We discuss various candidate explanations, among them forward model and retrieval algorithm errors, interchannel calibration inconsistencies and spectroscopic data inconsistencies. We show that forward-modelling errors as well as errors in the retrieval algorithm can be ruled out as an explanation because the bias can be reproduced with an entirely independent retrieval algorithm (GEOFIT), relying on a different forward radiative transfer model. Instrumental and calibration issues can also be refuted as an explanation because ozone retrievals based on balloon-borne measurements with a different instrument (MIPAS-B) and an independent level-1 data processing scheme produce a rather similar interchannel bias. Thus, spectroscopic inconsistencies in the MIPAS database used for ozone retrieval are practically the only reason left. To further investigate this issue, we performed retrievals using additional spectroscopic databases. Various versions of the HITRAN database generally produced rather similar channel AB–A differences. Use of a different database, namely GEISA-2015, led to similar results in channel AB, but to even higher ozone volume mixing ratios for channel A retrievals, i.e. to a reversal of the bias. We show that the differences in MIPAS channel A retrievals result from about 13 % lower air-broadening coefficients of the strongest lines in the GEISA-2015 database. Since the errors in line intensity of the major lines used in MIPAS channels A and AB are reported to be considerably lower than the observed bias, we posit that a major part of the channel AB–A differences can be attributed to inconsistent air-broadening coefficients as well. To corroborate this assumption we show some clearly inconsistent air-broadening coefficients in the HITRAN-2008 database. The interchannel bias in retrieved ozone amounts can be reduced by increasing the air-broadening coefficients of the lines in MIPAS channel AB in the HITRAN-2008 database by 6 %–8 %.


2018 ◽  
Vol 146 (4) ◽  
pp. 1197-1218
Author(s):  
Michèle De La Chevrotière ◽  
John Harlim

This paper demonstrates the efficacy of data-driven localization mappings for assimilating satellite-like observations in a dynamical system of intermediate complexity. In particular, a sparse network of synthetic brightness temperature measurements is simulated using an idealized radiative transfer model and assimilated to the monsoon–Hadley multicloud model, a nonlinear stochastic model containing several thousands of model coordinates. A serial ensemble Kalman filter is implemented in which the empirical correlation statistics are improved using localization maps obtained from a supervised learning algorithm. The impact of the localization mappings is assessed in perfect-model observing system simulation experiments (OSSEs) as well as in the presence of model errors resulting from the misspecification of key convective closure parameters. In perfect-model OSSEs, the localization mappings that use adjacent correlations to improve the correlation estimated from small ensemble sizes produce robust accurate analysis estimates. In the presence of model error, the filter skills of the localization maps trained on perfect- and imperfect-model data are comparable.


2005 ◽  
Vol 6 (4) ◽  
pp. 460-475 ◽  
Author(s):  
Jason A. Otkin ◽  
Martha C. Anderson ◽  
John R. Mecikalski ◽  
George R. Diak

Abstract Reliable procedures that accurately map surface insolation over large domains at high spatial and temporal resolution are a great benefit for making the predictions of potential and actual evapotranspiration that are required by a variety of hydrological and agricultural applications. Here, estimates of hourly and daily integrated insolation at 20-km resolution, based on Geostationary Operational Environmental Satellite (GOES) visible imagery are compared to pyranometer measurements made at 11 sites in the U.S. Climate Reference Network (USCRN) over a continuous 15-month period. Such a comprehensive survey is necessary in order to examine the accuracy of the satellite insolation estimates over a diverse range of seasons and land surface types. The relatively simple physical model of insolation that is tested here yields good results, with seasonally averaged model errors of 62 (19%) and 15 (10%) W m−2 for hourly and daily-averaged insolation, respectively, including both clear- and cloudy-sky conditions. This level of accuracy is comparable, or superior, to results that have been obtained with more complex models of atmospheric radiative transfer. Model performance can be improved in the future by addressing a small elevation-related bias in the physical model, which is likely the result of inaccurate model precipitable water inputs or cloud-height assessments.


Sign in / Sign up

Export Citation Format

Share Document