High Frequency(HF) radar current data is assimilated into a shelf sea circulation model based on optimal interpolation(OI) method. The purpose of this work is to develop a real-time computationally highly efficient as...High Frequency(HF) radar current data is assimilated into a shelf sea circulation model based on optimal interpolation(OI) method. The purpose of this work is to develop a real-time computationally highly efficient assimilation method to improve the forecast of shelf current. Since the true state of the ocean is not known, the specification of background error covariance is arduous. Usually, it is assumed or calculated from an ensemble of model states and is kept in constant. In our method, the spatial covariances of model forecast errors are derived from differences between the adjacent model forecast fields, which serve as the forecast tendencies. The assumption behind this is that forecast errors can resemble forecast tendencies, since variances are large when fields change quickly and small when fields change slowly. The implementation of HF radar data assimilation is found to yield good information for analyses. After assimilation, the root-mean-square error of model decreases significantly. Besides, three assimilation runs with variational observation density are implemented. The comparison of them indicates that the pattern described by observations is much more important than the amount of observations. It is more useful to expand the scope of observations than to increase the spatial interval. From our tests, the spatial interval of observation can be 5 times bigger than that of model grid.展开更多
Surface reconstruction from unorganized data points is a challenging problem in Computer Aided Design and Geometric Modeling. In this paper, we extend the mathematical model proposed by Juttler and Felis （Adv. Comput...Surface reconstruction from unorganized data points is a challenging problem in Computer Aided Design and Geometric Modeling. In this paper, we extend the mathematical model proposed by Juttler and Felis （Adv. Comput. Math., 17 （2002）, pp. 135-152） based on tensor product algebraic spline surfaces from fixed meshes to adaptive meshes. We start with a tensor product algebraic B-spline surface defined on an initial mesh to fit the given data based on an optimization approach. By measuring the fitting errors over each cell of the mesh, we recursively insert new knots in cells over which the errors are larger than some given threshold, and construct a new algebraic spline surface to better fit the given data locally. The algorithm terminates when the error over each cell is less than the threshold. We provide some examples to demonstrate our algorithm and compare it with Juttler＇s method. Examples suggest that our method is effective and is able to produce reconstruction surfaces of high quality.展开更多
基金supported by the State Oceanic Administration Young Marine Science Foundation (No. 2013201)the Shandong Provincial Key Laboratory of Marine Ecology and Environment & Disaster Prevention and Mitigation Foundation (No. 2012007)+1 种基金the Marine Public Foundation (No. 201005018)the North China Sea Branch Scientific Foundation (No. 2014B10)
文摘High Frequency(HF) radar current data is assimilated into a shelf sea circulation model based on optimal interpolation(OI) method. The purpose of this work is to develop a real-time computationally highly efficient assimilation method to improve the forecast of shelf current. Since the true state of the ocean is not known, the specification of background error covariance is arduous. Usually, it is assumed or calculated from an ensemble of model states and is kept in constant. In our method, the spatial covariances of model forecast errors are derived from differences between the adjacent model forecast fields, which serve as the forecast tendencies. The assumption behind this is that forecast errors can resemble forecast tendencies, since variances are large when fields change quickly and small when fields change slowly. The implementation of HF radar data assimilation is found to yield good information for analyses. After assimilation, the root-mean-square error of model decreases significantly. Besides, three assimilation runs with variational observation density are implemented. The comparison of them indicates that the pattern described by observations is much more important than the amount of observations. It is more useful to expand the scope of observations than to increase the spatial interval. From our tests, the spatial interval of observation can be 5 times bigger than that of model grid.
基金supported by the National Key Basic Research Project of China(No.2004CB318000)One Hundred Talent Project of the Chinese Academy of Sciences,the NSF of China(No.60225002,No.60533060)Doctorial Program of MOE of China and the 111 Project(No.B07033).
文摘Surface reconstruction from unorganized data points is a challenging problem in Computer Aided Design and Geometric Modeling. In this paper, we extend the mathematical model proposed by Juttler and Felis （Adv. Comput. Math., 17 （2002）, pp. 135-152） based on tensor product algebraic spline surfaces from fixed meshes to adaptive meshes. We start with a tensor product algebraic B-spline surface defined on an initial mesh to fit the given data based on an optimization approach. By measuring the fitting errors over each cell of the mesh, we recursively insert new knots in cells over which the errors are larger than some given threshold, and construct a new algebraic spline surface to better fit the given data locally. The algorithm terminates when the error over each cell is less than the threshold. We provide some examples to demonstrate our algorithm and compare it with Juttler＇s method. Examples suggest that our method is effective and is able to produce reconstruction surfaces of high quality.
