1 / 44

Cheas 2006 Meeting Marek Uliasz: Estimation of regional fluxes of CO 2 …

Cheas 2006 Meeting Marek Uliasz: Estimation of regional fluxes of CO 2 …. = 40m Sylvania flux tower with high-quality standard gases. = LI-820 sampling from 75m above ground on communication towers. = 447m WLEF tower. LI-820, CMDL in situ and flask measurements.

shiro
Download Presentation

Cheas 2006 Meeting Marek Uliasz: Estimation of regional fluxes of CO 2 …

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Cheas 2006 Meeting Marek Uliasz: Estimation of regional fluxes of CO2 …

  2. = 40m Sylvania flux tower with high-quality standard gases. = LI-820 sampling from 75m above ground on communication towers. = 447m WLEF tower. LI-820, CMDL in situ and flask measurements. • Problems of regional scale CO2 flux estimations by inversions: • limited domain • domain coverage by tower data

  3. ASSUMPTION: SiB-RAMS is capable to realistically reproduce diurnal cycle and spatial distribution of CO2 (assimilation and respiration) fluxes. Therefore, observation data are used to correct those fluxes for errors in atmospheric transport.

  4. time independent corrections to be estimated from concentration data for each inversion cycle CO2 flux respiration & assimilation fluxes simulated by SiB-RAMS

  5. time independent corrections to be estimated from concentration data for each inversion cycle CO2 flux respiration & assimilation fluxes simulated by SiB-RAMS

  6. MODELING FRAMEWORK SiB-RAMS typically run with several nested grids covering a continental scale meteo fields CO2 fields and fluxes LPDM run on any subdomain extracted from SiB-RAMS influence functions CO2 observations inversion techniques Bayesian MLEF corrected CO2 fluxes corrected within each inversion cycle

  7. Representation of atmospheric concentration sample with the aid of influence function, C*, derived from a backward in time run of the LPDM. concentration sample surface fluxes initial concentration inflow fluxes

  8. [10-10 x sm-3 ]

  9. [10-10 x sm-3 ]

  10. [10-10 x sm-3 ]

  11. [10-10 x sm-3 ]

  12. [10-10 x sm-3 ]

  13. [10-10 x sm-3 ]

  14. [10-10 x sm-3 ]

  15. [10-10 x sm-3 ]

  16. [10-10 x sm-3 ] sunrise

  17. [10-10 x sm-3 ]

  18. [10-10 x sm-3 ]

  19. [10-10 x sm-3 ]

  20. [10-10 x sm-3 ]

  21. influence functions for 396m WLEF tower integrated over unit flux for 7x10 day inversion cycles

  22. Implementation for a given inversion cycle C – observed concentration k – index over observations (sampling times and towers) i – index over source grid cell (both respiration & assimilation fluxes) C*R.A – influence function integrated with respiration & assimilation fluxes CIN – background concentration combining effect of the flow across lateral boundaries and initial concentration at the cycle start “beta’s” – corrections to be estimated

  23. INVERSION EXPERIMENTS SiB-RAMS simulation: 75 days starting on April 25th, 2004 on two nested grids (10 km grid spacing on the finer grid)

  24. INVERSION EXPERIMENTS SiB-RAMS simulation: 75 days starting on April 25th, 2004 on two nested grids (10 km grid spacing on the finer grid) LPDM and influence function domain: 600x600km centered at WLEF tower

  25. INVERSION EXPERIMENTS SiB-RAMS simulation: 75 days starting on April 25th, 2004 on two nested grids (10 km grid spacing on the finer grid) LPDM and influence function domain: 600x600km centered at WLEF tower Concentration pseudo-data were generated for WLEF and the ring of towers from SiB-RAMS assimilation and respiration fluxes using correction values of 1

  26. INVERSION EXPERIMENTS SiB-RAMS simulation: 75 days starting on April 25th, 2004 on two nested grids (10 km grid spacing on the finer grid) LPDM and influence function domain: 600x600km centered at WLEF tower Concentration pseudo-data were generated for WLEF and the ring of towers from SiB-RAMS assimilation and respiration fluxes using correction values of 1 Model-data mismatch error was assumed to be higher for lower towers: 1 ppm for towers>100m, 1.5 ppm for towers > 50m, and 3 ppm for towers < 50m and very high values for short towers during nighttime

  27. INVERSION EXPERIMENTS SiB-RAMS simulation: 75 days starting on April 25th, 2004 on two nested grids (10 km grid spacing on the finer grid) LPDM and influence function domain: 600x600km centered at WLEF tower Concentration pseudo-data were generated for WLEF and the ring of towers from SiB-RAMS assimilation and respiration fluxes using correction values of 1 Model-data mismatch error was assumed to be higher for lower towers: 1 ppm for towers>100m, 1.5 ppm for towers > 50m, and 3 ppm for towers < 50m and very high values for short towers during nighttime 7 x 10 day inversion cycles were performed using Bayesian inversion technique with concentration pseudo data (initial corrections = 0.75 and their standard deviations = 0.1)

  28. source area: 20x20 km NW of WLEF 10 day (cycle) average

  29. source area: 20x20 km NW of WLEF 24 hour average

  30. source area: 20x20 km NW of WLEF hourly average

  31. NEE uncertainty reduction [umol/m2/s] cycle #1

  32. NEE uncertainty reduction [umol/m2/s] cycle #2

  33. NEE uncertainty reduction [umol/m2/s] cycle #3

  34. NEE uncertainty reduction [umol/m2/s] cycle #4

  35. NEE uncertainty reduction [umol/m2/s] cycle #5

  36. NEE uncertainty reduction [umol/m2/s] cycle #6

  37. NEE uncertainty reduction [umol/m2/s] cycle #7

  38. NEE UNCERTAINTY: INITIAL, WLEF, RING aggregation of source areas

  39. NEE UNCERTAINTY: INITIAL, WLEF, RING aggregation of source areas

  40. NEE UNCERTAINTY: INITIAL, WLEF, RING aggregation of source areas

  41. Implementation for a given inversion cycle C – observed concentration k – index over observations (sampling times and towers) i – index over source grid cell (both respiration & assimilation fluxes) C*R.A – influence function integrated with respiration & assimilation fluxes CIN – background concentration combining effect of the flow across lateral boundaries and initial concentration at the cycle start “beta’s” – corrections to be estimated

  42. Implementation for a given inversion cycle C – observed concentration k – index over observations (sampling times and towers) i – index over source grid cell (both respiration & assimilation fluxes) l - index over time intervals C*R.A – influence function integrated with respiration & assimilation fluxes CIN – background concentration combining effect of the flow across lateral boundaries and initial concentration at the cycle start “beta’s” – corrections to be estimated

  43. Inversion experiments: • Pseudo-data experiments • The ring of towers (Bayesian, MLEF) • US continental scale (MLEF) • Real data experiments • The ring of towers (Bayesian) new SiB-RAMS simulations

  44. RUC-LPDM: • Influence functions to be integrated with user provided CO2 fluxes

More Related