1 / 25

Takashi MAKI ,Kazumi KAMIDE and Yukitomo TSUTSUMI Japan Meteorological Agency

OBSERVATIONAL DATA SCREENING TECHNIQUE USING TRANSPORT MODEL AND INVERSE MODEL IN ESTIMATING CO2 FLUX HISTORY. Takashi MAKI ,Kazumi KAMIDE and Yukitomo TSUTSUMI Japan Meteorological Agency TransCom 3 workshop in Paris, 15 th June 2005. Outline.

nodin
Download Presentation

Takashi MAKI ,Kazumi KAMIDE and Yukitomo TSUTSUMI Japan Meteorological Agency

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. OBSERVATIONAL DATA SCREENING TECHNIQUE USING TRANSPORT MODEL AND INVERSE MODEL IN ESTIMATING CO2 FLUX HISTORY Takashi MAKI ,Kazumi KAMIDE and Yukitomo TSUTSUMI Japan Meteorological Agency TransCom 3 workshop in Paris, 15th June 2005

  2. Outline 1. We have developed a new data selection technique which can treat raw (not smoothed) observational data in inter-annual inversion. 2. We have determined proper threshold value to data selection in inversion. 3. We have performed several sensitivity tests to evaluate this technique.

  3. Background However, the representatives of each observational data (not site) is different and the observational network changes year by year. We should use raw (not smoothed) observational data when we try to estimate CO2 flux. We should treat each observational data in accordance with their representatives. Monthly mean CO2 concentration reported to World Data Centre for Greenhouse Gases (WDCGG)

  4. 1. The latest data are not reported Problems in using raw observational data 2. The lack of observational data 3. Affected by local sources?? 4. We have to treat observational data as their representative.

  5. The latest data are not reported, the lack of observational data. How to resolve these problems Creating dummy data(given large uncertainties and reject from analysis in actual) by smoothing, interpolating and extrapolating. We have to treat observational data as their representative. Data uncertainty (representative) are determined by the difference between raw data and smoothed data (as in TC3L2). Affected by local sources?? Each observational data are selected by using inverse model.

  6. Inversion setup *We did not use fossil fuel correction as Baker et al 2005.

  7. Principle 1. The majority of raw observational data has a sub-continental scale representative. 2. CO2 fluxes are determined by the decision by majority of raw data. 3. Therefore, we can select raw observational data using transport model and inter-annual inversion (IAI). Data selection by inversion • Merit of this method • We can treat raw data without developing a new data selection technique. • We could adopt this method to several kind of observational data (surface, upper air and column data). • We do not need huge computational resources (we need only a few iteration of IAI).

  8. Method 1. We use all raw observational data in IAI at first. 2. We reject observational data which the data mismatch is larger than the threshold value. 3. We use selected observational data in inversion with same condition(prior flux and prior flux uncertainty). 4. We repeat process 2 and 3 until we have no rejected data. Data selection procedure

  9. 1. In smaller threshold value, we reject many observational data including back ground data. 2. In larger threshold value, we use many observational data including local representative data. How do we determine the threshold value? We have determined the proper threshold value by checking data selection rate of Global Atmosphere Watch (GAW) global stations (back ground site).

  10. Data selection rate in back ground sites To select 95% raw data in back ground sites, the threshold value (TV) should be larger than 2ppm(without Zugspitze). We choose this TV as the control case.

  11. Data selection rate for all raw data About 82% of the all input data are selected in control case. Recently, the data selection rate tend to smaller?

  12. Data availability* rate (1990-2002) The site which continues observation longer period could contribute well in this analysis. *Availability rate = selected data / analysis period (12*years)

  13. Data selection* rate (1990-2002) Remote sites show better selection rate than continental sites. *Selection rate = selected data / input raw data

  14. Iteration of inversion (TV=2ppm) Fluxes are remarkably changed from 1st to 2nd inversion.

  15. Sensitivity tests We have performed several sensitivity tests 1. Threshold value (TV) sensitivity We have changed TV from 1.0 to 3.0ppm. 2. Observational network sensitivity We have rejected specific site (Crozet, Izana, Mt. Waliguan) which have been shown by K. R. Gurney at last TransCom workshop from observational network.

  16. Sensitivity to the global CO2 flux 1990 1995 2000 Baker et al., 2005 In global scale, each result shows good agreement with each other except restrict selection (threshold value = 1ppm) case.

  17. Sub-continental case (1) In Asian region (relativelyless constrained region), Mt. Waliguan has a large role.

  18. Sub-continental case (2) Relatively many number of observational data region (North America and Europe), loose selection (TV=2.5ppm and 3ppm) tend to enlarge inter-annual CO2 flux variability. Izana has some roles in Temperate America.

  19. Sub-continental case (3) Less constrained region, loose selection tend to enlarge inter-annual CO2 flux variability. Crozet has some roles in South Africa.

  20. Threshold value = 2ppm Averaged flux uncertainty (1990-2002) This work could reduce estimated flux uncertainty especially in formerly less-constrained regions. This selection method could contribute more homogeneous uncertainty in global carbon cycle than T3L2.

  21. We have tried to use raw observational data in inter-anual inversion (IAI). We have developed data selection technique using IAI. In this selection method, we can use 80% of the raw observational data (from WDCGG) in IAI. In global scale, the difference in threshold value and observational network is small. In sub-continental scale, some specific site can effect less constrained region. Therefore, we need more observational data especially less constrained regions. We can reduce estimated uncertainty in using raw data. We can use raw observational data consistently in IAI! Conclusions

  22. We should have good model to data selection • Upgrade model and region spatial and temporal resolution. • From off-line model to on-line model? • We should determine threshold value by statistical analysis • Calculate chi-square (selected data). • 2ppm is too large (comparing with Baker et al) ?? Future plans • We should know the reason which makes the difference between this work and T3L2. • Inversion set-up (prior flux uncertainties)? • Observational data network?

  23. Acknowledgement We thank to TransCom project for providing us the protocol, analysis, inversion codes and many products. We also thank to WDCGG data contributors. (http://gaw.kishou.go.jp:/wdcgg/pdf_contributors.php). References Gurney, K.R., R.M. Law, A.S. Denning, P.J. Rayner, B. Pak, and the TransCom 3 L2 modelers, “Transcom 3 Inversion Intercomparison: Control results for the estimation of seasonal carbon sources and sinks”, Global Biogeochemical Cycles,18, GB1010, doi:10.1029/2003GB002111, 2004. Baker, D. F., R. M. Law, K. R. Gurney, P. Rayner, P. Peylin, A. S. Denning, et al., “TransCom 3 inversion intercomparison: Impact of transport model errors on the interannual variability of regional CO2 fluxes, 1988-2003”, Global Biogeochemical Cycles, 2005, in review. Acknowledgement and References

  24. Appendix 1 : WDC* data policyFrom Minutes of GAW World Data Centre Managers Meeting, March 2005. WDCGG (JMA) proposal - Requirement to acknowledge data contributors with co-authorship is dis-incentive where large scale use of data made in modeling study. Therefore either acknowledge data centre, or have online acknowledgement board. Acknowledgment of data originator is too variable. We forward option 2 of the JMA proposal to AREP and CAS-EPAC together with the recommendation from Albany (WDCPC) for further action. “The WDC managers recommended that the secretariat investigate means of ensuring that Journals require data originators and data centres be properly acknowledged in publications using their data”, with the request for advice and/or a standard acknowledgment/citation for use of WDC data. *WDC: World Data Centre, AREP: Atmospheric Research and Environment Program, CAS: Commission for Atmospheric Sciences, EPAC: Environmental Pollution and Atmospheric Chemistry,PC:Precipitation Chemistry

  25. Appendix 2 : TransCom continuous experiment (JMA experiment) • We have submitted our result at Apr. 2005. • The futures of our model (JMACDTM) • The model calculate vertical mixing coefficient (cumulus convection, turbulent diffusion and shallow convection) making use of our GCM code. • The model does not use PBL height and cloud cover rate explicitly. • The model has artificial parameter (minimum diffusion coefficient) in near surface layers. Thanks to Wouter, Christain, Rachel and Lori for coordinating!!

More Related