An estimation of tropospheric corrections using GPS and synoptic data: Improving Urmia Lake water level time series from Jason-2 and SARAL/AltiKa satellite altimetry

Tropospheric correction is one of the most important corrections in satellite altimetry measurements. Tropospheric wet and dry path delays have strong dependence on temperature, pressure and humidity. Tropospheric layer has particularly high variability over coastal regions due to humidity, wind and temperature gradients. Depending on the extent of water body and wind conditions over an inland water, Wet Tropospheric Correction (WTC) is within the ranges from a few centimeters to tens of centimeters. Therefore, an extra care is needed to estimate tropospheric corrections on the altimetric measurements over inland waters. This study assesses the role of tropospheric correction on the altimetric measurements over the Urmia Lake in Iran. For this purpose, four types of tropospheric corrections have been used: (i) microwave radiometer (MWR) observations, (ii) tropospheric corrections computed from meteorological models, (iii) GPS observations and (iv) synoptic station data. They have been applied to Jason-2 track no. 133 and SARAL/AltiKa track no. 741 and 356 corresponding to 117–153 and the 23–34 cycles, respectively. In addition, the corresponding measurements of PISTACH and PEACHI, include new retracking method and an innovative wet tropospheric correction, have also been used. Our results show that GPS observation leads to the most accurate tropospheric correction. The results obtained from the PISTACH and PEACHI projects confirm those obtained with the standard SGDR, i.e., the role of GPS in improving the tropospheric corrections. It is inferred that the MWR data from Jason-2 mission is appropriate for the tropospheric corrections, however the SARAL/AltiKa one is not proper because Jason-2 possesses an enhanced WTC near the coast. Furthermore, virtual stations are defined for assessment of the results in terms of time series of Water Level Height (WLH). The results show that GPS tropospheric corrections lead to the most accurate WLH estimation for the selected virtual stations, which improves the accuracy of the obtained WLH time series by about 5%.     

A DORIS determination of the absolute velocities of the Sorsdal and Mellor glaciers in Antarctica

The Lambert–Amery System is the largest glacier–ice shelf system in East Antarctica, draining a significant portion of the ice sheet. Variation in ice sheet discharge from Antarctica or Greenland has an impact on the rate of change in global mean sea level; which is a manifestation of climate change. In conjunction with a measure of ice thickness change, ice sheet discharge can be monitored by determining the absolute velocities of these glaciers.     

Integrating Non-Tidal Sea Level data from altimetry and tide gauges for coastal sea level prediction

The main objective of this paper is to integrate Non-Tidal Sea Level (NSL) from the joint TOPEX, Jason-1 and Jason-2 satellite altimetry with tide gauge data at the west and north coast of the United Kingdom for coastal sea level prediction. The temporal correlation coefficient between altimetric NSLs and tide gauge data reaches a maximum higher than 90% for each gauge. The results show that the multivariate regression approach can efficiently integrate the two types of data in the coastal waters of the area. The Multivariate Regression Model is established by integrating the along-track NSL from the joint TOPEX/Jason-1/Jason-2 altimeters with that from eleven tide gauges. The model results give a maximum hindcast skill of 0.95, which means maximum 95% of NSL variance can be explained by the model. The minimum Root Mean Square Error (RMSe) between altimetric observations and model predictions is 4.99 cm in the area. The validation of the model using Envisat satellite altimetric data gives a maximum temporal correlation coefficient of 0.96 and a minimum RMSe of 4.39 cm between altimetric observations and model predictions, respectively. The model is furthermore used to predict high frequency NSL variation (i.e., every 15 min) during a storm surge event at an independent tide gauge station at the Northeast of the UK (Aberdeen).     

Characterizing slope correction methods applied to satellite radar altimetry data: A case study around Dome Argus in East Antarctica

Slope correction is important to improve the accuracy of satellite radar elevation measurements by mitigating the slope-induced error (SE), especially over uneven ground surfaces. Although several slope correction methods have been proposed, guidance in the form of stepwise algorithm on how to implement these methods in processing radar altimetric data at the coding level, and the differences among these methods need to be presented and discussed systematically. In this paper, three existing types of slope correction methods—the direct method (DM), intermediate method (IM), and relocation method (RM, further divided into RM1 and RM2)—are described in detail. In addition, their main differences and features for various scientific applications are analyzed. We conduct a systematic experiment with CryoSat-2 Low Resolution Mode (LRM) data in a physically stable area around Dome Argus in East Antarctica, where in-situ measurements were available for comparison. The slope correction is implemented separately using the three methods, with the latest high-accuracy Reference Elevation Model of Antarctica (REMA) as the a-priori topography model. The bias and precision of the slope-corrected CryoSat-2 data results from the RM2 is ?0.18 ± 0.86 m based on the comparison with the field Global Navigation Satellite System (GNSS) data. The results from the RM2 indicate higher precision compared to those from the RM1. According to the correlation analysis of the slope-corrected CryoSat-2 data results (RM1 and RM2), the bias enlarges and the precision becomes worse when the surface slope increases from 0 to 0.85°. After a comprehensively comparative analysis, we find that the results from the RM1 and RM2 are superior in precision (0.93 m and 0.86 m) with respect to the GNSS data. The relatively low precision (1.22 m) from the IM is due to the potential error from the a-priori digital elevation model (DEM). The DM has the lowest precision (2.66 m). Another experiment over rough topography in West Antarctica is carried out for comparison, especially between the RM1 (precision of 15.27 m) and RM2 (precision of 16.25 m). In general, the RM is recommended for the SE elimination among the three methods. Moreover, the RM2 is firstly considered over smooth topography due to the superior performance in bias and precision, while the RM1 is more suggested over the rough topography because of the slightly smaller bias and better precision. The IM relies much on the accuracy of the a-prior DEM and is not usually recommended, because of the strict requirement in the sampling time between the radar altimetry data and the a-priori DEM to avoid any surface change over time.     

Preliminary study of differential force modelling on ERS-1 and ERS-2

The two near identical satellites ERS-1 and ERS-2 are orbiting the Earth in the same orbit and same orbital plane with ERS-1 some thirty-two minutes in the lead. If their altimetric datasets are differenced along the repeating ground-track any fixed geographically correlated signature will cancel. The residual signature of this differenced dataset will therefore arise from imprecise corrections to the altimetric dataset but more significantly from temporal variations in the forces affecting both satellites; in particular luni-solar attraction and aerodynamic force modelling. This study presents preliminary results of differential force modelling on the satellites with an emphasis on the aerodynamic effects.     

Periods of quasi-stationary conditions in interplanetary space according to sequences of SEP events

The values of the characteristic decay time of particle fluxes in SEP events vary, as a rule, considerably from event to event. We point out, however, that at times sequences of events having similar decay times were observed over long time intervals (up to one month, and even longer in a few cases). The values of the decay times, however, differed among different sequences. The constancy of the decay phase in each consecutive event of these series suggests that the interplanetary medium was in steady state during the event series, and, because of solar rotation, its uniformity within sectors extended to 90–180° in heliolongitude. The very rarely observed long series (up to 2–3 solar rotations) indicate the steadiness and homogeneity of the plasma and the interplanetary magnetic field (IMF) in the entire inner solar system in the course of this time span. It is pointed out that the neutral current sheet of the IMF does not represent a substantial obstacle for energetic charged particles. Both hemispheres are (above and below the current sheet), at least during the series of solar events, invariant with time, uniform and alike from the viewpoint of the propagation of charged particles. The investigation of such sequences of events can also be useful for forecasting characteristics of SEP events.     

Quality Assessment of the IDS Contribution to ITRF2008

Doppler Orbitography Radiopositionning Integrated by Satellite (DORIS) is one of the four fundamental techniques contributing to the ITRF. The optimal coverage over the globe of the DORIS observing sites and sites co-located with GPS, allow a strong embedding of DORIS within the ITRF network. DORIS contributes to the access to ITRF through precise orbit determination of altimetric satellites with onboard DORIS receivers. The DORIS contribution to the ITRF2008 is enhanced by the fact that the solutions of seven analysis centers were included in the submitted combined time series of weekly station positions and daily polar motion. We evaluate the quality of the DORIS combined solution in terms of its agreement with the other techniques (VLBI, SLR, GPS) contributing to the ITRF2008 combination. We show in particular that the precisions of the current IDS products range between 1.5 to 2.6 mm for station positions (at the epochs of minimum variances); better than 1 mm/yr in velocities and between 170 and 260 micro-arc-seconds for polar motion, a significant improvement by a factor of three to five, compared to past data used in the ITRF2005 combination. This improvement is certainly due to improved analysis strategies employed by the seven IDS analysis centers that contributed to the combined weekly submitted solutions of station positions and polar motion. A spectral analysis of DORIS station height time series indicates that annual and semi-annual signals are dominant. However, TOPEX draconitic period of about 118 days is still detected in about 20% of the station position power spectra. DORIS height annual signals correlate well with GPS annual signal estimated at some co-located stations, which show that DORIS technique is able to detect loading signals.     

Exploring the potential of Sentinel-3 delay Doppler altimetry for enhanced detection of coastal currents along the Northwest Atlantic shelf

In this study, we evaluate Sentinel-3A satellite synthetic aperture radar (SAR) altimeter observations along the Northwest Atlantic coast, spanning the Nova Scotian Shelf, Gulf of Maine, and Mid-Atlantic Bight. Comparisons are made of altimeter sea surface height (SSH) measurements from three different altimeter data processing approaches: fully-focused synthetic aperture radar (FFSAR), un-focused SAR (UFSAR), and conventional low-resolution mode (LRM). Results show that fully-focused SAR data always outperform LRM data and are comparable or slightly better than the nominal un-focused SAR product. SSH measurement noise in both SAR-mode datasets is significantly reduced compared to LRM. FFSAR SSH 20-Hz noise levels, derived from 80-Hz FFSAR data, are lower than 20-Hz UFSAR SSH with 25% noise reduction offshore of 5 km, and 55–70% within 5 km of the coast. The offshore noise improvement is most likely due to the higher native along-track data posting rate (80 Hz for FFSAR, and 20 Hz for UFSAR), while the large coastal improvement indicates an apparent FFSAR data processing advantage approaching the coastlines. FFSAR-derived geostrophic ocean current estimates exhibit the lowest bias and noise when compared to in situ buoy-measured currents. Assessment at short spatial scales of 5–20 km reveals that Sentinel-3A SAR data provide sharper and more realistic measurement of small-scale sea surface slopes associated with expected nearshore coastal currents and small-scale gyre features that are much less well resolved in conventional altimetric LRM data.     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

机翼前缘积冰对大飞机操稳特性的影响

机翼结冰影响飞机的操稳特性和飞行性能,对飞行安全造成危害。基于实验数据构建了典型的不同结冰严重程度的机翼前缘积冰冰形,采用高精度数值模拟方法得到背景飞机机翼前缘积冰的气动数据,建立了飞机六自由度非线性动力学模型,在此基础上设计了俯仰角保持、滚转角保持及高度保持模式的自动驾驶仪闭环仿真系统。通过开环仿真,分析了不同程度积冰对飞机配平特性、纵向长短周期模态及横航向模态的影响,比较了不同程度积冰情形下飞机动态响应的差异。通过闭环仿真,研究了积冰对3种模式下自动驾驶性能的影响。仿真结果表明:积冰对飞机配平特性、模态特性及开环动态响应特性均会造成一定的不良影响,威胁飞行安全。      

差分 GPS 鉴定地基雷达测量系统精度的方法研究

提出了利用差分ＧＰＳ分析地面雷测系统跟踪误差的思想,详细分析了各种误差源,建立了一组简洁的差分ＧＰＳ模型,构造了测量数据的比对分析算法以及具有良好抗野值能力的误差偏差最优统计算法,并利用模拟数据进行了仿真计算     

AE指数表现出的磁层的浑沌行为

时间序列的非线性预测误差可用来辨别产生这个时间序列的动力学系统是浑沌的还是随机的。据此,本文分析了自1982年1月30日12点开始的分辨率为1min的AE指数构成的时间序列,共15000个数据点。首先用其中的一部分数据构造了一个从非线性确定性的到线性统计的范围很宽的预测模型,然后利用此模型对另一部分数据进行预测。预测误差的结果表明,AE指数所表征的磁层具有浑沌行为。用同样的数据,利用滞后时间相空间重构法得到的系统吸引子的关联维为2.5.      

Cycle slip detection using multi-frequency GPS carrier phase observations: A simulation study

The detection and repair of the cycle slip or gross error is a key step for high precision global positioning system (GPS) carrier phase navigation and positioning due to interruption or unlocking of GPS signal. A number of methods have been developed to detect and repair cycle slips in the last two decades through cycle slip linear combinations of available GPS observations, but such approaches are subject to the changing GPS sampling and complex algorithms. Furthermore, the small cycle slip and gross error cannot be completely repaired or detected if the sampling is quite longer under some special observation conditions, such as Real Time Kinematic (RTK) positioning. With the development of the GPS modernization or Galileo system with three frequencies signals, it may be able to better detect and repair the cycle slip and gross error in the future. In this paper, the cycle slip and gross error of GPS carrier phase data are detected and repaired by using a new combination of the simulated multi-frequency GPS carrier phase data in different conditions. Results show that various real-time cycle slips are completely repaired with a gross error of up to 0.2 cycles.     

Evaluation and utilization of CloudSat and CALIPSO data to analyze the impact of dust aerosol on the microphysical properties of cirrus over the Tibetan Plateau

The present study elucidates on the evaluation of two versions (V3 and V4.10) of vertical feature mask (VFM) and aerosol sub-types data derived from the Cloud-Aerosol LiDAR and Infrared Pathfinder Satellite Observations (CALIPSO), and its utilization to analyze the impact of dust aerosol on the microphysical properties of cirrus over the Tibetan Plateau (TP). In conjunction to the CALIPSO, we have also used the CloudSat data to study the same during the summer season for the years 2007–2010 over the study area 25–40°N and 75–100°E. Compared to V3 of CALIPSO, V4.10 was found to have undergone substantial changes in the code, algorithm, and data products. Intercomparison of both versions of data products in the selected grid between 30–31°N and 83–84°E within the study area during 2007–2017 revealed that the VFM and aerosol sub-types are in good agreement of ～95.27% and ～82.80%, respectively. Dusty cirrus is defined as the clouds mixed with dust aerosols or existing in dust aerosol conditions, while the pure cirrus is that in a dust-free environment. The obtained results illustrated that the various microphysical properties of cirrus, namely ice water content (IWC), ice water path (IWP), ice distribution width (IDW), ice effective radius (IER), and ice number concentration (INC) noticed a decrease of 17%, 18%, 4%, 19%, and 10%, respectively due to the existence of dust aerosol, consistent with the classical “Twomey effect” for liquid clouds. Moreover, the aerosol optical depth (AOD) showed moderate negative correlations between ?0.4 and ?0.6 with the microphysical characteristics of cirrus. As our future studies, in addition to the present work undertaken, we planned to gain knowledge and interested to explore the impact of a variety of aerosols apart from the dust aerosol on the microphysical properties of cirrus in different regions of China.     

Searching for life in extreme environments relevant to Jovian’s Europa: Lessons from subglacial ice studies at Lake Vostok (East Antarctica)

The objective was to estimate the genuine microbial content of ice samples from refrozen water (accretion ice) from the subglacial Lake Vostok (Antarctica) buried beneath the 4-km thick East Antarctic ice sheet. The samples were extracted by heavy deep ice drilling from 3659 m below the surface. High pressure, a low carbon and chemical content, isolation, complete darkness and the probable excess of oxygen in water for millions of years characterize this extreme environment. A decontamination protocol was first applied to samples selected for the absence of cracks to remove the outer part contaminated by handling and drilling fluid. Preliminary indications showed the accretion ice samples to be almost gas free with a low impurity content. Flow cytometry showed the very low unevenly distributed biomass while repeated microscopic observations were unsuccessful.     

单星测频静态目标无源定位研究

提出了一种基于非线性最优化技术的单星测频无源定位算法,该方法利用单颗卫星单个通道在不同位置上测得的信号多普勒频率,实现对地面固定辐射源的无源定位,具有有效载荷简单、对卫星姿态无特殊要求、定位收敛快、精度较高等优点。同时,探讨了多普勒频率测量误差预处理技术,以及基于PAR模型的系统误差分离技术,通过对测量误差的分离与处理,可以大幅度提高频率测量数据的精度,从而保证后续测频定位精度。最后,通过仿真试验和工程实际,验证了单星测频定位技术的有效性和准确性。     

Characteristics and accuracies of the GRACE inter-satellite pointing

For almost 10 years, the Gravity Recovery and Climate Experiment (GRACE) has provided information about the Earth gravity field with unprecedented accuracy. Efforts are ongoing to approach the GRACE baseline accuracy as there still remains an order of magnitude between the present error level of the gravity field solutions and the GRACE baseline. At the current level of accuracy, thorough investigation of sensor related effects is necessary as they are one of the potential contributors to the error budget. In the science mode operations, the twin satellites are kept precisely pointed with their KBR antennas towards each other. It is the task of the onboard attitude and orbit control system (AOCS) to keep the satellites in the required formation. We analyzed long time series of the inter-satellite pointing variations as they reflect the AOCS performance and characteristics. We present significant systematic effects in the inter-satellite pointing and discuss their possible sources. Prominent features are especially related to the magnetic torquer characteristics, star cameras’ performance and KBR antenna calibration parameters. The relation between the magnetic torquer attitude control and the Earth magnetic field, impact of the different performance of the two star camera heads on the attitude control and the features due to uncertainties in the calibration parameters relating the star camera frame to K-frame are discussed in detail. Proper understanding of these effects will help to reduce their impact on the science data and subsequently increase the accuracy of the gravity field solutions. Moreover, understanding the complexity of the onboard system is essential not only for increasing the accuracy of the GRACE data but also for the development of the future gravity field satellite missions.     

奇异谱分析在故障时间序列分析中的应用

由于日益增长的飞行安全和飞机维护质量需求,飞机使用可靠性已经成为一个重要的研究领域。从某航空公司波音737飞机使用过程中现场所记录的18年的故障数据出发,应用奇异谱分析（SSA）方法,对故障时间序列进行了建模和预测,进一步以预测结果的均方根误差（RMSE）最小为优化目标对SSA模型参数进行了优选。在此基础上,提出了一种更为广泛的模型组合方法和实现算法,这种方法采用不同的时间序列模型来构造SSA分解出的趋势、周期和残差等成分。通过与三次指数平滑（Holt-Winters）、自回归移动平均（ARIMA）2种时间序列模型的实验结果对比,SSA及其参数优选和模型组合方法在故障时间序列分析中具有更好的拟合和预测精度。