共查询到20条相似文献，搜索用时 15 毫秒
1.
利用伴随法优化非线性潮汐模型的开边界条件Ⅱ.黄海、东海潮汐资料的同化试验 总被引：12，自引：4，他引：12
本研究基于最优控制理论,采用变分数据同化法,通过建立伴随模型,把观测资料同化到陆架海域潮汐数值模型中去,优化开边界条件,以便提高数值预报的精度.潮汐模型的控制方程为考虑平流项、非线性底摩擦和侧向涡动粘性项的非线性浅水方程组.在第Ⅰ部分建立伴随模型和进行“孪生”数值试验的基础上,给出利用验潮站的水位资料以及TOPEX/Poseidon卫星测高数据在黄海、东海进行变分数据同化试验的数值结果.试验表明利用上述资料对模型进行变分同化校正是可行的. 相似文献
2.
《Estuarine, Coastal and Shelf Science》2003,57(5-6):873-885
In the present work we explore the impact of assimilating local tide-gauge and altimetric data on the quality of predicting the major Adriatic tides (M2 and K1). To that end we compute optimal tidal open boundary conditions for a 3D high-resolution finite-element model by using an incremental assimilation formalism. The essence of the method is the use of two dynamical models where the solution in the complex 3D high-resolution model is sought via assimilation of prediction errors into the simpler 2D model with explicit inverse. In the central numerical experiment, harmonic constants from 12 tide gauges are assimilated and the results are analysed at 31 locations, hence 19 independent ones. The data assimilation contributes to the reduction of maximum amplitude error from 5.6 to 0.5 cm for M2 and from 3.9 to 0.1 cm for K1. The assimilation procedure is repeated by assimilating suitably processed Topex/Poseidon altimeter data, again validating the outcome at 31 tide gauge locations. The result was very similar to the gauge-data assimilation outcome. The model output is also validated with the current data, not used in the assimilation. At two locations and at three depths the model was able to reproduce the major and the minor semi-axes of tidal ellipses, as well as their orientations very well. 相似文献
3.
The variational assimilation method has been examined for ability of reconstructing mesoscale features in altimeter data using a simple dynamic model. A one-dimensional, two-layer Rossby wave model in a cross-track channel has been chosen. The simulated data are constructed from a theoretical solution, which is composed of any combination of two normal vertical (barotropic and baroclinic) modes. The data are collected along tracks and with repeat periods similar to those of the Geosat altimeter. The phase space of control variables is composed of initial and boundary conditions. A cost function is defined to measure differences between the simulated data and the model solution. Regularization (smoothing) terms are also included in the cost function in the form of secon-order spatial and time derivatives of the solution. In this paper, two potential problems existing in the altimeter data assimilation are addressed: one is low cross-track resolution, and the other is vertical projection of the data measured at the sea surface. A succesful metho is developed for reconstructing Rossby waves with wavelengths as short as twice the track intervals for any combination of two vertical modes. A key component to efficient assimilation is a preparation step prior to the actual variational assimilation: a uniform ratio of pressure amplitudes in the two layers is included as an optimization parameter. Starting with the first guess from the preparation step, the variational method is carried out based on adjoint equations without such constraint. Separation of the control variables into the two subsets of the initial and the boundary conditions is found useful. Characteristics of the Hessian matrix are related to the performance of this technique. The method developed for the linear system implies steps to be included in data assimilation for nonlinear meanders and eddies in a major current system as well. 相似文献
4.
Gennady A. Kivman 《Progress in Oceanography》1997,40(1-4)
Estimation of the open-boundary inputs solving a weak constraint variational formulation for an Arctic tide model is considered as an ill-posed problem in the sense that the solution is very sensitive to the data noise and to grid size. Mathematically, spatial discretization of a cost function to be minimized and penalization of normal flow through the open boundary act as regularization of the problem. An heuristic choosing rule for the regularization parameter is applied to assess a suitable spatial resolution and the weight referred to the open boundary penalty. It is shown that these provide a better fit of the solution to a control data set compared with a finer grid, the value of the energy flux through the open boundary being in agreement with other model estimates. The M2 solution obtained is much closer to the control data than other modern solutions while the accuracy of the simulated K1 constituent is within the same error level. The tidal maps for these waves exhibit certain distinctions in comparison with other charts. 相似文献
5.
黄海、渤海TOPEX/Poseidon高度计资料潮汐伴随同化 总被引：1，自引：0，他引：1
首先将大约10a的TOPEX/Poseidon(T/P)高度计资料沿星下轨迹点做潮汐调和分析,提取得到各分潮的调和常数,利用伴随同化方法,同化到二维非线性潮汐数值模式中,模拟了黄海、渤海区域M2,S2,O1,K1等4个潮汐分潮,并根据计算结果给出了各分潮的同潮图.将计算值与观测值的进行偏差统计,结果表明计算值与验潮站资料符合良好.研究过程中做了两类试验:一类试验是针对不同的参数进行优化,一类试验是针对不同的资料进行同化.第一类试验表明:将开边界条件和底摩擦系数同时作为模型优化的控制参数,其结果明显优于单独优化开边界条件;第二类试验表明:同时同化高度计资料与验潮站资料,比单独同化其中任一种资料,对模式计算结果都有较好的改进.研究结果表明,采用伴随同化方法,利用T/P高度计资料和验潮站资料作为同化数据能有效改进模拟结果,用来反演黄海、渤海的潮波系统是可行的. 相似文献
6.
《Ocean Modelling》2003,5(1):37-63
A stabilized finite-element (FE) algorithm for the solution of oceanic large scale circulation equations and optimization of the solutions is presented. Pseudo-residual-free bubble function (RFBF) stabilization technique is utilized to enforce robustness of the numerics and override limitations imposed by the Babuška–Brezzi condition on the choice of functional spaces. The numerical scheme is formulated on an unstructured tetrahedral 3d grid in velocity–pressure variables defined as piecewise linear continuous functions. The model is equipped with a standard variational data assimilation scheme, capable to perform optimization of the solutions with respect to open lateral boundary conditions and external forcing imposed at the ocean surface. We demonstrate the model performance in applications to idealized and realistic basin-scale flows. Using the adjoint method, the code is tested against a synthetic climatological data set for the South Atlantic ocean which includes hydrology, fluxes at the ocean surface and satellite altimetry. The optimized solution proves to be consistent with all these data sets, fitting them within the error bars.The presented diagnostic tool retains the advantages of existing FE ocean circulation models and in addition (1) improves resolution of the bottom boundary layer due to employment of the 3d tetrahedral elements; (2) enforces numerical robustness through utilization of the RFBF stabilization, and (3) provides an opportunity to optimize the solutions by means of 3d variational data assimilation. Numerical efficiency of the code makes this a desirable tool for dynamically constrained analyses of large datasets. 相似文献
7.
Correcting the errors in the initial conditions and wind stress in storm surge simulation using an adjoint optimal technique 总被引：1，自引：0，他引：1
An adjoint data assimilation methodology is applied to the Princeton Ocean Model and is evaluated by obtaining “optimal” initial conditions, sea surface forcing conditions, or both for coastal storm surge modelling. By prescribing different error sources and setting the corresponding control variables, we performed several sets of identical twin experiments by assimilating model-generated water levels. The experiment results show that, when the forecasting errors are caused by the initial (or surface boundary) conditions, adjusting initial (or surface boundary) conditions accordingly can significantly improve the storm surge simulation. However, when the forecasting errors are caused by surface boundary (or initial) conditions, data assimilation targeting improving the initial (or surface boundary) conditions is ineffective. When the forecasting errors are caused by both the initial and surface boundary conditions, adjusting both the initial and surface boundary conditions leads to the best results. In practice, we do not know whether the errors are caused by initial conditions or surface boundary conditions, therefore it is better to adjust both initial and surface boundary conditions in adjoint data assimilation. 相似文献
8.
9.
海洋中尺度涡在本质上是属于满足准地转平衡的大尺度运动,因此理论上,其在短时间内的运动将主要受到准地转平衡关系的约束,而外部强迫场的影响在短期内不会明显改变其运动特征。基于上述思想,我们提出了一种基于四维变分同化初始场的中尺度涡旋预报方案。为了检验该方案的可行性,本文使用区域海洋模式(regional ocean modeling system, ROMS)和其内建的增量强约束四维变分同化(incremental strong constraint four dimensional variational, I4D-Var)模块,建立了一个南海海洋同化模拟系统。首先,通过I4D-Var方法将AVISO卫星高度计资料同化到海洋数值模拟中,获得了理想的中尺度涡同化模拟结果。同化、模式模拟和观测三者的中尺度涡统计结果表明,该同化系统模拟的南海中尺度涡的路径、半径、海表高度异常和振幅等特征信息与AVISO(Archiving ValidationandInterpolationofSatelliteOceanographicData)观测结果高度吻合,同时在深度上的分析表明,涡旋对应的温度、盐度和密度均得到有效的调整。然后,将该同化系统的模拟结果做为初始场,对某一特定时段的南海中尺度涡进行了后报模拟和结果的定量化分析。通过比较后报模拟与观测资料中对应涡旋的海表面高度异常(sea surface height anomalies, SSHA)相关系数、涡心差距和半径绝对误差,证明该方案的中尺度涡后报时效至少可达10 d以上。后报实验结果验证了该中尺度涡预报方案的可行性,从而为中尺度涡的预报提供一定的理论基础和可行性方案。 相似文献
10.
A. A. Zelenko R. M. Vil’fand Yu. D. Resnyanskii B. S. Strukov M. D. Tsyrulnikov P. I. Svirenko 《Izvestiya Atmospheric and Oceanic Physics》2016,52(4):443-454
A new version of the ocean data assimilation system (ODAS) developed at the Hydrometcentre of Russia is presented. The assimilation is performed following the sequential scheme analysis–forecast–analysis. The main components of the ODAS are procedures for operational observation data processing, a variational analysis scheme, and an ocean general circulation model used to estimate the first guess fields involved in the analysis. In situ observations of temperature and salinity in the upper 1400-m ocean layer obtained from various observational platforms are used as input data. In the new ODAS version, the horizontal resolution of the assimilating model and of the output products is increased, the previous 2D-Var analysis scheme is replaced by a more general 3D-Var scheme, and a more flexible incremental analysis updating procedure is introduced to correct the model calculations. A reanalysis of the main World Ocean hydrophysical fields over the 2005–2015 period has been performed using the updated ODAS. The reanalysis results are compared with data from independent sources. 相似文献
11.
12.
13.
利用热带太平洋TAO浮标阵列评估全球高分辨率浪-潮-流耦合模式同化结果 总被引：3，自引：0，他引：3
In order to evaluate the assimilation results from a global high resolution ocean model, the buoy observations from tropical atmosphere ocean(TAO) during August 2014 to July 2015 are employed. The horizontal resolution of wave-tide-circulation coupled ocean model developed by The First Institute of Oceanography(FIOCOM model) is 0.1°×0.1°, and ensemble adjustment Kalman filter is used to assimilate the sea surface temperature(SST), sea level anomaly(SLA) and Argo temperature/salinity profiles. The simulation results with and without data assimilation are examined. First, the overall statistic errors of model results are analyzed. The scatter diagrams of model simulations versus observations and corresponding error probability density distribution show that the errors of all the observed variables, including the temperature, isotherm depth of 20°C(D20), salinity and two horizontal component of velocity are reduced to some extent with a maximum improvement of 54% after assimilation. Second, time-averaged variables are used to investigate the horizontal and vertical structures of the model results. Owing to the data assimilation, the biases of the time-averaged distribution are reduced more than70% for the temperature and D20 especially in the eastern Pacific. The obvious improvement of D20 which represents the upper mixed layer depth indicates that the structure of the temperature after the data assimilation becomes more close to the reality and the vertical structure of the upper ocean becomes more reasonable. At last,the physical processes of time series are compared with observations. The time evolution processes of all variables after the data assimilation are more consistent with the observations. The temperature bias and RMSE of D20 are reduced by 76% and 56% respectively with the data assimilation. More events during this period are also reproduced after the data assimilation. Under the condition of strong 2014/2016 El Ni?o, the Equatorial Undercurrent(EUC) from the TAO is gradually increased during August to November in 2014, and followed by a decreasing process. Since the improvement of the structure in the upper ocean, these events of the EUC can be clearly found in the assimilation results. In conclusion, the data assimilation in this global high resolution model has successfully reduced the model biases and improved the structures of the upper ocean, and the physical processes in reality can be well produced. 相似文献
14.
An ensemble optimal interpolation (EnOI) data assimilation method is applied in the BCC_CSM1.1 to investigate the impact of ocean data assimilations on seasonal forecasts in an idealized twin experiment framework. Pseudo-observations of sea surface temperature (SST), sea surface height (SSH), sea surface salinity (SSS), temperature and salinity (T/S) profiles were first generated in a free model run. Then, a series of sensitivity tests initialized with predefined bias were conducted for a one-year period; this involved a free run (CTR) and seven assimilation runs. These tests allowed us to check the analysis field accuracy against the “truth”. As expected, data assimilation improved all investigated quantities; the joint assimilation of all variables gave more improved results than assimilating them separately. One-year predictions initialized from the seven runs and CTR were then conducted and compared. The forecasts initialized from joint assimilation of surface data produced comparable SST root mean square errors to that from assimilation of T/S profiles, but the assimilation of T/S profiles is crucial to reduce subsurface deficiencies. The ocean surface currents in the tropics were better predicted when initial conditions produced by assimilating T/S profiles, while surface data assimilation became more important at higher latitudes, particularly near the western boundary currents. The predictions of ocean heat content and mixed layer depth are significantly improved initialized from the joint assimilation of all the variables. Finally, a central Pacific El Ni?o was well predicted from the joint assimilation of surface data, indicating the importance of joint assimilation of SST, SSH, and SSS for ENSO predictions. 相似文献
15.
Variational data analysis with control of the forecast bias 总被引：1，自引：0，他引：1
We propose a methodology for the treatment of the systematic model error in variational data assimilation. The principle of the method is to add a systematic error correction term in the model equations and to include it in the variational assimilation control vector.
This method is applied to a simplified ocean circulation model in an identical twin experiment framework. It shows a noticeable improvement compared to the result of a classical variational assimilation scheme in which the systematic error is not corrected. The estimated systematic error correction term is sufficiently consistent with that needed by the model that it allows improvements not just to the analysis, but also during the forecast phase. 相似文献
This method is applied to a simplified ocean circulation model in an identical twin experiment framework. It shows a noticeable improvement compared to the result of a classical variational assimilation scheme in which the systematic error is not corrected. The estimated systematic error correction term is sufficiently consistent with that needed by the model that it allows improvements not just to the analysis, but also during the forecast phase. 相似文献
16.
This study deals with the assimilation of TOPEX altimeter-derived sea level variability in a reduced-gravity model of the northwestern Indian Ocean. The assimilation has been done using variational assimilation method with adjoint technique. A cost function representing the misfit between the model and the data is minimized with the model equations acting as constraints. The initial conditions of the model are used as control parameters and the best-fit initial conditions are determined as a result of minimization carried out using a variable-storage quasi-Newton method. Experiment has been done with 10 days' cycle-averaged data of TOPEX. Assimilation for 10 days and 20 days has been performed and it has been seen that 20 days' assimilation of satellite data provides better results. 相似文献
17.
一个概化的潮汐河口羽状流动力学的初步研究 总被引：1，自引：0，他引：1
基于区域陆架海水动力-生态耦合模型(COHERENS)中的三维水动力模型,对一个概化的潮汐河口羽状流动力学进行了初步研究。在均一水深和无风的理想条件下,模拟得到了一个潮周期内概化的潮汐河口表层水体:a)羽状流流场和盐度场的平面分布;b)3种不同河流径流量作用下的羽状流流场和盐度场的平面分布;c)4种不同底床粗糙长度下的羽状流流场和盐度场的平面分布。结果分析表明:a)羽状流的长度在潮汐的作用下递增,而其最大宽度近似呈周期性的变化;b)羽状流的长度和最大宽度都随着河流径流量的增大而分别变长和变宽;c)羽状流的长度和最大宽度都随底床粗糙长度的增大而变小;d)无潮情况下羽状流的凸起明显并存在沿岸流,然而在有潮情况下凸起不明显且无沿岸流;e)水平扩散可能限制了沿岸羽状流的发展。 相似文献
18.
The impact of assimilating Argo data into an initial field on the short-term forecasting accuracy of temper- ature and salinity is quantitatively estimated by using a forecasting system of the western North Pacific, on the base of the Princeton ocean model with a generalized coordinate system （POMgcs）. This system uses a sequential multigrid three-dimensional variational （3DVAR） analysis scheme to assimilate observation da- ta. Two numerical experiments were conducted with and without Argo temperature and salinity profile data besides conventional temperature and salinity profile data and sea surface height anomaly （SSHa） and sea surface temperature （SST） in the process of assimilating data into the initial fields. The forecast errors are estimated by using independent temperature and salinity profiles during the forecasting period, including the vertical distributions of the horizontally averaged root mean square errors （H-RMSEs） and the horizontal distributions of the vertically averaged mean errors （MEs） and the temporal variation of spatially averaged root mean square errors （S-RMSEs）. Comparison between the two experiments shows that the assimila- tion of Argo data significantly improves the forecast accuracy, with 24% reduction of H-RMSE maximum for the temperature, and the salinity forecasts are improved more obviously, averagely dropping of 50% for H-RMSEs in depth shallower than 300 m. Such improvement is caused by relatively uniform sampling of both temperature and salinity from the Argo drifters in time and space. 相似文献
19.
中国海及邻近海域卫星观测资料同化试验 总被引：4，自引：0，他引：4
利用1个基于POMgcs海洋模式和多重网格三维变分同化方法建立的中国海及邻近海域海面高与三维温盐流数值预报模型,通过一系列数值试验,研究了同化卫星测高和卫星遥感海面温度观测资料对该模型预报能力的影响。试验结果表明,同化卫星测高资料可明显改善海面高度与三维温度和盐度的分析预报效果,使1 200 m以上的温度预报误差减小0.16℃,并能有效提高对海洋中尺度现象的预报能力;同化卫星遥感海面温度对100 m以上的温度和盐度的预报效果有所改善,可使海面温度的预报误差减小10%。 相似文献
20.
A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important. 相似文献