Precise orbit determination for the GRACE mission using only GPS data

The GRACE (gravity recovery and climate experiment) satellites, launched in March 2002, are each equipped with a BlackJack GPS onboard receiver for precise orbit determination and gravity field recovery. Since launch, there have been significant improvements in the background force models used for satellite orbit determination, most notably the model for the geopotential. This has resulted in significant improvements to orbit accuracy for very low altitude satellites. The purpose of this paper is to investigate how well the orbits of the GRACE satellites (about 470 km in altitude) can currently be determined using only GPS data and based on the current models and methods. The orbit accuracy is assessed using a number of tests, which include analysis of orbit fits, orbit overlaps, orbit connecting points, satellite Laser ranging residuals and K-band ranging (KBR) residuals. We show that 1-cm radial orbit accuracy for the GRACE satellites has probably been achieved. These precise GRACE orbits can be used for such purposes as improving gravity recovery from the GRACE KBR data and for atmospheric profiling, and they demonstrate the quality of the background force models being used.     

Efficient satellite orbit modelling using pseudo-stochastic parameters

If the force field acting on an artificial Earth satellite is not known a priori with sufficient accuracy to represent its observations on their accuracy level, one may introduce so-called pseudo-stochastic parameters into an orbit determination process, e.g. instantaneous velocity changes at user-defined epochs or piecewise constant accelerations in user-defined adjacent time subintervals or piecewise linear and continuous accelerations in adjacent time subintervals. The procedures, based on standard least-squares, associated with such parameterizations are well established, but they become inefficient (slow) if the number of pseudo-stochastic parameters becomes large. We develop two efficient methods to solve the orbit determination problem in the presence of pseudo-stochastic parameters. The results of the methods are identical to those obtained with conventional least-squares algorithms. The first efficient algorithm also provides the full variance–covariance matrix; the second, even more efficient algorithm, only parts of it.     

Research on Autonomous Orbit Determination of Navigation Satellite Based on Crosslink Range and Orientation Parameters Constraining

Introduction Thepotentialvulnerabilityofsatellitenaviga tionsystemthatreliesongroundstationsisthat thesystemwouldbreakdownifgroundstations weredestroyed,whichcannotmeettherequire mentofnavigationwarfare[1].Withthedevelop mentofsuchspace basedsystemsasgrou…     

PM10 and heavy metal measurements in an industrial area of southern Italy

Atmospheric particulate concentrations and heavy metal content are measured from March to July 2001 at an industrial site located in a rural zone of the southern Italy. PM₁₀ samples are collected by a low-volume sampler and each sample is analysed by AAS techniques for its content of Cd, Cr, Cu, Fe, Mn, Ni, Pb and Zn. We measure also temperature, atmospheric pressure and relative humidity, and we collect anemometric data. The study purpose is the investigation of pollutant levels in an industrial area located in a rather unpolluted region and the characterization of the correlation structure among particulate concentrations, heavy metal content and local meteorological parameters. Data analysis is carried out by means of univariate and multivariate statistical methods. In the investigated period, the average value of PM₁₀ daily concentrations (24 μg Nm^− 3) does not exceed the national standard of 40 μg Nm^− 3 and only nine values are higher than the European daily limit value of 50 μg Nm^− 3. Particularly, the occurrence of two anomalous values (183 μg Nm^− 3 in 3 March and 94 μg Nm^− 3 in 22 June) seems to be related to no-local events as confirmed both by in situ data measured in the AQM network of Potenza city (about 10 km far from the study area) and by remote measurements performed in the same days. Regarding the heavy metal levels, we observe high levels of Cr (34 ng Nm^− 3), Ni (85 ng Nm^− 3) and Zn (214 ng Nm^− 3) in agreement with the local industrial source pattern. The multivariate analysis, carried out using meteorological parameters as exogenous variables, allow to evaluate the role of the different variables as driving factors of the correlation structure among the metals.     

Genetic Programming Approach for Predicting Surface Subsidence Induced by Mining

INTRODUCTIONThe study of surface subsidence resulting fromunderground mining holds significance in view ofbeneficial effects on national economy and people’sliving conditions.As is known,a large quantity ofminerals need to be exploited for modernization con-struction.Mining leads to surface subsidence,someunderground engineering and production facilities aredestroyed accordingly,which blocks the developmentof production and retards people’s living standard.To solve this contradiction,su…     

Towards urbanisation of the non-hydrostatic numerical weather prediction model Lokalmodell (LM)

Numerical weather prediction models are increasingly employed for providing meteorological data for urban air quality applications. Model resolution, physiographic parameters and surface-layer parameterisations need to be adapted to the requirements of the urban boundary layer. The Lokalmodell of the German Weather Service was triple-nested down to a horizontal grid resolution of 1.1 km, urbanised physiographic parameters were implemented, and an additional anthropogenic heat source was introduced. Results of a sensitivity study for a spring dust episode in Helsinki show a clear urbanisation effect of these measures on temperature, humidity and the partitioning of surface fluxes, leading to an increased Bowen ratio and heat storage and an urban heat island effect.     

单纯形—模拟退火混合算法反求水文地质参数及其并行求解

本文利用单纯形法局部搜索速度快和模拟退火算法全局寻优的特点,同时为了克服各自算法的弊病,提出采用单纯形—模拟退火混合算法(SMSA)进行水文地质逆问题的求解。论文详细描述了SMSA算法的具体操作算子的实现,并将该算法应用于一个大型水源地的地下水模拟反演。计算结果表明,SMSA算法在水文地质参数反演计算具有求解速度快,精度高的特点,而且易于实现并行运算。     

单纯形一模拟退火混合算法反求水文地质参数及其并行求解

本文利用单纯形法局部搜索速度快和模拟退火算法全局寻优的特点，同时为了克服各自算法的弊病，提出采用单纯形一模拟退火混合算法（SMSA）进行水文地质逆问题的求解。论文详细描述了SMSA算法的具体操作算子的实现，并将该算法应用于一个大型水源地的地下水模拟反演。计算结果表明，SMSA算法在水文地质参数反演计算具有求解速度快，精度高的特点，而且易于实现并行运算。     

三维地震勘探激发参数选择与分析

文中介绍了山区进行三维地震勘探的重要性及其方法原理。以山西晋城煤业集团某勘探区为实例进行实地勘探和试验,并对试验结果进行了详细的分析,进一步验证了该方法的可行性。     

Source parameters of intermediate-depth Vrancea (Romania) earthquakes from empirical Green's functions modeling

Several source parameters (source dimensions, slip, particle velocity, static and dynamic stress drop) are determined for the moderate-size October 27th, 2004 (M_W = 5.8), and the large August 30th, 1986 (M_W = 7.1) and March 4th, 1977 (M_W = 7.4) Vrancea (Romania) intermediate-depth earthquakes. For this purpose, the empirical Green's functions method of Irikura [e.g. Irikura, K. (1983). Semi-Empirical Estimation of Strong Ground Motions during Large Earthquakes. Bull. Dis. Prev. Res. Inst., Kyoto Univ., 33, Part 2, No. 298, 63–104., Irikura, K. (1986). Prediction of strong acceleration motions using empirical Green's function, in Proceedings of the 7th Japan earthquake engineering symposium, 151–156., Irikura, K. (1999). Techniques for the simulation of strong ground motion and deterministic seismic hazard analysis, in Proceedings of the advanced study course seismotectonic and microzonation techniques in earthquake engineering: integrated training in earthquake risk reduction practices, Kefallinia, 453–554.] is used to generate synthetic time series from recordings of smaller events (with 4 ≤ M_W ≤ 5) in order to estimate several parameters characterizing the so-called strong motion generation area, which is defined as an extended area with homogeneous slip and rise time and, for crustal earthquakes, corresponds to an asperity of about 100 bar stress release [Miyake, H., T. Iwata and K. Irikura (2003). Source characterization for broadband ground-motion simulation: Kinematic heterogeneous source model and strong motion generation area. Bull. Seism. Soc. Am., 93, 2531–2545.] The parameters are obtained by acceleration envelope and displacement waveform inversion for the 2004 and 1986 events and MSK intensity pattern inversion for the 1977 event using a genetic algorithm. The strong motion recordings of the analyzed Vrancea earthquakes as well as the MSK intensity pattern of the 1977 earthquake can be well reproduced using relatively small strong motion generation areas, which corresponds to small asperities with high stress drops (300–1200 bar) and high particle velocities (3–5 m/s). These results imply a very efficient high-frequency radiation, which has to be taken into account for strong ground motion prediction, and indicate that the intermediate-depth Vrancea earthquakes are inherently different from crustal events.