Present the temperature of the Earth is 288 Kelvin. Without Earth would be 255 Kelvin. So the energy balance of the Earth add 33 Kelvin.

We optimize the variable GHG-effect and the optimal number is 0.29625 in this model.

Develop a daily time step simulation model that consists of a reservoir with a single inflow and single outflow (release)

Use units of million m3 and include any necessary parameters (e.g., capacity k) as separate adjustable variables

Implement the standard linear operating policy (SLOP). 

Assume the reservoir is 500 mcm (k=500). 

Develop yield-reliability results for a target (T) delivery values of 1, 3, 5, and 7 mcm/day.3 

(The mean inflow for the time series is 34.8 m3/s, or 3.0 million m3/day.)

The standard linear operating policy provides a basic rule for reservoir release. 

STEM-SM combines a simple ecosystem model (modified version of VSEM; Hartig et al. 2019) with a soil moisture model (Guswa et al. (2002) leaky bucket model). Outputs from the soil moisture model influence ecosystem dynamics in three ways. 

(1) The ratio of actual transpiration to maximum evapotranspiration (T/ETmax) modifies gross primary productivity (GPP).

(2) Degree of saturation of the soil (Sd) modifies the rate of soil heterotrophic respiration.

(3) Water limitation of GPP (by T/ETmax) and of soil nutrient availability (approximated by Sd) combine with leaf area limitation (approximated by fraction of incident photosynthetically-active radiation that is absorbed) to modify the allocation of net primary productivity to aboveground and belowground parts of the vegetation.

Ecosystem dynamics in turn influence flows of water in to and out of the soil moisture stock. The size of the aboveground biomass stock determines fractional vegetation cover, which modifies interception, soil evaporation and transpiration by plants.

STEM is a modified implementation of Hartig et al.'s (2019) Very Simple Ecosystem Model (VSEM). The vegetation part of the model has two stocks of biomass carbon (C): aboveground C and belowground C.  The soil part of the model has a single stock of soil organic C. Carbon flows into the biomass C stocks via net primary productivity (NPP). Carbon flows out of these stocks and into the soil organic C stock via the loss of aboveground/belowground C through senescence (i.e., abscission of dead leaves and roots). SOC loss is due to heterotrophic respiration of the soil organic matter.

A 'leaky bucket' soil moisture model, based on Guswa et al. (2002). Rain falls as discrete events. The mean depth and frequency of rainfall events are determined by total monthly rainfall and number of rain days. A portion of the rainfall is intercepted by vegetation and evaporates before reaching the soil. The remaining rainfall (throughflow) either infiltrates the soil or, if the soil has insufficient capacity, runs off immediately. Soil water exceeding the field capacity is lost by sub-surface leakage, at a rate determined by the degree of soil saturation. Degree of soil saturation also limits rates of soil evaporation and vegetation transpiration. The partitioning between evaporation and transpiration is influenced by fractional area covered by vegetation.

Reference:

Guswa, A.J., Celia, M.A., Rodriguez-Iturbe, I. (2002) Models of soil moisture dynamics in ecohydrology: a comparative study. Water Resources Research 38, 5-1 - 5-15.