![Page 1: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/1.jpg)
Virtual Tall Towers and Inversionsor
How to Make Productive Useof Continental CO2 Measurements
in Global Inversions
Martha ButlerThe Pennsylvania State University
Department of MeteorologyApril 24, 2007
TransCom Meeting
![Page 2: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/2.jpg)
Outline
• The Problem with Continental CO2 data
• Virtual Tall Tower Overview
• Inversion Setup
• Perfect Data Results
• Continuing Work
![Page 3: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/3.jpg)
Problem Definition
• Continental data are messy: big signal with lots of noise.
• Confounded by boundary layer processes (diurnal, synoptic, seasonal)
• Transport models have problems representing the continental boundary layer, especially near the land surface.
• Most continental observations are made in this surface layer.
• Observation site density does not permit global estimation of fluxes at the resolution we would like, unless we use continental data.
![Page 4: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/4.jpg)
What About Flux Towers?
• CO2 is measured continuously at every continental flux tower, but…
• These data are not required to be either accurate or precise to support CO2 flux calculations.
• What if?• Flux tower CO2 data could be high precision and calibrated
with the global standards (at reasonable cost).• The data could be “adjusted” (using meteorological and flux
data available at the site) to represent a mid-day CO2 mixing ratio value above the surface layer.
• This adjustment adds to the uncertainty of the observation; is it worth it?
![Page 5: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/5.jpg)
What Is This Adjustment?
Following the mixed layer similarity theory of Wyngaard & Brost [1984] and Moeng & Wyngaard [1989], the vertical gradient of a scalar in the boundary layer:
0
* *
izb t
i i i i
wcwcC z zg gz z w z z w z
where
gb and gt are bottom-up and top-down gradient functions scaled by boundary layer depth zi
w* is the convective velocity scale
wc0 and wczi are the surface and entrainment fluxes of the scalar C
![Page 6: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/6.jpg)
Example Adjusted Time Series
Hourly (mid-day hours) adjusted CO2 mixing ratios for:
Harvard Forest
Northern Old Black Spruce
WLEF
Courtesy of David Tyndall
See also www.amerifluxco2.psu.edu for high-precision, calibrated measurements at flux towers and www.eol.ucar.edu/~stephens/RACCOON/ for the Rocky Mountains sites.
![Page 7: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/7.jpg)
Inversion Setup
• Forward transport model: PCTM, 2.5° longitude x 2° latitude using GEOS-4 met data (3-hr) for 2000, 2001, 2002
• Background fluxes:– Seasonal fossil with 1995 spatial pattern (ORNL)– SiB3 hourly flux for 2000, 2001, 2002 (CSU)– Takahashi 2002 ocean
• Matrix inversion in the TransCom3 style solving for monthly fluxes in 36 land regions and 11 ocean regions in 2002
• Forward model sampling strategy: – Entire atmospheric column sampled hourly in 2002 at each observation
location.– Ten-day means of the atmospheric grid saved for analysis.
![Page 8: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/8.jpg)
Networks Tested (so far)
• 85 station network: NOAA GMD sites active in 2002• 159 station network: NOAA GMD + GAW• 145 station network:
– Remove some profiles and added VTT sites active in 2002• 165 station network:
– Add high-precision, well-calibrated sites active after 2002• Flux tower/VTT locations• Britt Stephen’s mountain sites• Bev Law’s Oregon transect• 2 established sites in Africa
• 174 station network:– Add a few other sites in under-represented regions
![Page 9: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/9.jpg)
![Page 10: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/10.jpg)
85 Station Network
![Page 11: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/11.jpg)
159 Station Network
![Page 12: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/12.jpg)
145 Station Network
![Page 13: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/13.jpg)
165 Station Network
![Page 14: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/14.jpg)
174 Station Network
![Page 15: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/15.jpg)
Tests of the Inversion Setup
• Sensitivity to Choice of Priors & Prior Uncertainty– Priors derived from TransCom3 Seasonal Mean Inversion– At least 0 ± 2 Gt C/yr for each region– 0 ± 5 Gt C/yr for each region– 0 ± 10 Gt C/yr for each region– At least 0 ± 3 Gt C/yr for each land region, TransCom for ocean regions
• Sensitivity of Posterior Flux Uncertainty to Choice of Networks– Boreal and Temperate North America– South America– Africa
![Page 16: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/16.jpg)
Global Aggregated Results
Thick colored lines show posterior flux (and the prior contained in the background flux).
Dashed lines show the posterior uncertainty for the aggregated area for the 85 station network.
Error bars show the posterior uncertainty for the 174 station network.
Note that flux scales differ for land and ocean.
![Page 17: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/17.jpg)
June Posterior Flux Uncertainties
Oregon transect Flux Towers
Flux Towers
![Page 18: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/18.jpg)
Boreal North America and Tropical America
Fraserdale Flux Towers
Tapajos
![Page 19: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/19.jpg)
Examples for the African Regions
Hombori
Skukuza
![Page 20: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/20.jpg)
Continuing Work
• Noisy data tests with simulated observations• Model sampling to match observation times• Preparation of observation data• “The real deal”• Daytime forward pulse runs from land regions
![Page 21: Martha Butler The Pennsylvania State University Department of Meteorology April 24, 2007](https://reader033.vdocument.in/reader033/viewer/2022051422/56814d04550346895dba322b/html5/thumbnails/21.jpg)
Acknowledgements
• People– Davis Group at Penn State– Denning Group at CSU– Randy Kawa at NASA GSFC– Everyone, everywhere in the CO2 measurement network
• Funding Agencies– NASA GSRP– NOAA, DOE, NSF– PSU, College of EMS Centennial Travel Grant
Comments and suggestions are welcome!
And thank you!