170 likes | 383 Views
Intro Page. Modeling Amazonian Carbon Release with Calibrated Soil-Vegetation-Atmosphere Transfer Models. NSF Center for Sustainability of semi-Arid Hydrology and Riparian Areas. Project Goal.
E N D
Intro Page Modeling Amazonian Carbon Release with Calibrated Soil-Vegetation-Atmosphere Transfer Models NSF Center for Sustainability of semi-Arid Hydrology and Riparian Areas
Project Goal • To investigate temporal and spatial variations and model-to-model differences in the calculated carbon exchange of the Amazônian forest ecosystem over the last 40-50 years using models of soil-vegetation-atmosphere interactions which have been calibrated against field data from the LBA field sites using modern multi-parameter estimation techniques.
Project Objectives • Obtain the available data from the LBA field sites relevant to the calibration of SVAT models and carry out a multi-parameter calibration of SiB2 and MOSES-TRIFFID using these data • Explore the variation in optimized parameters obtained by calibrating SiB2 and MOSES-TRIFFID against LBA data, to determine if and how these parameters are related to site-specific seasonal climate, disturbance regimes, underlying soil, and appropriate remotely sensed geophysical variables • Obtain the time series of near-surface forcing variables available from the re-analysis data sets from ECMWF and/or NCEP and validate these time-series against climate records for Amazônia and data from past and ongoing Amazonian field studies (e.g., LBA, ABRACOS, ARME, etc.) • Investigate the temporal and spatial variations and model-to-model differences in the calculated carbon exchange of the Amazônian forest ecosystem over the last 40-50 years by using the time series of meteorological variables [validated in (3)], to force two-dimensional arrays of calibrated SVAT models [specified from (1) and (2)]
Project Progress: Models • The models for which multi-parameter optimization is being (or will be) made: • BATS2: our existing “tried and tested” optimization model, currently being used as our “pathfinder” optimization to explore LBA data set availability/reliability • SiB2: optimization procedures are now largely developed, but refinements are still being made: problems include • unstable iteration methods in the original SiB2 pre-processor package; • some still poorly understood “timing” issues with modeled CO2 fluxes (implicit time of day in code?) • MOSES: we have not yet obtained a reliable source code for this model or begun setting up an optimization package • Simple-SiB: we are making exploratory optimization of this model for possible future use in CPTEC Eta model
Project Progress: LBA Data Sets Currently trying to use data from Beija-flor and trying to analyze for: • Tapajos national forest, Santarem, Para • Km 67 • Km 83 (logged after 1 year) • Reserva Biologica do Cuieiras, Manaus, Amazonas • ZF2 km 34 • ZF2 km 14 (EC data only) (km 14 not examined yet) • Reserva Boiologica Jaru (RBJ), Ji Parana, Rondonia • Floresta Nacional de Caxiuana, near Belem, Para • we know data exists but have not found it in Beja-flor yet (Andreae et al, JGR, 2002)
Data: Tapajos, Santarem, km 67 Currently optimizing with: • Measured EC fluxes of sensible and latent heat, and Net Ecosystem Exchange (NEE) calculated from storage flux and EC CO2 flux • Gaps in model forcing data were filled using km 83 data, where available, or mean data from pervious and next available time periods • Data period shown is that used for optimization • Eleanor’s comments: • Well documented data with good quality control eleanor@hwr.arizona.edu Approx. time EC tower installed
Data: Tapajos, Santarem, km 87 Approx. time EC tower installed Currently optimizing with: • Measured EC fluxes of sensible and latent heat, and Net Ecosystem Exchange (NEE) calculated from storage flux and EC CO2 flux • Forest logged after this period therefore not comparable with other primary forests • Data period shown is that used for optimization • Eleanor’s comments: • Well documented data with good quality control eleanor@hwr.arizona.edu
Data: Tapajos, Santarem, km 67 BATS2 Optimization: all data (including nighttime) included
Data: Tapajos, Santarem, km 67 BATS2 Optimization: no nighttime data included
Data: Tapajos, Santarem, km 83 BATS2 Optimization: all data (including nighttime) included
Data: Tapajos, Santarem, km 83 BATS2 Optimization: no nighttime data included
Data: Cuieiras, Manaus, km 34, 14 Currently working with: • km 34: measured EC fluxes and forcing and storage data accessible 1999-2003 • km 14: measured EC fluxes accessible 1999-2003 (close enough to use same forcing and CO2 storage data?) • Significant forcing data missing in 2001 and 2002 (is this correct?) Interpolated forcing data in red. • Currently optimizing using DOY 165-700
Data: Cuieiras, Manaus km 34 BATS2 Optimization: all data included (small spread in Pareto set?)
Data: Cuieiras, Manaus km 34 BATS2 Optimization: no nighttime data included
Data: Jaru, Rondonia BATS2 Optimization: all data included
Data: Some Comments Eleanor’s (eleanor@hwr.arizona.edu) comments: • Data in Beja-flor from: • Reserva Biologica do Cuieiras, Manaus, Amazonas • Reserva Boiologica Jaru (RBJ), Ji Parana, Rondonia Would be easier to use if: • Documentation specifying data were of same standard as that for Tapajos, Santarem • Link given in Beja-flor delivered data for Cuieiras and Jaru (?) • Time steps were included in data series for periods when no data were taken, but flagged as missing data (for continuity in models) • Net ecosystem were calculated and included (as at Santarem 67 km) • Can’t find (in Beja-flor) the data from: • Floresta Nacional de Caxiuana, near Belem, Para
Data: Some More Comments Jim’s comments (on the basis of discussion this week): • It will not be possible to complete this project in a timely and effective way using data obtained from Beja-flor. • Hence, we will have to go one-on-one with individual groups to obtain data • I believe this will be a common problem for all “synthesis” projects that active over the next 18-24 months • We need to define a list of primary (and secondary) data interface contacts for each flux tower group to service this intermediate need to distribute data rapidly to “synthesis teams” • Recent improvements in the understanding of the origin of flux loss in eddy correlation measurements mean it is possible that all LBA flux data will require reanalysis. • This has the potential to substantially delay progress in this and other “synthesis” projects angle of attack dependent calibration of anemometer and longer time periods for rotation analysis