Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Developing a Marine Information System by Integrating Existing Ocean Models Using Object-Oriented Technology

Information systems developed for different applications within the environmental domain have common characteristics, which can potentially be abstracted for sharing and reuse of design and software modules. This article presents an approach to designing for reuse by abstracting commonalities in the design of a Marine Information System (MIS), facilitating data management in a prototype operational monitoring, forecasting, and management system for the North Atlantic and the Nordic Seas. A detailed study of the requirements and data analysis was carried out, and Object-Oriented Technology (OOT) is employed to encapsulate abstractions and to promote reuse of code and design. This article identifies the Object-Oriented Frameworks (OOFW) required to build the MIS. It also provides guidelines to environmental scientists for restructuring legacy software and employing modern programming techniques.     

The impact of the M2 internal tide on data-assimilative model estimates of the surface tide

We use a synthetic data experiment to assess the accuracy of ocean tides estimated from satellite altimetry data, with emphasis on the impact of the phase-locked internal tide, which has a surface expression of several centimeters near its sites of genesis. Previous tidal estimates have regarded this signal as a random measurement error; however, it is deterministic and not scale-separated from the barotropic (surface) tide around complex bathymetric features. The synthetic data experiments show that the internal tide has a negligible impact on the barotropic tidal fields inferred under these circumstances, and the barotropic dissipation (a quadratic functional of the tidal fields) is in good agreement with the energetics of the three-dimensional regional primitive equations model which is the source of the synthetic data.     

A Reconstruction of Observed Profiles in the Sea East of Japan Using Vertical Coupled Temperature-Salinity EOF Modes

It is important to estimate hard-to-observe parameters in the ocean interior from easy-to-observe parameters. This study therefore demostrates a reconstruction of observed temperature and salinity profiles of the sea east of Japan (30°≈40°N, 140°≈150°E). The reconstruction was done by estimating suboptimal state from several values of the observed profiles and/or sea surface dynamic height (SDH) calculated from the profiles. The estimation used a variational method with vertical coupled temperature-salinity empirical orthogonal function (EOF) modes. Profiles of temperature and salinity in the subtropical region are effectively reconstructed from in situ temperature profile data, or sea surface temperature (SST) and SDH. For example, the analyzed temperature field from SST and SDH has an accuracy to within 1°C in the subtropical region. Salinity in the sea north of Kuroshio, however, is difficult to estimate because of its complex variability which is less correlated with temperature than in the subtropical region. Sea surface salinity is useful to estimate the subsurface structure. We also show the possibility that the estimation is improved by considering nonlinearity in the equation calculating SDH from temperature and salinity analysis values in order to examine the misfit between analysis and observation. Analysis using TOPEX/POSEIDON altimetry data instead of SDH was also performed. This revised version was published online in July 2006 with corrections to the Cover Date.     

哈尔滨市用地查询分析系统的设计与实现

随着3S技术的不断发展,客户要求的不断提高,原始的服务方式已经不能满足客户对多元信息获取的需要,在这样的背景之下,哈尔滨市用地查询分析系统应运而生,它的成功研制有效解决了客户需求,并对其他项目数据整合的方式与方法提供了有益的解决方案。     

GPS-RTK联合全站仪的野外数据采集

对于地形测绘和其它环境条件复杂的地区,用常规方法进行数据采集比较困难,传统的测绘方法越来越不能满足测量效率和成本的要求,但采用GPS RTK与全站仪联合作业,使RTK与全站仪优势互补,在保证测量精度的前提下,可大幅度提高测量速度和效率。     

GPS CORS系统中NTRIP传输技术的应用研究

NTRIP是一种播发实时差分定位改正信息的IP协议。本文讨论了NTRIP的组成、工作原理,并给出一个应用实例。     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.