A document by Chang Liao. Click on the document to view its contents.

Flow direction modeling consists of (1) an accurate representation of the river network and (2) digital elevation model (DEM) processing to preserve characteristics with hydrological significance. In part 1 of our study, we presented a mesh-independent approach to representing river networks on different types of meshes. This follow-up part 2 study presents a novel DEM processing approach for flow direction modeling. This approach consists of (1) a topological relationship-based hybrid breaching-filling method to conduct stream burning for the river network and (2) a modified depression removal method for rivers and hillslopes. Our methods minimize modifications to surface elevations and provide a robust two-step procedure to remove local depressions in DEM. They are mesh-independent and can be applied to both structured and unstructured meshes. We applied our new methods to the Susquehanna River Basin with different model configurations. The results show that topological relationship-based stream burning and depression-filling methods can reproduce the correct river networks, providing high-quality flow direction and other characteristics for hydrologic and Earth system models.

This work documents version two of the Department of Energy’s Energy Exascale Earth System Model (E3SM). E3SM version 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid-latitudes and 30 km at the equator and poles. The model performance is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima (DECK) simulations augmented with historical simulations as well as simulations to evaluate impact of different forcing agents. The simulated climate is generally realistic, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Programme (WCRP) assessment. However, E3SMv2 significantly underestimates the global mean temperature in the second half of the historical record. An analysis of single-forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol-related forcing.

The Delaware River is a major freshwater supplier of New York City (NYC). Nearly half of NYC drinking water is supplied by inter-basin transfer of surface water stored in reservoirs within the upper reaches of the Delaware River. In its lower reaches, the Delaware River is a tidal estuary, and upstream freshwater discharge provides a critical control on estuary salinity. During the record 1950–1960’s drought, NYC water withdrawals exacerbated low flows. Estuary salinity reached levels that threatened freshwater intakes and groundwater recharge, resulting in legal action and Supreme Court decrees. We revisit this classic case study in coupled human and natural systems using the Energy Exascale Earth System Model (E3SM). The E3SM water management sub-model is updated to include inter-basin water transfer and reservoir-specific operating rules. Model simulations are developed to investigate competition between NYC water demand and in-stream flow targets needed to maintain estuary salinity within regulatory guidelines under historic and future climate. To our knowledge, this is a first demonstration of an Earth System Model simulation with inter-basin water transfer, which, in this study area, provides water for nearly five million people living outside the Delaware River basin in New York City and New Jersey.

River networks are important features in surface hydrology. However, accurately representing river networks in spatially distributed hydrologic and Earth system models is often sensitive to the model’s spatial resolution. Specifically, river networks are often misrepresented because of the mismatch between the model’s spatial resolution and river network details, resulting in significant uncertainty in the projected flow direction. In this study, we developed a topological relationships-based river network representation method for spatially distributed hydrologic models. This novel method uses (1) graph theory algorithms to simplify real-world vector-based river networks and assist in mesh generation; and (2) a topological relationship-based method to reconstruct conceptual river networks. The main advantages of our method are that (1) it combines the strengths of vector-based and DEM raster-based river network extraction methods; and (2) it is mesh-independent and can be applied to both structured and unstructured meshes. This method paves a path for advanced terrain analysis and hydrologic modeling across different scales.