# Property:Extended model description

From CSDMS

This is a property of type Text.

A

ANUGA is a hydrodynamic model for simulating depth-averaged flows over 2D surfaces. This package adds two new modules (operators) to ANUGA. These are appropriate for reach-scale simulations of flows on mobile-bed streams with spatially extensive floodplain vegetation.
The mathematical framework for the sediment transport operator is described in Simpson and Castelltort (2006) and Davy and Lague (2009). This operator calculates an explicit sediment mass balance within the water column at every cell in order to handle the local disequilibria between entrainment and deposition that arise due to strong spatial variability in shear stress in complex flows.
The vegetation drag operator uses the mathematical approach of Nepf (1999) and Kean and Smith (2006), treating vegetation as arrays of objects (cylinders) that the flow must go around. Compared to methods that simulate the increased roughness of vegetation with a modified Manning's n, this method better accounts for the effects of drag on the body of the flow and the quantifiable differences between vegetation types and densities (as stem diameter and stem spacing). This operator can simulate uniform vegetation as well as spatially-varied vegetation across the domain. The vegetation drag module also accounts for the effects of vegetation on turbulent and mechanical diffusivity, following the equations in Nepf (1997, 1999). +

ANUGA is a hydrodynamic modelling tool that allows users to model realistic flow problems in complex 2D geometries. Examples include dam breaks or the effects of natural hazards such as riverine flooding, storm surges and tsunami. The user must specify a study area represented by a mesh of triangular cells, the topography and bathymetry, frictional resistance, initial values for water level (called stage within ANUGA), boundary conditions and forces such as rainfall, stream flows, windstress or pressure gradients if applicable.
ANUGA tracks the evolution of water depth and horizontal momentum within each cell over time by solving the shallow water wave governing equation using a finite-volume method.
ANUGA also incorporates a mesh generator that allows the user to set up the geometry of the problem interactively as well as tools for interpolation and surface fitting, and a number of auxiliary tools for visualising and interrogating the model output.
Most ANUGA components are written in the object-oriented programming language Python and most users will interact with ANUGA by writing small Python scripts based on the ANUGA library functions. Computationally intensive components are written for efficiency in C routines working directly with Python numpy structures. +

Acronym1D is an add on to Acronym1R in that it adds a flow duration curve to Acronym1R, which computes the volume bedload transport rate per unit width and bedload grain size distribution from a specified surface grain size distribution (with sand removed). +

Acronym1R computes the volume bedload transport rate per unit width and bedload grain size distribution from a specified surface grain size distribution (with sand removed). +

AeoLiS is a process-based model for simulating aeolian sediment transport in situations where supply-limiting factors are important, like in coastal environments. Supply-limitations currently supported are soil moisture contents, sediment sorting and armouring, bed slope effects, air humidity and roughness elements. +

F

Allow for quick estimation of water depths within a flooded domain using only the flood extent layer (polygon) and a DEM of the area. Useful for near-real-time flood analysis, especially from remote sensing mapping.
Version 2.0 offers improved capabilities in coastal areas. +

A

Alpine3D is a model for high resolution simulation of alpine surface processes, in particular snow processes. The model can be forced by measurements from automatic weather stations or by meteorological model outputs (this is handled by the MeteoIO pre-processing library). The core three-dimensional Alpine3D modules consist of a radiation balance model (which uses a view factor approach and includes shortwave scattering and longwave emission from terrain and tall vegetation) and a drifting snow model solving a diffusion equation for suspended snow and a saltation transport equation. The processes in the atmosphere are thus treated in three dimensions and coupled to a distributed one dimensional model of vegetation, snow and soil model (Snowpack) using the assumption that lateral exchange is small in these media. The model can be used to force a distributed catchment hydrology model (AlpineFlow). The model modules can be run in a parallel mode, using either OpenMP and/or MPI. Finally, the Inishell tool provides a GUI for configuring and running Alpine3D.
Alpine3D is a valuable tool to investigate surface dynamics in mountains and is currently used to investigate snow cover dynamics for avalanche warning and permafrost development and vegetation changes under climate change scenarios. It could also be used to create accurate soil moisture assessments for meteorological and flood forecasting. +

W

An extension of the WBMplus (WBM/WTM) model. Introduce a riverine sediment flux component based on the BQART and Psi models. +

An open-source Python package for flexible and customizable simulations of the water cycle that treats the physical components of the water cycle as nodes connected by arcs that convey water and pollutant fluxes between them. +

G

Another derivative of the original SEDSIM, completely rewritten from scratch. It uses finite differences (in addition to the original particle-cell method) to speed up steady flow calculations. It also incorporates compaction algorithms. A general description has been published. +

A

AquaTellUs models fluvial-dominated delta sedimentation. AquaTellUS uses a nested model approach; a 2D longitudinal profiles, embedded as a dynamical flowpath in a 3D grid-based space. A main channel belt is modeled as a 2D longitudinal profile that responds dynamically to changes in discharge, sediment load and sea level. Sediment flux is described by separate erosion and sedimentation components. Multiple grain-size classes are independently tracked. Erosion flux depends on discharge and slope, similar to process descriptions used in hill-slope models and is independent of grain-size. Offshore, where we assume unconfined flow, the erosion capacity decreases with increasing water depth. The erosion flux is a proxy for gravity flows in submarine channels close to the coast and for down-slope diffusion over the entire slope due to waves, tides and creep. Erosion is restricted to the main flowpath. This appears to be valid for the river-channel belt, but underestimates the spatial extent and variability of marine erosion processes.
Deposition flux depends on the stream velocity and on a travel-distance factor, which depends on grain size (i.e. settling velocity). The travel-distance factor is different in the fluvial and marine domains, which results in a sharp increase of the settling rate at the river mouth, mimicking bedload dumping.
Dynamic boundary conditions such as climatic changes over time are incorporated by increasing or decreasing discharge and sediment load for each time step. +

B

BATTRI does the mesh editing, bathymetry incorporation and interpolation, provides the grid generation and refinement properties, prepares the input file to Triangle and visualizes and saves the created grid. +

BIT Model aims to simulate the dynamics of the principal processes that govern the formation and evolution of a barrier island. The model includes sea-level oscillations and sediment distribution operated by waves and currents. Each process determines the deposition of a distinct sediment facies, separately schematized in the spatial domain. Therefore, at any temporal step, it is possible to recognize six different stratigraphic units: bedrock, transitional, overwash, shoreface aeolian and lagoonal. +

BRaKE is a 1-D bedrock channel profile evolution model. It calculates bedrock erosion in addition to treating the delivery, transport, degradation, and erosion-inhibiting effects of large, hillslope-derived blocks of rock. It uses a shear-stress bedrock erosion formulation with additional complexity related to flow resistance, block transport and erosion, and delivery of blocks from the hillslopes. +

Barrier3D is an exploratory model that resolves cross-shore and alongshore topographic variations to simulate the morphological evolution of a barrier segment over time scales of years to centuries. Barrier3D tackles the scale separation between event-based and long-term models by explicitly yet efficiently simulating dune evolution, storm overwash, and a dynamically evolving shoreface in response to individual storm events and sea-level rise. Ecological-geomorphological couplings of the barrier interior can be simulated with a shrub expansion and mortality module. +

BarrierBMFT is a coupled model framework for exploring morphodynamic interactions across components of the entire coastal barrier system, from the ocean shoreface to the mainland forest. The model framework couples Barrier3D (Reeves et al., 2021), a spatially explicit model of barrier evolution, with the Python version of the Coastal Landscape Transect model (CoLT; Valentine et al., 2023), known as PyBMFT-C (Bay-Marsh-Forest Transect Model with Carbon). In the BarrierBMFT coupled model framework, two PyBMFT-C simulations drive evolution of back-barrier marsh, bay, mainland marsh, and forest ecosystems, and a Barrier3D simulation drives evolution of barrier and back-barrier marsh ecosystems. As these model components simultaneously advance, they dynamically evolve together by sharing information annually to capture the effects of key cross-landscape couplings. BarrierBMFT contains no new governing equations or parameterizations itself, but rather is a framework for trading information between Barrier3D and PyBMFT-C.
The use of this coupled model framework requires Barrier3D v2.0 (https://doi.org/10.5281/zenodo.7604068)
and PyBMFT-C v1.0 (https://doi.org/10.5281/zenodo.7853803). +

R

Based on the publication:
Brown, RA, Pasternack, GB, Wallender, WW. 2013. Synthetic River Valleys: Creating Prescribed Topography for Form-Process Inquiry and River Rehabilitation Design. Geomorphology 214: 40–55. http://dx.doi.org/10.1016/j.geomorph.2014.02.025 +

B

Basin and Landscape Dynamics (Badlands) is a parallel TIN-based landscape evolution model, built to simulate topography development at various space and time scales. The model is presently capable of simulating hillslope processes (linear diffusion), fluvial incision ('modified' SPL: erosion/transport/deposition), spatially and temporally varying geodynamic (horizontal + vertical displacements) and climatic forces which can be used to simulate changes in base level, as well as effects of climate changes or sea-level fluctuations. +

Bifurcation is a morphodynamic model of a river delta bifurcation. Model outputs include flux partitioning and 1D bed elevation profiles, all of which can evolve through time. Interaction between the two branches occurs in the reach just upstream of the bifurcation, due to the development of a transverse bed slope. Aside from this interaction, the individual branches are modeled in 1D.
The model generates ongoing avulsion dynamics automatically, arising from the interaction between an upstream positive feedback and the negative feedback from branch progradation and/or aggradation. Depending on the choice of parameters, the model generates symmetry, soft avulsion, or full avulsion. Additionally, the model can include differential subsidence. It can also be run under bypass conditions, simulating the effect of an offshore sink, in which case ongoing avulsion dynamics do not occur.
Possible uses of the model include the study of avulsion, bifurcation stability, and the morphodynamic response of bifurcations to external changes. +

Blocklab treats landscape evolution in landscapes where surface rock may be released as large blocks of rock. The motion, degradation, and effects of large blocks do not play nicely with standard continuum sediment transport theory. BlockLab is intended to incorporate the effects of these large grains in a realistic way. +

C

CAESAR is a cellular landscape evolution model, with an emphasis on fluvial processes, including flow routing, multi grainsize sediment transport. It models morphological change in river catchments. +

CASCADE combines elements of two exploratory morphodynamic models of barrier evolution -- barrier3d (Reeves et al., 2021) and the BarrierR Inlet Environment (brie) model (Nienhuis & Lorenzo-Trueba, 2019) -- into a single model framework. Barrier3d, a spatially-explicit cellular exploratory model, is the core of CASCADE. It is used within the CASCADE framework to simulate the effects of individual storm events and SLR on shoreface evolution; dune dynamics, including dune growth, erosion, and migration; and overwash deposition by individual storms. BRIE is used to simulate large-scale coastline evolution arising from alongshore sediment transport processes; this is accomplished by connecting individual Barrier3d models through diffusive alongshore sediment transport. Human dynamics are incorporated in cascade in two separate modules. The first module simulates strategies for preventing roadway pavement damage during overwashing events, including rebuilding roadways at sufficiently low elevations to allow for burial by overwash, constructing large dunes, and relocating the road into the barrier interior. The second module incorporates management strategies for maintaining a coastal community, including beach nourishment, dune construction, and overwash removal. +

CHILD computes the time evolution of a topographic surface z(x,y,t) by fluvial and hillslope erosion and sediment transport. +

CICE is a computationally efficient model for simulating the growth, melting, and movement of polar sea ice. Designed as one component of coupled atmosphere-ocean-land-ice global climate models, today’s CICE model is the outcome of more than two decades of community collaboration in building a sea ice model suitable for multiple uses including process studies, operational forecasting, and climate simulation. +

CLUMondo is based on the land systems approach. Land systems are socio-ecological systems that reflect land use in a spatial unit in terms of land cover composition, spatial configuration, and the management activities employed. The precise definition of land systems depends on the scale of analysis, the purpose of modelling, and the case study region. In contrast to land cover classifications the role of land use intensity and livestock systems are explicitly addressed. Each land system can be characterized in terms of the fractional land covers.<br>Land systems are characterized based on the amount of forest in the landscape mosaic and the management type ranging from swidden cultivation to permanent cultivation and plantations. +

Caesar Lisflood is a geomorphological / Landscape evolution model that combines the Lisflood-FP 2d hydrodynamic flow model (Bates et al, 2010) with the CAESAR geomorphic model to simulate erosion and deposition in river catchments and reaches over time scales from hours to 1000's of years.
Featuring:
Landscape evolution model simulating erosion and deposition across river reaches and catchments
A hydrodynamic 2D flow model (based on the Lisflood FP code) that conserves mass and partial momentum. (model can be run as flow model alone)
designed to operate on multiple core processors (parallel processing of core functions)
Operates over a wide range to spatial and time scales (1km2 to 1000km2, <1year to 1000+ years)
Easy to use GUI +

P

Calculate the hypsometric integral for each pixel at the catchment. Each pixel is considered a local outlet and the hypsometric integral is calculated according to the characteristics of its contributing area. +

O

Calculate wave-generated bottom orbital velocities from measured surface wave parameters. Also permits calculation of surface wave spectra from wind conditions, from which bottom orbital velocities can be determined. +

S

Calculates non-equilibrium suspended load transport rates of various size-density fractions in the bed +

B

Calculates the bedload transport rates and weights per unit area for each size-density. NB. Bedload transport of different size-densities is proportioned according to the volumes in the bed. +

S

Calculates the constant terminal settling velocity of each size-density fraction's median size from Dietrich's equation. +

E

Calculates the critical Shields Theta for the median size of a distribution and then calculates the critical shear stress of the ith, jth fraction using a hiding function +

Calculates the critical shear stress for entrainment of the median size of each size-density fraction of a bed using Yalin and Karahan formulation, assuming no hiding +

T

Calculates the gaussian or log-gaussian distribution of instantaneous shear stresses on the bed, given a mean and coefficient of variation. +

L

Y

Calculates the total sediment transport rate in an open channel assuming a median bed grain size +

S

Calculation of Density Stratification Effects Associated with Suspended Sediment in Open Channels.
This program calculates the effect of sediment self-stratification on the streamwise velocity and suspended sediment concentration profiles in open-channel flow.
Two options are given. Either the near-bed reference concentration Cr can be specified by the user, or the user can specify a shear velocity due to skin friction u*s and compute Cr from the Garcia-Parker sediment entrainment relation. +

Calculation of Sediment Deposition in a Fan-Shaped Basin, undergoing Piston-Style Subsidence +

D

Calculator for 1D Subaerial Fluvial Fan-Delta with Channel of Constant Width. This model assumes a narrowly channelized 1D fan-delta prograding into standing water. The model uses a single grain size D, a generic total bed material load relation and a constant bed resistance coefficient. The channel is assumed to have a constant width. Water and sediment discharge are specified per unit width. The fan builds outward by forming a prograding delta front with an assigned foreset slope. The code employs a full backwater calculation. +

Calculator for 1D Subaerial Fluvial Fan-Delta with Channel of Constant Width. This model assumes a narrowly channelized 1D fan-delta prograding into standing water. The model uses a single grain size D, a generic total bed material load relation and a constant bed resistance coefficient. The channel is assumed to have a constant width. Water and sediment discharge are specified per unit width. The fan builds outward by forming a prograding delta front with an assigned foreset slope. The code employs the normal flow approximation rather than a full backwater calculation. +

C

CarboCAT uses a cellular automata to model horizontal and vertical distributions of carbonate lithofacies +

ChesROMS is a community ocean modeling system for the Chesapeake Bay region being developed by scientists in NOAA, University of Maryland, CRC (Chesapeake Research Consortium) and MD DNR (Maryland Department of Natural Resources) supported by the NOAA MERHAB program. The model is built based on the Rutgers Regional Ocean Modeling System (ROMS, http://www.myroms.org/) with significant adaptations for the Chesapeake Bay.
The model is developed to provide a community modeling system for nowcast and forecast of 3D hydrodynamic circulation, temperature and salinity, sediment transport, biogeochemical and ecosystem states with applications to ecosystem and human health in the bay. Model validation is based on bay wide satellite remote sensing, real-time in situ measurements and historical data provided by Chesapeake Bay Program.
http://ches.communitymodeling.org/models/ChesROMS/index.php +

Cliffs features:
Shallow-Water approximation;
Use of Cartesian or spherical (lon/lat) coordinates;
1D and 2D configurations;
Structured co-located grid with (optionally) varying spacing;
Run-up on land;
Initial conditions or boundary forcing;
Grid nesting with one-way coupling;
Parallelized with OpenMP;
NetCDF format of input/output data.
Cliffs utilizes VTCS-2 finite-difference scheme and dimensional splitting as in (Titov and Synolakis, 1998), and reflection and inundation computations as in (Tolkova, 2014).
References:
Titov, V.V., and C.E. Synolakis. Numerical modeling of tidal wave runup. J. Waterw. Port Coast. Ocean Eng., 124(4), 157–171 (1998)
Tolkova E. Land-Water Boundary Treatment for a Tsunami Model With Dimensional Splitting.
Pure and Applied Geophysics, 171(9), 2289-2314 (2014) +

B

Coastal barrier model that simulates storm overwash and tidal inlets and estimates coastal barrier transgression resulting from sea-level rise. +

D

Code for estimating long-term exhumation histories and spatial patterns of short-term erosion from the detrital thermochronometric data. +

M

Code functionality and purpose may be found in the following references:
References
# Zhang L., Parker, G., Stark, C.P., Inoue, T., Viparelli, V., Fu, X.D., and Izumi, N. 2015, "Macro-roughness model of bedrock–alluvial river morphodynamics", Earth Surface Dynamics, 3, 113–138.
# Zhang, L., Stark, C.P., Schumer, R., Kwang, J., Li, T.J., Fu, X.D., Wang, G.Q., and Parker, G. 2017, "The advective-diffusive morphodynamics of mixed bedrock-alluvial rivers subjected to spatiotemporally varying sediment supply" (submitted to JGR) +

G

Computes transient (semi-implicit numerical) and steady-state (analytical and numerical) solutions for the long-profile evolution of transport-limited gravel-bed rivers. Such rivers are assumed to have an equilibrium width (following Parker, 1978), experience flow resistance that is proportional to grain size, evolve primarily in response to a single dominant "channel-forming" or "geomorphically-effective" discharge (see Blom et al., 2017, for a recent study and justification of this assumption and how it can be applied), and transport gravel following the Meyer-Peter and Müller (1948) equation. This combination of variables results in a stream-power-like relationship for bed-material sediment discharge, which is then inserted into a valley-resolving Exner equation to compute long-profile evolution. +

C

CruAKtemp is a python 2.7 package that is a data component which serves to provide onthly temperature data over the 20th century for permafrost modeling. The original dataset at higher resolution can be found here:
http://ckan.snap.uaf.edu/dataset/historical-monthly-and-derived-temperature-products-771m-cru-ts
The geographical extent of this CRUAKtemp dataset has been reduced to greatly reduce the number of ocean or Canadian pixels. Also, the spatial resolution has been reduced by a factor of 13 in each direction, resulting in an effective pixel resolution of about 10km.
The data are monthly average temperatures for each month from January 1901 through December 2009. +

D

DFMFON stands for Delft3D-Flexible Mesh (DFM), and MesoFON (MFON) is an open-source software written in Python to simulate the Mangrove and Hydromorphology development mechanistically. We achieve that by coupling the multi-paradigm of the individual-based mangrove model MFON and process-based hydromorphodynamic model DFM. +

DHSVM is a distributed hydrology model that was developed at the University of Washington more than ten years ago. It has been applied both operationally, for streamflow prediction, and in a research capacity, to examine the effects of forest management on peak streamflow, among other things. +

DR3M is a watershed model for routing storm runoff through a Branched system of pipes and (or) natural channels using rainfall as input. DR3M provides detailed simulation of storm-runoff periods selected by the user. There is daily soil-moisture accounting between storms. A drainage basin is represented as a set of overland-flow, channel, and reservoir segments, which jointly describe the drainage features of the basin. This model is usually used to simulate small urban basins. Interflow and base flow are not simulated. Snow accumulation and snowmelt are not simulated. +

DROG3D tracks passive drogues with given harmonic velocity field(s) in a 3-D finite element mesh +

Dakota is a software toolkit, developed at Sandia National Laboratories, that provides an interface between models and a library of analysis methods, including support for sensitivity analysis, uncertainty quantification, optimization, and calibration techniques. Dakotathon is a Python package that wraps and extends Dakota’s file-based user interface. It simplifies the process of configuring and running a Dakota experiment, and it allows a Dakota experiment to be scripted. Any model written in Python that exposes a Basic Model Interface (BMI), as well as any model componentized in the CSDMS modeling framework, automatically works with Dakotathon. Currently, six Dakota analysis methods have been implemented from the much larger Dakota library:
* vector parameter study,
* centered parameter study,
* multidim parameter study,
* sampling,
* polynomial chaos, and
* stochastic collocation. +

C

Data component processed from the CRU-NCEP Climate Model Intercomparison Project - 5, also called CMIP 5. Data presented include the mean annual temperature for each gridcell, mean July temperature and mean January temperature over the period 1902 -2100. This dataset presents the mean of the CMIP5 models, and the original climate models were run for the representative concentration pathway RCP 8.5. +

D

DeltaRCM is a parcel-based cellular flux routing and sediment transport model for the formation of river deltas, which belongs to the broad category of rule-based exploratory models. It has the ability to resolve emergent channel behaviors including channel bifurcation, avulsion and migration. Sediment transport distinguishes two types of sediment: sand and mud, which have different transport and deposition/erosion rules. Stratigraphy is recorded as the sand fraction in layers.
Best usage of DeltaRCM is the investigation of autogenic processes in response to external forcings. +

Demeter is an open source Python package that was built to disaggregate projections of future land allocations generated by an integrated assessment model (IAM). Projected land allocation from IAMs is traditionally transferred to Earth System Models (ESMs) in a variety of gridded formats and spatial resolutions as inputs for simulating biophysical and biogeochemical fluxes. Existing tools for performing this translation generally require a number of manual steps which introduces error and is inefficient. Demeter makes this process seamless and repeatable by providing gridded land use and land cover change (LULCC) products derived directly from an IAM—in this case, the Global Change Assessment Model (GCAM)—in a variety of formats and resolutions commonly used by ESMs. +

W

Depth-Discharge and Bedload Calculator, uses:
# Wright-Parker formulation for flow resistance (without stratification correction)
# Ashida-Michiue formulation for bedload transport. +

D

Depth-Discharge and Total Load Calculator, uses:
# Wright-Parker formulation for flow resistance,
# Ashida-Michiue formulation for bedload transport,
# Wright-Parker formulation (without stratification) for suspended load. +

M

Derived from MOSART-WM (Model for Scale Adaptive River Transport with Water Management), mosasrtwmpy is a large-scale river-routing Python model used to study riverine dynamics of water, energy, and biogeochemistry cycles across local, regional, and global scales. The water management component represents river regulation through reservoir storage and release operations, diversions from reservoir releases, and allocation to sectoral water demands. The model allows an evaluation of the impact of water management over multiple river basins at once (global and continental scales) with consistent representation of human operations over the full domain. +

D

F

Directs flow by the D infinity method (Tarboton, 1997). Each node is assigned two flow directions, toward the two neighboring nodes that are on the steepest subtriangle. Partitioning of flow is done based on the aspect of the subtriangle. +

Directs flow by the multiple flow direction method. Each node is assigned multiple flow directions, toward all of the N neighboring nodes that are lower than it. If none of the neighboring nodes are lower, the location is identified as a pit. Flow proportions can be calculated as proportional to slope or proportional to the square root of slope, which is the solution to a steady kinematic wave. +

D

Dorado is a Python package for simulating passive Lagrangian particle transport over flow-fields from any 2D shallow-water hydrodynamic model using a weighted random walk methodology. +

DynEarthSol3D (Dynamic Earth Solver in Three Dimensions) is a flexible, open-source finite element code that solves the momentum balance and the heat transfer in Lagrangian form using unstructured meshes. It can be used to study the long-term deformation of Earth's lithosphere and problems alike. +

DynQual is a high-spatio-temporal-resolution surface water quality model, which can be used to simulate water temperature; concentrations of total dissolved solids to represent salinity pollution; biological oxygen demand to represent organic pollution; and fecal coliform as a coarse indicator for pathogen pollution. +

E

ECSimpleSnow is a simple snow model that employs an empirical algorithm to melt or accumulate snow based on surface temperature and precipitation that has fallen since the previous analysis step. +

EF5 was created by the Hydrometeorology and Remote Sensing Laboratory at the University of Oklahoma. The goal of EF5 is to have a framework for distributed hydrologic modeling that is user friendly, adaptable, expandable, all while being suitable for large scale (e.g. continental scale) modeling of flash floods with rapid forecast updates. Currently EF5 incorporates 3 water balance models including the Sacramento Soil Moisture Accouning Model (SAC-SMA), Coupled Routing and Excess Storage (CREST), and hydrophobic (HP). These water balance models can be coupled with either linear reservoir or kinematic wave routing. +

ELCIRC is an unstructured-grid model designed for the effective simulation of 3D baroclinic circulation across river-to-ocean scales. It uses a finite-volume/finite-difference Eulerian-Lagrangian algorithm to solve the shallow water equations, written to realistically address a wide range of physical processes and of atmospheric, ocean and river forcings. The numerical algorithm is low-order, but volume conservative, stable and computationally efficient. It also naturally incorporates wetting and drying of tidal flats. ELCIRC has been extensively tested against standard ocean/coastal benchmarks, and is starting to be applied to estuaries and continental shelves around the world. +

Ecopath with Ecosim (EwE) is an ecological modeling software suite for personal computers. EwE has three main components: Ecopath – a static, mass-balanced snapshot of the system; Ecosim – a time dynamic simulation module for policy exploration; and Ecospace – a spatial and temporal dynamic module primarily designed for exploring impact and placement of protected areas. The Ecopath software package can be used to:
*Address ecological questions;
*Evaluate ecosystem effects of fishing;
*Explore management policy options;
*Evaluate impact and placement of marine protected areas;
*Evaluate effect of environmental changes. +

Erode is a raster-based, fluvial landscape evolution model. The newest version (3.0) is written in Python and contains html help pages when running the program through the CSDMS Modeling Tool CMT (https://csdms.colorado.edu/wiki/Help:Ccaffeine_GUI). +

Erode-D8-Global is a raster, D8-based fluvial landscape evolution model (LEM) +

L

Exposures to heat and sunlight can be simulated and the resulting signals shown. For a detailed description of the underlying luminescence rate equations, or to cite your use of LuSS, please use Brown, (2020). +

S

Extended description for SINUOUS - Meander Evolution Model. The basic model simulates planform evolution of a meandering river starting from X,Y coordinates of centerline nodes, with specification of cross-sectional and flow parameters. If the model is intended to simulate evolution of an existing river, the success of the model can be evaluated by the included area between the simulated and the river centerline. In addition, topographic evolution of the surrounding floodplain can be simulated as a function of existing elevation, distance from the nearest channel, and time since the channel migrated through that location. Profile evolution of the channel can also be modeled by backwater flow routing and bed sediment transport relationships. +

F

FACET is a Python tool that uses open source modules to map the floodplain extent and derive reach-scale summaries of stream and floodplain geomorphic measurements from high-resolution digital elevation models (DEMs). Geomorphic measurements include channel width, stream bank height, floodplain width, and stream slope.<br>Current tool functionality is only meant to process DEMs within the Chesapeake Bay and Delaware River watersheds. FACET was developed to batch process 3-m resolution DEMs in the Chesapeake Bay and Delaware River watersheds. Future updates to FACET will allow users to process DEMs outside of the Chesapeake and Delaware basins.<br>FACET allows the user to hydrologically condition the DEM, generate the stream network, select one of two options for stream bank identification, map the floodplain extent using a Height Above Nearest Drainage (HAND) approach, and calculate stream and floodplain metrics using three approaches. +

FUNWAVE is a phase-resolving, time-stepping Boussinesq model for ocean surface wave propagation in the nearshore. +

FVCOM is a prognostic, unstructured-grid, finite-volume, free-surface, 3-D primitive equation coastal ocean circulation model developed by UMASSD-WHOI joint efforts. The model consists of momentum, continuity, temperature, salinity and density equations and is closed physically and mathematically using turbulence closure submodels. The horizontal grid is comprised of unstructured triangular cells and the irregular bottom is preseented using generalized terrain-following coordinates. The General Ocean Turbulent Model (GOTM) developed by Burchard’s research group in Germany (Burchard, 2002) has been added to FVCOM to provide optional vertical turbulent closure schemes. FVCOM is solved numerically by a second-order accurate discrete flux calculation in the integral form of the governing equations over an unstructured triangular grid. This approach combines the best features of finite-element methods (grid flexibility) and finite-difference methods (numerical efficiency and code simplicity) and provides a much better numerical representation of both local and global momentum, mass, salt, heat, and tracer conservation. The ability of FVCOM to accurately solve scalar conservation equations in addition to the topological flexibility provided by unstructured meshes and the simplicity of the coding structure has make FVCOM ideally suited for many coastal and interdisciplinary scientific applications. +

Fall velocity for spheres. Uses formulation of Dietrich (1982) +

Z

Finite difference approximations are great for modeling the erosion of landscapes. A paper by Densmore, Ellis, and Anderson provides details on application of landscape evolution models to the Basin and Range (USA) using complex rulesets that include landslides, tectonic displacements, and physically-based algorithms for hillslope sediment transport and fluvial transport. The solution given here is greatly simplified, only including the 1D approximation of the diffusion equation. The parallel development of the code is meant to be used as a class exercise +

G

Finite difference solution allows for calculations of flexural response in regions of variable elastic thickness / flexural rigidity. The direct solution technique means that it takes time to populate a cofactor matrix, but that once this has been done, flexural solutions may be obtained rapidly via a Thomas algorithm. This makes it less good for an individual solution where an iterative approach may be more computationally efficient, but better for modeling where elastic thickness does not change (meaning that you do not need to create a new cofactor matrix) but loads do. +

S

Finite element process based simulation model for fluid flow, clastic, carbonate and evaporate sedimentation. +

For each time step, this component calculates an infiltration rate for a given model location and updates surface water depths. Based on the Green-Ampt method, it follows the form of Julien et al., 1995. +

M

Fortran 95 routines to model the ocean carbonate system (mocsy). Mocsy take as input dissolved inorganic carbon CT and total alkalinity AT, the only two tracers of the ocean carbonate system that are unaffected by changes in temperature and salinity and conservative with respect to mixing, properties that make them ideally suited for ocean carbon models. With basic thermodynamic equilibria, mocsy compute surface-ocean pCO2 in order to simulate air-sea CO2 fluxes. The mocsy package goes beyond the OCMIP code by computing all other carbonate system variables (e.g., pH, CO32-, and CaCO3 saturation states) and by doing so throughout the water column. +

F

FuzzyReef is a three-dimensional (3D) numerical stratigraphic model that simulates the development of microbial reefs using fuzzy logic (multi-valued logic) modeling methods. The flexibility of the model allows for the examination of a large number of variables. This model has been used to examine the importance of local environmental conditions and global changes on the frequency of reef development relative to the temporal and spatial constraints from Upper Jurassic (Oxfordian) Smackover reef datasets from two Alabama oil fields.
The fuzzy model simulates the deposition of reefs and carbonate facies through integration of local and global variables. Local-scale factors include basement relief, sea-level change, climate, latitude, water energy, water depth, background sedimentation rate, and substrate conditions. Regional and global-scale changes include relative sea-level change, climate, and latitude. +

G

GENESIS calculates shoreline change produced by statial and temporal differences in longshore sand transport produced by breaking waves. The shoreline evolution portion of the numerical modeling system is based on one-line shoreline change theory, which assumes that the beach profile shape remains unchanged, allowing shoreline change to be described uniquely in terms of the translation of a single point (for example, Mean High Water shoreline) on the profile. +

GEOMBEST is a morphological-behaviour model that simulates the evolution of coastal morphology and stratigraphy resulting from changes in sea level and sediment volume within the shoreface, barrier, and estuary. +

GEOMBEST++ is a morphological-behaviour model that simulates the evolution of coastal morphology and stratigraphy resulting from changes in sea level and sediment volume within the shoreface, barrier, and estuary. GEOMBEST++ builds on previous iterations (i.e. GEOMBEST+) by incorporating the effects of waves into the backbarrier, providing a more physical basis for the evolution of the bay bottom and introducing wave erosion of marsh edges. +

GEOMBEST++Seagrass is a morphological-behaviour model that simulates the evolution of coastal morphology and stratigraphy resulting from changes in sea level and sediment volume within the shoreface, barrier, and estuary. GEOMBEST++Seagrass builds on previous iterations (i.e. GEOMBEST, GEOMBEST+, and GEOMBEST++) by incorporating seagrass dynamics into the back-barrier bay. +

GEOtop accommodates very complex topography and, besides the water balance integrates all the terms in the surface energy balance equation. For saturated and unsaturated subsurface flow, it uses the 3D Richards’ equation. An accurate treatment of radiation inputs is implemented in order to be able to return surface temperature.
The model GEOtop simulates the complete hydrological balance in a continuous way, during a whole year, inside a basin and combines the main features of the modern land surfaces models with the distributed rainfall-runoff models.
The new 0.875 version of GEOtop introduces the snow accumulation and melt module and describes sub-surface flows in an unsaturated media more accurately. With respect to the version 0.750 the updates are fundamental: the codex is completely eviewed, the energy and mass parametrizations are rewritten, the input/output file set is redifined.
GEOtop makes it possible to know the outgoing discharge at the basin's closing section, to estimate the local values at the ground of humidity, of soil temperature, of sensible and latent heat fluxes, of heat flux in the soil and of net radiation, together with other hydrometeorlogical distributed variables. Furthermore it describes the distributed snow water equivalent and surface snow temperature.
GEOtop is a model based on the use of Digital Elevation Models (DEMs). It makes also use of meteorological measurements obtained thought traditional instruments on the ground. Yet, it can also assimilate distributed data like those coming from radar measurements, from satellite terrain sensing or from micrometeorological models. +

GIPL(Geophysical Institute Permafrost Laboratory) is an implicit finite difference one-dimensional heat flow numerical model. The GIPL model uses the effect of snow layer and subsurface soil thermal properties to simulate ground temperatures and active layer thickness (ALT) by solving the 1D heat diffusion equation with phase change. The phase change associated with freezing and thawing process occurs within a range of temperatures below 0 degree centigrade, and is represented by the unfrozen water curve (Romanovsky and Osterkamp 2000). The model employs finite difference numerical scheme over a specified domain. The soil column is divided into several layers, each with distinct thermo-physical properties. The GIPL model has been successfully used to map permafrost dynamics in Alaska and validated using ground temperature measurements in shallow boreholes across Alaska (Nicolsky et al. 2009, Jafarov et al. 2012, Jafarov et al. 2013, Jafarov et al. 2014). +

GSFLOW was a coupled model based on the integration of the U.S. Geological Survey Precipitation-Runoff Modeling System (PRMS, Leavesley and others, 1983) and the U.S. Geological Survey Modular Groundwater Flow Model(MODFLOW-2005, Harbaugh, 2005). It was developed to simulate coupled groundwater/surface-water flow in one or more watersheds by simultaneously simulating flow across the land surface, within subsurface saturated and unsaturated materials, and within streams and lakes. +

A

Generates alluvial stratigraphy by channel migration and avulsion. Channel migration is handled via a random walk. Avulsions occur when the channel superelevates. Channels can create levees. Post-avulsion channel locations chosen at random, or based on topography. +

G

GeoFlood, a new open-source software package for solving shallow water equations (SWE) on a quadtree hierarchy of mapped, logically Cartesian grids managed by the parallel, adaptive library ForestClaw +

Glimmer is an open source (GPL) three-dimensional thermomechanical ice sheet model, designed to be interfaced to a range of global climate models. It can also be run in stand-alone mode. Glimmer was developed as part of the NERC GENIE project (www.genie.ac.uk). It's development follows the theoretical basis found in Payne (1999) and Payne (2001). Glimmer's structure contains numerous software design strategies that make it maintainable, extensible, and well documented. +

Grain Size Distribution Statistics Calculator +

Gridded Surface Subsurface Hydrologic Analysis (GSSHA) is a grid-based two-dimensional hydrologic model. Features include 2D overland flow, 1D stream flow, 1D infiltration, 2D groundwater, and full coupling between the groundwater, vadoze zone, streams, and overland flow. GSSHA can run in both single event and long-term modes. The fully coupled groundwater to surfacewater interaction allows GSSHA to model both Hortonian and Non-Hortonian basins.
New features of version 2.0 include support for small lakes and detention basins, wetlands, improved sediment transport, and an improved stream flow model.
GSSHA has been successfully used to predict soil moistures as well as runoff and flooding. +

W

Gridded water balance model using climate input forcings that calculate surface and subsurface runoff and ground water recharge for each grid cell. The surface and subsurface runoff is propagated horizontally along a prescribed gridded network using Musking type horizontal transport. +

H

HYPE is a semi-distributed hydrological model for water and water quality. It simulates water and nutrient concentrations in the landscape at the catchment scale. Its spatial division is related to catchments and sub-catchments, land use or land cover, soil type and elevation. Within a catchment the model will simulate different compartments; soil including shallow groundwater, rivers and lakes. It is a dynamical model forced with time series of precipitation and air temperature, typically on a daily time step. Forcing in the form of nutrient loads is not dynamical. Example includes atmospheric deposition, fertilizers and waste water. +

E

Here, we present a Python tool that includes a comprehensive set of relations that predicts the hydrodynamics, bed elevation and the patterns of channels and bars in mere seconds. Predictions are based on a combination of empirical relations derived from natural estuaries, including a novel predictor for cross-sectional depth distributions, which is dependent on the along-channel width profile. Flow velocity, an important habitat characteristic, is calculated with a new correlation between depth below high water level and peak tidal flow velocity, which was based on spatial numerical modelling. Salinity is calculated from estuarine geometry and flow conditions. The tool only requires an along-channel width profile and tidal amplitude, making it useful for quick assessments, for example of potential habitat in ecology, when only remotely-sensed imagery is available. +