Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 5.240
IF5.240
IF 5-year value: 5.768
IF 5-year
5.768
CiteScore value: 8.9
CiteScore
8.9
SNIP value: 1.713
SNIP1.713
IPP value: 5.53
IPP5.53
SJR value: 3.18
SJR3.18
Scimago H <br class='widget-line-break'>index value: 71
Scimago H
index
71
h5-index value: 51
h5-index51
Volume 5, issue 5
Geosci. Model Dev., 5, 1127–1136, 2012
https://doi.org/10.5194/gmd-5-1127-2012
© Author(s) 2012. This work is distributed under
the Creative Commons Attribution 3.0 License.
Geosci. Model Dev., 5, 1127–1136, 2012
https://doi.org/10.5194/gmd-5-1127-2012
© Author(s) 2012. This work is distributed under
the Creative Commons Attribution 3.0 License.

Development and technical paper 17 Sep 2012

Development and technical paper | 17 Sep 2012

Numerical uncertainty at mesoscale in a Lagrangian model in complex terrain

J. Brioude1,2, W. M. Angevine1,2, S. A. McKeen1,2, and E.-Y. Hsie1,2 J. Brioude et al.
  • 1Cooperative Institute for Research in Environmental Sciences, University of Colorado at Boulder, Boulder, Colorado, USA
  • 2Chemical Sciences Division, Earth System Research Laboratory, NOAA, Boulder, Colorado, USA

Abstract. Recently, it has been shown that mass conservation in Lagrangian models is improved by using time-average winds out of Eulerian models. In this study, we evaluate the mass conservation and trajectory uncertainties in complex terrain at mesoscale using the FLEXPART Lagrangian particle dispersion model coupled with the WRF mesoscale model. The specific form of vertical wind used is found to have a large effect. Time average wind with time average sigma dot (σ·), instantaneous wind with geometric cartesian vertical wind (w) and instantaneous wind with σ· are used to simulate mixing ratios of a passive tracer in forward and backward runs using different time interval outputs and horizontal resolutions in California. Mass conservation in the FLEXPART model was not an issue when using time-average wind or instantaneous wind with σ·. However, mass was poorly conserved using instantaneous wind with w, with a typical variation of 25% within 24 h.

Uncertainties in surface residence time (a backtrajectory product commonly used in source receptor studies or inverse modeling) calculated for each backtrajectory run were also analyzed. The smallest uncertainties were systematically found when using time-average wind. Uncertainties using instantaneous wind with σ· were slightly larger, as long as the time interval of output was sufficiently small. The largest uncertainties were found when using instantaneous wind with w. Those uncertainties were found to be linearly correlated with the local average gradient of orography. Differences in uncertainty were much smaller when trajectories were calculated over flat terrain. For a typical run at mesoscale in complex terrain, 4 km horizontal resolution and 1 h time interval output, the average uncertainty and bias in surface residence time is, respectively, 8.4% and −2.5% using time-average wind, and 13% and −3.7% using instantaneous wind with σ· in complex terrain. The corresponding values for instantaneous wind with cartesian w were 24% and −11%.

While the use of time-average wind systematically improves uncertainty in FLEXPART, the improvements are small, and therfore a systematic use of time-average wind in Lagrangian models is not necessarily required. Use of cartesian vertical wind in complex terrain, however, should be avoided.

Publications Copernicus
Download
Citation