17.6 Configuring STILT options in HYSPLIT




An earlier version of the HYSPLIT code was modified by Lin et al. (2003 - JGR, VOL. 108, NO. D16, 4493, doi:10.1029/2002JD003161) to optimize the application of the model to determine the CO2 fluxes over a large area based upon a backward advection-dispersion calculation from a few measurement points. This HYSPLIT modification is called the Stochastic Time-Inverted Lagrangian Transport (STILT) model and it is used to calculate a footprint field [ppm/(micro-mole/m2/s)] to represent the flux:

where m is the molar mass of air, ρ is the average density below h, and h is half the mixed layer depth. The summation represents the time that N particles spends in the lower half of the mixed layer. The use of STILT for inverse footprint calculations will not be addressed in this tutorial because more detailed instructions are available elsewhere.

Configuration: When HYSPLIT is configured in the STILT mode (ichem=8) the mass summation is divided by air density resulting in a mixing ratio output field and the lowest concentration summation layer (concentration layer top depth) is permitted to vary with the mixed layer depth. Although any option is permitted, if not already defined in SETUP.CFG, several other namelist parameters are automatically set as follows:

  • Mass consistent dispersion,
  • IDSP=2
  • Hanna turbulence equations,
  • KBLT=5
  • Modified Richardson number,
  • KMIXD=3
  • Variable Lagrangian time scale,
  • VSCALES=-1.0
  • Grell convective mixing,
  • CAPEMIN=-2.0
Convective mixing is not discussed in this tutorial because it is covered as on of the topics in the advanced tutorial. Briefly, a value of -1 for CAPEMIN, the default, skips the CAPE computation. If CAPE is not available in the input meteorological data file then it is not computed. When the value is greater than zero, CAPE enhanced mixing is turned on and if CAPE is not available in the input data file, it is computed for each grid point, which may substantially slow down the overall calculation. CAPE enhanced mixing results in particles in the cloud-layer being randomly redistributed within the cloud-layer if they reside in a grid cell where CAPE exceeds CAPEMIN. A initial value of 500 J/kg is suggested for this computation. Setting CAPEMIN to -2 will utilize the Grell convective parameterization, if Grell convective fluxes are in the meteorological input file. Setting CAPEMIN to -3 will utilize the extreme convection parameterization. The Grell convective scheme redistributes the particles differently that the default HYSPLIT scheme.

Output: When running in STILT mode, the source surface flux values at each time step (or interval OUTDT) will be output to files PARTICLE.DAT and PARTICLE_STILT.DAT. The file PARTICLE.DAT is an earlier version with the output variables fixed. While the output variables in PARTICLE_STILT.DAT are defined by the user through the namelist variables IVMAX and VARSIWANT. The surface flux output of PARTICLE.DAT represents particles that were below 50% of the mixed layer height while the output of PARTICLE_STILT.DAT represent particles that were below a user defined height (VEGHT). Some additional STILT specfic namelist parameters are discussed in the model documentation provided with the code:

  • fraction of the Lagrangian timescale
  • NTURB,
  • turbulence on/off
  • VEGHT,
  • fraction of PBL for particle tally
  • height error flag for PBL
  • wind error flag
  • IVMAX,
  • number variables to print
  • variable list to print