Search

High resolution (< 1 km) atmospheric modeling is increasingly used to study precipitation distributions in complex terrain and cryosphere-atmospheric processes. While this approach has yielded insightful results, studies over annual time-scales or at the spatial extents of watersheds remain unrealistic due to the computational costs of running most atmospheric models. In this paper we introduce a High-resolution variant of the Intermediate Complexity Atmospheric Research (ICAR) model, HICAR. We detail the model development that enabled HICAR simulations at the hectometer scale, including changes to the advection scheme and the wind solver. The latter uses near surface terrain parameters which allow HICAR to simulate complex topographic flow features. These model improvements clearly influence precipitation distributions at the ridge scale (50 m), suggesting that HICAR can approximate processes dependent on particle-flow interactions such as preferential deposition. A 250 m HICAR simulation over most of the Swiss Alps also shows monthly precipitation patterns similar to two different gridded precipitation products which assimilate available observations. Benchmarking runs show that HICAR uses 118x fewer computational resources than the WRF atmospheric model. This gain in efficiency makes dynamic downscaling accessible to ecohydrological research, where downscaled data is often required at hectometer resolution for whole basins at seasonal time scales. These results motivate further development of HICAR, including refinement of parameterizations used in the wind solver, and coupling of the model with an intermediate complexity snow model.

Rainwater is a vital resource and dynamic driver of terrestrial ecosystems. Yet, processes controlling precipitation inputs and interactions during storms are often poorly seen, and poorly sensed when direct observations are substituted with technological ones. We discuss how human observations complement technological ones, and the benefits of scientists spending more time in the storm. Human observation can reveal ephemeral storm-related phenomena such as biogeochemical ‘hot moments’, organismal responses, and sedimentary processes which can then be explored in greater resolution using sensors and virtual experimentation. Storm-related phenomena trigger lasting, oversized impacts on hydrologic and biogeochemical processes, organismal traits/functions, and ecosystem services. We provide examples of phenomena in forests, across disciplines and scales, to inspire mindful, holistic observation of ecosystems during storms. We conclude that technological observations alone are insufficient to trace the process complexity and unpredictability of fleeting biogeochemical or ecological events without the “shower thoughts” produced by scientists’ human sensory and cognitive systems during storms.

The snowpack over the Mediterranean mountains constitutes a key water resource for the downstream populations. However, its dynamics have not been studied in detail yet in many areas, mostly because of the scarcity of snowpack observations. In this work, we present a characterization of the snowpack over the two mountain ranges of Lebanon. To obtain the necessary snowpack information, we have developed a 1 km regional-scale snow reanalysis (ICAR_assim) covering the period 2010–2017. ICAR_assim was developed by means of an ensemble-based data assimilation of Moderate Resolution Imaging Spectroradiometer (MODIS) fractional snow-covered area (fSCA) through an energy and mass snow balance model, the Flexible Snow Model (FSM2), using the particle batch smoother (PBS). The meteorological forcing data were obtained by a regional atmospheric simulation from the Intermediate Complexity Atmospheric Research model (ICAR) nested inside a coarser regional simulation from the Weather Research and Forecasting model (WRF). The boundary and initial conditions of WRF were provided by the ERA5 atmospheric reanalysis. ICAR_assim showed very good agreement with MODIS gap-filled snow products, with a spatial correlation of R=0.98 in the snow probability (P(snow)) and a temporal correlation of R=0.88 on the day of peak snow water equivalent (SWE). Similarly, ICAR_assim has shown a correlation with the seasonal mean SWE of R=0.75 compared with in situ observations from automatic weather stations (AWSs). The results highlight the high temporal variability in the snowpack in the Lebanese mountain ranges, with the differences between Mount Lebanon and the Anti-Lebanon Mountains that cannot only be explained by hypsography as the Anti-Lebanon Mountains are in the rain shadow of Mount Lebanon. The maximum fresh water stored in the snowpack is in the middle elevations, approximately between 2200 and 2500 m a.s.l. (above sea level). Thus, the resilience to further warming is low for the snow water resources of Lebanon due to the proximity of the snowpack to the zero isotherm.

This study examines current and future western U.S. snowfall and snowpack through current and future climate simulations with a 4-km horizontal grid spacing cloud permitting regional climate model over the entire CONtinental U.S. for a 13-year period between 2001 and 2013. At this horizontal resolution, the spatiotemporal distribution of the orographic snowfall and snowpack is well captured partly due to the ability of the model to realistically simulate mesoscale and microphysical features such as orographically induced updrafts driving clouds and precipitation. The historical simulation well captures the observed snowfall and snowpack amounts and pattern in the western U.S. The future climate simulation uses the Pseudo-Global Warming approach, taking the climate change signal from CMIP5 multi-model ensemble-mean difference between 2070–2099 and 1976–2005. The results show that the thermodynamic impacts of climate change in the western U.S. can be characterized considering mountain ranges in two distinct geographic regions: the mountain ranges close to the Pacific Ocean (coastal ranges) and those in the inter-mountain west. Climate change out to 2100 significantly impacts all aspects of the water cycle, with pronounced climate change response in the coastal ranges. A notable result is that the snowpack in the Pacific Northwest is predicted to decrease by ~ 70% by 2100. Trends of this magnitude have already been observed in the historical data and in previous studies. The current Pseudo Global Warming future climate simulation and previous global climate simulations all suggest that these trends will continue to the point that most snowpack will be gone by 2100 in the Pacific Northwest for the most aggressive RCP8.5 climate scenario, even if annual precipitation increases by 10%. Future work will focus on extending the current convective permitting results to a full climate change simulation allowing for dynamical changes in the flow.

Snow models rely on accurate meteorological input data at the spatial scales at which they operate. However, even the highest resolution operational atmospheric models often run at horizontal resolutions at least an order of magnitude coarser than most snow models. Different downscaling techniques can be employed to bridge this scale gap, typically being sorted into either statistical or dynamical techniques. Statistical techniques often rely on temporally invariant spatial patterns or simplistic conceptual relationships, making them computationally cheap but prone to errors. Dynamical downscaling generally offers a counterpoint to this tradeoff: stronger physical basis but more computational demand. Efforts have been made to optimize this tradeoff of dynamic downscaling, reducing computational demand while maintaining physical accuracy of predicted variables as well as the interdependency of downscaled variables such as winds and precipitation. The Intermediate Complexity Atmospheric Research (ICAR) model recently demonstrated an ability to match precipitation patterns from WRF, but with computational costs at least two orders of magnitude lower. While promising, these results from a 4km comparison did not translate to finer spatial resolutions often needed as input to snow models.Thus, we introduce the High-resolution Intermediate Complexity Atmospheric Research Model (HICAR), a new variant of the ICAR model developed for spatial resolutions as high as 50m. Relative to a traditional atmospheric model like WRF, HICAR maintains the orders-of-magnitude reduction in computational demand which ICAR displayed, while resolving terrain-induced effects on the wind field not seen in ICAR. This is achieved through a novel combination of adjustments to a background wind field based on terrain descriptors with a wind solver. The solver enforces a mass-conservation constraint on the 3D wind field. These modifications successfully mimic dynamic effects such as flow blocking, ridge-crest speed up, and lee-side recirculation to be captured in the resulting wind field. These features are of particular importance for resolving snow deposition patterns, where the snow particles are particularly susceptible to advection by the near-surface flow field. We validate the accuracy of HICAR’s flow features using a wind LiDAR deployed in complex terrain and show a comparison between flow fields from HICAR and WRF at a horizontal resolution of 50 m. These comparisons demonstrate HICAR’s ability to resolve terrain-induced modifications to the flow field which result in increased heterogeneity of ridge-scale snowfall patterns. To this point, preliminary comparisons of snow deposition patterns in complex terrain between the HICAR and WRF models are offered. With this new model, physically-based downscaling of precipitation and other atmospheric variables which preserves their interdependencies is made available for high-resolutions (100m) and large-spatial extents (10,000 km2) which are often demanded by operational land-surface models.

Accurate knowledge of snow depth distributions in mountain catchments is critical for applications in hydrology and ecology. Recently, a method was proposed to map snow depth at meter-scale resolution from very-high-resolution stereo satellite imagery (e.g., Pléiades) with an accuracy close to 0.5 m. However, the validation was limited to probe measurements and unmanned aircraft vehicle (UAV) photogrammetry, which sampled a limited fraction of the topographic and snow depth variability. We improve upon this evaluation using accurate maps of the snow depth derived from Airborne Snow Observatory laser-scanning measurements in the Tuolumne river basin, USA. We find a good agreement between both datasets over a snow-covered area of 138 km2 on a 3 m grid, with a positive bias for a Pléiades snow depth of 0.08 m, a root mean square error of 0.80 m and a normalized median absolute deviation (NMAD) of 0.69 m. Satellite data capture the relationship between snow depth and elevation at the catchment scale and also small-scale features like snow drifts and avalanche deposits at a typical scale of tens of meters. The random error at the pixel level is lower in snow-free areas than in snow-covered areas, but it is reduced by a factor of 2 (NMAD of approximately 0.40 m for snow depth) when averaged to a 36 m grid. We conclude that satellite photogrammetry stands out as a convenient method to estimate the spatial distribution of snow depth in high mountain catchments., The Cryosphere, 14 (9), ISSN:1994-0416, ISSN:1994-0424

In mountain terrain, well-configured high-resolution atmospheric models are able to simulate total annual rain and snowfall better than spatial estimates derived from in situ observational networks of precipitation gauges, and significantly better than radar or satellite-derived estimates. This conclusion is primarily based on comparisons with streamflow and snow in basins across the western United States and in Iceland, Europe, and Asia. Even though they outperform gridded datasets based on gauge networks, atmospheric models still disagree with each other on annual average precipitation and often disagree more on their representation of individual storms. Research to address these difficulties must make use of a wide range of observations (snow, streamflow, ecology, radar, satellite) and bring together scientists from different disciplines and a wide range of communities.

Research in meteorological prediction on sub-seasonal to seasonal (S2S) timescales has seen growth in recent years. Concurrent with this, demand for seasonal drought forecasting has risen. While there is obvious synergy between these fields, S2S meteorological forecasting has typically focused on low resolution global models, while the development of drought can be sensitive to the local expression of weather anomalies and their interaction with local surface properties and processes. This suggests that downscaling might play an important role in the application of meteorological S2S forecasts to skillful forecasting of drought. Here, we apply the Generalized Analog Regression Downscaling (GARD) algorithm to downscale meteorological hindcasts from the NASA Goddard Earth Observing System (GEOS) global S2S forecast system. Downscaled meteorological fields are then applied to drive offline simulations with the Catchment Land Surface Model (CLSM) to forecast United States Drought Monitor (USDM) style drought indicators derived from simulated surface hydrology variables. We compare the representation of drought in these downscaled hindcasts to hindcasts that are not downscaled, using the North American Land Data Assimilation System Phase 2 (NLDAS-2) dataset as an observational reference. We find that downscaling using GARD improves hindcasts of temperature and temperature anomalies, but the results for precipitation are mixed and generally small. Overall, GARD downscaling led to improved hindcast skill for total drought across the Contiguous United States (CONUS), and improvements were greatest for extreme (D3) and exceptional (D4) drought categories.

Water managers are actively incorporating climate change information into their long- and short-term planning processes. This is generally seen as a step in the right direction because it supplements traditional methods, providing new insights that can help in planning for a non-stationary climate. However, the continuous evolution of climate change information can make it challenging to use available information appropriately. Advice on how to use the information is not always straightforward and typically requires extended dialogue between information producers and users, which is not always feasible. To help navigate better the ever-changing climate science landscape, this review is organized as a set of nine guidelines for water managers and planners that highlight better practices for incorporating climate change information into water resource planning and management. Each DOs and DON'Ts recommendation is given with context on why certain strategies are preferable and addresses frequently asked questions by exploring past studies and documents that provide guidance, including real-world examples mainly, though not exclusively, from the United States. This paper is intended to provide a foundation that can expand through continued dialogue within and between the climate science and application communities worldwide, a two-way information sharing that can increase the actionable nature of the information produced and promote greater utility and appropriate use.

The evaluation of models in general is a nontrivial task and can, due to epistemological and practical reasons, never be considered complete. Due to this incompleteness, a model may yield correct results for the wrong reasons, i.e., via a different chain of processes than found in observations. While guidelines and strategies exist in the atmospheric sciences to maximize the chances that models are correct for the right reasons, these are mostly applicable to full physics models, such as numerical weather prediction models. The Intermediate Complexity Atmospheric Research (ICAR) model is an atmospheric model employing linear mountain wave theory to represent the wind field. In this wind field, atmospheric quantities such as temperature and moisture are advected and a microphysics scheme is applied to represent the formation of clouds and precipitation. This study conducts an in-depth process-based evaluation of ICAR, employing idealized simulations to increase the understanding of the model and develop recommendations to maximize the probability that its results are correct for the right reasons. To contrast the obtained results from the linear-theory-based ICAR model to a full physics model, idealized simulations with the Weather Research and Forecasting (WRF) model are conducted. The impact of the developed recommendations is then demonstrated with a case study for the South Island of New Zealand. The results of this investigation suggest three modifications to improve different aspects of ICAR simulations. The representation of the wind field within the domain improves when the dry and the moist Brunt–Väisälä frequencies are calculated in accordance with linear mountain wave theory from the unperturbed base state rather than from the time-dependent perturbed atmosphere. Imposing boundary conditions at the upper boundary that are different to the standard zero-gradient boundary condition is shown to reduce errors in the potential temperature and water vapor fields. Furthermore, the results show that there is a lowest possible model top elevation that should not be undercut to avoid influences of the model top on cloud and precipitation processes within the domain. The method to determine the lowest model top elevation is applied to both the idealized simulations and the real terrain case study. Notable differences between the ICAR and WRF simulations are observed across all investigated quantities such as the wind field, water vapor and hydrometeor distributions, and the distribution of precipitation. The case study indicates that the precipitation maximum calculated by the ICAR simulation employing the developed recommendations is spatially shifted upwind in comparison to an unmodified version of ICAR. The cause for the shift is found in influences of the model top on cloud formation and precipitation processes in the ICAR simulations. Furthermore, the results show that when model skill is evaluated from statistical metrics based on comparisons to surface observations only, such an analysis may not reflect the skill of the model in capturing atmospheric processes like gravity waves and cloud formation.

Snow deposition patterns in complex terrain are heavily dependent on the underlying topography. This topography affects precipitating clouds at the kilometer-scale and causes changes to the wind field at the sub-kilometer scale, resulting in altered advection of falling hydrometeors. Snow particles are particularly sensitive to changes in the near-surface flow field due to their low density. Atmospheric models which run at the kilometer scale cannot resolve the actual heterogeneity of the underlying terrain, resulting in precipitation maps which do not capture terrain-affected precipitation patterns. Thus, snow-atmosphere interactions such as preferential deposition are often not resolved in precipitation data used as input to snow models. To bridge this spatial gap and resolve snow-atmosphere interactions at the sub-kilometer scale, we couple an intermediate complexity atmospheric model (ICAR) to the COSMO NWP model. Applying this model to sub-kilometer terrain (horizontal resolution of 50 and 250 m) required changes to ICAR’s computational grid, atmospheric dynamics, and boundary layer flow. As a result, the near-surface flow now accounts for surface roughness and topographically induced speed up. This has been achieved by using terrain descriptors calculated once at initialization which consider a point’s exposure or sheltering relative to surrounding terrain. In particular, the use of a 3-dimensional Sx parameter allows us to simulate areas of stagnation and recirculation on the lee of terrain features. Our approach maintains the accurate large-scale precipitation patterns from COSMO but resolves the dynamics induced by terrain at the sub-kilometer scale without adding additional computational burden. We find that solid precipitation patterns at the ridge scale, such as preferential deposition of snow, are better resolved in the high-resolution version of ICAR than the current ICAR or COSMO models. This updated version of ICAR presents a new tool to dynamically downscale NWP output for snow models and enables future studies of snow-atmosphere interactions at domain scales of 100’s of kilometers.

The Mississippi River Basin (MRB) is a flash flood hotspot in the United States, receiving the most frequent floods and highest rainfall accumulations across the country. In a future warmer climate, this region exhibits some of the greatest increases in rainfall associated with storms that produce flash floods. In order to better understand these future changes, convection-permitting simulations of a current and future climate are utilized to study changes to storm dynamics and precipitation in these convectively-driven flash flood-producing storms. First, nearly 500 flash flood-producing storms in the MRB are examined under a pseudo-global warming framework to examine the role of vertical velocity in modulating future rainfall changes. Three different categories of storms are designated based on their vertical velocity magnitude in the current climate–weak, moderate, and strong. While all storm categories display an increase in future rainfall accumulation, the amount of increase varies by the storm’s vertical velocity magnitude, which also changes in the future. Second, idealized WRF simulations are run based on a composite sounding of the flash flood-producing storms in the MRB that occurred during the warm season. Future temperature, moisture, and horizontal wind perturbations are added to the initial sounding using the CESM Large Ensemble Data Set under the RCP 8.5 emissions scenario. In these idealized simulations, the contribution of different precipitation modes to future changes in rainfall are examined. The relationship between changes in future precipitation mode and storm dynamics provides a better understanding of how storm processes influence future changes in rainfall in a flash flood prone region in the United States.

When formulating a hydrologic model, scientists rely on parameterizations of multiple processes based on field data, but literature review suggests that more frequently people select parameterizations that were included in pre-existing models rather than re-evaluating the underlying field experiments. Problems arise when limited field data exist, when “trusted” approaches do not get reevaluated, and when processes fundamentally change in different environments. The physics and dynamics of snow interception by conifers, including both loading and unloading of snow, is just such a case. The most commonly used interception parameterization is based on data from four trees from one site, but field study results are not directly transferable between environments. The process varies dramatically between locations with relatively warmer versus colder winters. Here, we combine a comprehensive literature review with a model to demonstrate essential improvements to model representations of snow interception. We recommend that, as a first and essential step, all models include increased loading due to increased adhesion and cohesion when temperatures rise from -3 and 0°C. The commonly used parameters of a fixed maximum value for loading and an e-folding time for unloading are not supported by observations or physical understanding and are not necessary to reproduce observations. In addition to unloading based on physical processes, such as wind or canopy warming, all models must represent melting of in-canopy snow so that it can be unloaded in liquid form. As a second step, we propose field experiments across climates and forest types to investigate: a) a representation of the force balance between adhesion and cohesion versus gravity for both interception efficiency and rates of unloading, b) wind effects during and between storms, and c) lubrication when snow melts. For greatest impact, this framework requires dedicated field measurements. These processes are essential for models to accurately represent the impacts of dynamically changing forest cover and snow cover on both global albedo and water supplies.

The Arctic has been warming faster than the global average during recent decades, and trends are projected to continue through the twenty-first century. Analysis of climate change impacts across the Arctic using dynamical models has almost exclusively been limited to outputs from global climate models or coarser regional climate models. Coarse resolution simulations limit the representation of physical processes, particularly in areas of complex topography and high land-surface heterogeneity. Here, current climate reference and future regional climate model simulations based on the RCP8.5 scenario over Alaska at 4 km grid spacing are compared to identify changes in snowfall and snowpack. In general, results show increases in total precipitation, large decreases in snowfall fractional contribution over 30% in some areas, decreases in snowpack season length by 50–100 days in lower elevations and along the southern Alaskan coastline, and decreases in snow water equivalent. However, increases in snowfall and snowpack of sometimes greater than 20% are evident for some colder northern areas and at the highest elevations in southern Alaska. The most significant changes in snow cover and snowfall fractional contributions occur during the spring and fall seasons. Finally, the spatial pattern of winter temperatures above freezing has small-scale spatial features tied to the topography. Such areas would not be resolved with coarser resolution regional or global climate model simulations.

This work added semi-Lagrangian convected air particles to the Intermediate Complexity Atmospheric Research (ICAR) model. The ICAR model is a simplified atmospheric model using quasi-dynamical downscaling to gain performance over more traditional atmospheric models. The ICAR model uses Fortran coarrays to split the domain amongst images and handle the halo region communication of the image’s boundary regions. The newly implemented convected air particles use trilinear interpolation to compute initial properties from the Eulerian domain and calculate humidity and buoyancy forces as the model runs. This paper investigated the performance cost and scaling attributes of executing unsaturated and saturated air particles versus the original particle-less model. An in-depth analysis was done on the communication patterns and performance of the semi-Lagrangian air particles, as well as the performance cost of a variety of initial conditions such as wind speed and saturation mixing ratios. This study found that given a linear increase in the number of particles communicated, there is an initial decrease in performance, but that it then levels out, indicating that over the runtime of the model, there is an initial cost of particle communication, but that the computational benefits quickly offset it. The study provided insight into the number of processors required to amortize the additional computational cost of the air particles.

The properties of hurricanes directly influence storm surges; however, the implications of projected changes to the climate are unclear. Here, we simulate the storm surges of historical storms under present day and end of century climate scenarios to assess the impact of climate change on storm surge inundation. We simulate 21 storms that impacted the Gulf of Mexico and Atlantic Coasts of the continental U.S. from 2000 to 2013. We find that the volume of inundation increases for 14 storms and the average change for all storms is +36%. The extent of inundation increases for 13 storms, and the average change for all storms is +25%. Notable increases in inundation occur near Texas, Louisiana, Mississippi, the west coast of Florida, the Carolinas, and New Jersey. Our calculations of inundation volume and extent suggest that at the end of the century, we can expect hurricanes to produce larger storm surge magnitudes in concentrated areas, as opposed to surges with lower magnitudes that are widespread. We examine changes in maximum wind speed, minimum central pressure, translation speed, and radius of the 33 ms−1 wind to assess the impacts of individual storm characteristics on storm surge. We find that there is no single storm characteristic that directly relates to storm surge inundation or its climate induced changes. Even when all the characteristics are considered together, the resulting influences are difficult to anticipate. This is likely due to the complexity of the hydrodynamics and interactions with local geography. This illustrates that even as climate change research advances and more is known about projected impacts to hurricanes, implications for storm surge will be difficult to predict without explicit numerical simulation.

The Snow Ensemble Uncertainty Project (SEUP) is an effort to establish a baseline characterization of snow water equivalent (SWE) uncertainty across North America with the goal of informing global snow observational needs. An ensemble-based modeling approach, encompassing a suite of current operational models, is used to assess the uncertainty in SWE and total snow storage (SWS) estimation over North America during the 2009&ndashl2017 period. The highest modeled SWE uncertainty is observed in mountainous regions, likely due to the relatively deep snow, forcing uncertainties, and variability between the different models in resolving the snow processes over complex terrain. This highlights a need for high-resolution observations in mountains to capture the high spatial SWE variability. The greatest SWS is found in Tundra regions where even though the spatiotemporal variability in modeled SWE is low, there is considerable uncertainty in the SWS estimates due to the large areal extent over which those estimates are spread. This highlights the need for high accuracy in snow estimations across the Tundra. In mid-latitude boreal forests, large uncertainties in both SWE and SWS indicate that vegetation-snow impacts are a critical area where focused improvements to modeled snow estimation efforts need to be made. Finally, the SEUP results indicate that SWE uncertainty is driving runoff uncertainty and measurements may be beneficial in reducing uncertainty in SWE and runoff, during the melt season at high latitudes (e.g., Tundra and Taiga regions) and in the Western mountain regions, whereas observations at (or near) peak SWE accumulation are more helpful over the mid-latitudes.

The snowpack over the Mediterranean mountains constitutes a key water resource for the downstream populations. However, its dynamics have not been studied in detail yet in many areas, mostly because of the scarcity of snowpack observations. In this work, we present a characterization of the snowpack over the two mountain ranges of Lebanon. To obtain the necessary snowpack information, we have developed a 1 km regional scale snow reanalysis (ICAR_assim) covering the period 2010–2017. ICAR_assim was developed by means of ensemble-based data assimilation of MODIS fractional snow-covered area (fSCA) through the energy and mass balance model the Flexible Snow Model (FSM2), using the Particle Batch Smoother (PBS). The meteorological forcing data was obtained by a regional atmospheric simulation developed through the Intermediate Complexity Atmospheric Research model (ICAR) nested inside a coarser regional simulation developed by the Weather Research and Forecasting model (WRF). The boundary and initial conditions of WRF were provided by the ERA5 atmospheric reanalysis. ICAR_assim showed very good agreement with MODIS gap-filled snow products, with a spatial correlation of R = 0.98 in the snow probability (P(snow)), and a temporal correlation of R = 0.88 in the day of peak snow water equivalent (SWE)Similarly, ICAR_assim has shown a correlation with the seasonal mean SWE of R = 0.75 compared with in-situ observations from Automatic Weather Stations (AWS). The results highlight the high temporal variability of the snowpack in the Lebanon ranges, with differences between Mount Lebanon and Anti-Lebanon that cannot be only explained by its hypsography been Anti-Lebanon in the rain shadow of Mount Lebanon. The maximum fresh water stored in the snowpack is in the middle elevations approximately between 2200 and 2500 m. a.s.l. Thus, the resilience to further warming is low for the snow water resources of Lebanon due to the proximity of the snowpack to the zero isotherm.

When is good enough, good enough? The spatio-temporal variability of precipitation makes measurements extremely challenging, particularly in the mountains. Simultaneously, the improvements in physical realism of atmospheric models makes them increasingly valuable for fields such as hydrology, particularly in the mountains. However, the computational cost of such models renders them impractical for many applications, in or out of the mountains. Here we describe an intermediate complexity atmospheric model (ICAR) capable of capturing around 90% of the variability in orographic precipitation for 1% of the computational cost of a state of the science non-hydrostatic atmospheric simulation. ICAR uses an analytical solution for flow perturbations created by topography and simulates the core processes responsible for orographic precipitation (e.g. orographic lifting, advection, cloud microphysical processes). We show that key aspects of orographic precipitation spatial patterns are well simulated in ICAR, including some that gridded observation based products are missing. We then show some early results when using ICAR to simulate regional climate changes forced by global models at higher spatial resolutions than it is currently practical to run traditional regional climate models. These simulations quantify plausible shifts in precipitation resulting in the transition from snow to rain, as well as elevation dependent warming caused by the snow albedo feedback. Further, the computational efficiency of ICAR permits us to run these simulations with many different physics configurations to better explore the sensitivity of these changes to assumptions in the microphysics and land surface model components.

We aim to discuss topics and questions raised by the recent book published under this presentation's title. The book presents research on precipitation partitioning processes in vegetated ecosystems, putting them into a global context. It describes the processes by which meteoric water comes into contact with the vegetation's canopy, typically the first surface contact of precipitation on land. It also discusses how precipitation partitioning by vegetation impacts the amount, patterning, and chemistry of water reaching the surface, as well as the amount and timing of evaporative return to the atmosphere. Although this process has been extensively studied, this is the first review of the global literature on the partitioning of precipitation by forests, shrubs, crops, grasslands and other less-studies plant types.

An accurate knowledge of snow depth distribution in mountain catchments is critical for applications in hydrology and ecology. A recent new method was proposed to map the snow depth at meter-scale resolution from very-high resolution stereo satellite imagery (e.g., Pléiades) with an accuracy close to 0.50 m. However, the validation was mainly done using probe measurements which sampled a limited fraction of the topographic and snow depth variability. We deepen this evaluation using accurate maps of the snow depth derived from ASO airborne lidar measurements in the Tuolumne river basin, USA. We find a good agreement between both datasets over a snow-covered area of 137 km2 on a 3 m grid with a positive bias for Pléiades snow depth of 0.08 m, a root-mean-square error of 0.80 m and a normalized median absolute deviation of 0.69 m. Satellite data capture the relationship between snow depth and elevation at the catchment scale, and also small-scale features like snow drifts and avalanche deposits. The random error on snow depth can be reduced by a factor two (up to approximately 0.40 m) when the snow depth map is spatially averaged to a ~ 20 m grid. The random error at the pixel level is lower on snow-free areas than on snow-covered areas, but errors on both terrain type converge at coarser resolutions, which is important for further applications of the method in areas without snow depth reference data. We conclude that satellite photogrammetry stands out as an efficient method to estimate the spatial distribution of snow depth in high mountain catchments.

The partitioning of precipitation by vegetation canopies into throughfall and evaporation can have a large effect on water availability for both ecosystems and human consumption. Canopies are often the first point of contact between precipitation and the land surface, yet the complexity and inaccessibility of this interface make measurements extremely difficult, and as a result our attempts to understand and model the relevant processes are only weakly constrained. Here, we describe common approaches to modeling canopy interception processes, particularly in global models. These models use a wide variety of parameterization and parameters internally, suggesting that we do not have a good understanding of how to model canopy interception on a global scale. To begin to quantify how big a problem this may be, we present a few ideal model experiments exploring a range of modeling approaches and assumptions to document what effect these choices have on the projection of changes in inputs to the eco-hydrologic system. Perhaps unsurprisingly, these effects vary with vegetation type and density, as well as precipitation type, intensity, and changes in precipitation. In many cases, these effects are likely dwarfed by the uncertainties in predicted changes in regional climate, but not accounting for our lack of certainty in canopy processes can lead to an overconfident, and in some cases likely incorrect, projection of future changes in water availability. This should serve as a call to action for those studying canopy interception processes to better document and consider how theories can be put into a numerical framework.

The interception of precipitation by vegetation has important consequences for climate and water resources. Although canopy interception has been studied for centuries, many fundamental unknowns remain. We present persistent questions that reflect challenges in measuring, representing, and understanding how terrestrial ecosystems intercept, partition, and transport precipitation—down to soils or back to the atmosphere. In summary of this book, we outline future needs and simultaneously provide a primer for those interested in precipitation interception processes.

Tropical cyclones have enormous costs to society through both loss of life and damage to infrastructure. There is good reason to believe that such storms will change in the future as a result of changes in the global climate system and that such changes may have important socioeconomic implications. Here a high-resolution regional climate modeling experiment is presented using the Weather Research and Forecasting (WRF) Model to investigate possible changes in tropical cyclones. These simulations were performed for the period 2001–13 using the ERA-Interim product for the boundary conditions, thus enabling a direct comparison between modeled and observed cyclone characteristics. The WRF simulation reproduced 30 of the 32 named storms that entered the model domain during this period. The model simulates the tropical cyclone tracks, storm radii, and translation speeds well, but the maximum wind speeds simulated were less than observed and the minimum central pressures were too large. This experiment is then repeated after imposing a future climate signal by adding changes in temperature, humidity, pressure, and wind speeds derived from phase 5 of the Coupled Model Intercomparison Project (CMIP5). In the current climate, 22 tracks were well simulated with little changes in future track locations. These simulations produced tropical cyclones with faster maximum winds, slower storm translation speeds, lower central pressures, and higher precipitation rates. Importantly, while these signals were statistically significant averaged across all 22 storms studied, changes varied substantially between individual storms. This illustrates the importance of using a large ensemble of storms to understand mean changes.

A growing body of work suggests that the extreme weather events that drive inland flooding are likely to increase in frequency and magnitude in a warming climate, thus potentially increasing flood damages in the future. We use hydrologic projections based on the Coupled Model Intercomparison Project Phase 5 (CMIP5) to estimate changes in the frequency of modeled 1 % annual exceedance probability (1 % AEP, or 100-year) flood events at 57 116 stream reaches across the contiguous United States (CONUS). We link these flood projections to a database of assets within mapped flood hazard zones to model changes in inland flooding damages throughout the CONUS over the remainder of the 21st century. Our model generates early 21st century flood damages that reasonably approximate the range of historical observations and trajectories of future damages that vary substantially depending on the greenhouse gas (GHG) emissions pathway. The difference in modeled flood damages between higher and lower emissions pathways approaches USD 4 billion per year by 2100 (in undiscounted 2014 dollars), suggesting that aggressive GHG emissions reductions could generate significant monetary benefits over the long term in terms of reduced flood damages. Although the downscaled hydrologic data we used have been applied to flood impacts studies elsewhere, this research expands on earlier work to quantify changes in flood risk by linking future flood exposure to assets and damages on a national scale. Our approach relies on a series of simplifications that could ultimately affect damage estimates (e.g., use of statistical downscaling, reliance on a nationwide hydrologic model, and linking damage estimates only to 1 % AEP floods). Although future work is needed to test the sensitivity of our results to these methodological choices, our results indicate that monetary damages from inland flooding could be significantly reduced through substantial GHG mitigation.

Estimating spatially distributed parameters remains one of the biggest challenges for large domain hydrologic modeling. Many large domain modeling efforts rely on spatially inconsistent parameter fields, e.g., patchwork patterns resulting from individual basin calibrations, parameter fields generated through default transfer functions that relate geophysical attributes to model parameters, or spatially constant, default parameter values. This paper provides an initial assessment of a multi-scale parameter regionalization (MPR) method over large geographical domains to derive seamless parameters in a spatially consistent manner. MPR applies transfer functions at the native scale of the geophysical data, and then scales these model parameters to the desired model resolution. We developed a stand-alone framework called MPR-flex for multi-model use and applied MPR-flex to the Variable Infiltration Capacity model to produce hydrologic simulations over the contiguous USA (CONUS). We first independently calibrate 531 basins across the CONUS to obtain a performance benchmark for each basin. To derive the CONUS parameter fields, we perform a joint MPR calibration using all but the poorest behaved basins to obtain a single set of transfer function parameters that are applied to the entire CONUS. Results show that the CONUS-wide calibration has similar performance compared to previous simulations using a patchwork quilt of partially calibrated parameter sets, but without the spatial discontinuities in parameters that characterize some previous CONUS-domain model simulations. Several avenues to improve CONUS-wide calibration remain, including selection of calibration basins, objective function formulation, as well as MPR-flex improvements including transfer function formations and scaling operator optimization.

Freezing temperatures can influence the material properties of trees. Understanding how forest biometeorological interactions respond to freezing temperatures is important for forest science and management, as its effects can cascade through coupled hydrological and ecological processes including limitations on tree growth, changes to canopy interception, and influences on atmospheric turbulence. This short communication details the effect of sub-freezing temperatures on the mechanical displacement of three in situ Pinus contorta Douglas (lodgepole pine) stems during a winter season (2013–2014) at Niwot Ridge Mountain Research Station (CO, USA). Although previous research on harvested trees suggests longitudinal stem expansion should occur as temperatures decrease below freezing, we observed longitudinal compression in live stems below −3 °C that linearl y correlated with sub-zero air temperatures ranging down to −22 °C. Freeze-related compression of stems frequently achieved displacement magnitudes comparable to applying 55 kN of compressive force (800 μm). Hypotheses are proposed to explain the significant mechanical displacement observed due to freezing: (1) internal gas release as freezing fronts propagate through stems; (2) short-term relocation of water via sap exudation or preferential ice crystal growth. As the observed stem compression response to freezing is substantial, future work on underlying processes and consequences is merited. The measurement technique used here may prove useful to others interested in the dynamics of stem freezing.

The coarse grid spacing of global circulation models necessitates the application of climate downscaling to investigate the local impact of a changing global climate. Difficulties arise for data sparse regions in complex topography which are computationally demanding for dynamic downscaling and often not suitable for statistical downscaling due to the lack of high quality observational data. The Intermediate Complexity Atmospheric Research Model (ICAR) is a physics-based model that can be applied without relying on measurements for training and is computationally more efficient than dynamic downscaling models. This study presents the first in-depth evaluation of multi-year precipitation time series generated with ICAR on a 4 x 4 km² grid for the South Island of New Zealand for the eleven-year period from 2007 to 2017. It focuses on complex topography and evaluates ICAR at 16 weather stations, eleven of which are situated in the Southern Alps between 700 m MSL and 2150 m MSL. ICAR is diagnosed with standard skill scores and the effect of model top elevation, topography, season, atmospheric background state and synoptic weather patterns on these scores are investigated. The results show a strong dependence of ICAR skill on the choice of the model top elevation, with the highest scores obtained for 4 km above topography. Furthermore, ICAR is found to provide added value over its ERA-Interim reanalysis forcing data set for alpine weather stations, improving mean squared errors (MSE) up to 53%. It performs similarly during all seasons with an MSE minimum during winter, while flow of higher linearity and atmospheric stability were found to increase skill scores. ICAR scores are highest during weather patterns associated with flow perpendicular to the Southern Alps and lowest for flow parallel to the alpine range. While measured precipitation is underestimated by ICAR, these results show the skill of ICAR in a real-world application, and may be improved upon by further observational tuning or bias correction techniques. Based on these findings ICAR shows the potential to generate downscaled fields for long term impact studies in data sparse regions with complex topography.

We examine the effects of regional climate model (RCM) horizontal resolution and forcing scaling (i.e., spatial aggregation of meteorological datasets) on the portrayal of climate change impacts. Specifically, we assess how the above decisions affect: (i) historical simulation of signature measures of hydrologic behavior, and (ii) projected changes in terms of annual water balance and hydrologic signature measures. To this end, we conduct our study in three catchments located in the headwaters of the Colorado River basin. Meteorological forcings for current and a future climate projection are obtained at three spatial resolutions (4-, 12- and 36-km) from dynamical downscaling with the Weather Research and Forecasting (WRF) regional climate model, and hydrologic changes are computed using four different hydrologic model structures. These projected changes are compared to those obtained from running hydrologic simulations with current and future 4-km WRF climate outputs re-scaled to 12- and 36-km. The results show that the horizontal resolution of WRF simulations heavily affects basin-averaged precipitation amounts, propagating into large differences in simulated signature measures across model structures. The implications of re-scaled forcing datasets on historical performance were primarily observed on simulated runoff seasonality. We also found that the effects of WRF grid resolution on projected changes in mean annual runoff and evapotranspiration may be larger than the effects of hydrologic model choice, which surpasses the effects from re-scaled forcings. Scaling effects on projected variations in hydrologic signature measures were found to be generally smaller than those coming from WRF resolution; however, forcing aggregation in many cases reversed the direction of projected changes in hydrologic behavior.

Orographic precipitation and snowpack provide a vital water resource for the western U.S., while convective precipitation accounts for a significant part of annual precipitation in the eastern U.S. As a result, water managers are keenly interested in their fate under climate change. However, previous studies of water cycle changes in the U.S. have been conducted with climate models of relatively coarse resolution, leading to potential misrepresentation of key physical processes. This paper presents results from a high-resolution climate change simulation that permits convection and resolves mesoscale orography at 4-km grid spacing over much of North America using the Weather Research and Forecasting (WRF) model. Two 13-year simulations were performed, consisting of a retrospective simulation (October 2000–September 2013) with initial and boundary conditions from ERA-interim and a future climate sensitivity simulation with modified reanalysis-derived initial and boundary conditions through adding the CMIP5 ensemble-mean high-end emission scenario climate change. The retrospective simulation is evaluated by validating against Snowpack Telemetry (SNOTEL) and an ensemble of gridded observational datasets. It shows overall good performance capturing the annual/seasonal/sub-seasonal precipitation and surface temperature climatology except for a summer dry and warm bias in the central U.S. In particular, the WRF seasonal precipitation agrees with SNOTEL observations within a few percent over the mountain ranges, providing confidence in the model’s estimation of western U.S. seasonal snowfall and snowpack. The future climate simulation forced with warmer and moister perturbed boundary conditions enhances annual and winter-spring-fall seasonal precipitation over most of the contiguous United States (CONUS), but suppresses summertime precipitation in the central U.S. The WRF-downscaled climate change simulations provide a high-resolution dataset (i.e., High-Resolution CONUS downscaling, HRCONUS) to the community for studying one possible scenario of regional climate changes and impacts.

A North American bark beetle infestation has killed billions of trees, increasing soil nitrogen and raising concern for N loss impacts on downstream ecosystems and water resources. There is surprisingly little evidence of stream N response in large basins, which may result from surviving vegetation uptake, gaseous loss, or dilution by streamflow from unimpacted stands. Observations are lacking along hydrologic flow paths connecting soils with streams, challenging our ability to determine where and how attenuation occurs. Here we quantified biogeochemical concentrations and fluxes at a lodgepole pine-dominated site where bark beetle infestation killed 50–60% of trees. We used nested observations along hydrologic flow paths connecting hillslope soils to streams of up to third order. We found soil water NO3 concentrations increased 100-fold compared to prior research at this and nearby southeast Wyoming sites. Nitrogen was lost below the major rooting zone to hillslope groundwater, where dissolved organic nitrogen (DON) increased by 3–10 times (mean 1.65 mg L−1) and NO3-N increased more than 100-fold (3.68 mg L−1) compared to preinfestation concentrations. Most of this N was removed as hillslope groundwater drained through riparian soils, and NO3 remained low in streams. DON entering the stream decreased 50% within 5 km downstream, to concentrations typical of unimpacted subalpine streams (~0.3 mg L−1). Although beetle outbreak caused hillslope N losses similar to other disturbances, up to 5.5 kg ha−1y−1, riparian and in-stream removal limited headwater catchment export to

With limited computational resources, there is a need for computationally frugal models. This is particularly the case for atmospheric sciences, which have long relied on either simplistic analytical solutions or computationally expensive numerical models. The simpler solutions are inadequate for many problems, while the cost of numerical models makes their use impossible for many problems, most notably high-resolution climate downscaling applications spanning large areas, long time periods, and many global climate projections. Here the Intermediate Complexity Atmospheric Research model (ICAR) is presented to provide a new step along the modeling complexity continuum. ICAR leverages an analytical solution for high-resolution perturbations to wind velocities, in conjunction with numerical physics schemes, that is, advection and cloud microphysics, to simulate the atmosphere. The focus of the initial development of ICAR is for predictions of precipitation, and eventually temperature, humidity, and radiation at the land surface. Comparisons between ICAR and the Weather Research and Forecasting (WRF) Model for simulations over an idealized mountain are presented, as well as among ICAR, WRF, and the Parameter-Elevation Regressions on Independent Slopes Model (PRISM) observation-based product for a year-long simulation over the Colorado Rockies. In the ideal simulations, ICAR matches WRF precipitation predictions across a range of environmental conditions with a coefficient of determination r2 of 0.92. In the Colorado Rockies, ICAR, WRF, and PRISM show very good agreement, with differences between ICAR and WRF comparable to the differences between WRF and PRISM in the cool season. For these simulations, WRF required 140–800 times more computational resources than ICAR.

Continental-domain assessments of climate change impacts on water resources typically rely on statistically downscaled climate model outputs to force hydrologic models at a finer spatial resolution. This study examines the effects of four statistical downscaling methods [bias-corrected constructed analog (BCCA), bias-corrected spatial disaggregation applied at daily (BCSDd) and monthly scales (BCSDm), and asynchronous regression (AR)] on retrospective hydrologic simulations using three hydrologic models with their default parameters (the Community Land Model, version 4.0; the Variable Infiltration Capacity model, version 4.1.2; and the Precipitation–Runoff Modeling System, version 3.0.4) over the contiguous United States (CONUS). Biases of hydrologic simulations forced by statistically downscaled climate data relative to the simulation with observation-based gridded data are presented. Each statistical downscaling method produces different meteorological portrayals including precipitation amount, wet-day frequency, and the energy input (i.e., shortwave radiation), and their interplay affects estimations of precipitation partitioning between evapotranspiration and runoff, extreme runoff, and hydrologic states (i.e., snow and soil moisture). The analyses show that BCCA underestimates annual precipitation by as much as −250 mm, leading to unreasonable hydrologic portrayals over the CONUS for all models. Although the other three statistical downscaling methods produce a comparable precipitation bias ranging from −10 to 8 mm across the CONUS, BCSDd severely overestimates the wet-day fraction by up to 0.25, leading to different precipitation partitioning compared to the simulations with other downscaled data. Overall, the choice of downscaling method contributes to less spread in runoff estimates (by a factor of 1.5–3) than the choice of hydrologic model with use of the default parameters if BCCA is excluded.

This paper reviews the benefits of a better understanding of atmospheric exchange processes over mountains. These processes affect weather and climate variables that are important in meteorological applications related to many scientific disciplines and sectors of the economy. We focus this review on examples of meteorological applications in hydrology, ecology, agriculture, urban planning, wind energy, transportation, air pollution, and climate change. These examples demonstrate the benefits of a more accurate knowledge of atmospheric exchange processes over mountains, including a better understanding of snow redistribution, microclimate, land-cover change, frost hazards, urban ventilation, wind gusts, road temperatures, air pollution, and the impacts of climate change. The examples show that continued research on atmospheric exchange processes over mountains is warranted, and that a recognition of the potential benefits can inspire new research directions. An awareness of the links between basic research topics and applications is important to the success and impact of new efforts that aim at better understanding atmospheric exchange processes over mountains. To maximize the benefits of future research for meteorological applications, coordinated international efforts involving scientists studying atmospheric exchange processes, as well as scientists and stakeholders representing many other scientific disciplines and economic sectors are required.

The rationale for using multi-model ensembles in climate change projections and impacts research is often based on the expectation that different models constitute independent estimates, so that a range of models allows a better characterisation of the uncertainties in the representation of the climate system than a single model. However, it is known that research groups share literature, ideas for representations of processes, parameterisations, evaluation data sets and even sections of model code. Thus, nominally different models might have similar biases because of similarities in the way they represent a subset of processes, or even be near duplicates of others, weakening the assumption that they constitute independent estimates. If there are near-replicates of some models, then treating all models equally is likely to bias the inferences made using these ensembles. The challenge is to establish the degree to which this might be true for any given application. While this issue is recognized by many in the community, quantifying and accounting for model dependence in anything other than an ad-hoc way is challenging. Here we present a synthesis of the range of disparate attempts to define, quantify and address model dependence in multi-model climate ensembles in a common conceptual framework, and provide guidance on how users can test the efficacy of approaches that move beyond the equally weighted ensemble. In the upcoming Coupled Model Intercomparison Project phase 6 (CMIP6), several new models that are closely related to existing models are anticipated, as well as large ensembles from some models. We argue that quantitatively accounting for dependence in addition to model performance, and thoroughly testing the effectiveness of the approach used will be key to a sound interpretation of the CMIP ensembles in future scientific studies.

Metrics based on streamflow and/or climate variables are used in water management for monitoring and evaluating available resources. To reflect future change in the hydrological regime, metrics are estimated using climate change information from Global Climate Models or from stochastic time series representing future climates. Whilst often simple to calculate, many metrics implicitly represent complex physical process. We evaluate the scientific validity of metrics used in a climate change context, demonstrating their use to reflect aspects of timing, magnitude, extreme values, variability, duration, state, system services, and performance. We raise awareness about the following generic issues (a) formulation: metrics often assume stationarity of the input data, which is invalid under climate change; and do not always consider potential changes to seasonality and the relevance of the temporal window used for analysis; (b) climate change input data: how well are the physical processes relevant to the metric represented in the climate change input data; what is the impact of bias correction on relevant spatial and temporal scale dependencies and relevant intervariable dependencies; how realistic are the data in representing sequencing of events and natural variability in large‐scale ocean–atmosphere systems; and (c) decision‐making context: are rules and values that frame the decision‐making process likely to remain constant or change in a future world.\ud \ud If critical climate or hydrological processes are not well represented by the metric constituents, these indices can be misleading about plausible future change. However, knowledge of how to construct a robust metric can safeguard against misleading interpretations about future change.

End users face a range of subjective decisions when evaluating climate change impacts on hydrology, but the importance of these decisions is rarely assessed. In this paper, we evaluate the implications of hydrologic modelling choices on projected changes in the annual water balance, monthly simulated processes, and signature measures (i.e. metrics that quantify characteristics of the hydrologic catchment response) under a future climate scenario. To this end, we compare hydrologic changes computed with four different model structures – whose parameters have been obtained using a common calibration strategy – with hydrologic changes computed with a single model structure and parameter sets from multiple options for different calibration decisions (objective function, local optima, and calibration forcing dataset). Results show that both model structure selection and the parameter estimation strategy affect the direction and magnitude of projected changes in the annual water balance, and that the relative effects of these decisions are basin dependent. The analysis of monthly changes illustrates that parameter estimation strategies can provide similar or larger uncertainties in simulations of some hydrologic processes when compared with uncertainties coming from model choice. We found that the relative effects of modelling decisions on projected changes in catchment behaviour depend on the signature measure analysed. Furthermore, parameter sets with similar performance, but located in different regions of the parameter space, provide very different projections for future catchment behaviour. More generally, the results obtained in this study prompt the need to incorporate parametric uncertainty in multi-model frameworks to avoid an over-confident portrayal of climate change impacts. Copyright © 2015 John Wiley & Sons, Ltd.

Gridded spatiotemporal maps of precipitation are essential for hydrometeorological and ecological analyses. In the United States, most of these datasets are developed using the Cooperative Observer (COOP) network of ground-based precipitation measurements, interpolation, and the Parameter–Elevation Regressions on Independent Slopes Model (PRISM) to map these measurements to places where data are not available. Here, we evaluate two daily datasets gridded at ° resolution against independent daily observations from over 100 snow pillows in California’s Sierra Nevada from 1990 to 2010. Over the entire period, the gridded datasets performed reasonably well, with median total water-year errors generally falling within ±10%. However, errors in individual storm events sometimes exceeded 50% for the median difference across all stations, and in many cases, the same underpredicted storms appear in both datasets. Synoptic analysis reveals that these underpredicted storms coincide with 700-hPa winds from the west or northwest, which are associated with post-cold-frontal flow and disproportionately small precipitation rates in low-elevation valley locations, where the COOP stations are primarily located. This atmospheric circulation leads to a stronger than normal valley-to-mountain precipitation gradient and underestimation of actual mountain precipitation. Because of the small average number of storms (

The assessment of climate change impacts on water resources involves several methodological decisions, including choices of global climate models (GCMs), emission scenarios, downscaling techniques, and hydrologic modeling approaches. Among these, hydrologic model structure selection and parameter calibration are particularly relevant and usually have a strong subjective component. The goal of this research is to improve understanding of the role of these decisions on the assessment of the effects of climate change on hydrologic processes. The study is conducted in three basins located in the Colorado headwaters region, using four different hydrologic model structures [PRMS, VIC, Noah LSM, and Noah LSM with multiparameterization options (Noah-MP)]. To better understand the role of parameter estimation, model performance and projected hydrologic changes (i.e., changes in the hydrology obtained from hydrologic models due to climate change) are compared before and after calibration with the University of Arizona shuffled complex evolution (SCE-UA) algorithm. Hydrologic changes are examined via a climate change scenario where the Community Climate System Model (CCSM) change signal is used to perturb the boundary conditions of the Weather Research and Forecasting (WRF) Model configured at 4-km resolution. Substantial intermodel differences (i.e., discrepancies between hydrologic models) in the portrayal of climate change impacts on water resources are demonstrated. Specifically, intermodel differences are larger than the mean signal from the CCSM–WRF climate scenario examined, even after the calibration process. Importantly, traditional single-objective calibration techniques aimed to reduce errors in runoff simulations do not necessarily improve intermodel agreement (i.e., same outputs from different hydrologic models) in projected changes of some hydrological processes such as evapotranspiration or snowpack.

This work advances a unified approach to process-based hydrologic modeling, which we term the “Structure for Unifying Multiple Modeling Alternatives (SUMMA).” The modeling framework, introduced in the companion paper, uses a general set of conservation equations with flexibility in the choice of process parameterizations (closure relationships) and spatial architecture. This second paper specifies the model equations and their spatial approximations, describes the hydrologic and biophysical process parameterizations currently supported within the framework, and illustrates how the framework can be used in conjunction with multivariate observations to identify model improvements and future research and data needs. The case studies illustrate the use of SUMMA to select among competing modeling approaches based on both observed data and theoretical considerations. Specific examples of preferable modeling approaches include the use of physiological methods to estimate stomatal resistance, careful specification of the shape of the within-canopy and below-canopy wind profile, explicitly accounting for dust concentrations within the snowpack, and explicitly representing distributed lateral flow processes. Results also demonstrate that changes in parameter values can make as much or more difference to the model predictions than changes in the process representation. This emphasizes that improvements in model fidelity require a sagacious choice of both process parameterizations and model parameters. In conclusion, we envisage that SUMMA can facilitate ongoing model development efforts, the diagnosis and correction of model structural errors, and improved characterization of model uncertainty.

This work advances a unified approach to process-based hydrologic modeling to enable controlled and systematic evaluation of multiple model representations (hypotheses) of hydrologic processes and scaling behavior. Our approach, which we term the Structure for Unifying Multiple Modeling Alternatives (SUMMA), formulates a general set of conservation equations, providing the flexibility to experiment with different spatial representations, different flux parameterizations, different model parameter values, and different time stepping schemes. In this paper, we introduce the general approach used in SUMMA, detailing the spatial organization and model simplifications, and how different representations of multiple physical processes can be combined within a single modeling framework. We discuss how SUMMA can be used to systematically pursue the method of multiple working hypotheses in hydrology. In particular, we discuss how SUMMA can help tackle major hydrologic modeling challenges, including defining the appropriate complexity of a model, selecting among competing flux parameterizations, representing spatial variability across a hierarchy of scales, identifying potential improvements in computational efficiency and numerical accuracy as part of the numerical solver, and improving understanding of the various sources of model uncertainty.

We examine the scalability and performance of an open-source, coarray Fortran (CAF) mini-application (mini-app) that implements the parallel, numerical algorithms that dominate the execution of The Intermediate Complexity Atmospheric Research (ICAR) [4] model developed at the the National Center for Atmospheric Research (NCAR). The Fortran 2008 mini-app includes one Fortran 2008 implementation of a collective subroutine defined in the Committee Draft of the upcoming Fortran 2018 standard. The ability of CAF to run atop various communication layers and the increasing CAF compiler availability facilitated evaluating several compilers, runtime libraries and hardware platforms. Results are presented for the GNU and Cray compilers, each of which offers different parallel runtime libraries employing one or more communication layers, including MPI, OpenSHMEM, and proprietary alternatives. We study performance on multi- and many-core processors in distributed memory. The results show promising scaling across a range of hardware, compiler, and runtime choices on up to ~100,000 cores.