共查询到20条相似文献，搜索用时 46 毫秒
1.
Alexander Barth Yajing Yan Aida Alvera-Azcárate Jean-Marie Beckers 《Ocean Dynamics》2016,66(12):1651-1664
Ensemble assimilation schemes applied in their original, global formulation respect linear conservation properties if the ensemble perturbations are set up accordingly. For realistic ocean systems, only a relatively small number of ensemble members can be calculated. A localization of the ensemble increment is therefore necessary to filter out spurious long-range correlations. The conservation of the global properties will be lost if the assimilation is performed locally, since the conservation requires a coupling between all model grid points which is removed by the localization. The distribution of ocean observations is often highly inhomogeneous. Systematic errors of the observed parts of the ocean state can lead to spurious adjustment of the non-observed parts via data assimilation and thus to a spurious increase or decrease in long-term simulations of global properties which should be conserved. In this paper, we propose a local assimilation scheme (with different variants and assumptions) which can satisfy global conservation properties. The proposed scheme can also be used for non-local observation operators. Different variants of the proposed scheme are tested in an idealized model and compared to the traditional covariance localization with an ad-hoc step enforcing conservation. It is shown that the inclusion of the conservation property reduces the total RMS error and that the presented stochastic and deterministic schemes avoiding error space rotation provide better results than the traditional covariance localization. 相似文献
2.
Li Xin Ma Hanqing Ran Youhua Wang Xufeng Zhu Gaofeng Liu Feng He Honglin Zhang Zhen Huang Chunlin 《中国科学:地球科学(英文版)》2021,64(10):1645-1657
The terrestrial carbon cycle is an important component of global biogeochemical cycling and is closely related to human well-being and sustainable development. However, large uncertainties exist in carbon cycle simulations and observations.Model-data fusion is a powerful technique that combines models and observational data to minimize the uncertainties in terrestrial carbon cycle estimation. In this paper, we comprehensively overview the sources and characteristics of the uncertainties in terrestrial carbon cycle models and observations. We present the mathematical principles of two model-data fusion methods, i.e., data assimilation and parameter estimation, both of which essentially achieve the optimal fusion of a model with observational data while considering the respective errors in the model and in the observations. Based upon reviewing the progress in carbon cycle models and observation techniques in recent years, we have highlighted the major challenges in terrestrial carbon cycle model-data fusion research, such as the "equifinality" of models, the identifiability of model parameters,the estimation of representativeness errors in surface fluxes and remote sensing observations, the potential role of the posterior probability distribution of parameters obtained from sensitivity analysis in determining the error covariance matrixes of the models, and opportunities that emerge by assimilating new remote sensing observations, such as solar-induced chlorophyll fluorescence. It is also noted that the synthesis of multisource observations into a coherent carbon data assimilation system is by no means an easy task, yet a breakthrough in this bottleneck is a prerequisite for the development of a new generation of global carbon data assimilation systems. This article also highlights the importance of carbon cycle data assimilation systems to generate reliable and physically consistent terrestrial carbon cycle reanalysis data products with high spatial resolution and longterm time series. These products are critical to the accurate estimation of carbon cycles at the global and regional scales and will help future carbon management strategies meet the goals of carbon neutrality. 相似文献
3.
Sujay V. Kumar Rolf H. Reichle Christa D. Peters-Lidard Randal D. Koster Xiwu Zhan Wade T. Crow John B. Eylander Paul R. Houser 《Advances in water resources》2008
The Land Information System (LIS) is an established land surface modeling framework that integrates various community land surface models, ground measurements, satellite-based observations, high performance computing and data management tools. The use of advanced software engineering principles in LIS allows interoperability of individual system components and thus enables assessment and prediction of hydrologic conditions at various spatial and temporal scales. In this work, we describe a sequential data assimilation extension of LIS that incorporates multiple observational sources, land surface models and assimilation algorithms. These capabilities are demonstrated here in a suite of experiments that use the ensemble Kalman filter (EnKF) and assimilation through direct insertion. In a soil moisture experiment, we discuss the impact of differences in modeling approaches on assimilation performance. Provided careful choice of model error parameters, we find that two entirely different hydrological modeling approaches offer comparable assimilation results. In a snow assimilation experiment, we investigate the relative merits of assimilating different types of observations (snow cover area and snow water equivalent). The experiments show that data assimilation enhancements in LIS are uniquely suited to compare the assimilation of various data types into different land surface models within a single framework. The high performance infrastructure provides adequate support for efficient data assimilation integrations of high computational granularity. 相似文献
4.
Local extreme rain usually resulted in disasters such as flash floods and landslides. Upon today, it is still one of the most difficult tasks for operational weather forecast centers to predict those events accurately. In this paper, we simulate an extreme precipitation event with ensemble Kalman filter (EnKF) assimilation of Doppler radial-velocity observations, and analyze the uncertainties of the assimilation. The results demonstrate that, without assimilation radar data, neither a single initialization of deterministic forecast nor an ensemble forecast with adding perturbations or multiple physical parameterizations can predict the location of strong precipitation. However, forecast was significantly improved with assimilation of radar data, especially the location of the precipitation. The direct cause of the improvement is the buildup of a deep mesoscale convection system with EnKF assimilation of radar data. Under a large scale background favorable for mesoscale convection, efficient perturbations of upstream mid-low level meridional wind and moisture are key factors for the assimilation and forecast. Uncertainty still exists for the forecast of this case due to its limited predictability. Both the difference of large scale initial fields and the difference of analysis obtained from EnKF assimilation due to small amplitude of initial perturbations could have critical influences to the event's prediction. Forecast could be improved through more cycles of EnKF assimilation. Sensitivity tests also support that more accurate forecasts are expected through improving numerical models and observations. 相似文献
5.
6.
The application of interferometric synthetic aperture radar (InSAR) has been increasingly used to improve capabilities to model land subsidence in hydrogeologic studies. A number of investigations over the last decade show how spatially detailed time‐lapse images of ground displacements could be utilized to advance our understanding for better predictions. In this work, we use simulated land subsidences as observed measurements, mimicking InSAR data to inversely infer inelastic specific storage in a stochastic framework. The inelastic specific storage is assumed as a random variable and modeled using a geostatistical method such that the detailed variations in space could be represented and also that the uncertainties of both characterization of specific storage and prediction of land subsidence can be assessed. The ensemble Kalman filter (EnKF), a real‐time data assimilation algorithm, is used to inversely calibrate a land subsidence model by matching simulated subsidences with InSAR data. The performance of the EnKF is demonstrated in a synthetic example in which simulated surface deformations using a reference field are assumed as InSAR data for inverse modeling. The results indicate: (1) the EnKF can be used successfully to calibrate a land subsidence model with InSAR data; the estimation of inelastic specific storage is improved, and uncertainty of prediction is reduced, when all the data are accounted for; and (2) if the same ensemble is used to estimate Kalman gain, the analysis errors could cause filter divergence; thus, it is essential to include localization in the EnKF for InSAR data assimilation. 相似文献
7.
Data assimilation as a method to predict variables, reduce uncertainties and explicitly handle various sources of uncertainties has recently received widespread attention and has been utilized to combine in situ and remotely sensed measurements with hydrological models. However, factors that significantly influence the capability of data assimilation still need testing and verifying. In this paper, synthetic surface soil moisture data are assimilated into the Soil and Water Assessment Tool (SWAT) model to evaluate their impact on other hydrological variables via the ensemble Kalman smoother (EnKS), using data from the Heihe River Basin, northwest China. The results show that the assimilation of surface soil moisture can moderately improve estimates of deep layer soil moisture, surface runoff and lateral flow, which reduces the negative influences of erroneous forcing and inaccurate parameters. The effects of the spatially heterogeneous input data (land cover and soil type) on the performance of the data assimilation technique are noteworthy. Moreover, the approaches including inflation and localization are specifically diagnosed to further extend the capability of the EnKS. 相似文献
8.
— Simulation of the complete earthquake generation process requires assimilation of observational data on long-term behavior of faults over multiple rupture cycles. Paleoseismology, the study of pre-instrumental earthquakes from geological field investigations, provides the only observations of long-term fault behavior. Paleoseismic data present challenges for assimilation into models because data sets are small, sparse, analog, and contain qualitative uncertainties. Observations can be categorized as primary measurements from field data, or non-primary extrapolations or interpretations of primary data. Assimilation of non-primary data could yield biased simulation results. We present a primary data set for the San Andreas and Imperial faults for comparison with non-primary data in existing databases, and propose standard formats and structure for paleoseismic data assimilation into numerical simulation models . 相似文献
9.
10.
We present a framework for design and deployment of decision support modeling based on metrics which have their roots in the scientific method. Application of these metrics to decision support modeling requires recognition of the importance of data assimilation and predictive uncertainty quantification in this type of modeling. The difficulties of implementing these procedures depend on the relationship between data that is available for assimilation and the nature of the prediction(s) that a decision support model is required to make. Three different data/prediction contexts are identified. Unfortunately, groundwater modeling is generally aligned with the most difficult of these. It is suggested that these difficulties can generally be ameliorated through appropriate model design. This design requires strategic abstraction of parameters and processes in a way that is optimal for the making of one particular prediction but is not necessarily optimal for the making of another. It is further suggested that the focus of decision support modeling should be on the ability of a model to provide receptacles for decision-pertinent information rather than on its purported ability to simulate environmental processes. While models are compromised in both of these roles, this view makes it clear that simulation should serve data assimilation and not the other way around. Data assimilation enables the uncertainties of decision-critical model predictions to be quantified and maybe reduced. Decision support modeling requires this. 相似文献
11.
A complex and highly dynamical ocean region, the Agulhas Current System plays an important role in the transfer of energy, nutrients and organic material from the Indian to the Atlantic Ocean. Its dynamics are not only important locally, but affect the global ocean-atmosphere system. In working towards improved ocean reanalysis and forecasting capabilities, it is important that numerical models simulate mesoscale variability accurately—especially given the scarcity of coherent observational platforms in the region. Data assimilation makes use of scarce observations, a dynamical model and their respective error statistics to estimate a new, improved model state that minimises the distance to the observations whilst preserving dynamical consistency. Qualitatively, it is unclear whether this minimisation directly translates to an improved representation of mesoscale dynamics. In this study, the impact of assimilating along-track sea-level anomaly (SLA) data into a regional Hybrid Coordinate Ocean Model (HYCOM) is investigated with regard to the simulation of mesoscale eddy characteristics. We use an eddy-tracking algorithm and compare the derived eddy characteristics of an assimilated (ASSIM) and an unassimilated (FREE) simulation experiment in HYCOM with gridded satellite altimetry-derived SLA data. Using an eddy tracking algorithm, we are able to quantitatively evaluate whether assimilation updates the model state estimate such that simulated mesoscale eddy characteristics are improved. Additionally, the analysis revealed limitations in the dynamical model and the data assimilation scheme, as well as artefacts introduced from the eddy tracking scheme. With some exceptions, ASSIM yields improvements over FREE in eddy density distribution and dynamics. Notably, it was found that FREE significantly underestimates the number of eddies south of Madagascar compared to gridded altimetry, with only slight improvements introduced through assimilation, highlighting the models’ limitation in sustaining mesoscale activity in this region. Interestingly, it was found that the threshold for the maximum eddy propagation velocity in the eddy detection scheme is often exceeded when data assimilation relocates an eddy, causing the algorithm to interpret the discontinuity as eddy genesis, which directly influences the eddy count, lifetime and propagation velocity, and indirectly influences other metrics such as non-linearity. Finally, the analysis allowed us to separate eddy kinetic energy into contributions from detected mesoscale eddies and meandering currents, revealing that the assimilation of SLA has a greater impact on mesoscale eddies than on meandering currents. 相似文献
12.
数据同化是提升复杂机理过程模型精度的关键技术之一,而湖泊藻类模型的敏感参数具有随时间动态变化的特征,导致数据同化过程中无法精准更新某一时段的敏感参数,影响数据同化的模型精度提升效果.针对上述问题,本研究耦合了参数敏感性分析与集合卡尔曼滤波,研发了一种能够实时识别模型敏感参数的新型数据同化算法;为验证研发算法的效率,依托巢湖的高频水质自动监测数据,测试算法对藻类动态模型的精度提升效果.测试结果表明：研发算法能够精准跟踪模型敏感参数的动态变化,并根据监测数据实时更新模型敏感参数,实现了水质高频自动监测数据与藻类动态模型的深度融合,藻类生物量模拟精度提升了55%,即纳什系数（NSE）从0.49提升到0.76,模拟精度提升效果也显著优于传统数据同化算法（NSE=0.63）.研发算法可应用于其它水生态环境模型的数据同化,为水生态环境相关要素的精准模拟预测提供关键技术支撑. 相似文献
13.
One-dimensional vertical and three-dimensional fine-resolution numerical models of sediment transport have been developed and applied to the Torres Strait region of northern Australia. The one-dimensional model, driven by measured waves and currents, was calibrated against measured suspended sediment concentrations using a sequential data assimilation algorithm. The algorithm produced a good match between model and data, but this was achieved only by allowing some temporal variability in parameter values, suggesting that there were underlying uncertainties in the model structure and forcing data. Implications of the assimilation results to the accuracy of the numerical modelling are discussed and the need for observational programmes having an extensive spatial and temporal coverage is highlighted. The three-dimensional sediment model, driven by modelled waves and currents, simulates sediment transport over the shelf during the monsoon and trade-wind seasons covering 1997–2000. The model predicts strong seasonal variability of the sediment transport on the shelf attributed to seasonally varying hydrodynamics, and illustrates significant inter-annual variability of the sediment fluxes driven by extreme events. The developed model provides a platform for testing scientific hypothesis. With additional calibration, including uncertainty analysis, it can also be used in a management context. 相似文献
14.
Land data assimilation (DA) has gradually developed into an important earth science research method because of its ability to combine model simulations and observations. Integrating new observations into a land surface model by the DA method can correct the predicted trajectory of the model and thus, improve the accuracy of state variables. It can also reduce uncertainties in the model by estimating some model parameters simultaneously. Among the various DA methods, the particle filter is free from the constraints of linear models and Gaussian error distributions, and can be applicable to any nonlinear and non-Gaussian state-space model; therefore, its importance in land data assimilation research has increased. In this study, a DA scheme was developed based on the residual resampling particle filter. Microwave brightness temperatures were assimilated into the macro-scale semi-distributed variance infiltration capacity model to estimate the surface soil moisture and three hydraulic parameters simultaneously. Finally, to verify the scheme, a series of comparative experiments was performed with experimental data obtained during the Soil Moisture Experiment of 2004 in Arizona. The results show that the scheme can improve the accuracy of soil moisture estimations significantly. In addition, the three hydraulic parameters were also well estimated, demonstrating the effectiveness of the DA scheme. 相似文献
15.
《Advances in water resources》2003,26(1):79-93
In order to improve the soil temperature profile predictions in land-surface models, an assimilation scheme using the extended Kalman filter is developed. This formulation is based on the discretized diffusion equation of heat transfer through the soil column. The scheme is designed to incorporate the knowledge of the uncertainties in both the model and the measurement. Model uncertainty is estimated by quantifying the model drift from observations when the model is initialized using the observed values. Furthermore, the initial error covariance has a significant influence on the performance of the assimilation scheme. It is shown that an inaccurate initial value for the error covariance can actually diminish the predictive capabilities of the model. When an appropriate initial error covariance is specified, using the top layer soil temperature observations in the assimilation scheme allows for improved predictive capabilities in lower layers. Observations at 30 min intervals have a significant effect on the model predictions in the lower layers. Assimilation of observations at 24 h intervals also has an effect on the lower layer predictive capability of the model, albeit more slowly than the 30 min assimilation scenario. 相似文献
16.
Ground shaking intensity varies spatially in earthquakes, and many studies have estimated correlations of intensity from past earthquake data. This paper presents a framework for quantifying uncertainty in the estimation of correlations and true variability in correlations from earthquake to earthquake. A procedure for evaluating estimation uncertainty is proposed and used to evaluate several methods that have been used in past studies to estimate correlations. The results indicate that a weighted least squares algorithm is most effective in estimating spatial correlation models and that earthquakes with at least 100 recordings are needed to produce informative earthquake-specific estimates of spatial correlations. The proposed procedure is also used to distinguish between estimation uncertainty and the true variability in model parameters that exist in a given data set. The estimation uncertainty is seen to vary between well-recorded and poorly recorded earthquakes, whereas the true variability is more stable. 相似文献
17.
Ensemble Kalman filter, EnKF, as a Monte Carlo sequential data assimilation method has emerged promisingly for subsurface media characterization during past decade. Due to high computational cost of large ensemble size, EnKF is limited to small ensemble set in practice. This results in appearance of spurious correlation in covariance structure leading to incorrect or probable divergence of updated realizations. In this paper, a universal/adaptive thresholding method is presented to remove and/or mitigate spurious correlation problem in the forecast covariance matrix. This method is, then, extended to regularize Kalman gain directly. Four different thresholding functions have been considered to threshold forecast covariance and gain matrices. These include hard, soft, lasso and Smoothly Clipped Absolute Deviation (SCAD) functions. Three benchmarks are used to evaluate the performances of these methods. These benchmarks include a small 1D linear model and two 2D water flooding (in petroleum reservoirs) cases whose levels of heterogeneity/nonlinearity are different. It should be noted that beside the adaptive thresholding, the standard distance dependant localization and bootstrap Kalman gain are also implemented for comparison purposes. We assessed each setup with different ensemble sets to investigate the sensitivity of each method on ensemble size. The results indicate that thresholding of forecast covariance yields more reliable performance than Kalman gain. Among thresholding function, SCAD is more robust for both covariance and gain estimation. Our analyses emphasize that not all assimilation cycles do require thresholding and it should be performed wisely during the early assimilation cycles. The proposed scheme of adaptive thresholding outperforms other methods for subsurface characterization of underlying benchmarks. 相似文献
18.
How Bayesian data assimilation can be used to estimate the mathematical structure of a model 总被引：1，自引：1，他引：0
Nataliya Bulygina Hoshin Gupta 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(6):925-937
In previous work, we presented a method for estimation and correction of non-linear mathematical model structures, within a Bayesian framework, by merging uncertain knowledge about process physics with uncertain and incomplete observations of dynamical input-state-output behavior. The resulting uncertainty in the model input-state-output mapping is expressed as a weighted combination of an uncertain conceptual model prior and a data-derived probability density function, with weights depending on the conditional data density. Our algorithm is based on the use of iterative data assimilation to update a conceptual model prior using observed system data, and thereby construct a posterior estimate of the model structure (the mathematical form of the equation itself, not just its parameters) that is consistent with both physically based prior knowledge and with the information in the data. An important aspect of the approach is that it facilitates a clear differentiation between the influences of different types of uncertainties (initial condition, input, and mapping structure) on the model prediction. Further, if some prior assumptions regarding the structural (mathematical) forms of the model equations exist, the procedure can help reveal errors in those forms and how they should be corrected. This paper examines the properties of the approach by investigating two case studies in considerable detail. The results show how, and to what degree, the structure of a dynamical hydrological model can be estimated without little or no prior knowledge (or under conditions of incorrect prior information) regarding the functional forms of the storage–streamflow and storage–evapotranspiration relationships. The importance and implications of careful specification of the model prior are illustrated and discussed. 相似文献
19.
Probabilistic eruption forecasting at short and long time scales 总被引：1，自引：1，他引：0
Any effective volcanic risk mitigation strategy requires a scientific assessment of the future evolution of a volcanic system and its eruptive behavior. Some consider the onus should be on volcanologists to provide simple but emphatic deterministic forecasts. This traditional way of thinking, however, does not deal with the implications of inherent uncertainties, both aleatoric and epistemic, that are inevitably present in observations, monitoring data, and interpretation of any natural system. In contrast to deterministic predictions, probabilistic eruption forecasting attempts to quantify these inherent uncertainties utilizing all available information to the extent that it can be relied upon and is informative. As with many other natural hazards, probabilistic eruption forecasting is becoming established as the primary scientific basis for planning rational risk mitigation actions: at short-term (hours to weeks or months), it allows decision-makers to prioritize actions in a crisis; and at long-term (years to decades), it is the basic component for land use and emergency planning. Probabilistic eruption forecasting consists of estimating the probability of an eruption event and where it sits in a complex multidimensional time–space–magnitude framework. In this review, we discuss the key developments and features of models that have been used to address the problem. 相似文献
20.
Remotely sensed observations of sea-level anomaly and sea-surface temperature have been assimilated into an implementation
of the Miami Isopycnic Coordinate Ocean Model (MICOM) for the Indian Ocean using the Ensemble Kalman Filter (EnKF). The system
has been applied in a hindcast validation experiment to examine the properties of the assimilation scheme when used with a
full ocean general circulation model and real observations. This work is considered as a first step towards an operational
ocean monitoring and forecasting system for the Indian Ocean. The assimilation of real data has demonstrated that the sequential
EnKF can efficiently control the model evolution in time. The use of data assimilation requires a significant amount of additional
processing and computational resources. However, we have tried to justify the cost of using a sophisticated assimilation scheme
by demonstrating strong regional and temporal dependencies of the covariance statistics, which include highly anisotropic
and flow-dependent correlation functions. In particular, we observed a marked difference between error statistics in the equatorial
region and at off-equatorial latitudes. We have also demonstrated how the assimilation of SLA and SST improves the model fields
with respect to real observations. Independent in situ temperature profiles have been used to examine the impact of assimilating
the remotely sensed observations. These intercomparisons have shown that the model temperature and salinity fields better
resemble in situ observations in the assimilation experiment than in a model free-run case. On the other hand, it is also
expected that assimilation of in situ profiles is needed to properly control the deep ocean circulation.
Received: 8 January 2002 / Accepted: 8 April 2002 相似文献