基于实时精密单点定位技术的暴雨短临预报

提出了一种将实时精密单点定位(Precise Point Positioning, PPP)技术用于暴雨短临预报的新方法.该方法首先基于GPS连续运行参考站网(Continuously Operating Reference Stations, CORS)实时估计的精密卫星钟差完成PPP解算,再以实时获取的对流层延迟(Zenith Tropospheric Delay, ZTD)及其增量变化为依据进行暴雨短临预报.研究结果表明:一般雷暴天气来临之前的2~6 h,ZTD增量表现为先后突破±5 mm/5 min,且后续记录到的实际降水量大小与ZTD维持在高水平阶段的时间长短有较好的对应关系;就热带气旋而言,在强风作用下,ZTD增量变化表现的异常活跃和复杂,规律性较弱,但对短临预报强降雨仍有一定的指示作用.     

Signature of atmospheric oscillations in GPS-measured tropospheric delay

Ground-based GPS finds potential applications in many atmospheric studies such as spatial distribution of columnar water vapor as well as the tidal oscillations in the atmosphere. The zenith tropospheric delay (ZTD) derived from GPS data at two Indian IGS stations are used to establish its potential for studying the atmospheric tidal, intra-seasonal and planetary oscillations. The major tidal oscillations observed in ZTD data are diurnal, semi-diurnal and their harmonics. Prominent intra-seasonal oscillations observed in ZTD are reported for the first time in this context. These intra-seasonal oscillations are Madden–Julian Oscillation (30–70 days, 60–90 days, 100–120 days) and planetary waves (like 27, 16 and 5–10 days periodicities). Quantification of these periodicities will provide a useful handle to improve the empirical models employed in the estimation of tropospheric delay.     

Assessment of GPS data for meteorological applications over Africa: Study of error sources and analysis of positioning accuracy

The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L₂ data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L₂ data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.     

顾及共性误差的南极半岛地区连续GPS站坐标时间序列分析

共性误差是区域连续GPS网络中存在的一种与时空相关的主要误差源.对区域网络使用空间滤波能够有效减小这种误差并改善坐标时间序列的精度.本文使用GAMIT/GLOBK 10.5软件解算南极半岛区域11个GPS测站2010—2014的数据,各测站独立估计的单天解组成坐标时间序列.利用堆栈法(stacking)、主分量分析法(principal component analysis,简称PCA)和Karhunen-Loeve展开法(KLE)对其进行空间滤波,在利用PCA和KLE方法进行空间滤波时,使用前两个主分量作为共有模式分量.结果表明:三种方法均能有效提取南极半岛区域共性误差,PCA方法的滤波效果明显优于stacking方法,略微优于KLE方法;空间滤波能够有效减小残差时间序列的振幅、功率和RMS;对坐标时间序列进行空间滤波能够有效降低其线性项和周期项误差,从而提高线性项和周期项估计的精度.共性误差的谱分析结果显示在U方向上存在9.4天、13.7天等短周期的信号,推测共性误差组成成分中可能含有与海潮相关的误差源.     

The influence of seasonal signals on the estimation of the tectonic motion in short continuous GPS time-series

Most GPS time-series exhibit a seasonal signal that can have an amplitude of a few millimetres. This seasonal signal can be removed by fitting an extra sinusoidal signal with a period of one year to the GPS data during the estimation of the linear trend.However, Blewitt and Lavallée (2002) showed that including an annual signal in the estimation process still can give a larger linear trend error than the trend error estimated from data from which the annual signal has been removed by other means. They assumed that the GPS data only contained white noise and we extend their result to the case of power-law plus white noise which is known to exist in most GPS observations. For the GPS stations CASC, LAGO, PDEL and TETN the difference in trend error between having or not having an annual signal in the data is around ten times larger when a power-law plus white noise model is used instead of a pure white noise model. Next, our methodology can be used to estimate for any station how much the accuracy of the linear trend will improve when one tries to subtract the annual signal from the GPS time-series by using a physical model.Finally, we demonstrate that for short time-series the trend error is more influenced by the fact that the noise properties also need to be estimated from the data. This causes on average an underestimation of the trend error.     

Temporal trends in δ18O composition of precipitation in Germany: insights from time series modelling and trend analysis

Temporal and spatial variations of stable oxygen (¹⁸O) and hydrogen (²H) isotope measurements in precipitation act as important proxies for changing hydro‐meteorological and regional and global climate patterns. Temporal trends in time series of the stable isotope composition in precipitation were rarely observed, and they are poorly understood. These might be a result of a lack of proper trend detection tools and effort for exploring trend processes. Here, we investigate temporal trends of δ¹⁸O in precipitation at 17 observation stations in Germany between 1978 and 2009. We test if significant trends in the isotope time series from different models can be observed. Mann–Kendall trend tests are applied on the isotope series, using general multiplicative seasonal autoregressive integrate moving average (ARIMA) models, which account for first and higher order serial correlations. Effects of temperature, precipitation, and geographic parameters on isotope trends are also investigated in the proposed models. To benchmark our proposed approach, the ARIMA results are compared with a trend‐free pre‐whitening procedure, the state of the art method for removing the first order autocorrelation in environmental trend studies. Moreover, we further explore whether higher order serial correlations in isotope series affects our trend results. Overall, three out of the 17 stations show significant changes when higher order autocorrelation are adjusted, and four show a significant trend when temperature and precipitation effects are considered. The significant trends in the isotope time series generally occur only at low elevation stations. Higher order autoregressive processes are shown to be important in the isotope time series analysis. Results suggest that the widely used trend analysis with only the first order autocorrelation adjustment may not adequately take account of the high order autocorrelated processes in the stable isotope series. The investigated time series analysis method including higher autocorrelation and external climate variable adjustments is shown to be a better alternative. Copyright © 2014 John Wiley & Sons, Ltd.     

一种新的全球对流层天顶延迟模型GZTD

对流层延迟是GNSS导航定位主要误差源之一,主要受气象参数(如总气压、温度和水汽压等)的影响,具有变化随机性强的特点.本文利用 GGOS Atmosphere提供的2002-2009年全球天顶对流层延迟格网时间序列研究了全球对流层天顶延迟的时空变化特征.并以此为基础对全球天顶对流层延迟(Zenith Troposphere Delay, ZTD)进行建模,提出了一种基于球谐函数的全球非气象参数对流层天顶延迟改正模型--GZTD模型.实验对比结果表明考虑ZTD经纬向变化的GZTD模型内符合精度全球统计结果(bias:0.2 cm,RMS:3.7 cm)优于只考虑ZTD纬向变化的UNB3m (bias:3.4 cm,RMS:6.0 cm)、UNB4 (bias:4.7 cm,RMS:7.4 cm)、UNB3 (bias:4.0 cm,RMS:7.0 cm)和EGNOS (bias:4.5 cm,RMS:6.9 cm)等模型.使用全球385个IGS站进行外符合检验,统计结果表明GZTD模型(bias:-0.02 cm,RMS:4.24 cm)同样优于其它模型.GZTD模型具有改正效果良好、使用简单、所需参数少等优点.     

川滇块体东边界主要断裂带现今运动特征分析

基于2009年以来的GPS观测数据,利用块体模型和GPS剖面方法分别计算川滇块体东边界主要断裂带的滑动速度,并结合跨断裂带的区域应变时间序列分析断裂带现今的运动特征。结果表明:从速度场变化来看,2013—2015期的速度场在川滇块体东北部有东向增加的微弱变化;从滑动速率结果来看,鲜水河北段的左旋走滑运动有所增强,拉张运动有所增加;小江断裂带的左旋走滑运动普遍有微弱的增强;从去掉线性的区域应变时间序列结果来看,小江断裂带南段主张应变在2014年底出现了趋势性转折,值得进一步关注。     

Long-term trend analysis using discrete wavelet components of annual precipitations measurements in Marmara region (Turkey)

The purpose of this study is to determine the possible trends in annual total precipitation series by using the non-parametric methods such as the wavelet analysis and Mann-Kendall test. The wavelet trend (W-T) analysis is for the first time presented in this study. Using discrete wavelet components of measurement series, we aimed to find which periodicities are mainly responsible for trend of the measurement series. We found that some periodic events clearly affect the trend of precipitation series. 16-yearly periodic component is the effective component on Bal?kesir annual precipitation data and is responsible for producing a real trend founded on the data. Also, global wavelet spectra and continuous wavelet transform were used for analysis to precipitation time series in order to clarify time-scale characteristics of the measured series. The effects of regional differences on W-T analysis are checked by using records of measurement stations located in different climatic areas. The data set spans from 1929 to 1993 and includes precipitation records from meteorological stations of Turkey. The trend analysis on DW components of the precipitation time series (W-T model) clearly explains the trend structure of data.     

南极半岛地区GPS坐标时间序列噪声分析及形变模式初探

有色噪声广泛存在于各种连续GPS站坐标时间序列,对GPS时间序列分析有重要影响.利用GAMIT/GLOBK软件解算了南极半岛地区8个GPS测站2010—2014年的实测数据,对坐标时间序列使用主分量分析法(PCA)进行了空间滤波,利用CATS软件估计了不同噪声模型下和空间滤波前后的噪声量级、站坐标时间序列参数及其不确定度,最后对南极半岛地区水平和垂向的形变模式进行了分析和讨论.结果表明,南极半岛地区GPS时间序列不仅存在白噪声,还存在较大量级的闪烁噪声,部分测站E方向在滤波前可能存在随机游走噪声;空间滤波能够有效降低这三种噪声的量级,从而有效减小线性项和周期项估计的不确定度;南极半岛地区在水平方向主要表现为板块运动,还可能存在局部性构造运动;在垂直方向上由冰川均衡调整(Glacial Isostatic Adjustment,GIA)因素引起的抬升较小,主要表现为现今冰雪质量损失引起的弹性抬升运动.     

Geodetic evidence on segmentation of Cascadia subduction zone based on episodic tremors and slips using multivariate harmonic analysis of GPS time series

Episodic tremor and slip (ETS) events with different recurrence intervals have been observed in abundance all along the Cascadia subduction zone margin. Analysis of seismic records as well as Global Positioning System (GPS) time series of the Pacific Northwest Geodetic Array (PANGA) has suggested three distinct coherent zones for the occurrence of these events. In this paper multivariate harmonic estimation has been deployed for further analysis of the segmentation in this area. Raw time series of 43 permanent GPS stations have been used for this purpose. The GPS stations have been geographically divided into three distinct groups including those in the northern, middle and southern parts of the study area. After the reduction of time series for the linear trend as well as annual and semiannual effects, the data series of each group has been analyzed using the multivariate harmonic estimation technique. Subsequently, different combinations of GPS stations including the stations located in the southern, northern and middle zones have been analyzed. Furthermore, the northern and middle, southern and middle as well as the northern and southern zone pair combinations have also been analyzed. The statistical measure devised for identifying the significant frequencies suggests common periods that are consistent with the recurrence intervals of the ETS events already reported for each of the above three geographic zones. Moreover, the method can provide geodetic evidence, in addition to geophysical ones, on the segmentation of ETSs, provided that the adopted time series are of a sufficient length. The geodetic evidence obtained in this research is consistent with the recurrence intervals as well as the boundaries obtained by the analysis of seismic records. Contrary to univariate harmonic estimation, multivariate approach using spatio-temporal correlation of the GPS time series is capable to detect those ETSs whose impacts on the time series are weak.     

Analysis of seasonal signals and long-term trends in the height time series of IGS sites in China

The seasonal signal and long-term trend in the height time series of 10 IGS sites in China are investigated in this paper. The offset detection and outlier removal as well as the removal of common mode error are performed on the raw GPS time-series data developed by the Scripps Orbit and Permanent Array Center (SOPAC). The seasonal-trend decomposition procedure based on LOESS (STL) is utilized to extract precise seasonal signals, followed by an estimation of the long-term trend with the application of maximum likelihood estimation (MLE) to the seasonally adjusted time series. The Up-components of all sites are featured by obvious seasonal variations, with significant phase and amplitude modulation on some sites. After Kendall’s tau test, a significant trend (99% confidence interval) for all sites is achieved. Furthermore, the trends at sites TCMS and TNML have significant changes at epochs 2009.5384 and 2009.1493 (95% confidence interval), respectively, using the Breaks For Additive Seasonal and Trend test. Finally, the velocities and their uncertainties for all sites are estimated using MLE with the white noise plus flicker noise model. And the results are analyzed and compared with those announced by SOPAC. The results obtained in this paper have a higher precision than the SOPAC results.     

亚洲地区ECMWF/NCEP资料计算 ZTD的精度分析

对流层延迟是卫星导航定位的主要误差源,气象观测的数值预报资料可用来计算对流层延迟改正量.本文通过分布于亚洲地区的49个GPS台站一年的实测ZTD资料,对利用欧洲中尺度天气预报中心(ECMWF)分析资料、美国国家环境预报中心(NCEP)再分析资料和NCEP预报资料,计算对流层天顶延迟(ZTD)改正的有效性和可能达到的精度进行了评估,分析了ECMWF和NCEP在亚洲地区的适用程度和其分辨率对计算ZTD精度的影响.研究结果表明:(1)相对于 GPS实测ZTD,用ECMWF资料计算ZTD的bias和rms分别为-1.0 cm 和2.7 cm,优于NCEP再分析资料,可用于高精度ZTD研究和应用;NCEP预报数据计算ZTD的bias和rms分别为2.4 cm 和 6.8 cm,足以满足广大GNSS实时导航定位用户对流层延迟改正的需要.(2)bias和rms呈现明显的季节性变化,总体上夏季大,冬季小;在空间分布上随着纬度的变化不明显,但随高度的增加rms总体上有递减趋势.另外还发现,亚洲东部地区夏季日平均bias和rms和南部热带地区冬季的日平均bias和rms变化相对较大.(3)ECMWF2.5°和0.5°的资料进行了对比分析,发现0.5°分辨率资料的rms比2.5°减小1~5 mm.这些结果,为在亚洲地区的空间大地测量、导航定位和INSAR等工作中,应用ECMWF/NCEP的资料进行对流层大气延迟改正的有效性和可能达到的精度提供了重要参考.     

Attribution analysis for the failure of CMIP5 climate models to simulate the recent global warming hiatus

The Coupled Model Inter-comparison Project Phase 5(CMIP5)contains a group of state-of-the-art climate models and represents the highest level of climate simulation thus far.However,these models significantly overestimated global mean surface temperature(GMST)during 2006-2014.Based on the ensemble empirical mode decomposition(EEMD)method,the long term change of the observed GMST time series of HadCRUT4 records during 1850-2014 was analyzed,then the simulated GMST by 33 CMIP5 climate models was assessed.The possible reason that climate models failed to project the recent global warming hiatus was revealed.Results show that during 1850-2014 the GMST on a centennial timescale rose with fluctuation,dominated by the secular trend and the multi-decadal variability(MDV).The secular trend was relatively steady beginning in the early 20th century,with an average warming rate of 0.0883℃/decade over the last 50 years.While the MDV(with a~65-year cycle)showed2.5 multi-decadal waves during 1850-2014,which deepened and steepened with time,the alarming warming over the last quarter of the 20th century was a result of the concurrence of the secular warming trend and the warming phase of the MDV,both of which accounted one third of the temperature increase during 1975-1998.Recently the slowdown of global warming emerged as the MDV approached its third peak,leading to a reduction in the warming rate.A comparative analysis between the GMST time series derived from HadCRUT4 records and 33 CMIP5 model outputs reveals that the GMSTs during the historical simulation period of 1850-2005 can be reproduced well by models,especially on the accelerated global warming over the last quarter of 20th century.However,the projected GMSTs and their linear trends during 2006-2014 under the RCP4.5 scenario were significantly higher than observed.This is because the CMIP5 models confused the MDV with secular trend underlying the GMST time series,which results in a fast secular trend and an improper MDV with irregular phases and small amplitudes.This implies that the role of atmospheric CO_2 in global warming may be overestimated,while the MDV which is an interior oscillation of the climate system may be underestimated,which should be related to insufficient understanding of key climatic internal dynamic processes.Our study puts forward an important criterion for the new generation of climate models:they should be able to simulate both the secular trend and the MDV of GMST.     

Standard normal homogeneity test as a tool to detect change points in climate-related river discharge variation: case study of the Kupa River Basin

ABSTRACT

The study addresses homogeneity testing of annual discharge time series for eight hydrological stations and five annual climate time series for one weather station in the Kupa River Basin, between Slovenia and Croatia, and global annual average surface temperature time series for the period 1961–2010. The standard normal homogeneity test (SNHT) was used to detect both abrupt and gradual linear trend homogeneity breaks. The results reveal natural change points at the beginning of the 1980s. Absolute homogeneity testing of average annual weather station-level air pressure, annual precipitation, differences between precipitation totals and potential evapotranspiration and surface runoff from the independent observation time series confirmed an abrupt shift, also at the beginning of the 1980s. The trend of local air temperature for 1985–2000, which partly coincides with global surface temperature trend for 1974–2005, strengthened the river discharge regime shift since the beginning of the 1980s. These results could improve climate variation monitoring and estimation of the impact of climate variation on the environment in the area. Generally, an indication of climate regime change points and an assessment of their duration could provide significant benefits for the society.     

基于多通道奇异谱分析的GNSS坐标时间序列共模误差的提取

共模误差是区域连续GNSS网中存在的一种与时空相关的主要误差源.为了有效的剔除共模误差,提高坐标时间序列的精度,本文提出了利用多通道奇异谱分析（Multi-channel Singular Spectrum Analysis,MSSA）提取共模误差的新思路,并利用实验区域18个测站9年（2002年到2010年）的GPS坐标时间序列进行实验,分析了共模误差对时间序列的影响和测站噪声特性的影响,并对共模误差序列进行周期探测,结果显示：通过MSSA能够有效的剔除共模误差,提高坐标时间序列的精度.     

The value of large-scale climate variables in climate change assessment: The case of Botswana’s rainfall

Climate change is expected to alter rainfall regimes across most parts of the world. The implications of this could be more severe in arid environments where rainfall is limited and highly variable in space and time. However, lack of good quality data, of sufficient record length and spatial coverage usually restricts model development and performance geared towards assessing the effects of climate change in these areas. This paper presents an analysis of rainfall and climate data in order to determine the time of change in rainfall series and identify possible correlations between rainfall and temperature. In addition, the paper aims to make predictions of future rainfall patterns in Botswana. This is achieved by using historical rainfall and climate data from rainfall stations spread across Botswana from 1965 to 2008. In addition, large scale reanalysis data from NCAR/NCEP and El Nino Southern Oscillation (ENSO) data were used to augment the limited observed spatial climate data series when developing a rainfall model. Temperature and ENSO indices have been used to predict rainfall regimes for the present climate. Based on these, the effects of climate change were quantified using a stochastic generalised linear rainfall model (GLM) driven by outputs of global climate models (GCMs). The results indicate that temperature is a significant rainfall predictor in Botswana compared to ENSO indices.     

Analysis of Abnormal Characteristics of Regional Crustal Deformation before the Menyuan MS6.4 Earthquake by GPS Continuous Data

In order to study the characteristics of crustal deformation around the epicenter before the 2016 M_S6.4 Menyuan earthquake, the GPS continuous stations of the period from 2010 to 2016 were selected according to the observation data of the tectonic environment monitoring network in Chinese Mainland. The deformation characteristics of the crust before the earthquake were discussed through inter-station baseline time series analysis and the strain time series analysis in the epicentral region. The results show that a trend turn of the baseline movement state around the epicenter region occurred after 2014, and the movement after 2014 reflects an obvious decreasing trend of compressional deformation. During this period, the stress field energy was in a certain accumulation state. Since the beginning of 2014, the EW-component linear strain and surface strain rate weakened gradually before the earthquake. It shows that there was an obvious deformation deficit at the epicentral area in the past two years, which indicates that the region accumulated a high degree of strain energy before the earthquake. Therefore, there was a significant background change in the area before the earthquake. The results of the study can provide basic research data for understanding the seismogenic process and mechanism of this earthquake.     

Diurnal and semidiurnal atmospheric tides observed by co-located GPS and VLBI measurements

The tropospheric zenith total delay (ZTD) derived from very long baseline interferometry (VLBI) is an important parameter of the atmosphere, reflecting various atmosphere-related processes and variations. In this paper, ZTD time series of the IVS rapid combined tropospheric product (2002–2006) with a 1-h resolution are used for the first time to investigate the diurnal and semidiurnal oscillations. Significant diurnal and semidiurnal variations of ZTD are found at all VLBI stations. The amplitude of the diurnal cycle S₁ is 0.6–1.2 mm at most of the VLBI stations, and the amplitude of the semidiurnal cycle S₂ is 0.2–1.9 mm, which nearly accord with the surface pressure tides S₁/S₂ and co-located GPS estimated S₁/S₂. The results indicate that the S₁ and S₂ behaviors are mainly dominated by the hydrostatic component, namely pressure tides. In general, the semidiurnal S₂ amplitudes are slightly larger than the diurnal S₁. While S₁ shows no clear dependency on site altitude, S₂ has a regular distribution with VLBI site altitude. The results are in accordance with predictions of the classic tidal theory [Chapman, S., Lindzen, R.S., 1970. Atmospheric Tides, Gordon and Breach, New York].     

Letters to the Editor

Abstract

This work investigates historical trends of meteorological drought in Taiwan by means of long-term precipitation records. Information on local climate change over the last century is also presented. Monthly and daily precipitation data for roughly 100 years, collected by 22 weather stations, were used as the study database. Meteorological droughts of different levels of severity are represented by the standardized precipitation index (SPI) at a three-monthly time scale. Additionally, change-point detection is used to identify meteorological drought trends in the SPI series. Results of the analysis indicate that the incidence of meteorological drought has decreased in northeastern Taiwan since around 1960, and increased in central and southern Taiwan. Long-term daily precipitation series show an increasing trend for dry days all over Taiwan. Finally, frequency analysis was performed to obtain further information on trends of return periods of drought characteristics.