[ Home  Homework  ]
Homework Assignments
Caution: Text in Red indicates assignments from 2 years ago. Please ignore them.
Due by 
Details 

. 
14 Sep 2021 
a) Read Chapter 19 Pollutant Dispersion from Stull, 2017: Practical Meteorology. b) Read the BC Primer on Air Quality Modeling. https://www2.gov.bc.ca/assets/gov/environment/airlandwater/air/reportspub/primerbcaqdispersionmodellingguideline.pdf 
Be prepared to discuss Chapter 19 of Stull in class. 
. 
16 Sep 
a) Read Chapter 18 Atmospheric Boundary Layers from Stull, 2017: Practical Meteorology. b) Skim the BC Modeling Guideline. https://www2.gov.bc.ca/assets/gov/environment/airlandwater/air/reportspub/bcdispersionmodellingguideline2015.pdf 
Be prepared to discuss Chapter 18 of Stull in class. 
21 Sep  a) Read Chapter 5 Atmospheric Stability from Stull, 2017: Practical Meteorology. b) Reinforce your understanding by reading p1 37 from the BC MoE 2017 talk. 
Be prepared to determine static stability from an atmospheric sounding, in class. We will also learn more about Gaussian plumes in class.  
. 
23 Sep 
# Gaussian Plume  simple # R. Stull, 8 Feb 2016, Modified 4 Jul 2018. # Givens ============== # meteorology m = 20. # wind speed (m/s) zi = 2000. # mixed layer depth (m) sigma_v = 1.3 # lateral velocity variance (m/s) sigma_w = 1.02 # vertical velocity variance (m/s) tl = 60. # Lagrangian time scale (s) # plume zcl = 150. # plume centerline height (m) q = 300. # emission rate (g/s) # domain, assuming origin (x, y, z) = (0,0,0) is a base of smoke stack xmax = 10000. # x downwind domain size for x = 0 to xmax (m) ymax = 500. # y crosswind domain size for y = ymax to +ymax (m) zmax = 500. # z vertical domain size for z = 0 to zmax (m) # spatial resolution for calculations of concentration delx = 200. # x increment (m) dely = 20. # y increment (m) delz = 20. # z increment(m) # Assignment ============== # Create contour plots of concentration (µg/m^3) # using Stull eq.(19.20) with eqs. (19.13) within the domain for: # a) horizontal (x,y) slice at earth's surface (z = 0) # b) horizontal (x,y) slice at height of plume centerline (z = zcl). # c) vertical alongwind (x,z) slice of thru plume centerline (y = 0) # d) vertical crosswind (y,z) slice at a downwind distance x # from the stack of xslice = 3000. # xlocation of crosswind (y,z) slice (m) # Hints: # Use any computer language; e.g., R, MatLab, python, fortran, excel. (I used R.) # Use a contour interval of 20 µg/m^3. # Check that your answer to part (a) is similar to the sample # application (Solved Example) in Stull p734 before you do the # other parts. (It won't be exactly equal, because of different # inputs in this HW assignment.) # Turn in your contour plots AND your code. # Please have your name on everything you turn in. # 
This is a
chance for you to learn about the simplest dispersion model   the
Gaussian plume model. It is useful for short range (< 50 km)
nearsource dispersion calculations for statically neutral or
statically stable conditions. You can write your code in any language; e.g., R, MatLab, python, fortran, excel, etc. (Note, although the assignment description at left is written using syntax for the R language, you may use any computer language you want. In R, the # symbol denotes the start of a comment.) 
. 
28 Sep inclass 
Install AERMOD on your computer, and run the test case. See details at https://www.eoas.ubc.ca/courses/atsc507/ADM/aermod/index.html and Install AERMET on your computer, and run the test case. See details at https://www.eoas.ubc.ca/courses/atsc507/ADM/aermod/index.html . 
We will do this in class on 28 Sep, so you might not need to do anything else for this HW. 
. 
5 Oct 
Install AERMAP on your computer, and run test case. See details at https://www.eoas.ubc.ca/courses/atsc507/ADM/aermod/index.html . Present a section of the AERMOD Model Formulation material to the class. a) Look over the EPA AERMOD home page Use the Slack chat space to communicate with me and other students to decide who wants to present which section, so that there are no duplicate presentations. One topic per student. Your topic choices are:

You don't need to prepare any powerpoint slides. Just be prepared to show the appropriate pages from the AERMOD Model Formulation manual as you explain your section to the class. If you don't understand every part of the section you present, that is OK   we can discuss it together during class. 
7 Oct  Finish student presentations (see above). Everyone, please go to the Hysplit page and register. It is free, but there could be a week delay between when you register and when you have access to the hysplit source code. So by starting now, you will be ready for the hysplit lessons in a couple weeks. https://www.ready.noaa.gov/HYSPLIT_register.php Also, Stull will lecture on:


6 
12 Oct 
Stull continues lecture on Lagrangian particle dispersion. Turn in your results and code for the following exercise: # Deardorff Plume Dispersion in Unstable PBL 
This
is a chance for you to learn about dispersion in a convective
(statically unstable) boundary layer   also known as a mixed
layer. It uses the discoveries by Deardoff and Willis, as
parameterized by me in my Practical Meteorology book. Read and follow the instructions in the order given   it will gently lead you towards a better understanding of dispersion in a convective pbl. You can write your code in any language; e.g., R, MatLab, python, fortran, excel, etc. 
14 Oct  Continue working on the AERMOD run, due by the next class meeting (see below). Note, you can work as a team and submit one report from the whole team.  
. 
19 Oct 
Due today is your run of AERMOD (and AIRMET and AIRMAP if needed) for a somewhat actual case: Battle River Generating Station Alberta, Canada about 150 km southeast of Edmonton, AB. Stack C (Unit 5) height of stacktop above ground: 161 m diameter of stack: 5.8 m fuel: coal SO2 emission rate: 22,961 tonnes/year Run your AERMOD for one year only. (You can pick the year, based on data availability.) Upper air data: Edmonton Stony Plain. Need to use only the morning (12 UTC) soundings. Select "Mandatory AND Significant Levels". Station identifier: WSE
Surface weather data. Nearby weather stations are: Coronation, or Red Deer, or Edmonton, AB. Note: if you are unable to find a surface data archive in the correct format for a Canadian station near Battle River, then for the purpose of this learningexercise, use Havre CityCounty Airport (KHVR) in Havre, Montana, USA. It is roughly the same distance east of the Rocky Mtns ad Battle River. Lat: 48.54278°NLon: 109.76333°WElev: 2589ft. , and pretend it is at Battle River. WMO Id: 72777. For roughness length, albedo, and Bowen ratio, use Google Maps street view to see the local surface near the power plant. (Hint you can see the smokestacks in the distance from the road, if you look in the right direction.) In your submitted report, justify your choice of these 3 variables. For tips on doing this, see the AERMOD implementation guide. https://www3.epa.gov/ttn/scram/models/aermod/aermod_implementation_guide.pdf But for this learningexercise, pretend the surrounding terrain is flat, and that all receptors are at the same elevation as the base of the smokestack. Put receptors on a cartesian grid, with 2 km grid spacing extending 50 km from the smokestack in each of the main compass directions (north, east, south, west). In addition to turning in the AERMOD output, please create a 1page summary/overview report, and indicate which receptor locations (if any) exceed the Canadian airquality standards for SO2. (Please see me if you need additional info or have questions for this exercise.) 
Bring your computers to class today, as Tim Chui will show us how to install hysplit. 
Other Info
[ Home  Homework  ]
http://www.eos.ubc.ca/courses/atsc507/ADM/
Copyright © 2016, 2018, 2021 by Roland Stull.