Abstract
We present a prototype of a Carbon Cycle Data Assimilation System (CCDAS), which is composed of a terrestrial biosphere model (BETHY) coupled to an atmospheric transport model (TM2), corresponding derivative codes and a derivative-based optimisation routine. In calibration mode, we use first and second derivatives to estimate model parameters and their uncertainties from atmospheric observations and their uncertainties. In prognostic mode, we use first derivatives to map model parameters and their uncertainties onto prognostic quantities and their uncertainties. For the initial version of BETHY the corresponding derivative codes have been generated automatically by FastOpt’s automatic differentiation (AD) tool Transformation of Algorithms in Fortran (TAF). From this point on, BETHY has been developed further within CCDAS, allowing immediate update of the derivative code by TAF. This yields, at each development step, both sensitivity information and systematic comparison with observational data meaning that CCDAS is supporting model development. The data assimilation activities, in turn, benefit from using the current model version. We describe generation and performance of the various derivative codes in CCDAS, i.e. reverse scalar (adjoint), forward over reverse (Hessian) as well as forward and reverse Jacobian plus detection of the Jacobian’s sparsity.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Giering, R., Kaminski, T.: Recipes for Adjoint Code Construction. ACM Trans. Math. Software 24 (1998) 437–474
Giering, R., Kaminski, T., Slawig, T.: Applying TAF to a Navier-Stokes solver that simulates an Euler flow around an airfoil. To appear in Future Generation Computer Systems (2003)
Stammer, D., Wunsch, C., Giering, R., Eckert, C., Heimbach, P., Marotzke, J., Adcroft, A., Hill, C.N., Marshall, J.: The global ocean circulation during 1992–1997, estimated from ocean observations and a general circulation model. J. Geophys. Res. 107 (doi:10.1029/2001JC000888, 2002)
Stammer, D., Wunsch, C., Giering, R., Eckert, C., Heimbach, P., Marotzke, J., Adcroft, A., Hill, C.N., Marshall, J.: Volume, heat and freshwater transports of the global ocean circulation 1992–1997, estimated from a general circulation model constrained by WOCE data. J. Geophys. Res. (doi:10.1029/2001JC001115, 2002)
Marshall, J., Adcroft, A., Hill, C., Perelman, L., Heisey, C.: A Finite-Volume, Incompressible Navier Stokes Model for Studies of the Ocean on Parallel Computers. Technical Report 36, Massachusetts Institut of Technology, Center for Global Change Science, Cambridge, MA 02139, USA (1995)
Adcroft, A., Campin, J.M., Heimbach, P., Hill, C., Marshall, J.: The MITgcm. Online documentation, Massachusetts Institute of Technology, USA (2002)
Griffies, S.M., Harrison, M.J., Pacanowski, R.C., Rosati, A.: The FMS MOM4-beta User Guide. Technical report, NOAA/Geophysical Fluid Dynamics Laboratory (2002)
Galanti, E., Tziperman, E., Harrison, M., Rosati, A., Giering, R., Sirkes, Z.: The equatorial thermocline outcropping-a seasonal control on the tropical pacific ocean-atmosphere instability. Journal of Climate 15 (2002) 2721–2739
Rayner, P., Knorr, W., Scholze, M., Giering, R., Kaminski, T., Heimann, M., Quere, C.L.: Inferring terrestrial biosphere carbon fluxes from combined inversions of atmospheric transport and process-based terrestrial ecosystem models. In: Proceedings of 6th Carbon dioxide conference at Sendai. (2001) 1015–1017
Knorr, W.: Annual and interannual CO2 exchanges of the terrestrial biosphere: process based simulations and uncertainties. Glob. Ecol. and Biogeogr. 9 (2000) 225–252
Heimann, M.: The global atmospheric tracer model TM2. Technical Report No.10, Max-Planck-Institut für Meteorologie, Hamburg, Germany (1995)
Knorr, W.: Satellitengestützte Fernerkundung und Modellierung des Globalen CO2-Austauschs der Landvegetation: Eine Synthese. PhD thesis, Max-Planck-Inst. für Meteorol., Hamburg, Germany (1997)
Wilson, M.F., Henderson-Sellers, A.: A global archive of land cover and soils data for use in general-circulation climate models. Journal of Climatology 5 (1985) 119–143
Nijssen, B., Schnur, R., Lettenmaier, D.: Retrospective estimation of soil moisture using the vic land surface model, 1980–1993. J. Climate (2001) 1790–1808
Schnur, R. (personal communication)
GLOBALVIEW-CO2: Cooperative Atmospheric Data Integration Project-Carbon Dioxide. CD-ROM, NOAA CMDL, Boulder, Colorado (2001) [Also available on Internet via anonymous FTP to ftp://ftp.cmdl.noaa.gov, Path: ftp://ccg/co2/GLOBALVIEW]
Kaminski, T., Heimann, M., Giering, R.: A coarse grid three dimensional global inverse model of the atmospheric transport, 1, Adjoint model and Jacobian matrix. J. Geophys. Res. 104 (1999) 18,535–18,553
Knorr, W., Heimann, M.: Impact of drought stress and other factors on seasonal land biosphere CO2 exchange studied through an atmospheric tracer transport model. Tellus, Ser. B 47 (1995) 471–489
Kaminski, T., Knorr, W., Rayner, P., Heimann, M.: Assimilating atmospheric data into a terrestrial biosphere model: A case study of the seasonal cycle. Global Biogeochemical Cycles 16 (2002) 14-1-14-16
Rayner et al.: The history of terrestrial carbon fluxes from 1980-2000: Results from a Data Assimilation System. Global Biogeochem. Cycles (in preparation 2003)
Takahashi, T., Wanninkhof, R.H., Feely, R.A., Weiss, R.F., Chipman, D.W., Bates, N., Olafsson, J., Sabine, C., Sutherland, S.C.: Net sea-air CO2flux over the global oceans: An improved estimate based on the sea-air pCO2 difference. In Nojiri, Y., ed.: Extended abstracts of the 2nd International CO2 in the Oceans Symposium, Tsukuba, Japan, January 18–22, 1999. (1999) 9–15
Le Quéré, C., Orr, J.C., Monfray, P., Aumont, O., Madec, G.: Interannual variability of the oceanic sink of CO2 from 1979 through 1997. Global Biogeochem. Cycles 14 (2000) 1247–1265
Houghton, R.A., Boone, R.D., Fruci, J.R., Hobbie, J., Melillo, J.M., Palm, C.A., Peterson, B.J., Shaver, G.R., Woodwell, G.M., Moore, B., Skole, D.L., Myers, N.: The flux of carbon from terrestrial ecosystems to the atmosphere in 1980 due to changes in land use: Geographic distribution of the global flux. Tellus, Ser. B 39 (1987) 122–139
Andres, R.J., Marland, G., Boden, T., Bischoff, S.: Carbon dioxide emissions from fossil fuel consumption and cement manufacture 1751 to 1991 and an estimate for their isotopic composition and latitudinal distribution. In Wigley, T.M.L., Schimel, D., eds.: The Carbon Cycle. Cambridge Univ., New York, in press (1999)
Marland, G., Boden, T.A., Andres, R.J.: Global, regional, and national CO2 emissions. In: Trends: A Compendium of Data on Global Change. Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, U.S. Department of Energy, Oak Ridge, Tenn. (2001)
Tarantola, A.: Inverse Problem Theory-Methods for Data Fitting and Model Parameter Estimation. Elsevier Sci., New York (1987)
Enting, I.G.: Inverse Problems in Atmospheric Constituent Transport. Cambridge University Press, Cambridge (2002)
Gilbert, J.C., Lemaréchal, C.: Some numerical experiments with variable-storage quasi-Newton algorithms. Mathematical Programming 45 (1989) 407–435
Scholze, M., Rayner, P., Knorr, W., Kaminski, T., Giering, R.: A prototype Carbon Cycle Data Assimilation System (CCDAS): Inferring interannual variations of vegetation-atmosphere CO2 fluxes. Abstract CG62A-05. Eos Trans. AGU 83 (2002)
Griewank, A.: Achieving logarithmic growth of temporal and spatial complexity in reverse automatic differentiation. Optimization Methods and Software 1 (1992) 35–54
Griewank, A.: Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation. Number 19 in Frontiers in Appl. Math. SIAM, Philadelphia (2000)
Giering, R., Kaminski, T.: Using TAMC to generate efficient adjoint code: Comparison of automatically generated code for evaluation of first and second order derivatives to hand written code from the Minpack-2 collection. In Faure, C., ed.: Automatic Differentiation for Adjoint Code Generation. INRIA, Sophia Antipolis, France (1998) 31–37
Curtis, A.R., Powell, M.J.D., Reid, J.K.: On the estimation of sparse Jacobian matrices. J. Inst. Math. Appl. 13 (1974) 117–119
Newsam, G.N., Ramsdell, J.D.: Estimation of sparse Jacobian matrices. SIAM J. Alg. Disc. Meth. 4 (1983) 404–417
Geitner, U., Utke, J., Griewank, A.: Automatic computation of sparse Jacobians by applying the method of Newsam and Ramsdell. In Berz, M., Bischof, C., Corliss, G., Griewank, A., eds.: Computational Differentiation: Techniques Applications, and Tools. SIAM, Philadelphia, Penn. (1996) 161–172
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2003 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Kaminski, T., Giering, R., Scholze, M., Rayner, P., Knorr, W. (2003). An Example of an Automatic Differentiation-Based Modelling System. In: Kumar, V., Gavrilova, M.L., Tan, C.J.K., L’Ecuyer, P. (eds) Computational Science and Its Applications — ICCSA 2003. ICCSA 2003. Lecture Notes in Computer Science, vol 2668. Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-44843-8_11
Download citation
DOI: https://doi.org/10.1007/3-540-44843-8_11
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-40161-2
Online ISBN: 978-3-540-44843-3
eBook Packages: Springer Book Archive