Volume 1 - The Atmospheric Studies Group at TRC

Volume 1 - The Atmospheric Studies Group at TRC Volume 1 - The Atmospheric Studies Group at TRC

29.01.2015 Views

(IMIXH = 1), or the newly implemented Batchvarova-Gryning model is selected (IMIXH = 2). CALMET must be started no later than the hour that ends at 0500 in the base time zone (CET) and it does not need to create meteorological fields for hours after the end of sampling, since no CALPUFF simulations are made after the sampling terminates. Therefore each experiment-day is simulated in a separate CALMET application. CALPUFF Configuration Cameron, Carpinteria, Pismo Beach, Ventura CALPUFF is applied to each of the OCD4 datasets in the same way, with a unit emission rate (1 g/s), starting one hour before the first sampling period. The peak simulated concentration each evaluation hour is retained for comparison with the corresponding peak observed concentration, scaled by the tracer emission rate (e.g., Χ/Q). Tracer emissions are characterized as non-reacting passive Gaussian puffs with no plume rise. Default plume-path-coefficient terrain adjustments are used, and ISClike downwash adjustments are applied to those releases made from boats and modeled with downwash in the OCD evaluations. None of these releases is from a platform with significant structures, so the new platform downwash module is not tested in these applications. Turbulence-based dispersion is selected, with the default transition to Heffter curves for σ y at σ y = 550m. For each of the five CALMET configurations, 16 CALPUFF configurations are run to identify the performance impact of using all combinations of: • Measured versus predicted lateral turbulence (I y ) • CALPUFF versus AERMOD turbulence profiling assumptions • Draxler (1976) F y curves for sigma-y growth with a fixed Lagrangian timescale versus a SCIPUFF-like computed Lagrangian timescale • Minimum σ v = 0.5 m/s (CALPUFF default) versus σ v = 0.37 m/s (OCD default) Oresund CALPUFF is applied to each experiment-day in the Oresund dataset with the actual emission rate, for the actual release period (in whole hours). The resulting hourly concentrations are postprocessed using CALPOST to select only those concentrations Final Report Vol.1 58

that are simulated during the one-hour sampling period. The peak simulated concentration on each sampling arc is retained for comparison with the corresponding peak observed concentration. Because the sampling period spans two simulation hours on May 16 th (1330-1430) and May 18 th (1220-1320), the appropriate simulation hour is ambiguous. For these two experiments, the peak one-hour concentration from the two-hour simulation period that contains the sampling hour is selected for each arc. Tracer emissions are characterized as non-reacting passive Gaussian puffs with no plume rise. Turbulence-based dispersion is selected, with the default transition to Heffter curves for σ y at σ y = 550m. The Draxler F y function is selected for lateral cloud growth with the standard CALPUFF turbulence profiling assumptions. Default plume-path-coefficient terrain adjustments are used. The subgrid TIBL option is used here to better resolve its onshore development because the grid resolution is 1 km. For each of the 4 CALMET configurations, two CALPUFF configurations are run to identify the performance impact of using a minimum σ v = 0.5 m/s (CALPUFF default) versus σ v = 0.37 m/s (OCD default). Evaluation Statistics The evaluation focuses on the peak observed concentrations each hour, paired in time but not in space. As in the OCD applications, the wind direction used in the CALMET input files for the OCD4 datasets (Carpinteria, Cameron, Pismo Beach and Ventura) is that which is directed from the source to the sampler with the peak concentration. Unlike OCD, this may not align the center of the CALPUFF footprint with a sampler location due to changes in the transport direction (CALPUFF trajectories are not straight lines). Therefore, additional receptors are placed along the sampling arcs in the OCD4 datasets to better resolve peak modeled concentrations within the sampling zone. Spacing for additional receptors along sampling lines is typically 100 m or less. For the Oresund simulations, wind directions are used as measured (adjusted for grid North), and no additional receptors are used. Model predictions are quantitatively compared against measured concentrations and the resulting performance measures are compared with those for OCD5 (Carpinteria, Cameron, Pismo Beach and Ventura). The original OCD4 input files are recast for OCD5, and the revised OCD model is rerun on these data. Using the peak observed (Co) and peak modeled (Cp) concentrations each hour, the evaluation measures are computed from ln(Co/Cp). This gives equal weight to overpredictions and underpredictions (e.g., ln(Co/Cp) for Co/Cp = ½ and two are equal in magnitude but of opposite sign). The performance measures are the geometric mean (MG), Final Report Vol.1 59

th<strong>at</strong> are simul<strong>at</strong>ed during the one-hour sampling period. <strong>The</strong> peak simul<strong>at</strong>ed<br />

concentr<strong>at</strong>ion on each sampling arc is retained for comparison with the corresponding<br />

peak observed concentr<strong>at</strong>ion. Because the sampling period spans two simul<strong>at</strong>ion<br />

hours on May 16 th (1330-1430) and May 18 th (1220-1320), the appropri<strong>at</strong>e simul<strong>at</strong>ion<br />

hour is ambiguous. For these two experiments, the peak one-hour concentr<strong>at</strong>ion from<br />

the two-hour simul<strong>at</strong>ion period th<strong>at</strong> contains the sampling hour is selected for each<br />

arc.<br />

Tracer emissions are characterized as non-reacting passive Gaussian puffs with no<br />

plume rise. Turbulence-based dispersion is selected, with the default transition to<br />

Heffter curves for σ y <strong>at</strong> σ y = 550m. <strong>The</strong> Draxler F y function is selected for l<strong>at</strong>eral<br />

cloud growth with the standard CALPUFF turbulence profiling assumptions. Default<br />

plume-p<strong>at</strong>h-coefficient terrain adjustments are used. <strong>The</strong> subgrid TIBL option is<br />

used here to better resolve its onshore development because the grid resolution is 1<br />

km.<br />

For each of the 4 CALMET configur<strong>at</strong>ions, two CALPUFF configur<strong>at</strong>ions are run to<br />

identify the performance impact of using a minimum σ v = 0.5 m/s (CALPUFF<br />

default) versus σ v = 0.37 m/s (OCD default).<br />

Evalu<strong>at</strong>ion St<strong>at</strong>istics<br />

<strong>The</strong> evalu<strong>at</strong>ion focuses on the peak observed concentr<strong>at</strong>ions each hour, paired in time<br />

but not in space. As in the OCD applic<strong>at</strong>ions, the wind direction used in the<br />

CALMET input files for the OCD4 d<strong>at</strong>asets (Carpinteria, Cameron, Pismo Beach and<br />

Ventura) is th<strong>at</strong> which is directed from the source to the sampler with the peak<br />

concentr<strong>at</strong>ion. Unlike OCD, this may not align the center of the CALPUFF footprint<br />

with a sampler loc<strong>at</strong>ion due to changes in the transport direction (CALPUFF<br />

trajectories are not straight lines). <strong>The</strong>refore, additional receptors are placed along<br />

the sampling arcs in the OCD4 d<strong>at</strong>asets to better resolve peak modeled<br />

concentr<strong>at</strong>ions within the sampling zone. Spacing for additional receptors along<br />

sampling lines is typically 100 m or less. For the Oresund simul<strong>at</strong>ions, wind<br />

directions are used as measured (adjusted for grid North), and no additional receptors<br />

are used.<br />

Model predictions are quantit<strong>at</strong>ively compared against measured concentr<strong>at</strong>ions and<br />

the resulting performance measures are compared with those for OCD5 (Carpinteria,<br />

Cameron, Pismo Beach and Ventura). <strong>The</strong> original OCD4 input files are recast for<br />

OCD5, and the revised OCD model is rerun on these d<strong>at</strong>a. Using the peak observed<br />

(Co) and peak modeled (Cp) concentr<strong>at</strong>ions each hour, the evalu<strong>at</strong>ion measures are<br />

computed from ln(Co/Cp). This gives equal weight to overpredictions and<br />

underpredictions (e.g., ln(Co/Cp) for Co/Cp = ½ and two are equal in magnitude but<br />

of opposite sign). <strong>The</strong> performance measures are the geometric mean (MG),<br />

Final Report Vol.1 59

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!