Aims Data assimilation is a useful tool to extract information from large datasets of the net ecosystem exchange(NEE)of CO_(2) obtained by eddy-flux measurements.However,the number of parameters in ecosystem models th...Aims Data assimilation is a useful tool to extract information from large datasets of the net ecosystem exchange(NEE)of CO_(2) obtained by eddy-flux measurements.However,the number of parameters in ecosystem models that can be constrained by eddy-flux data is limited by conventional inverse analysis that estimates parameter values based on one-time inversion.This study aimed to improve data assimilation to increase the number of constrained parameters.Methods In this study,we developed conditional Bayesian inversion to maximize the number of parameters to be constrained by NEE data in several steps.In each step,we conducted a Bayesian inversion to constrain parameters.The maximum likelihood estimates of the constrained parameters were then used as prior to fix parameter values in the next step of inversion.The conditional inversion was repeated until there were no more parameters that could be further constrained.We applied the conditional inversion to hourly NEE data from Harvard Forest with a physiologically based ecosystem model.Important Findings Results showed that the conventional inversion method constrained 6 of 16 parameters in the model while the conditional inversion method constrained 13 parameters after six steps.The cost function that indicates mismatch between the modeled and observed data decreased with each step of conditional Bayesian inversion.The Bayesian information criterion also decreased,suggesting reduced information loss with each step of conditional Bayesian inversion.A wavelet analysis reflected that model performance under conditional Bayesian inversion was better than that under conventional inversion at multiple time scales,except for seasonal and half-yearly scales.In addition,our analysis also demonstrated that parameter convergence in a subsequent step of the conditional inversion depended on correlations with the parameters constrained in a previous step.Overall,the conditional Bayesian inversion substantially increased the number of parameters to be constrained by NEE data and can be a powerful tool to be used in data assimilation in ecology.展开更多
Time dependent carbon transfer coefficients are estimated using ecosystem exchange databy minimizing over variable observational intervals,Kalman filter,and variational minimization techniques.Transfer coefficients ar...Time dependent carbon transfer coefficients are estimated using ecosystem exchange databy minimizing over variable observational intervals,Kalman filter,and variational minimization techniques.Transfer coefficients are determined by application of estimation procedures to subintervalsfrom a partition of the observational time period,minimizing the variance of analyzed errors withoutthe imposition of a priori transfer coefficient models in Kalman filters,and minimization with respectto transfer coefficients in variational fit-to-data functionals.Results are compared between methodsand seasonal variability is observed in the transfer coefficients.展开更多
基金National Science Foundation(DEB 0444518,DEB 0743778)Office of Science(BER),Department of Energy(DE-FG02-006ER64319)Midwestern Regional Center of the National Institute for Climatic Change Research at Michigan Technological University(Award Number DE-FC02-06ER64158).
文摘Aims Data assimilation is a useful tool to extract information from large datasets of the net ecosystem exchange(NEE)of CO_(2) obtained by eddy-flux measurements.However,the number of parameters in ecosystem models that can be constrained by eddy-flux data is limited by conventional inverse analysis that estimates parameter values based on one-time inversion.This study aimed to improve data assimilation to increase the number of constrained parameters.Methods In this study,we developed conditional Bayesian inversion to maximize the number of parameters to be constrained by NEE data in several steps.In each step,we conducted a Bayesian inversion to constrain parameters.The maximum likelihood estimates of the constrained parameters were then used as prior to fix parameter values in the next step of inversion.The conditional inversion was repeated until there were no more parameters that could be further constrained.We applied the conditional inversion to hourly NEE data from Harvard Forest with a physiologically based ecosystem model.Important Findings Results showed that the conventional inversion method constrained 6 of 16 parameters in the model while the conditional inversion method constrained 13 parameters after six steps.The cost function that indicates mismatch between the modeled and observed data decreased with each step of conditional Bayesian inversion.The Bayesian information criterion also decreased,suggesting reduced information loss with each step of conditional Bayesian inversion.A wavelet analysis reflected that model performance under conditional Bayesian inversion was better than that under conventional inversion at multiple time scales,except for seasonal and half-yearly scales.In addition,our analysis also demonstrated that parameter convergence in a subsequent step of the conditional inversion depended on correlations with the parameters constrained in a previous step.Overall,the conditional Bayesian inversion substantially increased the number of parameters to be constrained by NEE data and can be a powerful tool to be used in data assimilation in ecology.
文摘Time dependent carbon transfer coefficients are estimated using ecosystem exchange databy minimizing over variable observational intervals,Kalman filter,and variational minimization techniques.Transfer coefficients are determined by application of estimation procedures to subintervalsfrom a partition of the observational time period,minimizing the variance of analyzed errors withoutthe imposition of a priori transfer coefficient models in Kalman filters,and minimization with respectto transfer coefficients in variational fit-to-data functionals.Results are compared between methodsand seasonal variability is observed in the transfer coefficients.