Search Results
You are looking at 21 - 25 of 25 items for
- Author or Editor: Bart Nijssen x
- Refine by Access: All Content x
Abstract
The intent of this paper is to encourage improved numerical implementation of land models. Our contributions in this paper are twofold. First, we present a unified framework to formulate and implement land model equations. We separate the representation of physical processes from their numerical solution, enabling the use of established robust numerical methods to solve the model equations. Second, we introduce a set of synthetic test cases (the laugh tests) to evaluate the numerical implementation of land models. The test cases include storage and transmission of water in soils, lateral subsurface flow, coupled hydrological and thermodynamic processes in snow, and cryosuction processes in soil. We consider synthetic test cases as “laugh tests” for land models because they provide the most rudimentary test of model capabilities. The laugh tests presented in this paper are all solved with the Structure for Unifying Multiple Modeling Alternatives (SUMMA) model implemented using the Suite of Nonlinear and Differential/Algebraic Equation Solvers (SUNDIALS). The numerical simulations from SUMMA/SUNDIALS are compared against 1) solutions to the synthetic test cases from other models documented in the peer-reviewed literature, 2) analytical solutions, and 3) observations made in laboratory experiments. In all cases, the numerical simulations are similar to the benchmarks, building confidence in the numerical model implementation. We posit that some land models may have difficulty in solving these benchmark problems. Dedicating more effort to solving synthetic test cases is critical in order to build confidence in the numerical implementation of land models.
Abstract
The intent of this paper is to encourage improved numerical implementation of land models. Our contributions in this paper are twofold. First, we present a unified framework to formulate and implement land model equations. We separate the representation of physical processes from their numerical solution, enabling the use of established robust numerical methods to solve the model equations. Second, we introduce a set of synthetic test cases (the laugh tests) to evaluate the numerical implementation of land models. The test cases include storage and transmission of water in soils, lateral subsurface flow, coupled hydrological and thermodynamic processes in snow, and cryosuction processes in soil. We consider synthetic test cases as “laugh tests” for land models because they provide the most rudimentary test of model capabilities. The laugh tests presented in this paper are all solved with the Structure for Unifying Multiple Modeling Alternatives (SUMMA) model implemented using the Suite of Nonlinear and Differential/Algebraic Equation Solvers (SUNDIALS). The numerical simulations from SUMMA/SUNDIALS are compared against 1) solutions to the synthetic test cases from other models documented in the peer-reviewed literature, 2) analytical solutions, and 3) observations made in laboratory experiments. In all cases, the numerical simulations are similar to the benchmarks, building confidence in the numerical model implementation. We posit that some land models may have difficulty in solving these benchmark problems. Dedicating more effort to solving synthetic test cases is critical in order to build confidence in the numerical implementation of land models.
Abstract
The near-surface climate, including the atmosphere, ocean, sea ice, and land state and fluxes, in the initial version of the Regional Arctic System Model (RASM) are presented. The sensitivity of the RASM near-surface climate to changes in atmosphere, ocean, and sea ice parameters and physics is evaluated in four simulations. The near-surface atmospheric circulation is well simulated in all four RASM simulations but biases in surface temperature are caused by biases in downward surface radiative fluxes. Errors in radiative fluxes are due to biases in simulated clouds with different versions of RASM simulating either too much or too little cloud radiative impact over open ocean regions and all versions simulating too little cloud radiative impact over land areas. Cold surface temperature biases in the central Arctic in winter are likely due to too few or too radiatively thin clouds. The precipitation simulated by RASM is sensitive to changes in evaporation that were linked to sea surface temperature biases. Future work will explore changes in model microphysics aimed at minimizing the cloud and radiation biases identified in this work.
Abstract
The near-surface climate, including the atmosphere, ocean, sea ice, and land state and fluxes, in the initial version of the Regional Arctic System Model (RASM) are presented. The sensitivity of the RASM near-surface climate to changes in atmosphere, ocean, and sea ice parameters and physics is evaluated in four simulations. The near-surface atmospheric circulation is well simulated in all four RASM simulations but biases in surface temperature are caused by biases in downward surface radiative fluxes. Errors in radiative fluxes are due to biases in simulated clouds with different versions of RASM simulating either too much or too little cloud radiative impact over open ocean regions and all versions simulating too little cloud radiative impact over land areas. Cold surface temperature biases in the central Arctic in winter are likely due to too few or too radiatively thin clouds. The precipitation simulated by RASM is sensitive to changes in evaporation that were linked to sea surface temperature biases. Future work will explore changes in model microphysics aimed at minimizing the cloud and radiation biases identified in this work.
Abstract
Surface meteorological analyses are an essential input (termed “forcing”) for hydrologic modeling. This study investigated the sensitivity of different hydrologic model configurations to temporal variations of seven forcing variables (precipitation rate, air temperature, longwave radiation, specific humidity, shortwave radiation, wind speed, and air pressure). Specifically, the effects of temporally aggregating hourly forcings to hourly daily average forcings were examined. The analysis was based on 14 hydrological outputs from the Structure for Unifying Multiple Modeling Alternatives (SUMMA) model for the 671 Catchment Attributes and Meteorology for Large-Sample Studies (CAMELS) basins across the contiguous United States (CONUS). Results demonstrated that the hydrologic model sensitivity to temporally aggregating the forcing inputs varies across model output variables and model locations. We used Latin hypercube sampling to sample model parameters from eight combinations of three influential model physics choices (three model decisions with two options for each decision, i.e., eight model configurations). Results showed that the choice of model physics can change the relative influence of forcing on model outputs and the forcing importance may not be dependent on the parameter space. This allows for model output sensitivity to forcing aggregation to be tested prior to parameter calibration. More generally, this work provides a comprehensive analysis of the dependence of modeled outcomes on input forcing behavior, providing insight into the regional variability of forcing variable dominance on modeled outputs across CONUS.
Abstract
Surface meteorological analyses are an essential input (termed “forcing”) for hydrologic modeling. This study investigated the sensitivity of different hydrologic model configurations to temporal variations of seven forcing variables (precipitation rate, air temperature, longwave radiation, specific humidity, shortwave radiation, wind speed, and air pressure). Specifically, the effects of temporally aggregating hourly forcings to hourly daily average forcings were examined. The analysis was based on 14 hydrological outputs from the Structure for Unifying Multiple Modeling Alternatives (SUMMA) model for the 671 Catchment Attributes and Meteorology for Large-Sample Studies (CAMELS) basins across the contiguous United States (CONUS). Results demonstrated that the hydrologic model sensitivity to temporally aggregating the forcing inputs varies across model output variables and model locations. We used Latin hypercube sampling to sample model parameters from eight combinations of three influential model physics choices (three model decisions with two options for each decision, i.e., eight model configurations). Results showed that the choice of model physics can change the relative influence of forcing on model outputs and the forcing importance may not be dependent on the parameter space. This allows for model output sensitivity to forcing aggregation to be tested prior to parameter calibration. More generally, this work provides a comprehensive analysis of the dependence of modeled outcomes on input forcing behavior, providing insight into the regional variability of forcing variable dominance on modeled outputs across CONUS.
Abstract
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.
Abstract
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.
Abstract
The Olympic Mountains Experiment (OLYMPEX) took place during the 2015/16 fall–winter season in the vicinity of the mountainous Olympic Peninsula of Washington State. The goals of OLYMPEX were to provide physical and hydrologic ground validation for the U.S.–Japan Global Precipitation Measurement (GPM) satellite mission and, more specifically, to study how precipitation in Pacific frontal systems is modified by passage over coastal mountains. Four transportable scanning dual-polarization Doppler radars of various wavelengths were installed. Surface stations were placed at various altitudes to measure precipitation rates, particle size distributions, and fall velocities. Autonomous recording cameras monitored and recorded snow accumulation. Four research aircraft supplied by NASA investigated precipitation processes and snow cover, and supplemental rawinsondes and dropsondes were deployed during precipitation events. Numerous Pacific frontal systems were sampled, including several reaching “atmospheric river” status, warm- and cold-frontal systems, and postfrontal convection.
Abstract
The Olympic Mountains Experiment (OLYMPEX) took place during the 2015/16 fall–winter season in the vicinity of the mountainous Olympic Peninsula of Washington State. The goals of OLYMPEX were to provide physical and hydrologic ground validation for the U.S.–Japan Global Precipitation Measurement (GPM) satellite mission and, more specifically, to study how precipitation in Pacific frontal systems is modified by passage over coastal mountains. Four transportable scanning dual-polarization Doppler radars of various wavelengths were installed. Surface stations were placed at various altitudes to measure precipitation rates, particle size distributions, and fall velocities. Autonomous recording cameras monitored and recorded snow accumulation. Four research aircraft supplied by NASA investigated precipitation processes and snow cover, and supplemental rawinsondes and dropsondes were deployed during precipitation events. Numerous Pacific frontal systems were sampled, including several reaching “atmospheric river” status, warm- and cold-frontal systems, and postfrontal convection.
