Title: 2.32. Land-Only Mode — ctsm CTSM master documentation

URL Source: https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Land-Only_Mode/CLM50_Tech_Note_Land-Only_Mode.html

Markdown Content:
In land-only mode (uncoupled to an atmospheric model), the atmospheric forcing required by CLM ([Table 2.2.4](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-atmospheric-input-to-land-model)) is supplied by observed datasets. The standard forcing provided with the model is a 110-year (1901-2010) dataset provided by the Global Soil Wetness Project (GSWP3; NEED A REFERENCE). The GSWP3 dataset has a spatial resolution of 0.5° X 0.5° and a temporal resolution of three hours.

An alternative forcing dataset is also available, CRUNCEP, a 110-year (1901-2010) dataset (CRUNCEP; [Viovy 2011](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#viovy2011)) that is a combination of two existing datasets; the CRU TS3.2 0.5° X 0.5° monthly data covering the period 1901 to 2002 ([Mitchell and Jones 2005](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#mitchelljones2005)) and the NCEP reanalysis 2.5° X 2.5° 6-hourly data covering the period 1948 to 2010. The CRUNCEP dataset has been used to force CLM for studies of vegetation growth, evapotranspiration, and gross primary production ([Mao et al. 2012](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#maoetal2012), [Mao et al. 2013](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#maoetal2013), [Shi et al. 2013](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#shietal2013)) and for the TRENDY (trends in net land-atmosphere carbon exchange over the period 1980-2010) project ([Piao et al. 2012](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#piaoetal2012)). Version 7 is available here ([Viovy 2011](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#viovy2011)).

Here, the GSWP3 dataset, which does not include data for particular fields over oceans, lakes, and Antarctica is modified. This missing data is filled with [Qian et al. (2006)](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#qianetal2006) data from 1948 that is interpolated by the data atmosphere model to the 0.5° GSWP3 grid. This allows the model to be run over Antarctica and ensures data is available along coastlines regardless of model resolution.

The forcing data is ingested into a data atmosphere model in three “streams”; precipitation (\\(P\\)) (mm s\-1), solar radiation (\\(S\_{atm}\\) ) (W m\-2), and four other fields \[atmospheric pressure \\(P\_{atm}\\) (Pa), atmospheric specific humidity \\(q\_{atm}\\) (kg kg\-1), atmospheric temperature \\(T\_{atm}\\) (K), and atmospheric wind \\(W\_{atm}\\) (m s\-1)\]. These are separate streams because they are handled differently according to the type of field. In the GSWP3 dataset, the precipitation stream is provided at three hour intervals and the data atmosphere model prescribes the same precipitation rate for each model time step within the three hour period. The four fields that are grouped together in another stream (pressure, humidity, temperature, and wind) are provided at three hour intervals and the data atmosphere model linearly interpolates these fields to the time step of the model.

The total solar radiation is also provided at three hour intervals. The data is fit to the model time step using a diurnal function that depends on the cosine of the solar zenith angle \\(\\mu\\) to provide a smoother diurnal cycle of solar radiation and to ensure that all of the solar radiation supplied by the three-hourly forcing data is actually used. The solar radiation at model time step \\(t\_{M}\\) is

(2.32.1)[¶](#equation-31-1 "Permalink to this equation")\\\[\\begin{split}\\begin{array}{lr} S\_{atm} \\left(t\_{M} \\right)=\\frac{\\frac{\\Delta t\_{FD} }{\\Delta t\_{M} } S\_{atm} \\left(t\_{FD} \\right)\\mu \\left(t\_{M} \\right)}{\\sum \_{i=1}^{\\frac{\\Delta t\_{FD} }{\\Delta t\_{M} } }\\mu \\left(t\_{M\_{i} } \\right) } & \\qquad {\\rm for\\; }\\mu \\left(t\_{M} \\right)>0.001 \\\\ S\_{atm} \\left(t\_{M} \\right)=0 & \\qquad {\\rm for\\; }\\mu \\left(t\_{M} \\right)\\le 0.001 \\end{array}\\end{split}\\\]

where \\(\\Delta t\_{FD}\\) is the time step of the forcing data (3 hours \\(\\times\\) 3600 seconds hour\-1 = 10800 seconds), \\(\\Delta t\_{M}\\) is the model time step (seconds), \\(S\_{atm} \\left(t\_{FD} \\right)\\) is the three-hourly solar radiation from the forcing data (W m\-2), and \\(\\mu \\left(t\_{M} \\right)\\) is the cosine of the solar zenith angle at model time step \\(t\_{M}\\) (section [2.3.3](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Surface_Albedos/CLM50_Tech_Note_Surface_Albedos.html#solar-zenith-angle)). The term in the denominator of equation [(2.32.1)](#equation-31-1) is the sum of the cosine of the solar zenith angle for each model time step falling within the three hour period. For numerical purposes, \\(\\mu \\left(t\_{M\_{i} } \\right)\\ge 0.001\\).

The total incident solar radiation \\(S\_{atm}\\) at the model time step \\(t\_{M}\\) is then split into near-infrared and visible radiation and partitioned into direct and diffuse according to factors derived from one year’s worth of hourly CAM output from CAM version cam3\_5\_55 as

(2.32.2)[¶](#equation-31-2 "Permalink to this equation")\\\[S\_{atm} \\, \\downarrow \_{vis}^{\\mu } =R\_{vis} \\left(\\alpha S\_{atm} \\right)\\\]

(2.32.3)[¶](#equation-31-3 "Permalink to this equation")\\\[S\_{atm} \\, \\downarrow \_{nir}^{\\mu } =R\_{nir} \\left\[\\left(1-\\alpha \\right)S\_{atm} \\right\]\\\]

(2.32.4)[¶](#equation-31-4 "Permalink to this equation")\\\[S\_{atm} \\, \\downarrow \_{vis} =\\left(1-R\_{vis} \\right)\\left(\\alpha S\_{atm} \\right)\\\]

(2.32.5)[¶](#equation-31-5 "Permalink to this equation")\\\[S\_{atm} \\, \\downarrow \_{nir} =\\left(1-R\_{nir} \\right)\\left\[\\left(1-\\alpha \\right)S\_{atm} \\right\].\\\]

where \\(\\alpha\\), the ratio of visible to total incident solar radiation, is assumed to be

(2.32.6)[¶](#equation-31-6 "Permalink to this equation")\\\[\\alpha =\\frac{S\_{atm} \\, \\downarrow \_{vis}^{\\mu } +S\_{atm} \\, \\downarrow \_{vis}^{} }{S\_{atm} } =0.5.\\\]

The ratio of direct to total incident radiation in the visible \\(R\_{vis}\\) is

(2.32.7)[¶](#equation-31-7 "Permalink to this equation")\\\[R\_{vis} =a\_{0} +a\_{1} \\times \\alpha S\_{atm} +a\_{2} \\times \\left(\\alpha S\_{atm} \\right)^{2} +a\_{3} \\times \\left(\\alpha S\_{atm} \\right)^{3} \\qquad 0.01\\le R\_{vis} \\le 0.99\\\]

and in the near-infrared \\(R\_{nir}\\) is

(2.32.8)[¶](#equation-31-8 "Permalink to this equation")\\\[R\_{nir} =b\_{0} +b\_{1} \\times \\left(1-\\alpha \\right)S\_{atm} +b\_{2} \\times \\left\[\\left(1-\\alpha \\right)S\_{atm} \\right\]^{2} +b\_{3} \\times \\left\[\\left(1-\\alpha \\right)S\_{atm} \\right\]^{3} \\qquad 0.01\\le R\_{nir} \\le 0.99\\\]

where \\(a\_{0} =0.17639,\\, a\_{1} =0.00380,\\, a\_{2} =-9.0039\\times 10^{-6},\\, a\_{3} =8.1351\\times 10^{-9}\\) and \\(b\_{0} =0.29548,b\_{1} =0.00504,b\_{2} =-1.4957\\times 10^{-5},b\_{3} =1.4881\\times 10^{-8}\\) are coefficients from polynomial fits to the CAM data.

The additional atmospheric forcing variables required by [Table 2.2.4](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-atmospheric-input-to-land-model) are derived as follows. The atmospheric reference height \\(z'\_{atm}\\) (m) is set to 30 m. The directional wind components are derived as \\(u\_{atm} =v\_{atm} ={W\_{atm} \\mathord{\\left/ {\\vphantom {W\_{atm} \\sqrt{2} }} \\right.} \\sqrt{2} }\\). The potential temperature \\(\\overline{\\theta \_{atm} }\\) (K) is set to the atmospheric temperature \\(T\_{atm}\\). The atmospheric longwave radiation \\(L\_{atm} \\, \\downarrow\\) (W m\-2) is derived from the atmospheric vapor pressure \\(e\_{atm}\\) and temperature \\(T\_{atm}\\) ([Idso 1981](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#idso1981)) as

(2.32.9)[¶](#equation-31-9 "Permalink to this equation")\\\[L\_{atm} \\, \\downarrow =\\left\[0.70+5.95\\times 10^{-5} \\times 0.01e\_{atm} \\exp \\left(\\frac{1500}{T\_{atm} } \\right)\\right\]\\sigma T\_{atm}^{4}\\\]

where

(2.32.10)[¶](#equation-31-10 "Permalink to this equation")\\\[e\_{atm} =\\frac{P\_{atm} q\_{atm} }{0.622+0.378q\_{atm} }\\\]

and \\(\\sigma\\) is the Stefan-Boltzmann constant (W m\-2 K\-4) ([Table 2.2.7](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-physical-constants)). The fraction of precipitation \\(P\\) (mm s\-1) falling as rain and/or snow is

(2.32.11)[¶](#equation-31-11 "Permalink to this equation")\\\[q\_{rain} =P\\left(f\_{P} \\right),\\\]

(2.32.12)[¶](#equation-31-12 "Permalink to this equation")\\\[q\_{snow} =P\\left(1-f\_{P} \\right)\\\]

where

(2.32.13)[¶](#equation-31-13 "Permalink to this equation")\\\[f\_{P} =0<0.5\\left(T\_{atm} -T\_{f} \\right)<1.\\\]

The aerosol deposition rates \\(D\_{sp}\\) (14 rates as described in [Table 2.2.4](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-atmospheric-input-to-land-model)) are provided by a time-varying, globally-gridded aerosol deposition file developed by [Lamarque et al. (2010)](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#lamarqueetal2010).

If the user wishes to provide atmospheric forcing data from another source, the data format outlined above will need to be followed with the following exceptions. The data atmosphere model will accept a user-supplied relative humidity \\(RH\\) (%) and derive specific humidity \\(q\_{atm}\\) (kg kg\-1) from

(2.32.14)[¶](#equation-31-14 "Permalink to this equation")\\\[q\_{atm} =\\frac{0.622e\_{atm} }{P\_{atm} -0.378e\_{atm} }\\\]

where the atmospheric vapor pressure \\(e\_{atm}\\) (Pa) is derived from the water (\\(T\_{atm} >T\_{f}\\) ) or ice (\\(T\_{atm} \\le T\_{f}\\) ) saturation vapor pressure \\(e\_{sat}^{T\_{atm} }\\) as \\(e\_{atm} =\\frac{RH}{100} e\_{sat}^{T\_{atm} }\\) where \\(T\_{f}\\) is the freezing temperature of water (K) ([Table 2.2.7](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/Ecosystem/CLM50_Tech_Note_Ecosystem.html#table-physical-constants)), and \\(P\_{atm}\\) is the pressure at height \\(z\_{atm}\\) (Pa). The data atmosphere model will also accept a user-supplied dew point temperature \\(T\_{dew}\\) (K) and derive specific humidity \\(q\_{atm}\\) from

(2.32.15)[¶](#equation-31-15 "Permalink to this equation")\\\[q\_{atm} = \\frac{0.622e\_{sat}^{T\_{dew} } }{P\_{atm} -0.378e\_{sat}^{T\_{dew} } } .\\\]

Here, \\(e\_{sat}^{T}\\), the saturation vapor pressure as a function of temperature, is derived from [Lowe’s (1977)](https://escomp.github.io/ctsm-docs/versions/master/html/tech_note/References/CLM50_Tech_Note_References.html#lowe1977) polynomials. If not provided by the user, the atmospheric pressure \\(P\_{atm}\\) (Pa) is set equal to the standard atmospheric pressure \\(P\_{std} =101325\\) Pa, and surface pressure \\(P\_{srf}\\) (Pa) is set equal to\\(P\_{atm}\\).

The user may provide the total direct and diffuse solar radiation, \\(S\_{atm} \\, \\downarrow ^{\\mu }\\) and \\(S\_{atm} \\, \\downarrow\\). These will be time-interpolated using the procedure described above and then each term equally apportioned into the visible and near-infrared wavebands (e.g., \\(S\_{atm} \\, \\downarrow \_{vis}^{\\mu } =0.5S\_{atm} \\, \\downarrow ^{\\mu }\\), \\(S\_{atm} \\, \\downarrow \_{nir}^{\\mu } =0.5S\_{atm} \\, \\downarrow ^{\\mu }\\) ).

