Preprints
https://doi.org/10.5194/npg-2021-9
https://doi.org/10.5194/npg-2021-9

  19 Mar 2021

19 Mar 2021

Review status: this preprint is currently under review for the journal NPG.

Identification of linear response functions from arbitrary perturbation experiments in the presence of noise – Part I. Method development and toy model demonstration

Guilherme L. Torres Mendonça1,2, Julia Pongratz2,3, and Christian H. Reick2 Guilherme L. Torres Mendonça et al.
  • 1International Max Planck Research School on Earth System Modelling, Hamburg, Germany
  • 2Max Planck Institute for Meteorology, Hamburg, Germany
  • 3Ludwig-Maxmillians-Universität München, Munich, Germany

Abstract. Existent methods to identify linear response functions from data require tailored perturbation experiments, e.g. impulse or step experiments. And if the system is noisy, these experiments need to be repeated several times to obtain a good statistics. In contrast, for the method developed here, data from only a single perturbation experiment at arbitrary perturbation is sufficient if in addition data from an unperturbed (control) experiment is available. To identify the linear response function for this ill-posed problem we invoke regularization theory. The main novelty of our method lies in the determination of the level of background noise needed for a proper estimation of the regularization parameter: This is achieved by comparing the frequency spectrum of the perturbation experiment with that of the additional control experiment. The resulting noise level estimate can be further improved for linear response functions known to be monotonic. The robustness of our method and its advantages are investigated by means of a toy model. We discuss in detail the dependence of the identified response function on the quality of the data (signal-to-noise ratio) and on possible nonlinear contributions to the response. The method development presented here prepares in particular for the identification of carbon-cycle response functions in Part II of this study. But the core of our method, namely our new approach to obtain the noise level for a proper estimation of the regularization parameter, may find applications in solving also other types of linear ill-posed problems.

Guilherme L. Torres Mendonça et al.

Status: open (until 14 May 2021)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on npg-2021-9', Anonymous Referee #1, 24 Mar 2021 reply
    • AC1: 'Reply on RC1', Guilherme Torres Mendonça, 31 Mar 2021 reply
  • RC2: 'Comment on npg-2021-9', Anonymous Referee #2, 28 Mar 2021 reply
    • AC2: 'Reply on RC2', Guilherme Torres Mendonça, 08 Apr 2021 reply
  • CC1: 'Comment on npg-2021-9', Valerio Lucarini, 10 Apr 2021 reply

Guilherme L. Torres Mendonça et al.

Model code and software

Supplementary material for “Identification of linear response functions from arbitrary perturbation experiments in the presence of noise - Part I. Method development and toy model demonstration” Torres Mendonca, G., Pongratz, J., and Reick, C. H. https://pure.mpg.de/pubman/faces/ViewItemOverviewPage.jsp?itemId=item_3288741

Guilherme L. Torres Mendonça et al.

Viewed

Total article views: 316 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
268 37 11 316 2 2
  • HTML: 268
  • PDF: 37
  • XML: 11
  • Total: 316
  • BibTeX: 2
  • EndNote: 2
Views and downloads (calculated since 19 Mar 2021)
Cumulative views and downloads (calculated since 19 Mar 2021)

Viewed (geographical distribution)

Total article views: 311 (including HTML, PDF, and XML) Thereof 311 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 10 Apr 2021