Search Results
You are looking at 21 - 30 of 69 items for
- Author or Editor: Witold F. Krajewski x
- Refine by Access: All Content x
Abstract
IMERG provides state-of-the-art satellite-based precipitation estimates that combine observations from multiple satellite platforms. This study evaluates IMERG products by examining hydrologic simulations of streamflow at a range of spatial scales. The main objective of this study is to assess the predictive utility of the near-real-time product (IMERG-Early). The assessment also includes the IMERG-Final product that is not available in real time. The authors used MRMS precipitation estimates and USGS streamflow observation data as references for the precipitation and streamflow evaluations during a 5-yr period (2016–20). The precipitation evaluation results show that IMERG-Early yields significant overestimations, particularly during warm months, with higher variability in its conditional distributions, whereas the performance of IMERG-Final seems unbiased. The authors performed hydrologic simulations using the Iowa Flood Center’s Hillslope Link Model with three precipitation forcing products, i.e., MRMS, IMERG-Early, and IMERG-Final. The simulation results reveal that IMERG-Early leads to high hit and false alarm rates due to its overestimation in precipitation and has almost no skill, as measured by the overall performance metric Kling–Gupta efficiency (KGE), in streamflow prediction regarding basin scales ranging from 10 to 30 000 km2. This indicates that the product requires a bias correction before it is useful for real-time flood prediction. The streamflow prediction performance of IMERG-Final seems comparable to that of MRMS at spatial scales greater than 100 km2. This scale limitation is attributable to the IMERG’s product spatial resolution that is inadequate to capture the small-scale variability of precipitation.
Abstract
IMERG provides state-of-the-art satellite-based precipitation estimates that combine observations from multiple satellite platforms. This study evaluates IMERG products by examining hydrologic simulations of streamflow at a range of spatial scales. The main objective of this study is to assess the predictive utility of the near-real-time product (IMERG-Early). The assessment also includes the IMERG-Final product that is not available in real time. The authors used MRMS precipitation estimates and USGS streamflow observation data as references for the precipitation and streamflow evaluations during a 5-yr period (2016–20). The precipitation evaluation results show that IMERG-Early yields significant overestimations, particularly during warm months, with higher variability in its conditional distributions, whereas the performance of IMERG-Final seems unbiased. The authors performed hydrologic simulations using the Iowa Flood Center’s Hillslope Link Model with three precipitation forcing products, i.e., MRMS, IMERG-Early, and IMERG-Final. The simulation results reveal that IMERG-Early leads to high hit and false alarm rates due to its overestimation in precipitation and has almost no skill, as measured by the overall performance metric Kling–Gupta efficiency (KGE), in streamflow prediction regarding basin scales ranging from 10 to 30 000 km2. This indicates that the product requires a bias correction before it is useful for real-time flood prediction. The streamflow prediction performance of IMERG-Final seems comparable to that of MRMS at spatial scales greater than 100 km2. This scale limitation is attributable to the IMERG’s product spatial resolution that is inadequate to capture the small-scale variability of precipitation.
Abstract
Although it is broadly acknowledged that the radar-rainfall (RR) estimates based on the U.S. national network of Weather Surveillance Radar-1988 Doppler (WSR-88D) stations contain a high degree of uncertainty, no methods currently exist to inform users about its quantitative characteristics. The most comprehensive characterization of this uncertainty can be achieved by delivering the products in a probabilistic rather than the traditional deterministic form. The authors are developing a methodology for probabilistic quantitative precipitation estimation (PQPE) based on weather radar data. In this study, they present the central element of this methodology: an empirically based error structure model for the RR products.
The authors apply a product-error-driven (PED) approach to obtain a realistic uncertainty model. It is based on the analyses of six years of data from the Oklahoma City, Oklahoma, WSR-88D radar (KTLX) processed with the Precipitation Processing System algorithm of the NEXRAD system. The modeled functional-statistical relationship between RR estimates and corresponding true rainfall consists of two components: a systematic distortion function and a stochastic factor quantifying remaining random errors. The two components are identified using a nonparametric functional estimation apparatus. The true rainfall is approximated with rain gauge data from the Oklahoma Mesonet and the U.S. Department of Agriculture (USDA) Agricultural Research Service Micronet networks. The RR uncertainty model presented here accounts for different time scales, synoptic regimes, and distances from the radar. In addition, this study marks the first time in which results on RR error correlation in space and time are presented.
Abstract
Although it is broadly acknowledged that the radar-rainfall (RR) estimates based on the U.S. national network of Weather Surveillance Radar-1988 Doppler (WSR-88D) stations contain a high degree of uncertainty, no methods currently exist to inform users about its quantitative characteristics. The most comprehensive characterization of this uncertainty can be achieved by delivering the products in a probabilistic rather than the traditional deterministic form. The authors are developing a methodology for probabilistic quantitative precipitation estimation (PQPE) based on weather radar data. In this study, they present the central element of this methodology: an empirically based error structure model for the RR products.
The authors apply a product-error-driven (PED) approach to obtain a realistic uncertainty model. It is based on the analyses of six years of data from the Oklahoma City, Oklahoma, WSR-88D radar (KTLX) processed with the Precipitation Processing System algorithm of the NEXRAD system. The modeled functional-statistical relationship between RR estimates and corresponding true rainfall consists of two components: a systematic distortion function and a stochastic factor quantifying remaining random errors. The two components are identified using a nonparametric functional estimation apparatus. The true rainfall is approximated with rain gauge data from the Oklahoma Mesonet and the U.S. Department of Agriculture (USDA) Agricultural Research Service Micronet networks. The RR uncertainty model presented here accounts for different time scales, synoptic regimes, and distances from the radar. In addition, this study marks the first time in which results on RR error correlation in space and time are presented.
Abstract
The climatology of thunderstorms and flash floods in the Baltimore, Maryland, metropolitan region is examined through analyses of cloud-to-ground (CG) lightning observations from the National Lightning Detection Network (NLDN) and discharge observations from 11 U.S. Geological Survey (USGS) stream gauging stations. A point process framework is used for analyses of CG lightning strikes and the occurrences of flash floods. Analyses of lightning strikes as a space–time point process focus on the mean intensity function, from which the seasonal, diurnal, and spatial variation in mean lightning frequency are examined. Important elements of the spatial variation of mean lightning frequency are 1) initiation of thunderstorms along the Blue Ridge, 2) large variability of lightning frequency around the urban cores of Baltimore and Washington D.C., and 3) decreased lightning frequency over the Chesapeake Bay and Atlantic Ocean. Lightning frequency has a sharp seasonal maximum around mid-July, and the diurnal cycle of lightning frequency peaks between 2100 and 2200 UTC with a frequency that is more than an order of magnitude larger than the minimum frequency at 1200 UTC. The seasonal and diurnal variation of flash flood occurrence in urban streams of Baltimore mimics the seasonal and diurnal variation of lightning. The peak of the diurnal frequency of flash floods in Moores Run, a 9.1-km2 urban watershed in Baltimore City, occurs at 2200 UTC. Analyses of the lightning and flood peak data also show a close link between the occurrence of major thunderstorms systems and flash flooding on a regional scale.
Abstract
The climatology of thunderstorms and flash floods in the Baltimore, Maryland, metropolitan region is examined through analyses of cloud-to-ground (CG) lightning observations from the National Lightning Detection Network (NLDN) and discharge observations from 11 U.S. Geological Survey (USGS) stream gauging stations. A point process framework is used for analyses of CG lightning strikes and the occurrences of flash floods. Analyses of lightning strikes as a space–time point process focus on the mean intensity function, from which the seasonal, diurnal, and spatial variation in mean lightning frequency are examined. Important elements of the spatial variation of mean lightning frequency are 1) initiation of thunderstorms along the Blue Ridge, 2) large variability of lightning frequency around the urban cores of Baltimore and Washington D.C., and 3) decreased lightning frequency over the Chesapeake Bay and Atlantic Ocean. Lightning frequency has a sharp seasonal maximum around mid-July, and the diurnal cycle of lightning frequency peaks between 2100 and 2200 UTC with a frequency that is more than an order of magnitude larger than the minimum frequency at 1200 UTC. The seasonal and diurnal variation of flash flood occurrence in urban streams of Baltimore mimics the seasonal and diurnal variation of lightning. The peak of the diurnal frequency of flash floods in Moores Run, a 9.1-km2 urban watershed in Baltimore City, occurs at 2200 UTC. Analyses of the lightning and flood peak data also show a close link between the occurrence of major thunderstorms systems and flash flooding on a regional scale.
Abstract
This study demonstrates an implementation of the prototype quantitative precipitation R estimation algorithm using specific attenuation A for S-band polarimetric radar. The performance of R(A) algorithm is assessed, compared to the conventional algorithm using radar reflectivity Z, at multiple temporal scales. Because the factor α, defined as the net ratio of A to specific differential phase, is a key parameter of the algorithm characterized by drop size distributions (e.g., differential reflectivity Z dr dependence on Z), the estimation equations of α and a proper number of Z dr–Z samples required for a reliable α estimation are examined. Based on the dynamic estimation of α, the event-based evaluation using hourly rain gauge observations reveals that the performance of R(A) is superior to that of R(Z), with better agreement and lower variability. Despite its superiority, the study finds that R(A) leads to quite consistent overestimations of about 10%–30%. It is demonstrated that the application of uniform α over the entire radar domain yields the observed uncertainty because of the heterogeneity of precipitation in the domain. A climatological range-dependent feature of R(A) and R(Z) is inspected in the multiyear evaluation at yearly scale using rain totals for April–October. While R(Z) exposes a systematic shift and overestimation, each of which arise from the radar miscalibration and bright band effects, R(A) combining with multiple R(Z) values for solid/mixed precipitation shows relatively robust performance without those effects. The immunity of R(A) to partial beam blockage (PBB) based on both qualitative and quantitative analyses is also verified. However, the capability of R(A) regarding PBB is limited by the presence of the melting layer and its application requirement for the total span of differential phase (e.g., 3°), which is another challenge for light rain.
Abstract
This study demonstrates an implementation of the prototype quantitative precipitation R estimation algorithm using specific attenuation A for S-band polarimetric radar. The performance of R(A) algorithm is assessed, compared to the conventional algorithm using radar reflectivity Z, at multiple temporal scales. Because the factor α, defined as the net ratio of A to specific differential phase, is a key parameter of the algorithm characterized by drop size distributions (e.g., differential reflectivity Z dr dependence on Z), the estimation equations of α and a proper number of Z dr–Z samples required for a reliable α estimation are examined. Based on the dynamic estimation of α, the event-based evaluation using hourly rain gauge observations reveals that the performance of R(A) is superior to that of R(Z), with better agreement and lower variability. Despite its superiority, the study finds that R(A) leads to quite consistent overestimations of about 10%–30%. It is demonstrated that the application of uniform α over the entire radar domain yields the observed uncertainty because of the heterogeneity of precipitation in the domain. A climatological range-dependent feature of R(A) and R(Z) is inspected in the multiyear evaluation at yearly scale using rain totals for April–October. While R(Z) exposes a systematic shift and overestimation, each of which arise from the radar miscalibration and bright band effects, R(A) combining with multiple R(Z) values for solid/mixed precipitation shows relatively robust performance without those effects. The immunity of R(A) to partial beam blockage (PBB) based on both qualitative and quantitative analyses is also verified. However, the capability of R(A) regarding PBB is limited by the presence of the melting layer and its application requirement for the total span of differential phase (e.g., 3°), which is another challenge for light rain.
Abstract
This study addresses the uncertainty of High-Resolution Rapid Refresh (HRRR) quantitative precipitation forecasts (QPFs), which were recently appended to the operational hydrologic forecasting framework. In this study, we examine the uncertainty features of HRRR QPFs for an Iowa flooding event that occurred in September 2016. Our evaluation of HRRR QPFs is based on the conventional approach of QPF verification and the analysis of mean areal precipitation (MAP) with respect to forecast lead time. The QPF verification results show that the precipitation forecast skill of HRRR significantly drops during short lead times and then gradually decreases for further lead times. The MAP analysis also demonstrates that the QPF error sharply increases during short lead times and starts decreasing slightly beyond 4-h lead time. We found that the variability of QPF error measured in terms of MAP decreases as basin scale and lead time become larger and longer, respectively. The effects of QPF uncertainty on hydrologic prediction are quantified through the hillslope-link model (HLM) simulations using hydrologic performance metrics (e.g., Kling–Gupta efficiency). The simulation results agree to some degree with those from the MAP analysis, finding that the performance achieved from the QPF forcing decreases during 1–3-h lead times and starts increasing with 4–6-h lead times. The best performance acquired at the 1-h lead time does not seem acceptable because of the large overestimation of the flood peak, along with an erroneous early peak that is not observed in streamflow observations. This study provides further evidence that HRRR contains a well-known weakness at short lead times, and the QPF uncertainty (e.g., bias) described as a function of forecast lead times should be corrected before its use in hydrologic prediction.
Abstract
This study addresses the uncertainty of High-Resolution Rapid Refresh (HRRR) quantitative precipitation forecasts (QPFs), which were recently appended to the operational hydrologic forecasting framework. In this study, we examine the uncertainty features of HRRR QPFs for an Iowa flooding event that occurred in September 2016. Our evaluation of HRRR QPFs is based on the conventional approach of QPF verification and the analysis of mean areal precipitation (MAP) with respect to forecast lead time. The QPF verification results show that the precipitation forecast skill of HRRR significantly drops during short lead times and then gradually decreases for further lead times. The MAP analysis also demonstrates that the QPF error sharply increases during short lead times and starts decreasing slightly beyond 4-h lead time. We found that the variability of QPF error measured in terms of MAP decreases as basin scale and lead time become larger and longer, respectively. The effects of QPF uncertainty on hydrologic prediction are quantified through the hillslope-link model (HLM) simulations using hydrologic performance metrics (e.g., Kling–Gupta efficiency). The simulation results agree to some degree with those from the MAP analysis, finding that the performance achieved from the QPF forcing decreases during 1–3-h lead times and starts increasing with 4–6-h lead times. The best performance acquired at the 1-h lead time does not seem acceptable because of the large overestimation of the flood peak, along with an erroneous early peak that is not observed in streamflow observations. This study provides further evidence that HRRR contains a well-known weakness at short lead times, and the QPF uncertainty (e.g., bias) described as a function of forecast lead times should be corrected before its use in hydrologic prediction.
Abstract
In view of the importance of tropical rainfall and the ubiquitous need for its estimates in climate modeling, the authors assess the ability of the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) to characterize the scaling characteristics of rainfall by comparing the derived results with those obtained from the ground-based radar (GR) data. The analysis is based on 59 months of PR and GR rain rates at three TRMM ground validation (GV) sites: Houston, Texas; Melbourne, Florida; and Kwajalein Atoll, Republic of the Marshall Islands. The authors consider spatial scales ranging from about 4 to 64 km at a fixed temporal scale corresponding to the sensor “instantaneous” snapshots (∼15 min). The focus is on the scaling of the marginal moments, which allows estimation of the scaling parameters from a single scene of data. The standard rainfall products of the PR and the GR are compared in terms of distributions of the scaling parameter estimates, the connection between the scaling parameters and the large-scale spatial average rain rate, and deviations from scale invariance. The five main results are as follows: 1) the PR yields values of the rain intermittence scaling parameter within 20% of the GR estimate; 2) both the PR and GR data show a one-to-one relationship between the intermittence scaling parameter and the large-scale spatial average rain rate that can be fit with the same functional form; 3) the PR underestimates the curvature of the scaling function from 20% to 50%, implying that high rain-rate extremes would be missed in a downscaling procedure; 4) the majority of the scenes (>85%) from both the PR and GR are scale invariant at the moment orders q = 0 and 2; and 5) the scale-invariance property tends to break down more likely over ocean than over land; the rainfall regimes that are not scale invariant are dominated by light storms covering large areas. Our results further show that for a sampling size of one year of data, the TRMM temporal sampling does not significantly affect the derived scaling characteristics. The authors conclude that the TRMM PR has the ability to characterize the basic scaling properties of rainfall, though the resulting parameters are subject to some degree of uncertainty.
Abstract
In view of the importance of tropical rainfall and the ubiquitous need for its estimates in climate modeling, the authors assess the ability of the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) to characterize the scaling characteristics of rainfall by comparing the derived results with those obtained from the ground-based radar (GR) data. The analysis is based on 59 months of PR and GR rain rates at three TRMM ground validation (GV) sites: Houston, Texas; Melbourne, Florida; and Kwajalein Atoll, Republic of the Marshall Islands. The authors consider spatial scales ranging from about 4 to 64 km at a fixed temporal scale corresponding to the sensor “instantaneous” snapshots (∼15 min). The focus is on the scaling of the marginal moments, which allows estimation of the scaling parameters from a single scene of data. The standard rainfall products of the PR and the GR are compared in terms of distributions of the scaling parameter estimates, the connection between the scaling parameters and the large-scale spatial average rain rate, and deviations from scale invariance. The five main results are as follows: 1) the PR yields values of the rain intermittence scaling parameter within 20% of the GR estimate; 2) both the PR and GR data show a one-to-one relationship between the intermittence scaling parameter and the large-scale spatial average rain rate that can be fit with the same functional form; 3) the PR underestimates the curvature of the scaling function from 20% to 50%, implying that high rain-rate extremes would be missed in a downscaling procedure; 4) the majority of the scenes (>85%) from both the PR and GR are scale invariant at the moment orders q = 0 and 2; and 5) the scale-invariance property tends to break down more likely over ocean than over land; the rainfall regimes that are not scale invariant are dominated by light storms covering large areas. Our results further show that for a sampling size of one year of data, the TRMM temporal sampling does not significantly affect the derived scaling characteristics. The authors conclude that the TRMM PR has the ability to characterize the basic scaling properties of rainfall, though the resulting parameters are subject to some degree of uncertainty.
Abstract
Quantifying uncertainty associated with flash flood warning or forecast systems is required to enable informed decision making by those responsible for operation and management of natural hazard protection systems. The current system used by the U.S. National Weather Service (NWS) to issue flash-flood warnings and watches over the Unites States is a purely deterministic system. The authors propose a simple approach to augment the Flash Flood Guidance System (FFGS) with uncertainty propagation components. The authors briefly discuss the main components of the system, propose changes to improve it, and allow accounting for several sources of uncertainty. They illustrate their discussion with examples of uncertainty quantification procedures for several small basins of the Illinois River basin in Oklahoma. As the current FFGS is tightly coupled with two technologies, that is, threshold-runoff mapping and the Sacramento Soil Moisture Accounting Hydrologic Model, the authors discuss both as sources of uncertainty. To quantify and propagate those sources of uncertainty throughout the system, they develop a simple version of the Sacramento model and use Monte Carlo simulation to study several uncertainty scenarios. The results point out the significance of the stream characteristics such as top width and the hydraulic depth on the overall uncertainty of the Flash Flood Guidance System. They also show that the overall flash flood guidance uncertainty is higher under drier initial soil moisture conditions. The results presented herein, although limited, are a necessary first step toward the development of probabilistic operational flash flood guidance forecast-response systems.
Abstract
Quantifying uncertainty associated with flash flood warning or forecast systems is required to enable informed decision making by those responsible for operation and management of natural hazard protection systems. The current system used by the U.S. National Weather Service (NWS) to issue flash-flood warnings and watches over the Unites States is a purely deterministic system. The authors propose a simple approach to augment the Flash Flood Guidance System (FFGS) with uncertainty propagation components. The authors briefly discuss the main components of the system, propose changes to improve it, and allow accounting for several sources of uncertainty. They illustrate their discussion with examples of uncertainty quantification procedures for several small basins of the Illinois River basin in Oklahoma. As the current FFGS is tightly coupled with two technologies, that is, threshold-runoff mapping and the Sacramento Soil Moisture Accounting Hydrologic Model, the authors discuss both as sources of uncertainty. To quantify and propagate those sources of uncertainty throughout the system, they develop a simple version of the Sacramento model and use Monte Carlo simulation to study several uncertainty scenarios. The results point out the significance of the stream characteristics such as top width and the hydraulic depth on the overall uncertainty of the Flash Flood Guidance System. They also show that the overall flash flood guidance uncertainty is higher under drier initial soil moisture conditions. The results presented herein, although limited, are a necessary first step toward the development of probabilistic operational flash flood guidance forecast-response systems.
Abstract
The relationship between the fractional time raining and tropical rainfall amount is investigated using raingage data and a point process model of tropical rainfall. Both the strength and the nature of the relationship are dependent upon the resolution of the data used to estimate the fractional time raining. It is found that highly accurate estimates of rainfall amounts over periods of one month or greater can be obtained from the fractional time raining so long as high-time-resolution data are used. It is demonstrated that the relationship between the fractional time raining and monthly atoll rainfall is quasi-homogeneous within the monsoon trough region of the equatorial western Pacific.
Abstract
The relationship between the fractional time raining and tropical rainfall amount is investigated using raingage data and a point process model of tropical rainfall. Both the strength and the nature of the relationship are dependent upon the resolution of the data used to estimate the fractional time raining. It is found that highly accurate estimates of rainfall amounts over periods of one month or greater can be obtained from the fractional time raining so long as high-time-resolution data are used. It is demonstrated that the relationship between the fractional time raining and monthly atoll rainfall is quasi-homogeneous within the monsoon trough region of the equatorial western Pacific.
Abstract
The goal of this study is to improve understanding of the optimization criteria for radar rainfall (RR) products. Conditional bias (CB) is formally defined and discussed. The CB is defined as the difference between a given rain rate and the conditional average of its estimates. A simple analytical model is used to study the behavior of CB and its effect on the relationship between the estimates and the truth. This study shows the measurement errors of near-surface radar reflectivity and the natural reflectivity–rainfall rate variability can affect CB. This RR estimation error component is also compared with the commonly used mean-square error (MSE). A dilemma between the minimization of these two errors is demonstrated. Removing CB from the estimates significantly increases MSE, but minimizing MSE results in a large CB that manifests itself in underestimation of strong rainfalls.
Abstract
The goal of this study is to improve understanding of the optimization criteria for radar rainfall (RR) products. Conditional bias (CB) is formally defined and discussed. The CB is defined as the difference between a given rain rate and the conditional average of its estimates. A simple analytical model is used to study the behavior of CB and its effect on the relationship between the estimates and the truth. This study shows the measurement errors of near-surface radar reflectivity and the natural reflectivity–rainfall rate variability can affect CB. This RR estimation error component is also compared with the commonly used mean-square error (MSE). A dilemma between the minimization of these two errors is demonstrated. Removing CB from the estimates significantly increases MSE, but minimizing MSE results in a large CB that manifests itself in underestimation of strong rainfalls.
Abstract
Brightband effects are one of the more important causes of vertical variability of reflectivity and severely affect the accuracy of rainfall estimates from ground-based radar. Monte Carlo simulation experiments are performed to investigate the efficiency of a procedure for the correction of errors related to the vertical variability of reflectivity. The simulation model generates three-dimensional radar reflectivity fields. Brightband effects are simulated through a physically based model of melting-layer reflectivity observations. Sensitivity of the correction procedure for a number of different precipitation scenarios and radar systems is analyzed. Overall, the identification method is found to be a robust procedure for correction of brightband effects. Results indicate a dependence of the effectiveness of the correction procedure on mean altitude and spatial variability of the melting layer.
Abstract
Brightband effects are one of the more important causes of vertical variability of reflectivity and severely affect the accuracy of rainfall estimates from ground-based radar. Monte Carlo simulation experiments are performed to investigate the efficiency of a procedure for the correction of errors related to the vertical variability of reflectivity. The simulation model generates three-dimensional radar reflectivity fields. Brightband effects are simulated through a physically based model of melting-layer reflectivity observations. Sensitivity of the correction procedure for a number of different precipitation scenarios and radar systems is analyzed. Overall, the identification method is found to be a robust procedure for correction of brightband effects. Results indicate a dependence of the effectiveness of the correction procedure on mean altitude and spatial variability of the melting layer.
