In the 1940s, KIRKENDALL showed that diffusion in binary solid solutions cannot be described by only one diffusion coefficient. Rather, one has to consider the diffusivity of both species. His findings changed the tre...In the 1940s, KIRKENDALL showed that diffusion in binary solid solutions cannot be described by only one diffusion coefficient. Rather, one has to consider the diffusivity of both species. His findings changed the treatment of diffusion data and the theory of diffusion itself. A diffusion-based framework was successfully employed to explain the behaviour of the Kirkendall plane. Nonetheless, the complexity of a multiphase diffusion zone and the morphological evolution during interdiffusion requires a physico-chemical approach. The interactions in binary and more complex systems are key issues from both the fundamental and technological points of view. This paper reviews the Kirkendall effect from the circumstances of its discovery to recent developments in its understanding, with broad applicability in materials science and engineering.展开更多
We study the random motion of a charged test particle with a normal classical constant velocity in a spacetime with a perfectly reflecting plane boundary and calculate both the velocity and position dispersions of the...We study the random motion of a charged test particle with a normal classical constant velocity in a spacetime with a perfectly reflecting plane boundary and calculate both the velocity and position dispersions of the test particle. Our results show that the dispersions in the normal direction are weakened while those in the parallel directions are strengthened as compared to the classical static case when the test particle classically moves away from the boundary. However, if the classical motion reverses its direction, then the dispersions in the normal direction are reinforced while those in the parallel directions get weakened.展开更多
Based on sounding data from 1975 to 2005 and TM/ETM+ remote sensing images in 1989, 2001 and 2007, the climate changes in Harbin City, Northeast China in recent 30 years were analyzed and forecasted. Results show that...Based on sounding data from 1975 to 2005 and TM/ETM+ remote sensing images in 1989, 2001 and 2007, the climate changes in Harbin City, Northeast China in recent 30 years were analyzed and forecasted. Results show that in the lower troposphere the meridional wind speed and mean annual wind speed decrease, and in the lower stratosphere the temperature decreases while the meridional wind speed increases significantly. In the study area, the climate is becoming warmer and wetter in the middle lower troposphere. The expansion of urban area has great effects on the surface air temperature and the wind speed, leading to the increase of the surface air temperature, the decrease of the surface wind speed, and the increase of the area of urban high temperature zone. The quantitative equations have been established among the surface air temperature, the carbon dioxide (CO2) concentration and the specific humidity (the water vapor content). It is predicted that the future increasing rate of the surface air temperature is 0.85℃/10yr if emission concentration of CO2 remains unchanged; if emission concentration of CO2 decreases to 75%, 50% and 25%, respectively, the surface air temperature will increase 0.65℃/10yr, 0.46℃/10yr and 0.27℃/10yr, respectively. The rise of the surface air temperature in the study area is higher than that of the global mean temperature forecasted by IPCC.展开更多
Ships which have large structures above water surface,such as pure car carriers(PCCs) and container vessels,have large speed reduction by wind pressure.In the present study,the running speed of a large PCC with two or...Ships which have large structures above water surface,such as pure car carriers(PCCs) and container vessels,have large speed reduction by wind pressure.In the present study,the running speed of a large PCC with two or more sails for using wind power is simulated.The simulated results demonstrate that the ship can keep a constant service speed even in winds of 20m/s except head and bow winds.This sail system can shorten annual average navigation time by about 4 hours per voyage.展开更多
The formation of shear bands for time and length scales appropriate for deformation processes in the upper Lithosphere is investigated in plane strain finite element simulations under predominantly uniaxial extension ...The formation of shear bands for time and length scales appropriate for deformation processes in the upper Lithosphere is investigated in plane strain finite element simulations under predominantly uniaxial extension and compression, respectively. The direction of gravity is assumed orthogonal to the extension/compression axis. Mathematically, the formation of shear zones may be explained as a consequence of changes in the type of the governing model equations. Such changes or bifurcations depend strongly on the details of the constitutive relationships such as strain softening, thermal or chemical effects, associated or non-associated--coaxial or non-coaxial flow rules. Here we focus on strain softening and coaxial and non-coaxial flow rules. In the simulations, we consider an initially rectangular domain with the dimensions Lo, Ho in the horizontal, vertical directions, respectively. The domain is extended or compressed by prescribing a uniform, horizontal velocity field along one of the vertical boundaries while keeping the opposite boundary fixed. An important global descriptor of the deformation process is the relationship between the horizontal stress resultant (average horizontal stress) and the strain ln(L/Lo), where L is the deformed length of the domain. The main goal of this paper is to investigate key factors influencing the phenomenology of the localization process such as flow rule, coaxial, non-coaxial and strain softening. Different origins of the mesh sensitivity of deformations involving localization are also investigated.展开更多
The sea level derived from TOPEX/Poseidon (T/P) altimetry data shows prominent long term trend and inter-annual variability. The global mean sea level rising rate during 1993-2003 was 2.9mm a^-1. The T/P sea level t...The sea level derived from TOPEX/Poseidon (T/P) altimetry data shows prominent long term trend and inter-annual variability. The global mean sea level rising rate during 1993-2003 was 2.9mm a^-1. The T/P sea level trend maps the geographical variability. In the Northern Hemisphere (15°-64°N), the sea level rise is very fast at the mid-latitude (20°-40°N) but much slower at the high-latitude, for example, only 0.5 mm a^-1 in the latitude band 40°-50°N. In the Southern Hemisphere, the sea level shows high rising rate both in mid-latitude and high-latitude areas, for example, 5.1 mm a^-1 in the band 40°- 50°S. The global thermosteric sea level (TSL) derived from Ishii temperature data was rising during 1993-2003 at a rate of 1.2 mm a^-1 and accounted for more than 40% of the global T/P sea level rise. The contributions of the TSL distribution are not spatially uniform; for instance, the percentage is 67% for the Northern Hemisphere and only 29% for the Southern Hemisphere (15°-64°S) and the maximum thermosteric contribution appears in the Pacific Ocean, which contributes more than 60% of the global TSL. The sea level change trend in tropical ocean is mainly caused by the thermosteric effect, which is different from the case of seasonal variability in this area. The TSL variability dominates the T/P sea level rise in the North Atlantic, but it is small in other areas, and shows negative trend at the high-latitude area (40°-60°N, and 50°-60°S). The global TSL during 1945-2003 showed obvious rising trend with the rate of about 0.3 mm a-l and striking inter-annual and decadal variability with period of 20 years. In the past 60 years, the Atlantic TSL was rising continuously and remarkably, contributing 38% to the global TSL rising. The TSL in the Pacific and Indian Ocean rose with significant in- ter-annual and decadal variability. The first EOF mode of the global TSL from Ishii temperature data was the ENSO mode in which the time series of the first mode showed steady rising trend. Among the three oceans, the first mode of the Pacific TSL presented the ENSO mode; there was relatively steady rising trend in the Atlantic Ocean, and no dominant mode in the Indian Ocean.展开更多
Pearl millet at various moisture content (6.2, 9.4 and 12.3%, d.b.) and feed rates (3, 6 and 9 kg/h) was ground using hammer mill and its physical properties viz. arithmetic mean diameter, geometric mean diameter,...Pearl millet at various moisture content (6.2, 9.4 and 12.3%, d.b.) and feed rates (3, 6 and 9 kg/h) was ground using hammer mill and its physical properties viz. arithmetic mean diameter, geometric mean diameter, thousand grain weight, aspect ratio, specific surface area, surface area and bulk density were studied. Sieve analysis results showed that the increase in moisture content produced more medium sized particles with decreased percent weight retained in pan. Bond's work index, Kick's constant and average particle size were increased with the decrease in total surface area at higher moisture levels. The highest energy (2.34 KWh/kg) was consumed for 12.3% moisture content. Various grinding characteristics were significantly affected by moisture content and feed rate either individually or in combination and correlated in terms of Bond's work index, Kick's constant, total surface area, average particle size, effectiveness of milling and bulk density. Milling loss was found to be higher at lower moisture level and decreased with the increase of moisture content as well as feed rate. The loose and compact bulk density was ranged between 46.8-199.5 kg m^-3and 53.5-254.1 kg m^-3, respectively among the entire sieve fractions. Water absorption capacity increased with the decrease in particle size,展开更多
Rupture directivity effect causes spatial variation in strong ground motion parameters. It causes difference between the strike- normal (V.) and strike-parallel (Vp) components of horizontal ground motion amplitud...Rupture directivity effect causes spatial variation in strong ground motion parameters. It causes difference between the strike- normal (V.) and strike-parallel (Vp) components of horizontal ground motion amplitudes. These variations become significant for strong ground motion velocity and the authors have developed a modification to define directivity effect factor to account for the effect of rupture directivity in empirical velocity attenuation relations which are based on modeling Silakhor earthquake, using finite element method by ANSYS. The ground motion parameters that are modified include ratio of Vn/Vp component of horizontal velocity and Vn component to average horizontal velocity (V). The ratio of Vn to Vp is large in both the forward directivity direction, where velocity is larger, and in the backward directivity direction, where velocity is smaller. Therefore the authors expected that the Vn/Vp was mainly controlled by directivity angle. Also the variation of fault normal velocity to average horizontal velocity ratio by directivity angle (0) is defined from earthquake modeling. It shows Vn/V is controlled by directivity angle, distance between the site, epicenter and rupture length. This ratio has the same trend in Silakhor earthquake strong ground velocity data. In this paper the equation for Vn/Vp variations by directivity angle is recommended. The authors used Somervill et al. (1997) directivity model parameters as (R/L) cos2 ~ to define directivity effect on Vn/V ratio and therefore directivity factor is determined to account in near field empirical strong ground velocity attenuation relationships.展开更多
On the basis of elastic rebound theory,using the horizontal velocity field of the Chinese mainland calculated from GPS data during three observation periods from 1999 to 2007,the velocity components that are parallel ...On the basis of elastic rebound theory,using the horizontal velocity field of the Chinese mainland calculated from GPS data during three observation periods from 1999 to 2007,the velocity components that are parallel and plumb to the fault zone are calculated respectively for different periods,and then relative ground movements of two sides of the fault zone are analyzed with power function fitting and graphics. The results show that the relative ground movement shows right-lateral shear deformation before the Wenchuan MS8. 0 earthquake,and at the same time the movement was hindered by the Longmenshan fault zone. Thus,this result has positive significance for distinguishing the elastic strain energy accumulation and deformation anomaly in an earthquake preparation process,and for conducting further research on earthquake prediction.展开更多
We present a bulk micromachined in-plane capacitive accelerometer fabricated with an improved process flow,by etching only one-fifth of the wafer thickness at the back of the silicon while forming the bar-structure el...We present a bulk micromachined in-plane capacitive accelerometer fabricated with an improved process flow,by etching only one-fifth of the wafer thickness at the back of the silicon while forming the bar-structure electrode for the sensing capacitor.The improved flow greatly lowers the footing effect during deep reactive ion etching(DRIE),and increases the proof mass by 54% compared to the traditional way,resulting in both improved device quality and a higher yield rate.Acceleration in the X direction is sensed capacitively by varying the overlapped area of a differential capacitor pair,which eliminates the nonlinear behavior by fixing the parallel-plate gap.The damping coefficient of the sensing motion is low due to the slide-film damping.A large proof mass is made using DRIE,which also ensures that dimensions of the spring beams in the Y and Z directions can be made large to lower cross axis coupling and increase the pull-in voltage.The theoretical Brownian noise floor is 0.47 μg/Hz1/2 at room temperature and atmospheric pressure.The tested frequency response of a prototype complies with the low damping design scheme.Output data for input acceleration from ?1 g to 1 g are recorded by a digital multimeter and show very good linearity.The tested random bias of the prototype is 130 μg at an averaging time of around 6 s.展开更多
A modified Gauss-Markov model with weighted constraints was constructed by combining satellite altimeter and tide gauge records. Vertical motion rates of nine tide gauge stations around the Bohal Sea and Yellow Sea ar...A modified Gauss-Markov model with weighted constraints was constructed by combining satellite altimeter and tide gauge records. Vertical motion rates of nine tide gauge stations around the Bohal Sea and Yellow Sea are estimated. This is the first time systematic estimates have been derived in this region. Downward trends were seen at the six tide gauge stations located at Tanggu, Longkou, Laohutan, Bayuquan, Xiaochangshan, and Yantai; with vertical motion rates of-1.82±0.50, -1.65±0.46, -0.88±0.42, -0.58±0.62, -0.13±0.43, and -0.01±0.43 mm/yr, respectively. Upward trends were seen at the three tide gauge stations located at Qinhuangdao, Huludao and Chengshantou; with vertical motion rates of 1.12±0.46, 0.55±0.49 and 0.26±0.44 mm/yr, respectively. There was significant subsidence in Tanggu and Longkou, and a rising trend in Qinhuangdao. According to our results, the rate of sea level rise calculated from these tide gauge records can be improved using a more accurate measurement of the land elevation accounting for lifting or subsidence. The model derived can be used to estimate vertical motions of tide gauge stations, and can be widely applied to revise the benchmark levels of tide gauges.展开更多
文摘In the 1940s, KIRKENDALL showed that diffusion in binary solid solutions cannot be described by only one diffusion coefficient. Rather, one has to consider the diffusivity of both species. His findings changed the treatment of diffusion data and the theory of diffusion itself. A diffusion-based framework was successfully employed to explain the behaviour of the Kirkendall plane. Nonetheless, the complexity of a multiphase diffusion zone and the morphological evolution during interdiffusion requires a physico-chemical approach. The interactions in binary and more complex systems are key issues from both the fundamental and technological points of view. This paper reviews the Kirkendall effect from the circumstances of its discovery to recent developments in its understanding, with broad applicability in materials science and engineering.
基金The project supported by the National Natural Science Foundation of China under Grant No. 10575035
文摘We study the random motion of a charged test particle with a normal classical constant velocity in a spacetime with a perfectly reflecting plane boundary and calculate both the velocity and position dispersions of the test particle. Our results show that the dispersions in the normal direction are weakened while those in the parallel directions are strengthened as compared to the classical static case when the test particle classically moves away from the boundary. However, if the classical motion reverses its direction, then the dispersions in the normal direction are reinforced while those in the parallel directions get weakened.
基金Under the auspices of Major State Basic Research Development Program of China (No. 2010CB950900)Technology Innovation Program of Harbin City (No. 2007RFXXS029)
文摘Based on sounding data from 1975 to 2005 and TM/ETM+ remote sensing images in 1989, 2001 and 2007, the climate changes in Harbin City, Northeast China in recent 30 years were analyzed and forecasted. Results show that in the lower troposphere the meridional wind speed and mean annual wind speed decrease, and in the lower stratosphere the temperature decreases while the meridional wind speed increases significantly. In the study area, the climate is becoming warmer and wetter in the middle lower troposphere. The expansion of urban area has great effects on the surface air temperature and the wind speed, leading to the increase of the surface air temperature, the decrease of the surface wind speed, and the increase of the area of urban high temperature zone. The quantitative equations have been established among the surface air temperature, the carbon dioxide (CO2) concentration and the specific humidity (the water vapor content). It is predicted that the future increasing rate of the surface air temperature is 0.85℃/10yr if emission concentration of CO2 remains unchanged; if emission concentration of CO2 decreases to 75%, 50% and 25%, respectively, the surface air temperature will increase 0.65℃/10yr, 0.46℃/10yr and 0.27℃/10yr, respectively. The rise of the surface air temperature in the study area is higher than that of the global mean temperature forecasted by IPCC.
文摘Ships which have large structures above water surface,such as pure car carriers(PCCs) and container vessels,have large speed reduction by wind pressure.In the present study,the running speed of a large PCC with two or more sails for using wind power is simulated.The simulated results demonstrate that the ship can keep a constant service speed even in winds of 20m/s except head and bow winds.This sail system can shorten annual average navigation time by about 4 hours per voyage.
文摘The formation of shear bands for time and length scales appropriate for deformation processes in the upper Lithosphere is investigated in plane strain finite element simulations under predominantly uniaxial extension and compression, respectively. The direction of gravity is assumed orthogonal to the extension/compression axis. Mathematically, the formation of shear zones may be explained as a consequence of changes in the type of the governing model equations. Such changes or bifurcations depend strongly on the details of the constitutive relationships such as strain softening, thermal or chemical effects, associated or non-associated--coaxial or non-coaxial flow rules. Here we focus on strain softening and coaxial and non-coaxial flow rules. In the simulations, we consider an initially rectangular domain with the dimensions Lo, Ho in the horizontal, vertical directions, respectively. The domain is extended or compressed by prescribing a uniform, horizontal velocity field along one of the vertical boundaries while keeping the opposite boundary fixed. An important global descriptor of the deformation process is the relationship between the horizontal stress resultant (average horizontal stress) and the strain ln(L/Lo), where L is the deformed length of the domain. The main goal of this paper is to investigate key factors influencing the phenomenology of the localization process such as flow rule, coaxial, non-coaxial and strain softening. Different origins of the mesh sensitivity of deformations involving localization are also investigated.
基金supported by the National Basic Research Program of China (No 2007CB411807)the NSFC project (Nos 40976006 and 40906002)+1 种基金the National Key Technology R&D Program (No 2007BAC03A06-06)the project of Key Laboratory of Coastal Disasters and Defence (No 200802)
文摘The sea level derived from TOPEX/Poseidon (T/P) altimetry data shows prominent long term trend and inter-annual variability. The global mean sea level rising rate during 1993-2003 was 2.9mm a^-1. The T/P sea level trend maps the geographical variability. In the Northern Hemisphere (15°-64°N), the sea level rise is very fast at the mid-latitude (20°-40°N) but much slower at the high-latitude, for example, only 0.5 mm a^-1 in the latitude band 40°-50°N. In the Southern Hemisphere, the sea level shows high rising rate both in mid-latitude and high-latitude areas, for example, 5.1 mm a^-1 in the band 40°- 50°S. The global thermosteric sea level (TSL) derived from Ishii temperature data was rising during 1993-2003 at a rate of 1.2 mm a^-1 and accounted for more than 40% of the global T/P sea level rise. The contributions of the TSL distribution are not spatially uniform; for instance, the percentage is 67% for the Northern Hemisphere and only 29% for the Southern Hemisphere (15°-64°S) and the maximum thermosteric contribution appears in the Pacific Ocean, which contributes more than 60% of the global TSL. The sea level change trend in tropical ocean is mainly caused by the thermosteric effect, which is different from the case of seasonal variability in this area. The TSL variability dominates the T/P sea level rise in the North Atlantic, but it is small in other areas, and shows negative trend at the high-latitude area (40°-60°N, and 50°-60°S). The global TSL during 1945-2003 showed obvious rising trend with the rate of about 0.3 mm a-l and striking inter-annual and decadal variability with period of 20 years. In the past 60 years, the Atlantic TSL was rising continuously and remarkably, contributing 38% to the global TSL rising. The TSL in the Pacific and Indian Ocean rose with significant in- ter-annual and decadal variability. The first EOF mode of the global TSL from Ishii temperature data was the ENSO mode in which the time series of the first mode showed steady rising trend. Among the three oceans, the first mode of the Pacific TSL presented the ENSO mode; there was relatively steady rising trend in the Atlantic Ocean, and no dominant mode in the Indian Ocean.
文摘Pearl millet at various moisture content (6.2, 9.4 and 12.3%, d.b.) and feed rates (3, 6 and 9 kg/h) was ground using hammer mill and its physical properties viz. arithmetic mean diameter, geometric mean diameter, thousand grain weight, aspect ratio, specific surface area, surface area and bulk density were studied. Sieve analysis results showed that the increase in moisture content produced more medium sized particles with decreased percent weight retained in pan. Bond's work index, Kick's constant and average particle size were increased with the decrease in total surface area at higher moisture levels. The highest energy (2.34 KWh/kg) was consumed for 12.3% moisture content. Various grinding characteristics were significantly affected by moisture content and feed rate either individually or in combination and correlated in terms of Bond's work index, Kick's constant, total surface area, average particle size, effectiveness of milling and bulk density. Milling loss was found to be higher at lower moisture level and decreased with the increase of moisture content as well as feed rate. The loose and compact bulk density was ranged between 46.8-199.5 kg m^-3and 53.5-254.1 kg m^-3, respectively among the entire sieve fractions. Water absorption capacity increased with the decrease in particle size,
文摘Rupture directivity effect causes spatial variation in strong ground motion parameters. It causes difference between the strike- normal (V.) and strike-parallel (Vp) components of horizontal ground motion amplitudes. These variations become significant for strong ground motion velocity and the authors have developed a modification to define directivity effect factor to account for the effect of rupture directivity in empirical velocity attenuation relations which are based on modeling Silakhor earthquake, using finite element method by ANSYS. The ground motion parameters that are modified include ratio of Vn/Vp component of horizontal velocity and Vn component to average horizontal velocity (V). The ratio of Vn to Vp is large in both the forward directivity direction, where velocity is larger, and in the backward directivity direction, where velocity is smaller. Therefore the authors expected that the Vn/Vp was mainly controlled by directivity angle. Also the variation of fault normal velocity to average horizontal velocity ratio by directivity angle (0) is defined from earthquake modeling. It shows Vn/V is controlled by directivity angle, distance between the site, epicenter and rupture length. This ratio has the same trend in Silakhor earthquake strong ground velocity data. In this paper the equation for Vn/Vp variations by directivity angle is recommended. The authors used Somervill et al. (1997) directivity model parameters as (R/L) cos2 ~ to define directivity effect on Vn/V ratio and therefore directivity factor is determined to account in near field empirical strong ground velocity attenuation relationships.
基金funded by the Special Subject of the National Key Technology R&D Program for the 11th "Five-year Plan" of China(2006BAC01B02-02-02)
文摘On the basis of elastic rebound theory,using the horizontal velocity field of the Chinese mainland calculated from GPS data during three observation periods from 1999 to 2007,the velocity components that are parallel and plumb to the fault zone are calculated respectively for different periods,and then relative ground movements of two sides of the fault zone are analyzed with power function fitting and graphics. The results show that the relative ground movement shows right-lateral shear deformation before the Wenchuan MS8. 0 earthquake,and at the same time the movement was hindered by the Longmenshan fault zone. Thus,this result has positive significance for distinguishing the elastic strain energy accumulation and deformation anomaly in an earthquake preparation process,and for conducting further research on earthquake prediction.
基金Project(No.NCET-06-0514)supported by the Program for New Century Excellent Talents in University,China
文摘We present a bulk micromachined in-plane capacitive accelerometer fabricated with an improved process flow,by etching only one-fifth of the wafer thickness at the back of the silicon while forming the bar-structure electrode for the sensing capacitor.The improved flow greatly lowers the footing effect during deep reactive ion etching(DRIE),and increases the proof mass by 54% compared to the traditional way,resulting in both improved device quality and a higher yield rate.Acceleration in the X direction is sensed capacitively by varying the overlapped area of a differential capacitor pair,which eliminates the nonlinear behavior by fixing the parallel-plate gap.The damping coefficient of the sensing motion is low due to the slide-film damping.A large proof mass is made using DRIE,which also ensures that dimensions of the spring beams in the Y and Z directions can be made large to lower cross axis coupling and increase the pull-in voltage.The theoretical Brownian noise floor is 0.47 μg/Hz1/2 at room temperature and atmospheric pressure.The tested frequency response of a prototype complies with the low damping design scheme.Output data for input acceleration from ?1 g to 1 g are recorded by a digital multimeter and show very good linearity.The tested random bias of the prototype is 130 μg at an averaging time of around 6 s.
基金supported by the National Natural Science Foundation of China(Grant Nos.4140603241206021&41376014)+4 种基金the National Basic Research Program of China(Grant No.2012CB955601)the National Key Technology R&D Program of China(Grant No.2014BAB12B02)the Key Technology R&D Program of Tianjin(Grant No.14ZCZDSF00012)the Open Fund of the State Key Laboratory of Satellite Ocean Environment Dynamics at the Second Institute of OceanographySOA(Grant No.SOED1305)
文摘A modified Gauss-Markov model with weighted constraints was constructed by combining satellite altimeter and tide gauge records. Vertical motion rates of nine tide gauge stations around the Bohal Sea and Yellow Sea are estimated. This is the first time systematic estimates have been derived in this region. Downward trends were seen at the six tide gauge stations located at Tanggu, Longkou, Laohutan, Bayuquan, Xiaochangshan, and Yantai; with vertical motion rates of-1.82±0.50, -1.65±0.46, -0.88±0.42, -0.58±0.62, -0.13±0.43, and -0.01±0.43 mm/yr, respectively. Upward trends were seen at the three tide gauge stations located at Qinhuangdao, Huludao and Chengshantou; with vertical motion rates of 1.12±0.46, 0.55±0.49 and 0.26±0.44 mm/yr, respectively. There was significant subsidence in Tanggu and Longkou, and a rising trend in Qinhuangdao. According to our results, the rate of sea level rise calculated from these tide gauge records can be improved using a more accurate measurement of the land elevation accounting for lifting or subsidence. The model derived can be used to estimate vertical motions of tide gauge stations, and can be widely applied to revise the benchmark levels of tide gauges.