Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

Calibration of JASON-1 Altimeter over Lake Erie

This article describes absolute calibration results for both JASON-1 and TOPEX Side B (TSB) altimeters obtained at the Lake Erie calibration site, Marblehead, Ohio, USA. Using 15 overflights, the estimated JASON altimeter bias at Marblehead is 58 ± 38 mm, with an uncertainty of 19 mm based on detailed error analysis. Assuming that the TSB bias is negligible, relative bias estimates using both data from the TSB-JASON formation flight period and data from 48 water level gauges around the entire Great Lakes confirmed the Marblehead results. Global analyses using both the formation flight data and dual-satellite (TSB and JASON) crossovers yield a similar relative bias estimate of 146 ± 59 mm, which agrees well with open ocean absolute calibration results obtained at Harvest, Corsica, and Bass Strait (e.g., Watson et al. 2003). We find that there is a strong dependence of bias estimates on the choice of sea state bias (SSB) models. Results indicate that the invariant JASON instrument bias estimated oceanwide is 71 mm, with additional biases of 76 mm or 28 mm contributed by the choice of Collecte Localisation Satellites (CLS) SSB or Center for Space Research (CSR) SSB model, respectively. Similar analysis in the Great Lakes yields the invariant JASON instrument bias at 19 mm, with the SSB contributed biases at 58 mm or 13 mm, respectively. The reason for the discrepancy is currently unknown and warrants further investigation. Finally, comparison of the TOPEX/POSEIDON mission (1992-2002) data with the Great Lakes water level gauge measurements yields a negligible TOPEX altimeter drift of 0.1 mm/yr.     

The gravity field and GGOS

The gravity field of the earth is a natural element of the Global Geodetic Observing System (GGOS). Gravity field quantities are like spatial geodetic observations of potential very high accuracy, with measurements, currently at part-per-billion (ppb) accuracy, but gravity field quantities are also unique as they can be globally represented by harmonic functions (long-wavelength geopotential model primarily from satellite gravity field missions), or based on point sampling (airborne and in situ absolute and superconducting gravimetry). From a GGOS global perspective, one of the main challenges is to ensure the consistency of the global and regional geopotential and geoid models, and the temporal changes of the gravity field at large spatial scales. The International Gravity Field Service, an umbrella “level-2” IAG service (incorporating the International Gravity Bureau, International Geoid Service, International Center for Earth Tides, International Center for Global Earth models, and other future new services for, e.g., digital terrain models), would be a natural key element contributing to GGOS. Major parts of the work of the services would, however, remain complementary to the GGOS contributions, which focus on the long-wavelength components of the geopotential and its temporal variations, the consistent procedures for regional data processing in a unified vertical datum and Terrestrial Reference Frame, and the ensuring validations of long-wavelength gravity field data products.     

Satellite altimetry for measuring river stages in remote regions

The advent in satellite altimetry with the most accurate satellite radar altimeter since 1992 and its successive missions have enabled the routine global monitoring of water-level (or stage) for surface waters and changes in the quantities of dammed water reservoirs. However, satellite altimeter measurements typically have spatial resolution capable of observing only large water bodies, such as major lakes and rivers. This paper addresses the challenges of how to investigate water levels in medium (~?1 km in width) to small (~?100 m and narrower) rivers. Comparisons between the ENVISAT altimetry ICE-1 waveform retracking height and standard water-level measurements for multiple sections of Ohio River, Columbia River, and Red River of the North in the United States (US) reveals that the satellite altimetry measured water levels agree well with those observed at nearby US Geological Survey gaging stations over the 10-year period starting from 2002. The significant results include those obtained at Thompson, North Dakota (ND, correlation coefficient or R value of 0.76 between satellite and in situ water-level measurements) and Fargo, ND (R?=?0.74), where the stream channels of Red River are merely?~?50 m and ~?40 m wide, respectively, under normal climatic conditions. In addition, demonstrations of the approach over largely inaccessible portions of Tigris–Euphrates Rivers and Helmand River in the Middle East aided in understanding hydrology in these systems. This study demonstrates the ability of satellite radar altimetry to characterize rivers in these study regions which are much narrower than 100 m in width.     

The role of horizontal impulses of the faulting continental slope in generating the 26 December 2004 tsunami

For a long time, people have believed that the vertical displacement of seafloor due to undersea earthquakes is the primary cause of tsunami genesis. However, seismically-inverted seafloor deformation of the 2004 Sumatra–Andaman earthquake shows that the total vertical displacement is not enough to have generated the powerful Indian Ocean tsunami. Based on the seismically-inverted data and a three-dimensional ocean general circulation model (OGCM), we show that the momentum force, transferred by the horizontal impulses of the faulting continental slope in that earthquake, has accounted for two thirds of the satellite-observed tsunami height and generated kinetic energy 5 times larger than the potential energy due to the vertical displacement. The asymmetric tsunami pattern, recorded by tide-gauges showing leading-elevation waves toward Sri Lanka and leading-depression waves toward Thailand, is best explained by the horizontally-forced mechanism. The same mechanism has also explained the March 2005 Nias earthquake and tsunami data, suggesting that the horizontal motions of faulting have played more important roles in tsunami genesis than previously thought.     

Triangulated spherical splines for geopotential reconstruction

We present an alternate mathematical technique than contemporary spherical harmonics to approximate the geopotential based on triangulated spherical spline functions, which are smooth piecewise spherical harmonic polynomials over spherical triangulations. The new method is capable of multi-spatial resolution modeling and could thus enhance spatial resolutions for regional gravity field inversion using data from space gravimetry missions such as CHAMP, GRACE or GOCE. First, we propose to use the minimal energy spherical spline interpolation to find a good approximation of the geopotential at the orbital altitude of the satellite. Then we explain how to solve Laplace’s equation on the Earth’s exterior to compute a spherical spline to approximate the geopotential at the Earth’s surface. We propose a domain decomposition technique, which can compute an approximation of the minimal energy spherical spline interpolation on the orbital altitude and a multiple star technique to compute the spherical spline approximation by the collocation method. We prove that the spherical spline constructed by means of the domain decomposition technique converges to the minimal energy spline interpolation. We also prove that the modeled spline geopotential is continuous from the satellite altitude down to the Earth’s surface. We have implemented the two computational algorithms and applied them in a numerical experiment using simulated CHAMP geopotential observations computed at satellite altitude (450 km) assuming EGM96 (n _max = 90) is the truth model. We then validate our approach by comparing the computed geopotential values using the resulting spherical spline model down to the Earth’s surface, with the truth EGM96 values over several study regions. Our numerical evidence demonstrates that the algorithms produce a viable alternative of regional gravity field solution potentially exploiting the full accuracy of data from space gravimetry missions. The major advantage of our method is that it allows us to compute the geopotential over the regions of interest as well as enhancing the spatial resolution commensurable with the characteristics of satellite coverage, which could not be done using a global spherical harmonic representation. The results in this paper are based on the research supported by the National Science Foundation under the grant no. 0327577.     

Effects of economic recession and local weather on climate change attitudes

How are large companies responding to the challenges of reducing their greenhouse gas (GHG) emissions? An analysis of the published climate change policies and performance of 125 large European companies is presented. The results suggest that most large European companies have now developed the management systems and processes necessary for them to effectively manage their GHG emissions and related business risks. However, there is a significant disconnect between the targets that companies set for themselves and the more ambitious targets being set by the European Union (which has committed to a 20% reduction in its emissions by 2020 against a 1990 baseline). Of the companies surveyed, just over one-third had stabilized or reduced their total GHG emissions over the period 2002–2007, and fewer than one-third expected their emissions to stabilize or reduce in the coming years. The relationship between the quality of corporate policies and performance outcomes (in terms of GHG emissions) suggests that while companies with stronger policies are likely to have relatively better performance, only a minority of those companies with the highest-quality policies are committing to absolute reductions in their GHG emissions.     

利用海洋重力场变化分析洋底板块运动

洋底板块运动是地球动力学和全球变化研究的重要内容.本文根据质量迁移与地球外部重力场变化的对应关系,利用不同时期测高资料推算的1995—2019全球海洋重力场变化结果,反演分析全球洋底板块运动特征.研究表明,板块汇聚边界、板块内无震海岭、海山群、断裂带等区域重力异常变化显著,而在板块离散边界无明显变化趋势;西南印度洋中脊、大西洋中脊、中印度洋中脊等地区重力异常垂直梯度变化显著,且在西太平洋俯冲带、部分海岭区域也存在明显变化,其空间分布与地形基本吻合.海洋重力场变化整体上准确反映了全球洋底板块构造运动.相较于重力异常变化反演结果,重力垂直梯度的变化能够更为准确地反映洋底板块运动特征,特别是在洋中脊区域,扩张速率越小,垂直重力梯度变化越显著.此外,详细讨论了测高海洋重力场不确定因素对洋底板块运动分析结果的影响,海面坡度改正是主要因素之一.     

Regional Temperature Forecast for the Next Day in Hong Kong

For a century or so, the Hong Kong Observatory (HKO) has been providing temperature forecast for the whole of Hong Kong with the HKO Headquarters as the reference location. In recent decades, due to spreading of population from the main urban center to satellite towns, there is an increasing demand for regional temperature forecasts. To support such provision, the HKO has developed a regression model to provide objective guidance to forecasters in formulating forecasts of maximum and minimum temperatures for the next day at various locations in Hong Kong. In this paper, the regression model is presented, together with the assessment of its performance. Based on the verification of one year of forecasts, it is found that the root mean square errors (RMSEs) of maximum (minimum) temperature forecasts are from about 1.3 to 2.1 (1.1 to 1.4) degrees, respectively. The regression model is shown to have generally out-performed the operational regional spectral model then operated by HKO. Regional temperature forecast methods of other meteorological or research centers are also surveyed. Equipped with the regression model, the HKO has launched an online regional temperature forecast service for the next day in Hong Kong since March 2008.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.