Numerical uncertainty at mesoscale in a Lagrangian model in complex terrain
Open Access
- 17 September 2012
- journal article
- Published by Copernicus GmbH in Geoscientific Model Development
- Vol. 5 (5), 1127-1136
- https://doi.org/10.5194/gmd-5-1127-2012
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.Keywords
This publication has 8 references indexed in Scilit:
- Meteorological Model Evaluation for CalNex 2010Monthly Weather Review, 2012
- A new inversion method to calculate emission inventories without a prior at mesoscale: Application to the anthropogenic CO2 emission from Houston, TexasJournal of Geophysical Research: Solid Earth, 2012
- Top-down estimate of anthropogenic emission inventories and their interannual variability in Houston using a mesoscale inverse modeling techniqueJournal of Geophysical Research, 2011
- Coupled weather research and forecasting–stochastic time-inverted lagrangian transport (WRF–STILT) modelArchiv für Meteorologie, Geophysik und Bioklimatologie Serie A, 2010
- Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2Atmospheric Chemistry and Physics, 2005
- Online versus offline air quality modeling on cloud‐resolving scalesGeophysical Research Letters, 2004
- A near‐field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time‐Inverted Lagrangian Transport (STILT) modelJournal of Geophysical Research: Solid Earth, 2003
- Computation, accuracy and applications of trajectories—A review and bibliographyAtmospheric Environment, 1998