Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantile regression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

What factors affect the prices of low-priced U.S. solar PV systems?

The price of solar PV systems has declined rapidly, yet there are some much lower-priced systems than others. This study explores the factors that determine prices in these low-priced (LP) systems. Using a data set of 42,611 residential-scale PV systems installed in the U.S. in 2013, we use quantile regressions to estimate the importance of factors affecting the installed prices for LP systems (those at the 10th percentile) in comparison to median-priced systems. We find that the value of solar to consumers-a variable that accounts for subsidies, electric rates, and PV generation levels-is associated with lower prices for LP systems but higher prices for median priced systems. Conversely, systems installed in new home construction are associated with lower prices at the median but higher prices for LP. Other variables have larger price-reducing effects on LP than on median priced systems: systems installed in Arizona and Florida, as well as commercial and thin film systems. In contrast, the following have a smaller effect on prices for LP systems than median priced systems: tracking systems, self-installations, systems installed in Massachusetts, the system size, and installer experience. Furthermore, these results highlight the complex factors at play that lead to LP systems and shed light into how such LP systems can come about.

14 SOLAR ENERGY↗
DOE OSTI2017

Constructing probabilistic scenarios for wide-area solar power generation

Optimizing thermal generation commitments and dispatch in the presence of high penetrations of renewable resources such as solar energy requires a characterization of their stochastic properties. In this study, we describe novel methods designed to create day-ahead, wide-area probabilistic solar power scenarios based only on historical forecasts and associated observations of solar power production. Each scenario represents a possible trajectory for solar power in next-day operations with an associated probability computed by algorithms that use historical forecast errors. Scenarios are created by segmentation of historic data, fitting non-parametric error distributions using epi-splines, and then computing specific quantiles from these distributions. Additionally, we address the challenge of establishing an upper bound on solar power output. Our specific application driver is for use in stochastic variants of core power systems operations optimization problems, e.g., unit commitment and economic dispatch. These problems require as input a range of possible future realizations of renewables production. However, the utility of such probabilistic scenarios extends to other contexts, e.g., operator and trader situational awareness. Finally, we compare the performance of our approach to a recently proposed method based on quantile regression, and demonstrate that our method performs comparably to this approach in terms of two widely used methods for assessing the quality of probabilistic scenarios: the Energy score and the Variogram score.

14 SOLAR ENERGY↗
DOE OSTI2017

Economic consequences of aviation system disruptions: A reduced-form computable general equilibrium analysis

The state of the art approach to economic consequence analysis (ECA) is computable general equilibrium (CGE) modeling. However, such models contain thousands of equations and cannot readily be incorporated into computerized systems used by policy analysts to yield estimates of economic impacts of various types of transportation system failures due to natural hazards, human related attacks or technological accidents. This paper presents a reduced-form approach to simplify the analytical content of CGE models to make them more transparent and enhance their utilization potential. The reduced-form CGE analysis is conducted by first running simulations one hundred times, varying key parameters, such as magnitude of the initial shock, duration, location, remediation, and resilience, according to a Latin Hypercube sampling procedure. Statistical analysis is then applied to the “synthetic data” results in the form of both ordinary least squares and quantile regression. The analysis yields linear equations that are incorporated into a computerized system and utilized along with Monte Carlo simulation methods for propagating uncertainties in economic consequences. Although our demonstration and discussion focuses on aviation system disruptions caused by terrorist attacks, the approach can be applied to a broad range of threat scenarios.

Aviation System Disruptions↗
DOE OSTI2019

Spatial Predictions and Associated Uncertainty of Annual Soil Respiration at the Global Scale

Abstract Soil respiration (Rs), the soil‐to‐atmosphere CO 2 flux produced by microbes and plant roots, is a critical but uncertain component of the global carbon cycle. Our current understanding of the variability and dynamics is limited by the coarse spatial resolution of existing estimates. We predicted annual Rs and associated uncertainty across the world at 1‐km resolution using a quantile regression forest algorithm trained with observations from the global Soil Respiration Database spanning from 1961 to 2011. This model yielded a global annual Rs estimate of 87.9 Pg C/year with an associated global uncertainty of 18.6 (mean absolute error) and 40.4 (root mean square error) Pg C/year. The estimated annual heterotrophic respiration (Rh), derived from empirical relationships with Rs, was 49.7 Pg C/year over the same period. Predicted Rs rates and associated uncertainty varied widely across vegetation types, with the greatest predicted rates of Rs in evergreen broadleaf forests (accounting for 20.9% of global Rs). The greatest prediction uncertainties were in northern latitudes and arid to semiarid ecosystems, suggesting that these areas should be targeted in future measurement campaigns. This study provides predictions of Rs (and associated prediction uncertainty) at unprecedentedly high spatial resolution across the globe that could help constrain local‐to‐global process‐based models. Furthermore, it provides insights into the large variability of Rs and Rh across vegetation classes and identifies regions and vegetation types with poor model performance that should be prioritized for future data collection.

Warner, D. L.↗
DOE OSTI2023

Uncertainty quantification of the virial black hole mass with conformal prediction

Precise measurements of the black hole mass are essential to gain insight on the black hole and host galaxy co-evolution. A direct measure of the black hole mass is often restricted to nearest galaxies and instead, an indirect method using the single-epoch virial black hole mass estimation is used for objects at high redshifts. However, this method is subjected to biases and uncertainties as it is reliant on the scaling relation from a small sample of local active galactic nuclei. In this study, we propose the application of conformalized quantile regression (CQR) to quantify the uncertainties of the black hole predictions in a machine learning setting. We compare CQR with various prediction interval techniques and demonstrated that CQR can provide a more useful prediction interval indicator. In contrast to baseline approaches for prediction interval estimation, we show that the CQR method provides prediction intervals that adjust to the black hole mass and its related properties. That is it yields a tighter constraint on the prediction interval (hence more certain) for a larger black hole mass, and accordingly, bright and broad spectral line width source. Using a combination of neural network model and CQR framework, the recovered virial black hole mass predictions and uncertainties are comparable to those measured from the Sloan Digital Sky Survey. The code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Deep Learning and Likelihood Approaches for Viral Phylogeography Converge on the Same Answers Whether the Inference Model Is Right or Wrong

Abstract Analysis of phylogenetic trees has become an essential tool in epidemiology. Likelihood-based methods fit models to phylogenies to draw inferences about the phylodynamics and history of viral transmission. However, these methods are often computationally expensive, which limits the complexity and realism of phylodynamic models and makes them ill-suited for informing policy decisions in real-time during rapidly developing outbreaks. Likelihood-free methods using deep learning are pushing the boundaries of inference beyond these constraints. In this paper, we extend, compare, and contrast a recently developed deep learning method for likelihood-free inference from trees. We trained multiple deep neural networks using phylogenies from simulated outbreaks that spread among 5 locations and found they achieve close to the same levels of accuracy as Bayesian inference under the true simulation model. We compared robustness to model misspecification of a trained neural network to that of a Bayesian method. We found that both models had comparable performance, converging on similar biases. We also implemented a method of uncertainty quantification called conformalized quantile regression that we demonstrate has similar patterns of sensitivity to model misspecification as Bayesian highest posterior density (HPD) and greatly overlap with HPDs, but have lower precision (more conservative). Finally, we trained and tested a neural network against phylogeographic data from a recent study of the SARS-Cov-2 pandemic in Europe and obtained similar estimates of region-specific epidemiological parameters and the location of the common ancestor in Europe. Along with being as accurate and robust as likelihood-based methods, our trained neural networks are on average over 3 orders of magnitude faster after training. Our results support the notion that neural networks can be trained with simulated data to accurately mimic the good and bad statistical properties of the likelihood functions of generative phylogenetic models.

Evolutionary Biology↗
DOE OSTI2018

Probabilistic Short-Term Wind Forecasting Based on Pinball Loss Optimization

Probabilistic wind power forecasts that quantify the uncertainty in wind output have the potential to aid in the economic grid integration of wind power at large penetration levels. In this paper, a novel probabilistic wind forecasting approach based on pinball loss optimization is proposed, in conjunction with a multi-model machine learning based ensemble deterministic forecasting framework. By assuming the point-forecasted value as the mean at each point, one unknown parameter (i.e., standard deviation) of a predictive distribution at each forecasting point is determined by minimizing the pinball loss. A surrogate model is developed to represent the unknown distribution parameter as a function of deterministic forecasts. This surrogate model can be used together with deterministic forecasts to predict the unknown distribution parameter and thereby generate probabilistic forecasts. Numerical results of case studies show that the proposed method has improved the pinball loss by up to 35% compared to a baseline quantile regression forecasting model.

17 WIND ENERGY↗
DOE OSTI2023

Bi-Level Adaptive Storage Expansion Strategy for Microgrids Using Deep Reinforcement Learning

Battery energy storage (BES) is a versatile resource for the secure and economic operation of microgrids (MGs). Prevailing stochastic optimization-based approaches for BES expansion planning for MGs are computationally complicated. This work proposes a data-driven bi-level multi-period BES expansion planning framework to determine the siting, sizing, and timing of BES installations. The proposed planning framework unifies deep reinforcement learning (DRL) and linear programming, thereby decoupling the determinations for the integer and continuous decision variables in two time scales, respectively. In the upper level, a rainbow DRL agent with quantile regression is trained to provide dynamic planning policies to accommodate stochastic renewable energy resources (RESs), load, and battery price changes efficiently. Further, the lower level computes the optimal operation of MGs with frequency constraints to hedge the islanding contingency. The two levels communicate with one another by exchanging storage configuration and operating expenses in order to accomplish the shared goal of minimizing investment and operation costs. Comparative case studies on an MG are carried out to demonstrate the superiority of the proposed DRL-based solution to the mixed-integer linear programming counterpart on efficiency, scalability, and adaptability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Spatial biases of information influence global estimates of soil respiration: how can we improve global predictions?

Soil respiration (Rs), the efflux of CO2 from soils to the atmosphere, is a major component of the terrestrial carbon cycle, but is poorly constrained from regional to global scales. The global soil respiration database (SRDB) is a compilation of in-situ Rs observations from around the globe that has been consistently updated with new measurements over the past decade. It is unclear whether the addition of data to new versions of the SRDB has produced better-constrained global Rs estimates. We compared two versions of the SRDB (v3.0 n=5173 and v5.0 n=10366) to determine how additional data influenced the global sum and uncertainty of Rs as predicted by a quantile regression forest model. A model parameterized using SRDBv3 yielded a global Rs sum of 88.6 Pg C yr-1, and associated uncertainty of 29.9 (mean absolute error) and 57.9 (standard deviation) Pg C yr-1, whereas parameterization using SRDBv5 yielded 96.5 Pg C yr-1 and associated uncertainty of 30.2 (mean average error) and 73.4 (standard deviation) Pg C yr-1. Empirically estimated global heterotrophic respiration (Rh) from v3 and v5 were 49.9-50.2 (mean 50.1) and 53.3-53.5 (mean 53.4) Pg C yr-1, respectively. Interestingly, SRDBv5’s inclusion of new data from underrepresented regions (e.g., Asia, Africa, South America) thus resulted in overall higher model uncertainty. The SRDBv5 is still biased towards northern latitudes and temperate zones, so we tested that an optimized global distribution of Rs measurements resulted in a global sum of 96.4 ± 21.4 Pg C yr-1 and lower model uncertainty. These results support current global estimates of Rs but highlight spatial inconsistencies that influence model parameterization and interpretation that should be addressed to improve regional- to global-scale predictions.

Stell, Emma↗
DOE OSTI2024

Methane fluxes in tidal marshes of the conterminous United States

Abstract Methane (CH 4 ) is a potent greenhouse gas (GHG) with atmospheric concentrations that have nearly tripled since pre‐industrial times. Wetlands account for a large share of global CH 4 emissions, yet the magnitude and factors controlling CH 4 fluxes in tidal wetlands remain uncertain. We synthesized CH 4 flux data from 100 chamber and 9 eddy covariance (EC) sites across tidal marshes in the conterminous United States to assess controlling factors and improve predictions of CH 4 emissions. This effort included creating an open‐source database of chamber‐based GHG fluxes ( https://doi.org/10.25573/serc.14227085 ). Annual fluxes across chamber and EC sites averaged 26 ± 53 g CH 4 m −2 year −1 , with a median of 3.9 g CH 4 m −2 year −1 , and only 25% of sites exceeding 18 g CH 4 m −2 year −1 . The highest fluxes were observed at fresh‐oligohaline sites with daily maximum temperature normals (MATmax) above 25.6°C. These were followed by frequently inundated low and mid‐fresh‐oligohaline marshes with MATmax ≤25.6°C, and mesohaline sites with MATmax >19°C. Quantile regressions of paired chamber CH 4 flux and porewater biogeochemistry revealed that the 90th percentile of fluxes fell below 5 ± 3 nmol m −2 s −1 at sulfate concentrations >4.7 ± 0.6 mM, porewater salinity >21 ± 2 psu, or surface water salinity >15 ± 3 psu. Across sites, salinity was the dominant predictor of annual CH 4 fluxes, while within sites, temperature, gross primary productivity (GPP), and tidal height controlled variability at diel and seasonal scales. At the diel scale, GPP preceded temperature in importance for predicting CH 4 flux changes, while the opposite was observed at the seasonal scale. Water levels influenced the timing and pathway of diel CH 4 fluxes, with pulsed releases of stored CH 4 at low to rising tide. This study provides data and methods to improve tidal marsh CH 4 emission estimates, support blue carbon assessments, and refine national and global GHG inventories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

What Factors Affect the Prices of Low-Priced U.S. Solar PV Systems?

The price of solar PV systems has declined rapidly, yet there are some much lower-priced systems than others. This study explores the factors leading some systems to be so much lower priced than others. Using a data set of 42,611 residential-scale PV systems installed in the U.S. in 2013, we use quantile regressions to estimate the importance of factors affecting the installed prices for low-priced (LP) systems (those at the 10th percentile) in comparison to median-priced systems. We find that the value of solar to consumers–a variable that accounts for subsidies, electric rates, and PV generation levels–is associated with lower prices for LP systems but higher prices for median priced systems. Conversely, systems installed in new home construction are associated with lower prices at the median but higher prices for LP. Other variables have larger cost-reducing effects on LP than on median priced systems: systems installed in Arizona and Florida, as well as commercial and thin film systems. In contrast, the following have a smaller effect on prices for LP systems than median priced systems: tracking systems, self-installations, systems installed in Massachusetts, the system size, and installer experience. These results highlight the complex factors at play that lead to LP systems and shed light into how such LP systems can come about.

14 SOLAR ENERGY↗
DOE OSTI2025

Advancing Multiscale Simulation of Plasma-Surface Interfaces

We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Examples of Mission-driven Data Science from Jefferson Lab and ACES

This presentation details mission-driven data science initiatives at Jefferson Lab and the Joint Institute for Advanced Computing on Environmental Studies (ACES). JLab, a U.S. Department of Energy Office of Science national laboratory, operates the Continuous Electron Beam Accelerator Facility (CEBAF), and is the lead institute for the new High Performance Data Facility (HPDF) Hub. The Joint Institute for ACES brings together interdisciplinary teams in health informatics, climate modeling, computer science, and physics to address environmental challenges, including flood modeling. The Hampton Roads region, particularly Norfolk and Virginia Beach, faces increasing flood risks, motivating the need for rapid, reliable, and risk-aware decision support. ACES’s flooding work has a focus on uncertainty quantification (UQ) and machine learning (ML) for coastal flood management. The work is motivated by the increasing vulnerability of communities such as Norfolk and Virginia Beach, Virginia, to frequent coastal flooding events, and the need for rapid, reliable decision support. The research develops computationally efficient ML surrogate models to forecast water levels and flooding risk. A central theme is the quantification and calibration of predictive uncertainty, especially for out-of-distribution (OOD) scenarios, using techniques such as Monte Carlo Dropout, Deep Ensembles, Gaussian Processes, and Deep Quantile Regression (DQR). The study demonstrates that distance-aware UQ is critical for reliable scientific AI, particularly in high-dimensional, safety-critical, and real-time applications.

McSpadden, Diana [Thomas Jefferson National Accele↗
DOE OSTI2025

UNCERTAINTY-AWARE DEEP LEARNING FRAMEWORK FOR FORECASTING COASTAL WATER LEVEL IN VIRGINIA BEACH

Coastal areas like Virginia Beach, USA, are increasingly vulnerable to flooding. To mitigate the impact of flooding, it is crucial for the City of Virginia Beach to have reliable 72-hour-ahead (3 days) forecasts of water levels at key gauge locations. To support this effort, several sensors have been installed throughout the city to monitor water levels and other environmental parameters such as wind speed, precipitation, and atmospheric pressure. Leveraging sensor data from one of these locations, we developed an uncertainty-aware deep learning model to forecast water levels. We employed deep quantile regression (DQR) to quantify variability in the predictions and examined the performance of three different model architectures. In addition to exclusively including historical data, we investigated the improvement wind forecasts provide to the accuracy of 72-hour-ahead water level predictions. The results show a twelvefold improvement in the flood forecast for a real flooding event.

Hasan, Mahmud [Thomas Jefferson National Accelerat↗
DOE OSTI2023

Patterns, drivers, and a predictive model of dam removal cost in the United States

Given the burgeoning dam removal movement and the large number of dams approaching obsolescence in the United States, cost estimating data and tools are needed for dam removal prioritization, planning, and execution. We used the list of removed dams compiled by American Rivers to search for publicly available reported costs for dam removal projects. Total cost information could include component costs related to project planning, dam deconstruction, monitoring, and several categories of mitigation activities. We compiled reported costs from 455 unique sources for 668 dams removed in the United States from 1965 to 2020. The dam removals occurred within 571 unique projects involving 1–18 dams. When adjusted for inflation into 2020 USD, cost of these projects totaled $\$1.522$ billion, with per-dam costs ranging from $\$1$ thousand (k) to $\$268.8$ million (M). The median cost for dam removals was $\$157$k, $\$823$k, and $\$6.2$M for dams that were< 5 m, between 5–10 m, and > 10 m in height, respectively. Geographic differences in total costs showed that northern states in general, and the Pacific Northwest in particular, spent the most on dam removal. The Midwest and the Northeast spent proportionally more on removal of dams less than 5 m in height, whereas the Northwest and Southwest spent the most on larger dam removals > 10 m tall. We used stochastic gradient boosting with quantile regression to model dam removal cost against potential predictor variables including dam characteristics (dam height and material), hydrography (average annual discharge and drainage area), project complexity (inferred from construction and sediment management, mitigation, and post-removal cost drivers), and geographic region. Dam height, annual average discharge at the dam site, and project complexity were the predominant drivers of removal cost. The final model had an R 2 of 57% and when applied to a test dataset model predictions had a root mean squared error of $\$5.09$M and a mean absolute error of $\$1.45$M, indicating its potential utility to predict estimated costs of dam removal. We developed a R shiny application for estimating dam removal costs using customized model inputs for exploratory analyses and potential dam removal planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2016

The star cluster mass–galactocentric radius relation: Implications for cluster formation

Whether or not the initial star cluster mass function is established through a universal, galactocentric-distance-independent stochastic process, on the scales of individual galaxies, remains an unsolved problem. This debate has recently gained new impetus through the publication of a study that concluded that the maximum cluster mass in a given population is not solely determined by size-of-sample effects. Here, we revisit the evidence in favor and against stochastic cluster formation by examining the young (≲ a few ×10{sup 8} year old) star cluster mass–galactocentric radius relation in M33, M51, M83, and the Large Magellanic Cloud. To eliminate size-of-sample effects, we first adopt radial bin sizes containing constant numbers of clusters, which we use to quantify the radial distribution of the first- to fifth-ranked most massive clusters using ordinary least-squares fitting. We supplement this analysis with an application of quantile regression, a binless approach to rank-based regression taking an absolute-value-distance penalty. Both methods yield, within the 1σ to 3σ uncertainties, near-zero slopes in the diagnostic plane, largely irrespective of the maximum age or minimum mass imposed on our sample selection, or of the radial bin size adopted. We conclude that, at least in our four well-studied sample galaxies, star cluster formation does not necessarily require an environment-dependent cluster formation scenario, which thus supports the notion of stochastic star cluster formation as the dominant star cluster-formation process within a given galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Experimental trials of species-specific bat flight responses to an ultrasonic deterrent

Unintended consequences of increasing wind energy production include bat mortalities from wind turbine blade strikes. Ultrasonic deterrents (UDs) have been developed to reduce bat mortalities at wind turbines. Our goal was to experimentally assess the species-specific effectiveness of three emission treatments from the UD developed by NRG Systems. We conducted trials in a flight cage measuring approximately 60 m × 10 m × 4.4 m (length × width × height) from July 2020 to May 2021 in San Marcos, Texas, USA. A single UD was placed at either end of the flight cage, and we randomly selected one for each night of field trials. Trials focused on a red bat species group (Lasiurus borealis and Lasiurus blossevillii; n = 46) and four species: cave myotis (Myotis velifer; n = 57), Brazilian free-tailed bats (Tadarida brasiliensis; n = 73), evening bats (Nycteceius humeralis; n = 53), and tricolored bats (Perimyotis subflavus; n = 17). The trials occurred during three treatment emissions: low (emissions from subarrays at 20, 26, and 32 kHz), high (emissions from subarrays at 38, 44, and 50 kHz), and combined (all six emission frequencies). We placed one wild-captured bat into the flight cage for each trial, which consisted of an acclimation period, a control period with the UD powered off, and the three emission treatments (order randomly selected), each interspersed with a control period. We tracked bat flight using four thermal cameras placed outside the flight cage. We quantified the effectiveness of each treatment by comparing the distances each bat flew from the UD during each treatment vs. the control period using quantile regression. Additionally, we conducted an exploratory analysis of differences between sex and season and sex within season using analysis of variance. Broadly, UDs were effective at altering the bats’ flight paths as they flew farther from the UD during treatments than during controls; however, results varied by species, sex, season, and sex within season. For the red bat group, bats flew farther from the UD during all treatments than during the control period at all percentiles (p < 0.001), and treatments were comparable in effectiveness. For cave myotis, all percentile distances were farther from the UD during each of the treatments than during the control, except the 90th percentile distance during high, and low was most effective. For evening bats and Brazilian free-tailed bats, results were inconsistent, but high and low were most effective, respectively. For tricolored bats, combined and low were significant at the 10th–75th percentiles, high was significant at all percentiles, and combined was most effective. Results suggest UDs may be an effective means of reducing bat mortalities due to wind turbine blade strikes. We recommend that continued research on UDs focus on low emission treatments, which have decreased sound attenuation and demonstrated effectiveness across the bat species evaluated in this study.

17 WIND ENERGY↗
NASA NTRS2017

Statistical Approaches for the Definition of Landslide Rainfall Thresholds and their Uncertainty Using Rain Gauge and Satellite Data

Models for forecasting rainfall-induced landslides are mostly based on the identification of empirical rainfall thresholds obtained exploiting rain gauge data. Despite their increased availability, satellite rainfall estimates are scarcely used for this purpose. Satellite data should be useful in ungauged and remote areas, or should provide a significant spatial and temporal reference in gauged areas. In this paper, the analysis of the reliability of rainfall thresholds based on rainfall remote sensed and rain gauge data for the prediction of landslide occurrence is carried out. To date, the estimation of the uncertainty associated with the empirical rainfall thresholds is mostly based on a bootstrap resampling of the rainfall duration and the cumulated event rainfall pairs (D,E) characterizing rainfall events responsible for past failures. This estimation does not consider the measurement uncertainty associated with D and E. In the paper, we propose (i) a new automated procedure to reconstruct ED conditions responsible for the landslide triggering and their uncertainties, and (ii) three new methods to identify rainfall threshold for the possible landslide occurrence, exploiting rain gauge and satellite data. In particular, the proposed methods are based on Least Square (LS), Quantile Regression (QR) and Nonlinear Least Square (NLS) statistical approaches. We applied the new procedure and methods to define empirical rainfall thresholds and their associated uncertainties in the Umbria region (central Italy) using both rain-gauge measurements and satellite estimates. We finally validated the thresholds and tested the effectiveness of the different threshold definition methods with independent landslide information. The NLS method among the others performed better in calculating thresholds in the full range of rainfall durations. We found that the thresholds obtained from satellite data are lower than those obtained from rain gauge measurements. This is in agreement with the literature, where satellite rainfall data underestimate the 'ground' rainfall registered by rain gauges.

landslide prediction↗