Rupture process of the 2011 Tohoku earthquake from the joint inversion of teleseismic and GPS data

Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.

A Rapid Optimization Algorithm for GPS Data Assimilation

Global Positioning System (GPS) meteorology data variational assimilation can be reduced to the problem of a large-scale unconstrained optimization. Because the dimension of this problem is too large, most optimal algorithms cannot be performed. In order to make GPS/MET data assimilation able to satisfy the demand of numerical weather prediction, finding an algorithm with a great convergence rate of iteration will be the most important thing. A new method is presented that dynamically combines the limited memory BFGS (L-BFGS) method with the Hessian-free Newton(HFN) method, and it has a good rate of convergence in iteration. The numerical tests indicate that the computational efficiency of the method is better than the L-BFGS and HFN methods.

Source of the 2019 M_(W)6.9 Banten Intraslab earthquake modelled with GPS data inversion

The 2019 M_(W)6.9 Banten Intraslab earthquake occurred at~100 km to the northeast of the Sunda Trench with two nodal plane models estimated by the Global Centroid-Moment-Tensor(GCMT)Project with a strike of 200°and a dip of 65°.Continuous GPS data from 11 GPS sites were used to model the source of the earthquake in three-components.The coseismic displacements and its uncertainties are obtained from the coordinates of these GPS sites from 7 days before to 7 days after the earthquake.The coseismic slip is the inversion result of those displacements based on the best fit in an elastic half-space.The maximum displacement is~5 cm with a large uncertainty that is comparable to the amplitude of displacement.A seismic moment of the best model(strike of 65°and dip of 54°)is 2.79×10^(19) Nm or equivalent to M_(W)6.89.The fault model of the earthquake is highly presumed as a continuation of Sumatran Fault Zone.

Relationship between regional stress field variation and earthquake activities from GPS data

The relationship between the earth stress field, earth surface displacement field and the gravity variation is deduced. Algorithms based on the boundary element method to compute the earth stress variation using the earth surface displacement is discussed. The stress field variation in Jiashi region, Xinjiang, China is obtained from the GPS data observed in 1997 and 1998, respectively, and the relationship among the local stress field variation, seismic activities and fault tectonic activities is discussed.

A Survey of Sediment Fineness and Moisture Content in the Soyang Lake Floodplain Using GPS Data

Soyang Lake is the largest lake in Republic of Korea bordering Chuncheon,Yanggu,and Inje in Gangwon Province.It is widely used as an environmental resource for hydropower,flood control,and water supply.Therefore,we conducted a survey of the floodplain of Soyang Lake to analyze the sediments in the area.We used global positioning system(GPS)data and aerial photography to monitor sediment deposits in the Soyang Lake floodplain.Data from three GPS units were compared to determine the accuracy of sampling location measurement.Sediment samples were collected at three sites:two in the eastern region of the floodplain and one in the western region.A total of eight samples were collected:Three samples were collected at 10 cm intervals to a depth of 30 cm from each site of the eastern sampling point,and two samples were collected at depths of 10 and 30 cm at the western sampling point.Samples were collected and analyzed for vertical and horizontal trends in particle size and moisture content.The sizes of the sediment samples ranged from coarse to very coarse sediments with a negative slope,which indicate eastward movement from the breach.The probability of a breach was indicated by the high water content at the eastern side of the floodplain,with the eastern sites showing a higher probability than the western sites.The results of this study indicate that analyses of grain fineness,moisture content,sediment deposits,and sediment removal rates can be used to understand and predict the direction of breach movement and sediment distribution in Soyang Lake.

Estimation of postseismic deformation parameters from continuous GPS data in northern Sumatra after the 2004 Sumatra-Andaman earthquake

Continuous global positioning system （GPS） in northern Sumatra detected signal of the ongoing physical process of postseismic deformation after the M9.2 2004 Sumatra-Andaman earthquake. We analyze the character- istics of postseismic deformation of the 2004 earthquake based on GPS networks operated by BIG, and the others named AGNeSS, and SuGAr networks located in northern Sumatra. We use a simple analytical logarithmic and exponential function to evaluate the postseismic deforma- tion parameters of the 2004 earthquake. We find that GPS data in northern Sumatra during time periods of 2005-2012 are fit better using the logarithmic function with dog of 104.2 ± 0.1 than using the exponential function. Our result clearly indicates that other physical mechanisms of post- seismic deformation should be taken into account rather than a single physical mechanism of afterslip only.

Co-seismic Slip Distribution of the 2001 Kokoxili Earthquake Inverted by the GPS Data

The authors analyze co_seismic displacement field derived from the Global Position System (GPS) observations collected before and after the 2001 Kokoxili earthquake, western China. Using the co_seismic displacement data, and constrained with surface rupture data, they invert co_seismic slip distribution along the seismic fault. Their result shows that the earthquake ruptured the upper crust down to a depth of 13.1～22km (at 70% certainty), with its optimal estimate at 16.5km. A 2～3m left_lateral strike slip is resolved between the Sun Lake segment and the west end of the main rupture zone, although surface rupture has not been observed there. The surface rupture of this earthquake is ended at the Sun Lake to the west, but left_lateral slip of 1.5～2.0m seems to exist beyond the east end of surface rupture observed from field geology. Seismic moment release estimated using GPS and surface rupture measurement is 6.0×10 20 N·m, which is in good agreement with the result obtained from seismic wave inversion.

Water storage changes in North America retrieved from GRACE gravity and GPS data

As global warming continues,the monitoring of changes in terrestrial water storage becomes increasingly important since it plays a critical role in understanding global change and water resource management.In North America as elsewhere in the world,changes in water resources strongly impact agriculture and animal husbandry.From a combination of Gravity Recovery and Climate Experiment（GRACE） gravity and Global Positioning System（GPS） data,it is recently found that water storage from August,2002 to March,2011 recovered after the extreme Canadian Prairies drought between 1999 and 2005.In this paper,we use GRACE monthly gravity data of Release 5 to track the water storage change from August,2002 to June,2014.In Canadian Prairies and the Great Lakes areas,the total water storage is found to have increased during the last decade by a rate of 73.8 ± 14.5 Gt/a,which is larger than that found in the previous study due to the longer time span of GRACE observations used and the reduction of the leakage error.We also find a long term decrease of water storage at a rate of-12.0 ± 4.2 Gt/a in Ungava Peninsula,possibly due to permafrost degradation and less snow accumulation during the winter in the region.In addition,the effect of total mass gain in the surveyed area,on present-day sea level,amounts to-0.18 mm/a,and thus should be taken into account in studies of global sea level change.

Estimation of the 2001 Kunlun earthquake fault slip from GPS coseismic data using Hori’s inverse method

The Hori＇s inverse method based on spectral decomposition was applied to estimate coseismic slip distribution on the rupture plane of the 14 November 2001 Ms8.1 Kunlun earthquake based on GPS survey results. The inversion result shows that the six sliding models can be constrained by the coseismic GPS data. The established slips mainly concentrated along the eastern segment of the fault rupture, and the maximum magnitude is about 7 m. Slip on the eastern segment of the fault rupture represents as purely left-lateral strike-slip. Slip on the western segment of the seismic rupture represents as mainly dip-stip with the maximum dip-slip about 1 m. Total predicted scalar seismic moment is 5.196× 10^2° N.m. Our results constrained by geodetic data are consistent with seismological results.

Atmospheric density determination using high-accuracy satellite GPS data

Atmospheric drag is the main source of error in the determination and prediction of the orbit of low Earth orbit （LEO） satellites; however, empirical models that are used to account for this often have density errors of around 15%-30%. Atmospheric density determination has thus become an important topic for researchers. Based on the relationship between file atmospheric drag force and the decay of the semi-major axis of the orbit, we derived atmospheric density along the trajectory of challenging mini-satellite payload （CHAMP） satellite with its rapid science orbit （RSO） data. Three primary parameters--the ratio of cross-sectional area to mass, the drag coefficient, and the decay of the semi-major axis caused by atmospheric drag--were calculated. We also analyse the source of the error and made a comparison between the GPS-derived and reference density. The result for December 2, 2008, showed that the mean error of the GPS-derived density could be decreased from 29.21% to 9.20%, if the time span adopted for the process of computation was increased from 10 min to 50 min. The result for the entire month of December indicated that a density precision of 10% could be achieved, when the time span meets the condition that the amplitude of the decay of the semi-major axis is much greater than its standard deviation.

On the Accuracy Assessment of the Latest Releases of GOCE Satellite-Based Geopotential Models with EGM2008 and Terrestrial GPS/Levelling and Gravity Data over Egypt

The Global Geopotential Models (GGMs) of GOCE (Gravity Recovery and steady- state Ocean Circulation Explorer) differ globally as well as regionally in their accuracy and resolution based on the maximum degree and order (d/o) of the fully normalized spherical harmonic (SH) coefficients, which express each GGM. The main idea of this study is to compare the free-air gravity anomalies and quasi geoid heights determined from several recent GOCE-based GGMs with the corresponding ones from the Earth Gravitational Model 2008 (EGM2008) over Egypt on the one hand and with ground-based measurements on the other hand. The results regarding to the comparison of GOCE-based GGMs with terrestrial gravity and GPS/levelling data provide better improvement with respect to EGM2008. The 4th release GOCE-based GGM developed with the use of space-wise solution strategy (SPW_R4) approximates the gravity field well over the Egyptian region. The SPW_R4 model is accordingly suggested as a reference model for recovering the long wavelength (up to SH d/o 200) components of quasi geoid heights when modelling the gravimetric quasi-geoid over the Egypt. Finally, three types of transformation models: Four-, Five- and Seven-parameter transformations have been applied to reduce the data biases and to provide a better fitting of quasi geoid heights obtained from the studied GOCE-based GGMs to those from GPS/levelling data. These models reveal that the standard deviation of vertical datum over Egypt is at the level of about 32 cm.

Integrated Data Processing Method for GPS and INS Field Test over Rocky Mountain

The method of integrated data processing for GPS and INS(inertial navigation system) field test over the Rocky Mountains using the adaptive Kalman filtering technique is presented. On the basis of the known GPS outputs and the offset of GPS and INS, state equations and observations are designed to perform the calculation and improve the navigation accuracy. An example shows that with the method the reliable navigation parameters have been obtained.

Global Statistical Study of Ionospheric Waves Based on COSMIC GPS Radio Occultation Data

Extracting and parameterizing ionospheric waves globally and statistically is a longstanding problem. Based on the multichannel maximum entropy method（MMEM） used for studying ionospheric waves by previous work, we calculate the parameters of ionospheric waves by applying the MMEM to numerously temporally approximate and spatially close global-positioning-system radio occultation total electron content profile triples provided by the unique clustered satellites flight between years 2006 and 2007 right after the constellation observing system for meteorology, ionosphere, and climate（COSMIC） mission launch. The results show that the amplitude of ionospheric waves increases at the low and high latitudes（~0.15 TECU） and decreases in the mid-latitudes（~0.05 TECU）. The vertical wavelength of the ionospheric waves increases in the mid-latitudes（e.g., ~50 km at altitudes of 200–250 km） and decreases at the low and high latitudes（e.g., ~35 km at altitudes of 200–250 km）.The horizontal wavelength shows a similar result（e.g., ~1400 km in the mid-latitudes and ~800 km at the low and high latitudes）.

Evaluation of driving behavior based on massive vehicle trajectory data

Based on the driver surveillance video data and controller area network(CAN)data,the methods of studying commercial vehicles’driving behavior is relatively advanced.However,these methods have difficulty in covering private vehicles.Naturalistic driving studies have disadvantages of small sample size and high cost,one new driving behavior evaluation method using massive vehicle trajectory data is put forward.An automatic encoding machine is used to reduce the noise of raw data,and then driving dynamics and self-organizing mapping(SOM)classification are used to give thresholds or the judgement method of overspeed,rapid speed change,rapid turning and rapid lane changing.The proportion of different driving behaviors and typical dangerous driving behaviors is calculated,then the temporal and spatial distribution of drivers’driving behavior and the driving behavior characteristics on typical roads are analyzed.Driving behaviors on accident-prone road sections and normal road sections are compared.Results show that in Shenzhen,frequent lane changing and overspeed are the most common unsafe driving behaviors;16.1%drivers have relatively aggressive driving behavior;the proportion of dangerous driving behavior is higher outside the original economic special zone;dangerous driving behavior is highly correlated with traffic accident frequency.

GBO algorithm for seismic source parameters inversion

The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.

Investigating distance halo effect of fixed automated speed camera based on taxi GPS trajectory data

Background:The deterrence effect of automated speed camera(ASC)is still inconclusive.Moreover,it is pointed out that ASC may have varying deterrence effects on different types of road users(e.g.,taxis).Objective:This study intends to investigate the distance halo effect of fixed ASC(hereafter called ASC)on taxis.Method:More than 1.34 million taxis’GPS trajectory data were collected.A novel indicator,the delta speed(defined as the difference between the traveling speed and the speed limit),was proposed to continuously describe the variations in traveling speeds.The upstream and downstream critical delta speeds during each time period on weekdays and weekends were obtained by using K-means clustering method,respectively.Based on the critical delta speeds,the ranges of upstream and downstream distance halo effects of ASC during different time periods on weekdays and weekends were determined separately and compared.Results:The downstream critical delta speed is smaller than the upstream one.The upstream and downstream distance halo effects of ASC on taxis are within a range of 8-2180 m and an area of 10-580 m to the ASC location,respectively.There are no obvious difference in the ranges of upstream and downstream distance halo effects of ASC on taxis between different time periods or between weekdays and weekends.Conclusion:The present study confirms that the upstream and downstream distance halo effects of ASC on taxis have different ranges and the stabilities of time-of-day and day-of-week.Practical application:The findings of this study can provide a basic reference for reasonably deploying ASCs within a region.

High Precision Orbit Determination of CHAMP Satellite

The precision orbit determination of challenging minisatellite payload(CHAMP) satellite was done based on position and navigation data analyst(PANDA) software which is developed in Wuhan University, using the onboard GPS data of year 2002 from day 126 to 131. The orbit accuracy was assessed by analyzing the difference from GFZ post-processed science orbits (PSO), the GPS carrier and pseudo-range data residuals and the satellite laser ranging (SLR) residuals.

Numerical Simulation of Annual Change Patterns of Contemporary Tectonic Stress-Strain Field of the Chinese Mainland

Based on the active crustal block structures, the Holocene active faults and the wave velocity structures with a resolution of 1°×1°, a two-dimensional finite element model for the tectonic stress-strain field of the Chinese mainland is constructed in the paper. Using GPS measurements, the velocity boundary conditions for the model are deduced, then, the annual change patterns of the present-day stress-strain field of the Chinese mainland are simulated. The results show that (1) the general pattern of the recent tectonic deformation in the Chinese mainland is governed by the interactions of its surrounding plates, of which, the Indian Plate plays a major role. There is a NNE-directed velocity distribution in the west of the Chinese mainland. The maximum slip rate appears at the collision boundary. The north-directed components decrease, while the east-directed components increase gradually from south to north and from west to east. In the east part, there is a general east-directed movement, with a certain amount of south-directed components. (2) The present-day tectonic stress field in the Chinese mainland has undergone the process of enhancement in recent years, and this process presents a general pattern of radiating eastwards from the Qinghai-Xizang (Tibet) Plateau as the center. The general pattern is similar to the ambient tectonic stress field, indicating the inheritance of contemporary tectonic deformation on the Chinese mainland. (3) The maximum principal strain presents an obvious pattern of being high in the west and low in the east. The tectonic movement in the west is stronger than that in the east. Large active faults are all located in the high-value zones of maximum principal strain. However, the magnitude of strain is smaller in the interior of the active crustal blocks, which are enclosed by these faults. (4) The stress-strain field of the Sichuan-Yunnan region is unique. It may not be governed by collision of plates alone but a combination of the movement of peripheral active blocks, material flow in the lower crust or upper mantle and special tectonic geometry (such as the eastern Himalayan syntax) as well.

Modeling Freight Truck Trips in Birmingham Using Tour-Based Approach

This research study explores the use of an innovative freight tour-based approach to model truck trips as an alternative to the conventional trip-based approach. The tour-based approach is more realistic as it captures the intermediate stops of each truck and reflects the implications of those stops on vehicle miles traveled (VMT). The paper describes the truck tour-based model concept, and presents the framework of a truck tour-based travel demand forecasting approach. As a case study, Global Positioning System (GPS) truck data are used to determine origin, destination, and truck stops for trucks moving within the Birmingham, Alabama region. Such information is then utilized to model truck movements within the study region as individual truck tours. The tour-based model is ran, and the resulting performance measures are contrasted to those obtained from the conventional trip-based planning model used by the Regional Planning Commission of Greater Birmingham (RPCGB). This case study demonstrates the feasibility of using a tour-based freight demand forecasting model as an alternative to the conventional 4-step process currently used to estimate truck trips in the Birmingham region. The results and lessons learned from the Birmingham case study are expected to improve truck movement modeling practices in the region and advance the accuracy of truck travel demand forecasting models at other locations in the future.

Earth orientation parameters from VLBI determined with a Kalman filter

This paper introduces the reader to our Kalman filter developed for geodetic VLBI（very long baseline interferometry） data analysis. The focus lies on the EOP（Earth Orientation Parameter） determination based on the Continuous VLBI Campaign 2014（CONT14） data, but also earlier CONT campaigns are analyzed. For validation and comparison purposes we use EOP determined with the classical LSM（least squares method） estimated from the same VLBI data set as the Kalman solution with a daily resolution. To gain higher resolved EOP from LSM we run solutions which yield hourly estimates for polar motion and dUTl = Universal Time（UT1）-Coordinated Universal Time（UTC）. As an external validation data set we use a GPS（Global Positioning System） solution providing hourly polar motion results.Further, we describe our approach for determining the noise driving the Kalman filter. It has to be chosen carefully, since it can lead to a significant degradation of the results. We illustrate this issue in context with the de-correlation of polar motion and nutation.Finally, we find that the agreement with respect to GPS can be improved by up to 50% using our filter compared to the LSM approach, reaching a similar precision than the GPS solution. Especially the power of erroneous high-frequency signals can be reduced dramatically, opening up new possibilities for highfrequency EOP studies and investigations of the models involved in VLBI data analysis.We prove that the Kalman filter is more than on par with the classical least squares method and that it is a valuable alternative, especially on the advent of the VLBI2010 Global Observing System and within the GGOS frame work.