Prelim. outline for Data Analysis

Monthly scale – refers to 5 weeks of data plotted 1-2-3-4…as days

Plot: CO2 continuous data from picarro, 5-min running mean

Analyses:

-       regression line

o   8 oct à 27 oct, fitted to 2 weeks of data from window

o   28 oct à 14 nov, fitted to the other 3 weeks of data (from the roof)

o   slope from regression line is the ‘quantifiable’ increase in [CO2]/time

-       periodicity – interpolate onto graph to see how data of individual days relates to ‘average’ periodicity of the dataset

o   are our hypothesized timescales of interest reflected in the periodicity of the data?

-       Fourier Transform – see how noisy the data is, look at frequency domain

o   Substantial differences between the roof and window data?

o   Is data useful or just a mess (extremely noisy)?

-       regression and correlation for each factor to [CO2] and see which is dominant (has the greatest ‘relationship’ to the pattern of CO2

o   wind speed and (direction), ppn, T

§  do you have to transform CO2 to the same timescale to analyze or do you interpolate the factor data

-       do Principle Component Analysis?

-       regression between δ13C month-long plot & [CO2] (negative)

-       Mass balance and keeling plot

Weekly scale (M, T, W, Th, F, Sa, Su stacked on top of each other – average)

x2, for 2 weeks of window measurements and 3 weeks of roof measurements

                                                                                               

Analyses:

Just for traffic and CO2

-       must test which days of the calendar week best represent our hypothesis for ‘weekday’ traffic pattern being different from ‘weekend’ traffic pattern

o   i.e. how do we define a weekday and a weekend? (aside from solely based on the calendar)

1^st test:

plot Mà F vs Sa, Su

                  2^nd test:

plot T à Th vs F, Sa, Su, M

  

From these tests:

- Is there a significant different between the traffic and CO2 patterns on weekdays (defined by tests 1 & 2) and weekends?

Analyzing all factors:

Now, take out a ‘test week’ from roof and window datasets (picarro) and look at these two datasets in relation to averaged out weekly CO2

-       again correlation, regression etc

o   do the values change as compared to those calculated for monthly?

o   If so, can we infer something about the degree of relationship (between each factor and CO2, or for traffic and biogenic - δ13C) changing over time or according to timescale being analysed

§  Is our method of looking at these relationships skewing the results or do they change over time

-       for CO2, look at variance/std deviation etc in relation to averaged out week to see if there is some sort of representativeness in the average values

-       pick weeks with 2+ days of traffic data from counting cars (‘real-time’ instead of averaged out google data)

Diurnal – Averaged out so that it represents more than a single day i.e. is representative of all days stacked into 1 24h period

-       there are many options for how we construct this 24 hour period

o   Based on weekly analysis of the CO2 pattern

1)   2 plots, one for a weekday (defined as Mà F or T à Th) and one for a weekend (Sa, Su or F, Sa, Su, M)

o   x2 for the 2 weeks of window data and the 3 weeks of roof data

2)   1 plot using data from weekends and weekdays

o   x2 for window and roof

3)   1 plot using data from weekends and weekdays for all 5 weeks

Analyses: will compare this ‘ideal’ or averaged out day to actual days where the influence of one factor is predominant (e.g. Oct 22, very heavy winds)

-       see if some factors have  a relatively constant value for correlation despite not being ‘dominant’ on that day

-       will want to look at amplitude of peaks, [CO2] values, variation (fluctuations), etc

o   basically, characterize or describe the CO2 pattern, δ13C as an ideal/representative pattern, and then see how it changes on a case-by-case basis

Case-by-case analysis – Has been  expanded and will be explained in class and then added.

Update on CO2 data

Here are some graphs of the recent Picarro data, which Gregor gave me last week. The fist graph shows part of the data we looked at on last Tuesday (this time only form 17th to 26th). The second graph is the most recent CO2 data that i have. It ranges from October 26th to November 4th.

I thought it is a good idea to show this now, such that everyone can get an idea what the recent data looks like.

There are a couple of gaps in the readings which corresponds to the time when flask measurements were taken. Except the large one between Tuesday 2^nd and Wednesday 3^rd, which is due to the power shut down.

The black crosses correspond to the raw 12-CO2 data. We can see that the CO2 signal becomes much more spread starting at mid -Thursday (28^th of October). I’m guessing that this corresponds to the event when the tube was set up on the roof.

I also attached a 3rd graph showing the standard deviation of the 2nd CO2 time series (the one between October 26th and November 4th). You can see that the error increases significantly after the 28th, and then even more after the 29th of October.

The red line is a running mean with a window of 100 time units, corresponding to about 15 minutes.  The running mean shows a single curve with a detectable trend, but after the 28th of October it has a much larger deviation associated with it.

In addition, it seems that the "spread of the Co2 signal" decreases during night time, but is prevalent through the day until late night.

That is what I can say now. Let me know if you have some questions and feel free to comment on this.

Even more things to look at

Comparison between our first Panorama (18th of Oct) and the one on Nov.1. 

-g

Something to look at!

Hi, This is the graph showing the assembled 12-CO2 Picarro measurements from the 9th to the 26th of October.

Enjoy!

- Jan-Erik

Transportation Study in Montreal

This is from a report (following a transportation conference) by Chapleau and Morency (from the Civil Engineering Department (Transportation) at Ecole Polytechnique de Montreal): "Generic in nature but specific to the methodological procedures undertaken in the Greater Montreal Area (GMA), the CATI (Computer-Assisted Telephone Interview) household survey is conducted about every five years over a 5% sample. Typically, it represents about 160,000 people belonging to 65,000 households declaring some 400,000 individual trip records for an average weekday. Individual trips are geo-referenced for the residence, trip origin and destination, modal junction points (kiss-and-ride and park-and-ride locations), and are described for their household and personal characteristics (age, gender, car license, car ownership, income) in addition to the trip attributes (purpose, mode, departure time, train-subway-bus routes taken if traveling by transit, bridges and highway taken if traveling by car)."

They have a lot of GIS analysis done on transportation ranging from buses, cars to subway trains. A nice map they made was the percentage of motorized trips in Montreal am rush hour:

We could contact them for the data if we find it is useful!

-Angela.

Mount-Royal Flasks

This morning at 7am (it was still dark outside), Pogo and myself took our first flask measurment! I placed a stick with a Yellow Flag in the ground and I tied an Orange tape across a branch in the air at the spot where i took the measurments. It is directly facing you when you reach the top of the stairs.

-Atleast it is warm today, Angela.

Principles of Urban Meteorology (re atmospheric CO2)

(by Andrew)

The generalization of fundamental meteorological principles over urban environments is more problematic than to other land surfaces. This is for three reasons. First, the heterogeneity of roughness elements within the complex 3-D geometry of the urban “canopy” creates spatial variability in turbulence patterns. Second, the multiple sources and sinks of momentum, heat, moisture, and emissions creates spatial variability in fluxes and concentrations. Third, the impact of human activities to continually re-shape and alter, in new and distinct ways, the urban environment itself limits the long-term validity of any findings. In short, there is considerable uncertainty concerning the dynamics of the unique microclimates found within previously understudied and increasingly complex urban areas.

Both across cities and between cities, significant spatial and temporal variability in CO2 concentrations and fluxes can be expected as a consequence of the distribution of anthropogenic sources (mobile and fixed), processes of urban vegetation (including irrigation, and patterns of atmospheric convection and advection. Previous studies have shown that there is a marked and distinct diurnal cycle in the concentration of CO₂ with a morning peak attributable to anthropogenic (largely traffic), biogenic (nocturnal respiration), and meteorological (atmospheric stability) factors. In contrast, a mid-afternoon minimum can be attributed to vegetative photosynthesis and strong convective turbulence; concentrations then begin to rise again during the evening “rush-hour” traffic-flow.

....

The main roughness elements in urban environments are trees and buildings. Other than being large, trees and buildings share few characteristics of meteorological significance. Aerodynamically, buildings are true “bluff” bodies because of their impermeability, inflexibility, and sharp edges. When exposed to airflow they create strong positive and negative pressure differences over their surface, leading to flow separation and vortex shedding. Trees are also good generators of mechanical turbulence but buildings have to be judged as more effective roughness elements. The effects of smaller roughness elements, such as cars or paved surfaces, are minimum in comparison.

Three spatial scales are commonly utilized for studying urban environments:

-       the micro-scale (10¹–10²m) involves spatial differences in response to individual roughness elements (variability in building/canyon dimensions, trees) and proximity to localized emissions sources (e.g. roads, vegetation);

-       the local-scale (10²–10⁴ m) represents the integrated response of an array of roughness elements with spatial variability reflecting the unique characteristics of different neighborhoods/land-uses;

-       and the meso-scale (10⁴–10⁵ m) considers the city in its entirely, and differentiated from its surroundings, areas of forest, agriculture, etc.

The urban canopy layer (UCL) is defined as being from the ground to the mean height of the roughness elements, usually just below roof –level, where micro-scale effects of the site characteristics are dominated. The UCL is most clearly delineated in areas of high building density; it may be discontinuous or absent in less densely developed suburban areas.

The layer extending from the top of the UCL, to a height where urban surface influences are no longer perceptible, is defined the urban boundary layer (UBL). It includes the roughness sub-layer immediately affected by the individual roughness elements, the turbulent surface layer (local-scale), and the outer mixed layer (meso-scale).

....

Of those studies employing atmospheric-based measurement methods to study CO2 concentrations in urban environments to date, virtually all, with a few exceptions, have focused on the micro-scale, considering processes and patterns within the UCL. Inadequate attention has as yet focused on how micro-scale results can be extrapolated to larger scales and on how to accurately study the local-scale using atmospheric-based measurement methods. 

In regards to the latter, current debates focus on determining the height (or “depth”) of the roughness sub-layer, in which the perturbations caused by individual roughness elements are “blended” together due to atmospheric turbulence. It is as this height that instruments are to be placed in order to study at the local-scale that is spatially representative of a distinct urban neighbourhood/land-use. To be sure, placing instruments are greater heights than this leads to increased risk of incurring errors due to advection from dissimilar upwind surfaces and storage changes below the measurement level due to vertical flux divergence.

It is known that the height of the roughness sub-layer is a function of both the length/height of roughness elements (z_H) and their horizontal spacing. More recent research suggests that the latter factor may in fact be the primary determinant. It has been estimated that, as a general “rule-of-thumb”, instruments must be mounted at a height at least twice the mean height of the roughness elements (approximately 20-90m) to ensure that they are above the influence of individual roughness elements and, therefore, that the measurements represent an integrated response at the local-scale.

....

References

Grimmond, C.S.B., et al. (2006) Progress in measuring and observing the urban atmosphere.

Theoretical and Applied Climatology 84, 3-22.

Grimmond, C.S.B., et al. (2002) Local-scale fluxes of carbon dioxide in urban environments:

methodological challenges and results from Chicago. Environmental Pollution 116,

243-254.

Kanada, Manabu. (2007) Progress in Urban Meteorology: A Review. Journal of the

Meteorological Society of Japan 85B, 363-383.

Koerner, B. and J. Klopatek. (2002) Anthropogenic and natural CO2 emission sources in an

arid urban environment. Environmental Pollution 116, 45-51.

Nemitz, E., K. J. Hargreaves, A. G. McDonald, J. R. Dorsey, and D. Fowler (2002) Micrometerological Measurements of the Urban Heat Budget and CO₂ Emissions on a City Scale. Environ. Sci. Technology 36, 3139-3146.

Oke, T.R., et al. (1988) The urban energy balance. Progress in Physical Geography 12, 471-

483.

Oke, T.R., et al. (1989) The Micrometeorology of the Urban Forest. Philosophical

Transactions of the Royal Society of London 324, 335-349.

Wentz, Elizabeth A., et al. (2002) Spatial Patterns and Determinants of Winter Atmospheric

Carbon Dioxide Concentrations in an Urban Environment. Annals of the Association of American Geographers 99(1), 15-28.

Random Picture and Maps

I posted some pictures and maps I have that may be usefull for our website!
- Angela.

Montreal Land Use Map: map showing land use in montreal from the City of Montreal website.

Drivers to Work: Here is a map showing the percentage of drivers for each bourough from The City of Montreal website (Source is Stats Canada though).

Studay Area: This is a Satelite image from Google maps.

Site Area: Map from google maps.

3D Site Area: From Google Earth.

Mc Tavish and Burnside Wind Speed and Direction (Oct 21-22)

I looked at wind speeds and directions from thursday 21st to friday 22nd. The attached graphs show you what I have found out. First of all, I realised that the Burnside data (for whatever reason) is shifted by 4 hours. So the wind data are recorded with a time tb = tt + 4. To correct for this difference I simply subtracted 4 hours in the time columns for burnside. That is why the burnside data actually starts at t = -4. Note that the datasets end at 8pm local time on the 22nd October (that means at t = 44h, 4 hours before midnight).
If we correct for this shift the two data records correlate quite well. However, this is only based on qualitative estimation. As you can see from the Windrose plots, during relatively large time intervals Burnside and McTavish have more or less the same wind direction.
However, I tried some regression analysis and the correlation coefficients for the wind direction is always very low. That mean we don't have a quantitative proof yet to confirm our hypothesis. Maybe we can ask one of the profs to find out how we can use statistics to proof our finding quantitatively.

- Jan-erik

Weekend CO2 Events

Over the weekend I witnessed a few potential sources of CO2. A building was on fire Friday afternoon and Saturday day and night University street was closed for construction on the sewers.

Friday, on René-Lévesque Blvd. near the intersection with de Lorimier Ave. at 2 p.m. fire crews were on the scene of a two-alarm blaze.

Saturday, morning till night, university street right next to our machine was emmitting this huge cloud of white smoke. I asked the workers what they were doing and they said repairing the sewers. (See pictures posted below, taken at 10pm Saturday night).

Next week I think it would be cool to see if we can detect these events with the use of piccaro and wind speed and directions!

- Angela

Old Weekly Data

Here the data from last week. -g

Point Sites on Mt. Royal

I've picked these sites based on accessibility, their elevation and how close they are to a road. All three of them are in a forested area, roughy the same elevation (Location C is higher, but also the most remote).

They can be accessed at any day and usually have people around them so they would be safe at night as well. Check out the Google Map for their exact location.

Location C

Location B

Location A

Gregor

Panorama

Here is the first Panorama. Monday 18-10-2010 at 11.30am. (click for high res)

-Gregor

http://www.quebec511.gouv.qc.ca/en/cameras/montreal/

http://www.radar.mcgill.ca/data-archive.html

sounding data

-rachelle

some few things

Teaaammm!!

Some few things...

(1) We've decided that we should all to try to meet twice a week. Once on Tuesday (during class), and then once on Friday. I propose that we meet in the ESYS 500 lounge (which now has a coffee maker + kettle and many chairs) or the GIC.. I've created a doddle poll for everyone to fill in their availabilities. Please fill it out ASAP.. http://doodle.com/participation.html?pollId=qzirrw45scsvxemu

(2) Once the work plan is finished (today by 5) and emailed around, it might be good idea for everyone to quickly read through it -- esp. the hypotheses section and the workplan section with the weekly timetable.... The next two weeks should be pretty hectic as we set-up data-collection procedures (which Gregor is coordinating) and begin data-analysis, and so we weren't planning to begin developing website content until these two weeks are done.

(3) Adam and I are going to be doing some more final, 'background' research this week on urban micro-meteorology and biogenic respiration.. Please let us know ASAP if there are any other specific areas that you think require more research so that we're prepared to interpret our data.

Thanks all,



Andrew

Thanksgiving Data

[Gregor]

Hey Team,

here is the outdoor data from this weekend. Ignore the bottom graph (only for reference).

- 12C
- Delta (5 min avg.)
- H2O
- Day

Burnside winds vs. McTavish winds

(Angela and Jan-Erik)

We have compared the wind directions and speeds at the McTavish weather station and the Burnside weather station to see if they are similar (if they were we could also assume FDA will have similar wind patterns).

Unfortunatly the 2 weather stations do not agree in their wind speeds (in the past hour McT is 6km/h while Burn is 16km/h) or their wind directions (in the past 24 hours the winds at Burnside were mainly SW and N while McT were mainly S and SW (see attached photos)).

So the questions I have in mind are:
1. Should we use wind data?
2. If yes, then from where? McT or Burn?
3. Should/Can we get our own anamomter?
4. Can it be installed on the roof and left unatended?

Traffic Data

(Gregor) Hey Team, a quick update on the traffic data. I carried out a quick experiment on my way home and it looks like the google traffic data is not full-proof. These photos all show the same location. Sherbrooke road, the blue dot is me and the photo is of that area looking SW !!! not NE.

Here is a link as to how they measure traffic in an urban area (crowdsourcing): http://googleblog.blogspot.com/2009/08/bright-side-of-sitting-in-traffic.html

workplan!!!!

Teammmmmm!!!!!!

Here is a preliminary outline for the 10-page workplan that we have to have done by next Tuesday.

Recall from syllabus what the workplan is supposed to include…

-       Statement of, and brief introduction to, the research question

-       Review of literature outlining critical findings and current gaps in knowledge

-       Proposed work including a description of how the research questions will be answered

-       Presentation of research methods and work plan outlining individuals responsibilities

-       Proposed timeline of research progress

Please have your final sections to the two of us by the end of Sunday at the latest so that we have time to make edits re coherence. (For some of you, this just means editing sections you already sent us to make them more concise.)

Please use the following format: 12 pt. font, Times New Roman, single-spaced, 2.5 margins, MLA citation style.  

INTRODUCTION – Urban Metabolism?

-       1 page

-       by Andrew and Rachelle

STUDY – What and Why?

-       2 pages

-       by Andrew and Rachelle

METHODS

-       Instrument / set-up

o   1 page

o    by Nadine and Gregor

-       Determining background CO2 / differentiating between sources

o   1 page

o   by Ariella and Adam

o   (summary of lit. review concerning methods for determining background, specific challenges with doing so for this project, and future research that will be done (e.g. Keeling Curve) -- perhaps include Appendix with more technical details)

-       Weather data and traffic data (including wind-rose and eddy mixing with height)

o   ¾ of a page

o   by Angela and Jan-Erik

o   (I am guessing you’ll have to do some research about what traffic data we can obtain.)

-       Data-Analysis (using matlab)

o   ¼ page

o   by Nadine/Gregor + Jan-Erik

o   (really brief few sentences about how we are going to convert data into MATLAB format so as to apply statistics)

HYPOTHESES –

-       2 pages

-       As with data interpretation, it would likely make sense if we all contributed to developing our research hypotheses… We’ve created a shared word document with the paragraphs that Adam and Arielle already wrote concerning temporal variation… Everyone should try put in their ideas by the end and Sunday, and then the two of us will finalize the section…

-       https://docs.google.com/document/edit?id=1p8ds2hxc_Q9zLJC0xeAMi_9vSAz1wQiIaPdbi1ZfKEk&hl=en&authkey=CI_Mk9sB

WORKPLAN

-       2 pages

-       by Andrew and Rachelle

-       outline of workload with timetable of progress

+ Appendices (!)

-       calendar with time-table

-       spreadsheet that divides up responsibilities among all of us

-       ….

We have also created a spreadsheet with all the tasks and responsibilities that we have to divide up. Please fill your name into things you’d be interested in doing, and also add things that you think we missed.

https://spreadsheets.google.com/ccc?key=0AkXrxolOozuddEMtVU5tbUI2S2Z5YVNWbW1TVmhBNVE&hl=en&authkey=CImI7dwE

Thanks all…

Andrew and Rachelle