Ethan T. Coon

Managing complexity in simulations of land surface and near-surface processes

Increasing computing power and the growing role of simulation in Earth systems science have led to an increase in the number and complexity of processes in modern simulators. We present a multiphysics framework that specifies interfaces for coupled processes and automates weak and strong coupling strategies to manage this complexity. Process management is enabled by viewing the system of equations as a tree, where individual equations are associated with leaf nodes and coupling strategies with internal nodes. A dynamically generated dependency graph connects a variable to its dependencies, streamlining and automating model evaluation, easing model development, and ensuring models are modular and flexible. Additionally, the dependency graph is used to ensure that data requirements are consistent between all processes in a given simulation. Here we discuss the design and implementation of these concepts within the Arcos framework, and demonstrate their use for verification testing and hypothesis evaluation in numerical experiments. We describe a conceptual model for managing complexity in multiphysics software.A process tree describes how individual equations are coupled.A dependency graph describes how variables are dependent upon each other.The model is implemented within the Arcos software framework.Examples of code and model runs are shown to demonstrate the idea.

The integrated hydrologic model intercomparison project, IH‐MIP2: A second set of benchmark results to diagnose integrated hydrology and feedbacks

Emphasizing the physical intricacies of integrated hydrology and feedbacks in simulating connected, variably saturated groundwater-surface water systems, the Integrated Hydrologic Model Intercomparison Project initiated a second phase (IH-MIP2), increasing the complexity of the benchmarks of the first phase. The models that took part in the intercomparison were ATS, Cast3M, CATHY, GEOtop, HydroGeoSphere, MIKE-SHE, and ParFlow. IH-MIP2 benchmarks included a tilted v-catchment with 3-D subsurface; a superslab case expanding the slab case of the first phase with an additional horizontal subsurface heterogeneity; and the Borden field rainfall-runoff experiment. The analyses encompassed time series of saturated, unsaturated, and ponded storages, as well as discharge. Vertical cross sections and profiles were also inspected in the superslab and Borden benchmarks. An analysis of agreement was performed including systematic and unsystematic deviations between the different models. Results show generally good agreement between the different models, which lends confidence in the fundamental physical and numerical implementation of the governing equations in the different models. Differences can be attributed to the varying level of detail in the mathematical and numerical representation or in the parameterization of physical processes, in particular with regard to ponded storage and friction slope in the calculation of overland flow. These differences may become important for specific applications such as detailed inundation modeling or when strong inhomogeneities are present in the simulation domain.

Taxila LBM: a parallel, modular lattice Boltzmann framework for simulating pore-scale flow in porous media

The lattice Boltzmann method is a popular tool for pore-scale simulation of flow. This is likely due to the ease of including complex geometries such as porous media and representing multiphase and multifluid flows. Many advancements, including multiple relaxation times, increased isotropy, and others have improved the accuracy and physical fidelity of the method. Additionally, the lattice Bolzmann method is computationally very efficient, thanks to the explicit nature of the algorithm and relatively large amount of local work. The combination of many algorithmic options and efficiency means that a software framework enabling the usage and comparison of these advancements on computers from laptops to large clusters has much to offer. In this paper, we introduce Taxila LBM, an open-source software framework for lattice Boltzmann simulations. We discuss the design of the framework and lay out the features available, including both methods in the literature and a few new enhancements which generalize methods to complex geometries. We discuss the trade-off of accuracy and performance in various methods, noting how the Taxila LBM makes it easy to perform these comparisons for real problems. And finally, we demonstrate a few common applications in pore-scale simulation, including the characterization of permeability of a Berea sandstone and analysis of multifluid flow in heterogenous micromodels.

Thermal effects of groundwater flow through subarctic fens: A case study based on field observations and numerical modeling

Modeling and observation of ground temperature dynamics are the main tools for understanding current permafrost thermal regimes and projecting future thaw. Until recently, most studies on permafrost have focused on vertical ground heat fluxes. Groundwater can transport heat in both lateral and vertical directions but its influence on ground temperatures at local scales in permafrost environments is not well understood. In this study we combine field observations from a subarctic fen in the sporadic permafrost zone with numerical simulations of coupled water and thermal fluxes. At the Tavvavuoma study site in northern Sweden, ground temperature profiles and groundwater levels were observed in boreholes. These observations were used to set up one- and two-dimensional simulations down to 2 m depth across a gradient of permafrost conditions within and surrounding the fen. Two-dimensional scenarios representing the fen under various hydraulic gradients were developed to quantify the influence of groundwater flow on ground temperature. Our observations suggest that lateral groundwater flow significantly affects ground temperatures. This is corroborated by modeling results that show seasonal ground ice melts 1 month earlier when a lateral groundwater flux is present. Further, although the thermal regime may be dominated by vertically conducted heat fluxes during most of the year, isolated high groundwater flow rate events such as the spring freshet are potentially important for ground temperatures. As sporadic permafrost environments often contain substantial portions of unfrozen ground with active groundwater flow paths, knowledge of this heat transport mechanism is important for understanding permafrost dynamics in these environments.

Influences and interactions of inundation, peat, and snow on active layer thickness

The effect of three environmental conditions: 1) thickness of organic soil, 2) snow depth, and 3) soil moisture content or water table height above and below the soil surface, on active layer thickness (ALT) are investigated using an ensemble of 1D thermal hydrology models. Sensitivity analyses of the ensemble exposed the isolated influence of each environmental condition on ALT and their multivariate interactions. The primary and interactive influences are illustrated in the form of color maps of ALT change. Results show that organic layer acts as a strong insulator, and its thickness is the dominant control of ALT, but the strength of the effect of organic layer thickness is dependent on the saturation state. Snow depth, subsurface saturation, and ponded water depth are strongly codependent and positively correlated to ALT.

From documentation to prediction: raising the bar for thermokarst research

Here we report that to date the majority of published research on thermokarst has been directed at documenting its form, occurrence, and rates of occurrence. The fundamental processes driving thermokarst have long been largely understood. However, the detailed physical couplings between, water, air, soil, and the thermal dynamics governing freeze-thaw and soil mechanics is less understood and not captured in models aimed at predicting the response of frozen soils to warming and thaw. As computational resources increase more sophisticated mechanistic models can be applied; these show great promise as predictive tools. These models will be capable of simulating the response of soil deformation to thawing/freezing cycles and the long-term, non-recoverable response of the land surface to the loss of ice. At the same time, advances in remote sensing of permafrost environments also show promise in providing detailed and spatially extensive estimates in the rates and patterns of subsidence. These datasets provide key constraints to calibrate and evaluate the predictive power of mechanistic models. In conclusion, in the coming decade, these emerging technologies will greatly increase our ability to predict when, where, and how thermokarst will occur in a changing climate.

Influences and interactions of inundation, peat, and snow on active layer thickness: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication currently in review [4/2016]. The Advanced Terrestrial Simulator (ATS) was used to simulate thermal hydrological conditions across varied environmental conditions for an ensemble of 1D models of Arctic permafrost. The thickness of organic soil is varied from 2 to 40cm, snow depth is varied from approximately 0 to 1.2 meters, water table depth was varied from -51cm below the soil surface to 31 cm above the soil surface. A total of 15,960 ensemble members are included. Data produced includes the third and fourth simulation year: active layer thickness, time of deepest thaw depth, temperature of the unfrozen soil, and unfrozen liquid saturation, for each ensemble member. Input files used to run the ensemble are also included.

Predicting Climate Feedbacks and Impacts in the Terrestrial Arctic: w14_terraarctic progress report

Regarding the Arctic Terrestrial Simulator (ATS), previous work solved integrated hydrology (coupled surface/subsurface flow) on multiple polygons, and surface flow over larger domains to guide landscape characterization. Solved thermal hydrology with freeze/thaw dynamics in three dimensions. Ongoing efforts apply state of the art thermal hydrology model to complex topography, and include mesh deformation processes.

Calibrated Hydrothermal Parameters, Barrow, Alaska, 2013

A model-observation-experiment process (ModEx) is used to generate three 1D models of characteristic micro-topographical land-formations, which are capable of simulating present active thaw layer (ALT) from current climate conditions. Each column was used in a coupled calibration to identify moss, peat and mineral soil hydrothermal properties to be used in up-scaled simulations. Observational soil temperature data from a tundra site located near Barrow, AK (Area C) is used to calibrate thermal properties of moss, peat, and sandy loam soil to be used in the multiphysics Advanced Terrestrial Simulator (ATS) models. Simulation results are a list of calibrated hydrothermal parameters for moss, peat, and mineral soil hydrothermal parameters.