Atmospheric rivers (ARs) are extratropical storms that produce extreme precipitation on the west coasts of the world’s major landmasses. In the United States, ARs cause significant flooding, yet their economic impacts have not been quantified. Here, using 40 years of data from the National Flood Insurance Program, we show that ARs are the primary drivers of flood damages in the western United States. Using a recently developed AR scale, which varies from category 1 to 5, we find that flood damages increase exponentially with AR intensity and duration: Each increase in category corresponds to a roughly 10-fold increase in damages. Category 4 and 5 ARs cause median damages in the tens and hundreds of millions of dollars, respectively. Rising population, increased development, and climate change are expected to worsen the risk of AR-driven flood damage in future decades.

[1] A new one-dimensional infiltration and redistribution method is proposed as an alternative to the Richards equation (RE) for coupled surface and subsurface models. The proposed method discretizes soil water content into hypothetical hydraulically interacting bins. The entry and propagation of displacement fronts in each bin are simulated by means of explicit infiltration and drainage approximations based on capillary and gravitational driving forces. Wetting front advances within bins create water deficits that are satisfied by capillary-driven interbin flow. The method inherently provides numerical stability by precluding the need to directly estimate nonlinear gradients through numerical schemes. Comparisons of the performance of this method against RE solutions for theoretical, laboratory, and field data in both well-drained and high water table conditions are presented. The new method produces infiltration flux estimates with errors less than 10% compared to a widely used RE solution in tests on multiple soil textures with and without high water table conditions while providing unconditionally guaranteed conservation of mass.Citation: Talbot, C. A., and F. L. Ogden (2008), A method for computing infiltration and redistribution in a discretized moisture content domain, Water Resour. Res., 44, W08453,

We have developed an alternative to the one-dimensional partial differential equation (PDE) attributed to Richards (1931) that describes unsaturated porous media flow in homogeneous soil layers. Our solution is a set of three ordinary differential equations (ODEs) derived from unsaturated flux and mass conservation principles. We used a hodograph transformation, the Method of Lines, and a finite water-content discretization to produce ODEs that accurately simulate infiltration, falling slugs, and groundwater table dynamic effects on vadose zone fluxes. This formulation, which we refer to as ''finite water-content'', simulates sharp fronts and is guaranteed to conserve mass using a finite-volume solution. Our ODE solution method is explicitly integrable, does not require iterations and therefore has no convergence limits and is computationally efficient. The method accepts boundary fluxes including arbitrary precipitation, bare soil evaporation, and evapotranspiration. The method can simulate heterogeneous soils using layers. Results are presented in terms of fluxes and water content profiles. Comparing our method against analytical solutions, laboratory data, and the Hydrus-1D solver, we find that predictive performance of our finite water-content ODE method is comparable to or in some cases exceeds that of the solution of Richards' equation, with or without a shallow water table. The presented ODE method is transformative in that it offers accuracy comparable to the Richards (1931) PDE numerical solution, without the numerical complexity, in a form that is robust, continuous, and suitable for use in large watershed and land-atmosphere simulation models, including regional-scale models of coupled climate and hydrology.