Browse
Abstract
This article describes a new precipitation analysis algorithm developed by NASA for time-sensitive operations at the United States Air Force. Implemented as part of the Land Information System—a land modeling and data assimilation software framework—this NASA–Air Force Precipitation Analysis (NAFPA) combines numerical weather prediction model outputs with rain gauge measurements and satellite estimates to produce global, gridded 3-h accumulated precipitation fields at approximately 10-km resolution. Input observations are subjected to quality control checks before being used by the Bratseth analysis algorithm that converges to optimal interpolation. NAFPA assimilates up to 3.5 million observations without artificial data thinning or selection. To evaluate this new approach, a multiyear reanalysis is generated and intercompared with eight alternative precipitation products across the contiguous United States, Africa, and the monsoon region of eastern Asia. NAFPA yields superior accuracy and correlation over low-latency (up to 14 h) alternatives (numerical weather prediction and satellite retrievals), and often outperforms high-latency (up to 3.5 months) products, although the details for the latter vary by region and product. The development of NAFPA offers a high-quality, near-real-time product for use in meteorological, land surface, and hydrological research and applications.
Significance Statement
Precipitation is a key input to land modeling systems due to effects on soil moisture and other parts of the hydrologic cycle. It is also of interest to government decision-makers due to impacts on human activities. Here we present a new precipitation analysis based on available near-real-time data. By running the program for prior years and comparing with alternative products, we demonstrate that our analysis provides better accuracy and usually less bias than near-real-time satellite data alone, and better accuracy and correlation than data provided by numerical weather models. Our analysis is also competitive with other products created months after the fact, justifying confidence in using our analysis in near-real-time operations.
Abstract
This article describes a new precipitation analysis algorithm developed by NASA for time-sensitive operations at the United States Air Force. Implemented as part of the Land Information System—a land modeling and data assimilation software framework—this NASA–Air Force Precipitation Analysis (NAFPA) combines numerical weather prediction model outputs with rain gauge measurements and satellite estimates to produce global, gridded 3-h accumulated precipitation fields at approximately 10-km resolution. Input observations are subjected to quality control checks before being used by the Bratseth analysis algorithm that converges to optimal interpolation. NAFPA assimilates up to 3.5 million observations without artificial data thinning or selection. To evaluate this new approach, a multiyear reanalysis is generated and intercompared with eight alternative precipitation products across the contiguous United States, Africa, and the monsoon region of eastern Asia. NAFPA yields superior accuracy and correlation over low-latency (up to 14 h) alternatives (numerical weather prediction and satellite retrievals), and often outperforms high-latency (up to 3.5 months) products, although the details for the latter vary by region and product. The development of NAFPA offers a high-quality, near-real-time product for use in meteorological, land surface, and hydrological research and applications.
Significance Statement
Precipitation is a key input to land modeling systems due to effects on soil moisture and other parts of the hydrologic cycle. It is also of interest to government decision-makers due to impacts on human activities. Here we present a new precipitation analysis based on available near-real-time data. By running the program for prior years and comparing with alternative products, we demonstrate that our analysis provides better accuracy and usually less bias than near-real-time satellite data alone, and better accuracy and correlation than data provided by numerical weather models. Our analysis is also competitive with other products created months after the fact, justifying confidence in using our analysis in near-real-time operations.
Abstract
The water resources of the western United States have enormous agricultural and municipal demands. At the same time, droughts like the one enveloping the West in the summer of 2021 have disrupted supply of this strained and precious resource. Historically, seasonal forecasts of cool-season (November–March) precipitation from dynamical models such as North American Multi-Model Ensemble (NMME) and the Seasonal Forecasting System 5 (SEAS5) from the European Centre for Medium-Range Weather Forecasts have lacked sufficient skill to aid in Western stakeholders’ and water managers’ decision-making. Here, we propose a new empirical–statistical framework to improve cool-season precipitation forecasts across the contiguous United States (CONUS). This newly developed framework is called the Statistical Climate Ensemble Forecast (SCEF) model. The SCEF framework applies a principal component regression model to predictors and predictands that have undergone dimensionality reduction, where the predictors are large-scale meteorological variables that have been prefiltered in space. The forecasts of the SCEF model captures 12.0% of the total CONUS-wide standardized observed variance over the period 1982/83–2019/20, whereas NMME captures 7.2%. Over the more recent period 2000/01–2019/20, the SCEF, NMME, and SEAS5 models respectively capture 11.8%, 4.0%, and 4.1% of the total CONUS-wide standardized observed variance. An important finding is that much of the improved skill in the SCEF, with respect to models such as NMME and SEAS5, can be attributed to better forecasts across most of the western United States.
Abstract
The water resources of the western United States have enormous agricultural and municipal demands. At the same time, droughts like the one enveloping the West in the summer of 2021 have disrupted supply of this strained and precious resource. Historically, seasonal forecasts of cool-season (November–March) precipitation from dynamical models such as North American Multi-Model Ensemble (NMME) and the Seasonal Forecasting System 5 (SEAS5) from the European Centre for Medium-Range Weather Forecasts have lacked sufficient skill to aid in Western stakeholders’ and water managers’ decision-making. Here, we propose a new empirical–statistical framework to improve cool-season precipitation forecasts across the contiguous United States (CONUS). This newly developed framework is called the Statistical Climate Ensemble Forecast (SCEF) model. The SCEF framework applies a principal component regression model to predictors and predictands that have undergone dimensionality reduction, where the predictors are large-scale meteorological variables that have been prefiltered in space. The forecasts of the SCEF model captures 12.0% of the total CONUS-wide standardized observed variance over the period 1982/83–2019/20, whereas NMME captures 7.2%. Over the more recent period 2000/01–2019/20, the SCEF, NMME, and SEAS5 models respectively capture 11.8%, 4.0%, and 4.1% of the total CONUS-wide standardized observed variance. An important finding is that much of the improved skill in the SCEF, with respect to models such as NMME and SEAS5, can be attributed to better forecasts across most of the western United States.
Abstract
Estimating lake evaporation is a challenge due to both practical considerations and theoretical assumptions embedded in indirect methods. For the first time, we evaluated measurements from an optical microwave scintillometer (OMS) system over an open-water body under arid conditions. The OMS is a line-of-sight remote sensing technique that can be used to measure the sensible and latent heat fluxes over horizontal areas with pathlengths ranging from 0.5 to 10 km. We installed an OMS at a saline lake surrounded by a wet-salt crust in the Salar del Huasco, a heterogeneous desert landscape in the Atacama Desert. As a reference, we used eddy covariance systems installed over the two main surfaces in the OMS footprint. We performed a footprint analysis to reconstruct the surface contribution to the OMS measured fluxes (80% water and 20% wet salt). Furthermore, we investigated the applicability of the Monin–Obukhov similarity theory (MOST), which was needed to infer fluxes from the OMS-derived structure parameters to the fluxes. The OMS structure parameters and MOST were compromised, which we mitigated by fitting MOST coefficients to the site conditions. We argue that the MOST deviation from values found in the literature is due to the effects of the surface heterogeneity and the nonlocal processes induced by regional circulation. With the available dataset we were not able to rule out instrument issues, such as additional fluctuations to the scintillation signal due to absorption or the effect of vibration in high-wind conditions. The adjusted MOST coefficients lowered by a factor of 1.64 compared to using standard MOST coefficients. For H and LυE, we obtained zero-intercept linear regressions with correlations, R 2, of 0.92 and 0.96, respectively. We conclude that advances in MOST are needed to successfully apply the OMS method in landscapes characterized by complex heterogeneity such as the Salar del Huasco.
Abstract
Estimating lake evaporation is a challenge due to both practical considerations and theoretical assumptions embedded in indirect methods. For the first time, we evaluated measurements from an optical microwave scintillometer (OMS) system over an open-water body under arid conditions. The OMS is a line-of-sight remote sensing technique that can be used to measure the sensible and latent heat fluxes over horizontal areas with pathlengths ranging from 0.5 to 10 km. We installed an OMS at a saline lake surrounded by a wet-salt crust in the Salar del Huasco, a heterogeneous desert landscape in the Atacama Desert. As a reference, we used eddy covariance systems installed over the two main surfaces in the OMS footprint. We performed a footprint analysis to reconstruct the surface contribution to the OMS measured fluxes (80% water and 20% wet salt). Furthermore, we investigated the applicability of the Monin–Obukhov similarity theory (MOST), which was needed to infer fluxes from the OMS-derived structure parameters to the fluxes. The OMS structure parameters and MOST were compromised, which we mitigated by fitting MOST coefficients to the site conditions. We argue that the MOST deviation from values found in the literature is due to the effects of the surface heterogeneity and the nonlocal processes induced by regional circulation. With the available dataset we were not able to rule out instrument issues, such as additional fluctuations to the scintillation signal due to absorption or the effect of vibration in high-wind conditions. The adjusted MOST coefficients lowered by a factor of 1.64 compared to using standard MOST coefficients. For H and LυE, we obtained zero-intercept linear regressions with correlations, R 2, of 0.92 and 0.96, respectively. We conclude that advances in MOST are needed to successfully apply the OMS method in landscapes characterized by complex heterogeneity such as the Salar del Huasco.
Abstract
Increased operational use of convection-allowing models and ensembles offers substantial improvements for some aspects of convective weather forecasting; however, errors in quantitative precipitation forecasts (QPFs) from these models, especially those related to incorrect placement of heavy rainfall systems, limit their usefulness as an input into hydrological models. To improve understanding of QPF location errors, this study quantifies the displacement errors for the centroids of both 0–18-h accumulated rainfall and rainfall in the first hour after initiation of precipitation systems in both the High-Resolution Rapid Refresh Ensemble (HRRRE) and the High-Resolution Ensemble Forecast (HREF) for 30 events in the 2018 warm season. Ensemble member QPFs are compared to quantitative precipitation estimates (QPEs) obtained from the North Central River Forecast Center (NCRFC). HRRRE is found to have less spread in centroid locations than HREF, and both HRRRE and HREF 0–18-h QPF accumulations have less spread than the 1-h QPF accumulation when the precipitation event initiates. Furthermore, QPF centroids are most often displaced to the west in HRRRE for both 0–18-h QPF accumulation and the 1-h QPF accumulation when the precipitation event initiates. The 0–18-h QPF accumulation displacement errors can be reduced when adjustments are made to the forecasted position based upon displacement errors present in the first hour of precipitation, but only when the adjustments are a function of the intercardinal quadrant in which the initial hour QPF centroid was displaced.
Abstract
Increased operational use of convection-allowing models and ensembles offers substantial improvements for some aspects of convective weather forecasting; however, errors in quantitative precipitation forecasts (QPFs) from these models, especially those related to incorrect placement of heavy rainfall systems, limit their usefulness as an input into hydrological models. To improve understanding of QPF location errors, this study quantifies the displacement errors for the centroids of both 0–18-h accumulated rainfall and rainfall in the first hour after initiation of precipitation systems in both the High-Resolution Rapid Refresh Ensemble (HRRRE) and the High-Resolution Ensemble Forecast (HREF) for 30 events in the 2018 warm season. Ensemble member QPFs are compared to quantitative precipitation estimates (QPEs) obtained from the North Central River Forecast Center (NCRFC). HRRRE is found to have less spread in centroid locations than HREF, and both HRRRE and HREF 0–18-h QPF accumulations have less spread than the 1-h QPF accumulation when the precipitation event initiates. Furthermore, QPF centroids are most often displaced to the west in HRRRE for both 0–18-h QPF accumulation and the 1-h QPF accumulation when the precipitation event initiates. The 0–18-h QPF accumulation displacement errors can be reduced when adjustments are made to the forecasted position based upon displacement errors present in the first hour of precipitation, but only when the adjustments are a function of the intercardinal quadrant in which the initial hour QPF centroid was displaced.
Abstract
In this study, we used hourly observations to investigate the cooling effect of summer rainfall on surface air temperature (Ta) in a subtropical area, Guangdong province, South China. Data were categorized step-by-step by rainfall system (convection, monsoon, and typhoon), daily rainfall amount, and relative humidity (RH) level. Moreover, the average hourly Ta variation due to solar radiation was removed from all observations before statistical analysis. The results showed that the linear relationship between hourly Ta variation and rainfall intensity did not exist. However, the cooling effect of rainfall on Ta variation was dominant. In addition, convective rainfall does cause a greater temperature drop than the other two rainfall systems. After further partitioning all samples by RH level preceding the rainfall, the relationship between hourly Ta variation and rainfall intensity became distinctive. When RH was below 70%, rainfall-induced cooling became more substantial and scaled linearly with event intensity, but when RH exceeded 70%, the rainfall cooling effect was generally restrained by the RH increase. A strong correlation between hourly Ta variation and RH level preceding the rainfall suggests the importance of RH on the rainfall cooling effect.
Abstract
In this study, we used hourly observations to investigate the cooling effect of summer rainfall on surface air temperature (Ta) in a subtropical area, Guangdong province, South China. Data were categorized step-by-step by rainfall system (convection, monsoon, and typhoon), daily rainfall amount, and relative humidity (RH) level. Moreover, the average hourly Ta variation due to solar radiation was removed from all observations before statistical analysis. The results showed that the linear relationship between hourly Ta variation and rainfall intensity did not exist. However, the cooling effect of rainfall on Ta variation was dominant. In addition, convective rainfall does cause a greater temperature drop than the other two rainfall systems. After further partitioning all samples by RH level preceding the rainfall, the relationship between hourly Ta variation and rainfall intensity became distinctive. When RH was below 70%, rainfall-induced cooling became more substantial and scaled linearly with event intensity, but when RH exceeded 70%, the rainfall cooling effect was generally restrained by the RH increase. A strong correlation between hourly Ta variation and RH level preceding the rainfall suggests the importance of RH on the rainfall cooling effect.
Abstract
Tropical cyclone (TC) rainfall hazard assessment is subject to the bias in TC climatology estimation from climate simulations or synthetic downscaling. In this study, we investigate the uncertainty in TC rainfall hazard assessment induced by this bias using both rain gauge and radar observations and synthetic-storm-model-coupled TC rainfall simulations. We identify the storm’s maximum intensity, impact duration, and minimal distance to the site to be the three most important storm parameters for TC rainfall hazard, and the relationship between the important storm parameters and TC rainfall can be well captured by a physics-based TC rainfall model. The uncertainty in the synthetic rainfall hazard induced by the bias in TC climatology can be largely explained by the bias in the important storm parameters simulated by the synthetic storm model. Correcting the distribution of the most biased parameter may significantly improve rainfall hazard estimation. Bias correction based on the joint distribution of the important parameters may render more accurate rainfall hazard estimations; however, the general technical difficulties in resampling from high-dimensional joint probability distributions prevent more accurate estimations in some cases. The results of the study also support future investigation of the impact of climate change on TC rainfall hazards through the lens of future changes in the identified important storm parameters.
Abstract
Tropical cyclone (TC) rainfall hazard assessment is subject to the bias in TC climatology estimation from climate simulations or synthetic downscaling. In this study, we investigate the uncertainty in TC rainfall hazard assessment induced by this bias using both rain gauge and radar observations and synthetic-storm-model-coupled TC rainfall simulations. We identify the storm’s maximum intensity, impact duration, and minimal distance to the site to be the three most important storm parameters for TC rainfall hazard, and the relationship between the important storm parameters and TC rainfall can be well captured by a physics-based TC rainfall model. The uncertainty in the synthetic rainfall hazard induced by the bias in TC climatology can be largely explained by the bias in the important storm parameters simulated by the synthetic storm model. Correcting the distribution of the most biased parameter may significantly improve rainfall hazard estimation. Bias correction based on the joint distribution of the important parameters may render more accurate rainfall hazard estimations; however, the general technical difficulties in resampling from high-dimensional joint probability distributions prevent more accurate estimations in some cases. The results of the study also support future investigation of the impact of climate change on TC rainfall hazards through the lens of future changes in the identified important storm parameters.
Abstract
The use of global climate model (GCM) precipitation simulations typically requires corrections for precipitation biases at subgrid spatial scales, typically at daily or monthly time scales. However, over many regions GCMs underestimate the magnitudes of multiyear precipitation extremes in the observed climate, resulting in a likely underestimation of the magnitudes of multiyear precipitation extremes in future scenarios. The objective of this study is to propose a method to extract from GCMs more realistic scenarios of multiyear precipitation extremes over time horizons of decades to one century. This proposed correction method is analogous to widely used bias correction methods, except that it is applied to variability at longer time scales than previous implementations (i.e., multiyear rather than daily or monthly). A case study of precipitation over a basin from the New York City water supply system demonstrates the potential magnitude of the underestimation of multiyear precipitation using uncorrected GCM scenarios, and the potential impact of the correction on multiyear hydrological extremes. Overall, it is a practical, conceptually simple approach meant for water supply system impact studies, but can be used for any impact studies that require more realistic multiyear extreme precipitation extreme scenarios.
Significance Statement
The purpose of this study is to present a practical method to address a particular difficulty that in some regions arises in climate change impact studies: global climate models tend to underestimate the multiyear variability of precipitation over some regions, resulting in an underestimation of the magnitudes and/or intensities of prolonged droughts as well as prolonged wet periods. The method is analogous to widely used bias correction methods, except it is applied to variability at longer time scales than previous implementations (i.e., multiyear rather than daily or monthly). It is designed to provide more realistic estimates of extreme hydrological scenarios during the twenty-first century. Our particular interest is for managers of water supply systems, but the method may be of interest to others for whom multiyear precipitation extremes are critical.
Abstract
The use of global climate model (GCM) precipitation simulations typically requires corrections for precipitation biases at subgrid spatial scales, typically at daily or monthly time scales. However, over many regions GCMs underestimate the magnitudes of multiyear precipitation extremes in the observed climate, resulting in a likely underestimation of the magnitudes of multiyear precipitation extremes in future scenarios. The objective of this study is to propose a method to extract from GCMs more realistic scenarios of multiyear precipitation extremes over time horizons of decades to one century. This proposed correction method is analogous to widely used bias correction methods, except that it is applied to variability at longer time scales than previous implementations (i.e., multiyear rather than daily or monthly). A case study of precipitation over a basin from the New York City water supply system demonstrates the potential magnitude of the underestimation of multiyear precipitation using uncorrected GCM scenarios, and the potential impact of the correction on multiyear hydrological extremes. Overall, it is a practical, conceptually simple approach meant for water supply system impact studies, but can be used for any impact studies that require more realistic multiyear extreme precipitation extreme scenarios.
Significance Statement
The purpose of this study is to present a practical method to address a particular difficulty that in some regions arises in climate change impact studies: global climate models tend to underestimate the multiyear variability of precipitation over some regions, resulting in an underestimation of the magnitudes and/or intensities of prolonged droughts as well as prolonged wet periods. The method is analogous to widely used bias correction methods, except it is applied to variability at longer time scales than previous implementations (i.e., multiyear rather than daily or monthly). It is designed to provide more realistic estimates of extreme hydrological scenarios during the twenty-first century. Our particular interest is for managers of water supply systems, but the method may be of interest to others for whom multiyear precipitation extremes are critical.
Abstract
Soil hydrophysical properties are necessary components in weather and climate simulation, yet the parameter inaccuracies may introduce considerable uncertainty in the representation of surface water and energy fluxes. This study uses seasonal coupled simulations to examine the uncertainties in the North American atmospheric water cycle that result from the use of different soil datasets. Two soil datasets are considered: the State Soil Geographic dataset (STATSGO) from the U.S. Department of Agriculture and the Global Soil Dataset for Earth System Modeling (GSDE) from Beijing Normal University. Two simulations are conducted from 1 June to 31 August 2016–18 using the Weather Research and Forecasting (WRF) Model coupled with the Community Land Model (CLM) version 4 and applying each soil dataset. It is found that changes in soil texture lead to statistically significant differences in daily mean surface water and energy fluxes. The boundary layer thermodynamic structure responds to these changes in surface fluxes resulting in differences in mean CAPE and CIN, leading to conditions that are less conducive for precipitation. The soil-texture-related surface fluxes instigate dynamic responses as well. Low-level wind fields are altered, resulting in differences in the associated vertically integrated moisture fluxes and in vertically integrated moisture flux convergence in the same regions. Through land–atmosphere interactions, it is shown that soil parameters can affect each component of the atmospheric water budget.
Abstract
Soil hydrophysical properties are necessary components in weather and climate simulation, yet the parameter inaccuracies may introduce considerable uncertainty in the representation of surface water and energy fluxes. This study uses seasonal coupled simulations to examine the uncertainties in the North American atmospheric water cycle that result from the use of different soil datasets. Two soil datasets are considered: the State Soil Geographic dataset (STATSGO) from the U.S. Department of Agriculture and the Global Soil Dataset for Earth System Modeling (GSDE) from Beijing Normal University. Two simulations are conducted from 1 June to 31 August 2016–18 using the Weather Research and Forecasting (WRF) Model coupled with the Community Land Model (CLM) version 4 and applying each soil dataset. It is found that changes in soil texture lead to statistically significant differences in daily mean surface water and energy fluxes. The boundary layer thermodynamic structure responds to these changes in surface fluxes resulting in differences in mean CAPE and CIN, leading to conditions that are less conducive for precipitation. The soil-texture-related surface fluxes instigate dynamic responses as well. Low-level wind fields are altered, resulting in differences in the associated vertically integrated moisture fluxes and in vertically integrated moisture flux convergence in the same regions. Through land–atmosphere interactions, it is shown that soil parameters can affect each component of the atmospheric water budget.
Abstract
The Rwenzori Mountains, in southwest Uganda, are prone to precipitation-related hazards such as flash floods and landslides. These natural hazards highly impact the lives and livelihoods of the people living in the region. However, our understanding of the precipitation patterns and their impact on related hazardous events and/or agricultural productivity is hampered by a dearth of in situ precipitation observations. Here, we propose an evaluation of gridded precipitation products as potential candidates filling this hiatus. We evaluate three state-of-the-art gridded products, the ERA5 reanalysis, IMERG satellite observations, and a simulation from the convection-permitting climate model (CPM), COSMO-CLM, for their ability to represent precipitation totals, timing, and precipitation probability density function. The evaluation is performed against observations from 11 gauge stations that provide at least 2.5 years of hourly and half-hourly data, recorded between 2011 and 2016. Results indicate a poor performance of ERA5 with a persistent wet bias, mostly for stations in the rain shadow of the mountains. IMERG gives the best representation of the precipitation totals as indicated by bias score comparisons. The CPM outperforms both ERA5 and IMERG in representing the probability density function, while both IMERG and the CPM have a good skill in capturing precipitation seasonal and diurnal cycles. The better performance of CPM is attributable to its higher resolution. This study highlights the potential of using IMERG and CPM precipitation estimates for hydrological and impact modeling over the Rwenzori Mountains, preferring IMERG for precipitation totals and CPM for precipitation extremes.
Abstract
The Rwenzori Mountains, in southwest Uganda, are prone to precipitation-related hazards such as flash floods and landslides. These natural hazards highly impact the lives and livelihoods of the people living in the region. However, our understanding of the precipitation patterns and their impact on related hazardous events and/or agricultural productivity is hampered by a dearth of in situ precipitation observations. Here, we propose an evaluation of gridded precipitation products as potential candidates filling this hiatus. We evaluate three state-of-the-art gridded products, the ERA5 reanalysis, IMERG satellite observations, and a simulation from the convection-permitting climate model (CPM), COSMO-CLM, for their ability to represent precipitation totals, timing, and precipitation probability density function. The evaluation is performed against observations from 11 gauge stations that provide at least 2.5 years of hourly and half-hourly data, recorded between 2011 and 2016. Results indicate a poor performance of ERA5 with a persistent wet bias, mostly for stations in the rain shadow of the mountains. IMERG gives the best representation of the precipitation totals as indicated by bias score comparisons. The CPM outperforms both ERA5 and IMERG in representing the probability density function, while both IMERG and the CPM have a good skill in capturing precipitation seasonal and diurnal cycles. The better performance of CPM is attributable to its higher resolution. This study highlights the potential of using IMERG and CPM precipitation estimates for hydrological and impact modeling over the Rwenzori Mountains, preferring IMERG for precipitation totals and CPM for precipitation extremes.
Abstract
Satellite-based and reanalysis precipitation estimates are an alternative and important supplement to rain gauge data. However, performance of China’s Fengyun (FY) satellite precipitation product and how it compares with other mainstream satellite and reanalysis precipitation products over China remain largely unknown. Here five satellite-based precipitation products (i.e., FY-2 precipitation product, IMERG, GSMaP, CMORPH, and PERSIANN-CDR) and one reanalysis product (i.e., ERA5) are intercompared and evaluated based on in situ daily precipitation measurements over mainland China during 2007–17. Results show that the performance of these precipitation products varies with regions and seasons, with better statistical metrics over wet regions and during warm seasons. The infrared–microwave combined precipitation [i.e., IMERG, GSMaP, and CMORPH, with median KGE (Kling–Gupta efficiency) values of 0.53, 0.52, 0.59, respectively] reveals better performance than the infrared-based only product (i.e., PERSIANN-CDR, with a median KGE of 0.31) and the reanalysis product (i.e., ERA5, with a median KGE of 0.43). IMERG performs well in retrieving precipitation intensity and occurrence over China, while GSMaP performs well in the middle to low reaches of the Yangtze River basin but poorly over sparsely gauged regions, e.g., Xinjiang in northwest China and the Tibetan Plateau. CMORPH performs well over most regions and has a greater ability to detect precipitation events than GSMaP. The FY-2 precipitation product can capture the overall spatial distribution of precipitation in terms of both precipitation intensity and occurrence (median KGE and CSI of 0.54 and 0.55), and shows better performance than other satellite precipitation products in winter and over sparsely gauged regions. Annual precipitation from different products is generally consistent, though underestimation exists in the FY-2 precipitation product during 2015–17.
Significance Statement
Intercomparison between the FY-2 precipitation product and mainstream precipitation products is valuable to guide applications of satellite precipitation products to China and its subregions. This study illustrates uncertainties in various satellite precipitation products, and could guide optimization of algorithms of precipitation retrieval and data fusion/merging to improve the accuracy and resolution of satellite precipitation products.
Abstract
Satellite-based and reanalysis precipitation estimates are an alternative and important supplement to rain gauge data. However, performance of China’s Fengyun (FY) satellite precipitation product and how it compares with other mainstream satellite and reanalysis precipitation products over China remain largely unknown. Here five satellite-based precipitation products (i.e., FY-2 precipitation product, IMERG, GSMaP, CMORPH, and PERSIANN-CDR) and one reanalysis product (i.e., ERA5) are intercompared and evaluated based on in situ daily precipitation measurements over mainland China during 2007–17. Results show that the performance of these precipitation products varies with regions and seasons, with better statistical metrics over wet regions and during warm seasons. The infrared–microwave combined precipitation [i.e., IMERG, GSMaP, and CMORPH, with median KGE (Kling–Gupta efficiency) values of 0.53, 0.52, 0.59, respectively] reveals better performance than the infrared-based only product (i.e., PERSIANN-CDR, with a median KGE of 0.31) and the reanalysis product (i.e., ERA5, with a median KGE of 0.43). IMERG performs well in retrieving precipitation intensity and occurrence over China, while GSMaP performs well in the middle to low reaches of the Yangtze River basin but poorly over sparsely gauged regions, e.g., Xinjiang in northwest China and the Tibetan Plateau. CMORPH performs well over most regions and has a greater ability to detect precipitation events than GSMaP. The FY-2 precipitation product can capture the overall spatial distribution of precipitation in terms of both precipitation intensity and occurrence (median KGE and CSI of 0.54 and 0.55), and shows better performance than other satellite precipitation products in winter and over sparsely gauged regions. Annual precipitation from different products is generally consistent, though underestimation exists in the FY-2 precipitation product during 2015–17.
Significance Statement
Intercomparison between the FY-2 precipitation product and mainstream precipitation products is valuable to guide applications of satellite precipitation products to China and its subregions. This study illustrates uncertainties in various satellite precipitation products, and could guide optimization of algorithms of precipitation retrieval and data fusion/merging to improve the accuracy and resolution of satellite precipitation products.
