## MOSART (Model for Scale Adaptive River Transport)

MOSART is a river transport model designed for applications across local, regional and global scales

river transportlocalregionalglobalscales
427

#### Authorship

Affiliation:
Pacific Northwest National Laboratory
Is authorship not correct? Feed back

#### Classification(s)

Application-focused categoriesNatural-perspectiveLand regions

#### Model Description

English {{currentDetailLanguage}} English

MOSART is a river transport model designed for applications across local, regional and global scales (Li et al., 2013b). A major purpose of MOSART is to provide freshwater input for the ocean model in coupled Earth System Models. MOSART also provides an effective way of evaluating and diagnosing the soil hydrology simulated by land surface models through direct comparison of the simulated river flow with observations of natural streamflow at gauging stations (Li et al., 2015a). Moreover, MOSART provides a modeling framework for representing riverine transport and transformation of energy and biogeochemical fluxes under both natural and human-influenced conditions ( (Li et al., 2015b).

## Routing Processes

MOSART divides each spatial unit such as a lat/lon grid or watershed into three categories of hydrologic units (as shown in Figure 2.14.1): hillslopes that convert both surface and subsurface runoff into tributaries, tributaries that discharge into a single main channel, and the main channel that connects the local spatial unit with upstream/downstream units through the river network. MOSART assumes that all the tributaries within a spatial unit can be treated as a single hypothetical sub-network channel with a transport capacity equivalent to all the tributaries combined. Correspondingly, three routing processes are represented in MOSART: 1) hillslope routing: in each spatial unit, surface runoff is routed as overland flow into the sub-network channel, while subsurface runoff generated in the spatial unit directly enters the sub-network channel; 2) sub-network channel routing: the sub-network channel receives water from the hillslopes, routes water through the channel and discharges it into the main channel; 3) main channel routing: the main channel receives water from the sub-network channel and/or inflow, if any, from the upstream spatial units, and discharges the water to its downstream spatial unit or the ocean.

MOSART only routes positive runoff, although negative runoff can be generated occasionally by the land model (e.g., ). Negative runoff in any runoff component including  is not routed through MOSART, but instead is mapped directly from the spatial unit where it is generated at any time step to the coupler.

In MOSART, the travel velocities of water across hillslopes, sub-network and main channel are all estimated using Manning’s equation with different levels of simplifications. Generally the Manning’s equation is in the form of

(2.14.1)

where  is the travel velocity (m s -1 ),  is the hydraulic radius (m).  is the friction slope that accounts for the effects of gravity, friction, inertia and other forces on the water. If the channel slope is steep enough, the gravity force dominates over the others so one can approximate  by the channel bed slope  , which is the key assumption underpinning the kinematic wave method.  is the Manning’s roughness coefficient, which is mainly controlled by surface roughness and sinuosity of the flow path.

If the water surface is sufficiently large or the water depth  is sufficiently shallow, the hydraulic radius can be approximated by the water depth. This is the case for both hillslope and sub-network channel routing.

(2.14.2)

Here  (m) and  (m) are hydraulic radius for hillslope and sub-network channel routing respectively, and  (m) and  (m) are water depth during hillslope and sub-network channel routing respectively.

For the main channel, the hydraulic radius is given by

(2.14.3)

where  (m 2 ) is the wetted area defined as the part of the channel cross-section area below the water surface,  (m) is the wetted perimeter, the perimeter confined in the wetted area.

For hillslopes, sub-network and main channels, a common continuity equation can be written as

(2.14.4)

where  (m 3 s -1 ) is the main channel flow from the upstream grid(s) into the main channel of the current grid, which is zero for hillslope and sub-network routing.  (m 3 s -1 ) is the outflow rate from hillslope into the sub-network, from the sub-network into the main channel, or from the current main channel to the main channel of its downstream grid (if not the outlet grid) or ocean (if the current grid is the basin outlet).  (m 3 s -1 ) is a source term, which could be the surface runoff generation rate for hillslopes, or lateral inflow (from hillslopes) into sub-network channel or water-atmosphere exchange fluxes such as precipitation and evaporation. It is assumed that surface runoff is generated uniformly across all the hillslopes. Currently, MOSART does not exchange water with the atmosphere or return water to the land model so its function is strictly to transport water from runoff generation through the hillslope, tributaries, and main channels to the basin outlets.

## Numerical Solution

The numerical implementation of MOSART is mainly based on a subcycling scheme and a local time-stepping algorithm. There are two levels of subcycling. For convenience, we denote  (s),  (s),  (s) and  (s) as the time steps of runoff inputs (from CLM to MOSART via the flux coupler), MOSART routing, hillslope routing, and channel routing, respectively. The first level of subcycling is between the runoff inputs and MOSART routing. If  is 10800s and  is 3600s, three MOSART time steps will be invoked each time the runoff inputs are updated. The second level of subcycling is between the hillslope routing and channel routing. This is to account for the fact that the travel velocity of water across hillslopes is usually much slower than that in the channels.  is usually set as the same as , but within each time step of hillslope routing there are a few time steps for channel routing, i.e., . The local time-stepping algorithm is to account for the fact that the travel velocity of water is much faster in some river channels (e.g., with steeper bed slope, narrower channel width) than others. That is, for each channel (either a sub-network or main channel), the final time step of local channel routing is given as  is currently estimated empirically as a function of local channel slope, width, length and upstream drainage area. If MOSART crashes due to a numerical issue, we recommend increasing  and, if the issue remains, reducing .

## Parameters and Input Data

MOSART is supported by a comprehensive, global hydrography dataset at 0.5 o resolution. As such, the fundamental spatial unit of MOSART is a 0.5 o lat/lon grid. The topographic parameters (such as flow direction, channel length, topographic and channel slopes, etc.) were derived using the Dominant River Tracing (DRT) algorithm (Wu et al., 2011 ; Wu et al. 2012). The DRT algorithm produces the topographic parameters in a scale-consistent way to preserve/upscale the key features of a baseline high-resolution hydrography dataset at multiple coarser spatial resolutions. Here the baseline high-resolution hydrography dataset is the 1km resolution Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS) (Lehner and Döll, 2004 ; Lehner et al., 2008). The channel geometry parameters, e.g., bankfull width and depth, were estimated from empirical hydraulic geometry relationships as functions of the mean annual discharge. The Manning roughness coefficients for overland and channel flow were calculated as functions of landcover and water depth. For more details on the methodology to derive channel geometry and the Manning’s roughness coefficients, please refer to Getirana et al. (2012) . The full list of parameters included in this global hydrography dataset is provided in Table 2.14.1. Evaluation of global simulations by MOSART using the aforementioned parameters is described in Li et al. (2015b) .

Table 2.14.1 List of parameters in the global hydrography dataset
Name Unit Description
- The D8 single flow direction for each coarse grid cell coded using 1 (E), 2 (SE), 4 (S), 8 (SW), 16 (W), 32 (NW), 64 (N), 128 (NE)
km 2 The upstream drainage area of each coarse grid cell
m The dominant river length for each coarse grid cell
- The average channel slope for each coarse grid cell
- The average topographic slope (for overland flow routing) for each coarse grid cell
km 2 The surface area for each coarse grid cell
-1 Drainage density, calculated as the total channel length within each coarse grid cell divided by the local cell area
m The bankfull depth of main channel
m The bankfull width of main channel
m The average bankfull depth of tributary channels
m The average bankfull width of tributary channels
- Manning’s roughness coefficient for channel flow routing
- Manning’s roughness coefficient for overland flow routing

Model Domain
{{domain}}

Principles
{{principle}}
Incorporated Models
{{incorporatedModel}}
Model part of larger framework: {{metadata.design.framework}}
Incorporated Models
{{process}}

Inputs
{{input}}
Outputs
{{output}}

#### How to Cite

Hongyi Li et al. (2019). MOSART (Model for Scale Adaptive River Transport), Model Item, OpenGMS, https://geomodeling.njnu.edu.cn/modelItem/04b97595-bbae-4198-bf66-2aee232b0b6e

#### Comment(s)

{{comment.date}}
{{comment.content}}
{{subComm.date}}
{{subComm.content}}

#### Authorship

Affiliation:
Pacific Northwest National Laboratory
Is authorship not correct? Feedback

#### QR Code

• {{curRelation.name}}
{{curRelation.name}}

{{curRelation.overview}}
{{curRelation.author.join('; ')}}
{{curRelation.journal}}

You can link related {{typeName}} from repository to this model item, or you can create a new {{typeName.toLowerCase()}}.

Related Items
{{ props.row.description }}
{{ props.row.description }}
Related Items
{{props.row.name}}

You can link resource from repository to this model item, or you can create a new {{typeName.toLowerCase()}}.

{{ props.row.description }}
{{ props.row.description }}
Drop the file here, orclick to upload.
File size should not exceed 10m.
Select From My Space

These authorship information will be submitted to the contributor to review.

Cancel Submit
Model Classifications
Cancel Submit
{{ item.label }} {{ item.value }}
{{props.row.localName}}
Model Name :
Cancel Submit
Name:
Version:
Model Type:
Model Domain:
Scale:
Purpose:
Principles:
Incorporated models:

Model part of

larger framework

Process:
Information:
Initialization:
Hardware Requirements:
Software Requirements:
Inputs:
Outputs:
Cancel Submit
Title Author Date Journal Volume(Issue) Pages Links Doi Operation
Cancel Submit

Yes, this is it Cancel

OK
Cancel Confirm
Model Classifications 1
Model Classifications 2
Title Author Date Journal Volume(Issue) Pages Links Doi Operation

#### NEW

Name:
Affiliation:
Email:
Homepage:

Yes, this is it Cancel

Confirm
path
:
/{{path.label}}
search results of '{{searchContentShown}}'

#### No content to show

{{item.label}}

.

{{item.suffix}}

.{{item.suffix}}

{{item.fileName}}.{{item.suffix}}

Copy
Delete
Rename
/{{path.label}}
Change
/{{path.label}}
Select File
Cancel Confirm
path
:
/{{path.label}}
/..
Cancel Confirm
{{ node.label }}
##### You have select  {{multipleSelection.length+multipleSelectionMyData.length}} data .
• Output Data
• {{item.computableName}}@{{formatDate(item.runTime)}}
{{scope.row.type}}
{{ scope.row.tag }}