First laser measurements to space debris in Poland

The Borowiec Satellite Laser Ranging station (BORL 7811, Borowiec) being a part of the Space Research Centre of the Polish Academy of Sciences (SRC PAS) went through modernization in 2014–2015. One of the main tasks of the modernization was the installation of a high-energy laser module dedicated to space debris tracking. Surelite III by Continuum is a Nd:YAG pulse laser with 10?Hz repetition rate, a pulse width of 3–5?ns and a pulse energy of 450?mJ for green (532?nm). This new laser unit was integrated with the SLR system at Borowiec performing standard satellite tracking. In 2016 BORL 7811 participated actively to the observational campaigns related to the space debris targets from LEO region managed by the Space Debris Study Group (SDSG) of the International Laser Ranging Service (ILRS).Currently, Borowiec station regularly tracks 36 space debris from the LEO regime, including typical rocket bodies (Russian/Chinese) and cooperative targets like the inactive TOPEX/Poseidon, ENVISAT, OICETS and others. In this paper the first results of space debris laser measurements obtained by the Borowiec station in period August 2016 – January 2017 are presented. The results gained by the SRC PAS Borowiec station confirm the rotation of the defunct TOPEX/Poseidon satellite which spins with a period of approximately 10?s. The novelty of this work is the presentation of the sample results of the Chinese CZ-2C R/B target (NORAD catalogue number 31114) which is equipped (probably) with retroreflectors. Laser measurements to space debris is a very desirable topic for the next years, especially in the context of the Space Surveillance and Tracking (SST) activity. Some targets are very easy to track like defunct ENVISAT or TOPEX/Poseidon. On the other hand, there is a big population of different LEO targets with different orbital and physical parameters, which are challenging for laser ranging like small irregular debris and rocket boosters.     

New improved orbit solutions for the ERS-1 and ERS-2 satellites

Improved orbit solutions of the European Remote Sensing Satellites ERS-1 and ERS-2 have been computed in the ITRF2005 terrestrial reference frame using the recent models based mainly on IERS Conventions 2003. These solutions cover the periods 3 August 1991 to 8 July 1996 for ERS-1, and 3 May 1995 to 4 July 2003 for ERS-2. For each satellite, the final orbit solution is based on a combination of three separate orbit solutions independently computed at the Delft Institute of Earth Observation and Space Systems (DEOS) of the Delft University of Technology (The Netherlands), the Navigation Support Office of the European Space Operations Centre (ESOC, Germany) and the Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences (Germany) using three different software packages for precise orbit determination, but using the same models in the same terrestrial reference frame within the European Space Agency (ESA) project ‘Reprocessing of Altimeter Products for ERS (REAPER)’. Validation using radar altimeter data indicates that the new combined orbits of ERS-1 and ERS-2 computed by us are significantly more accurate, approaching the 2–3 cm level in radial direction, than previously available orbit solutions.     

Earth and ocean tides from long-term analysis of satellite orbits

Precision laser tracking data on three close-earth satellites (Geos 1, Geos 2, and Lageos) are analyzed to obtain measurements of ocean tides, core-mantle resonance, and other mass displacements.Mean elements are computed for 6.0 years for Geos 1, 4.9 years for Geos 2, and 2.1 years for Lageos. Ocean tides T2, S2, K2, P1, K1, and φ1 are measured, along with the core-mantle resonance at the K1 and P1 frequencies, as predicted by Molodensky. Also measured are the Sa and Ssa tides and radiation-pressure area-to-mass ratios.     

High accuracy satellite drag model (HASDM)

The dominant error source in force models used to predict low-perigee satellite trajectories is atmospheric drag. Errors in operational thermospheric density models cause significant errors in predicted satellite positions, since these models do not account for dynamic changes in atmospheric drag for orbit predictions. The Air Force Space Battlelab’s High Accuracy Satellite Drag Model (HASDM) estimates and predicts (out three days) a dynamically varying global density field. HASDM includes the Dynamic Calibration Atmosphere (DCA) algorithm that solves for the phases and amplitudes of the diurnal and semidiurnal variations of thermospheric density near real-time from the observed drag effects on a set of Low Earth Orbit (LEO) calibration satellites. The density correction is expressed as a function of latitude, local solar time and altitude. In HASDM, a time series prediction filter relates the extreme ultraviolet (EUV) energy index E_10.7 and the geomagnetic storm index a_p, to the DCA density correction parameters. The E_10.7 index is generated by the SOLAR2000 model, the first full spectrum model of solar irradiance. The estimated and predicted density fields will be used operationally to significantly improve the accuracy of predicted trajectories for all low-perigee satellites.     

Inconsistency of magnetic field and plasma velocity variations in the distant plasma sheet: Violation of the “frozen-in” criterion?

Measurements of the magnetic field and low energy plasma by the GEOTAIL spacecraft have been used to study the relationship between variations of the plasma velocity and of the magnetic field in the distant (100–200 R_E) and middle (40–80 R_E) tail. The analysis was carried out separately for the tail lobes and the plasma sheet. It is shown that the absolute values of the magnetic field and plasma velocity, as well as their corresponding components (V_X and B_X, V_Y and B_Y, V_Z and B_Z), are linearly connected in the tail lobes. In the plasma sheet, however, the plasma velocity and the magnetic field do not seem to be related to one another. The distant plasma sheet seems to be in a regime of turbulence. The diffusion coefficients estimated from our data set of the velocity parameters in the plasma sheet are in good agreement with the theoretical predictions of Antonova and Ovchinnikov (1996, 1999).     

INTERBALL statistical study of ion flow fluctuations in the plasma sheet

We use ion distribution measurements with CORALL instrument on-board the INTERBALL/Tail spacecraft to study plasma flows in the mid-tail (−9> X> −27 R_E) plasma sheet. Three velocity components computed every 2 minutes exhibit two types of velocity variations: Earthward bursty bulk flows (BBFs) and random flow fluctuations. Their properties are in a good agreement with the observations of the ISEE-2 spacecraft (Borovsky et al., 1997). The INTERBALL/Tail spacecraft configuration favors measurements of V_z component, in contrast to previous experiments in which only V_x and V_y were measured reliably. In the outer part of the plasma sheet V_y and V_z fluctuations were close to each other (variances σ(V_y) and σ(V_z) were about 160 and 110 km/s, respectively), but in the inner part at the dusk flank amplitude of V_y fluctuations increased and was 2 times higher than that of V_z component. This asymmetry of fluctuations should be taken into account during modern theoretical analysis and simulations.     

GRACE gravity field modeling with an investigation on correlation between nuisance parameters and gravity field coefficients

The GRACE (Gravity Recovery And Climate Experiment) monthly gravity models have been independently produced and published by several research institutions, such as Center for Space Research (CSR), GeoForschungsZentrum (GFZ), Jet Propulsion Laboratory (JPL), Centre National d’Etudes Spatiales (CNES) and Delft Institute of Earth Observation and Space Systems (DEOS). According to their processing standards, above institutions use the traditional variational approach except that the DEOS exploits the acceleration approach. The background force models employed are rather similar. The produced gravity field models generally agree with one another in the spatial pattern. However, there are some discrepancies in the gravity signal amplitude between solutions produced by different institutions. In particular, 10%–30% signal amplitude differences in some river basins can be observed. In this paper, we implemented a variant of the traditional variational approach and computed two sets of monthly gravity field solutions using the data from January 2005 to December 2006. The input data are K-band range-rates (KBRR) and kinematic orbits of GRACE satellites. The main difference in the production of our two types of models is how to deal with nuisance parameters. This type of parameters is necessary to absorb low-frequency errors in the data, which are mainly the aliasing and instrument errors. One way is to remove the nuisance parameters before estimating the geopotential coefficients, called NPARB approach in the paper. The other way is to estimate the nuisance parameters and geopotential coefficients simultaneously, called NPESS approach. These two types of solutions mainly differ in geopotential coefficients from degree 2 to 5. This can be explained by the fact that the nuisance parameters and the gravity field coefficients are highly correlated, particularly at low degrees. We compare these solutions with the official and published ones by means of spectral analysis. It is found that our solutions are, in general, consistent with others in the spatial pattern. The water storage variations of the Amazon, Chari and Ganges river basins have also been computed. The variations computed with the NPARB approach are closer to those produced by JPL and DEOS solutions, while the variations produced with the NPESS approach are in good agreement with those produced by the CSR and GFZ solutions. A simulation study is implemented with considering realistic noise and low-frequency error. The two approaches are used to recover the true model. The NPESS solution appears closer to the true one. Therefore we are inclined to estimate the nuisance parameters simultaneously with the geopential coefficients.     

两种用于NS-DSMC耦合方法的连续失效参数对比

NS-DSMC（Navier Stokes-Direct Simulation Monte Carlo）耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参数,对目前广泛使用的两种连续失效参数Kn_Q和B参数进行了理论上的分析和比较,表明虽然他们建立的出发点不同,但在数学形式上具有一定相似性.通过圆柱绕流问题的数值试验进一步验证了两种连续失效参数具有一定对应关系.Kn_Q和B两种连续失效参数在NS-DSMC耦合方法中应用效果相差不大,但Kn_Q 的阈值0.05适用范围较广,而B参数的阈值对不同流动问题会有变化.      

Plasma convection in the low-latitude mantle of the near-earth magnetosphere

A detailed analysis of Geotail observations on June 1, 1997 provides evidence for the appearance of the plasma mantle at low-latitudes, Z_GSM 3 R_E for southward interplanetary magnetic field (IMF), and clarifies the characteristics of the plasma flow in the low-latitude mantle. The tailward flow speed in the mantle is 80 160 km/s, smaller than the flow speed (300 km/s) in the adjacent magnetosheath. The component of the velocity perpendicular to the magnetic field lines is predominantly poleward at speeds up to 60 km/s, and the convection tends to be stronger for |;B_Y|<|;B_Z| than for IMF |B_Y|>|B_Z|;. This polewardconvection can be interpreted as being driven by a J × B force operating on the field lines.     

Numerical simulation of three-dimensional reconnection due to the instability of collisionless current sheets

Based on analytical calculations we have currently argued that spontaneous reconnection through thin collisionless current sheets is an essentially three-dimensional (3 D) process (Büchner, 1996 a, b). Since 3 D kinetic PIC codes have become available, the three dimensional nature of the collisionless current sheet decay are now illustrated by numerical simulations (Büchner and Kuska, 1996; Pritchett and Coroniti, 1996; Zhu and Winglee, 1996). While the latter two claim a coupling to a longer wavelength kink mode as a main factor, destabilizing thin current sheets in 3 D, our simulations have revealed that even shorter scale perturbations in the current direction suffice to destabilize thin sheets very quickly. Since past simulation runs, however, were limited to mass ratios near unity, the influence of the electrons was not treated adequately. We have now investigated the stability of thin collisionless current sheets including 64 times lighter negatively charged particles. We can now show that while the two-dimensional tearing instability slows down for M = M_p/m_e = 64, the three-dimensional current sheet decay is a much faster process — practically as fast as the mass ratio M = 1 3 D sheet decay, even without kinking the sheet. We further conclude that, unlike the two-dimensional tearing instability, the three-dimensional decay of thin current sheets is not controlled by the electrons. For a sheet width comparable with the ion inertial length, we also recovered signatures of the Hall effect as predicted by Vasyliunas (1975) in the mass ratio M = 64 case. The ion inertial length seems to be the critical scale at which the sheet starts to decay.     

利用"火星快车"三程多普勒跟踪数据限定局部洛伦兹不变性

目前航天器的三程多普勒跟踪技术已经在深空探测的控制与导航领域起到了重要作用。利用包含了对局部洛伦兹不变性(LLI)以及局部位置不变性(LPI)原理有破坏的三程多普勒跟踪理论,研究分析了"火星快车"(MEX)三程多普勒跟踪数据的残差。这些多普勒观测于2009年8月7日和8日进行,利用了欧洲航天局(ESA)在澳大利亚新诺舍(New Norcia)的上行站和三个分别在中国上海、昆明以及乌鲁木齐的下行站。我们发现,这些观测结果给出的LLI上限在10-2的量级。但由于各观测站本身对频率测量的精度有限,这些数据并不适合于检验LPI。     

子午工程首次火箭探空数据准单色惯性重力波特性分析

利用2010年6月3日子午工程首次气象火箭探测的温度和风场数据,采用矢端曲线法分别从平流层（20~50km）和对流层（0~15km）廓线提取了海南火箭发射场上空准单色惯性重力波参数.火箭探测的平流层和对流层两个准单色惯性重力波分别向上和向下逆风传播,固有周期为20.1h和22.4h,垂直波长为9.5km和4.0km,水平波长为2900km和753km,垂直群速度c_gz为0.0887m·-1和0.0298m·-1,水平群速度c_gh为12.7m·-1和3.65m·-1,λ_h/λ_z为305:1和188:1,c_gh/c_gz为143:1和122:1.      

The variability and predictability of the IRI B0, B1 parameters over Grahamstown,South Africa

The International Reference Ionosphere (IRI) parameters B₀ and B₁ provide a representation of the thickness and shape, respectively, of the F2 layer of the bottomside ionosphere. These parameters can be derived from electron density profiles that are determined from vertical incidence ionograms. This paper aims to illustrate the variability of these parameters for a single mid latitude station and demonstrate the ability of the Neural Network (NN) modeling technique for developing a predictive model for these parameters. Grahamstown, South Africa (33.3°S, 26.5°E) was chosen as the mid latitude station used in this study and the B₀ and B₁ parameters for an 11 year period were determined from electron density profiles recorded at that station with a University of Massachusetts Lowell Center for Atmospheric Research (UMLCAR) Digisonde. A preliminary single station NN model was then developed using the Grahamstown data from 1996 to 2005 as a training database, and input parameters known to affect the behaviour of the F2 layer, such as day number, hour, solar and magnetic indices. An analysis of the diurnal, seasonal and solar variations of these parameters was undertaken for the years 2000, 2005 and 2006 using hourly monthly median values. Comparisons between the values derived from measured data and those predicted using the two available IRI-2001 methods (IRI tables and Gulyaeva, T. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7(6), 39–48, 1987.) and the newly developed NN model are also shown in this paper. The preliminary NN model showed that it is feasible to use the NN technique to develop a prediction tool for the IRI thickness and shape parameters and first results from this model reveal that for the mid latitude location used in this study the NN model provides a more accurate prediction than the current IRI model options.     

Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

LMXB luminosity function and connection to globular cluster in early type galaxies

We derive bias-corrected X-ray luminosity functions (XLFs) of LMXBs detected in 14 E and S0 galaxies observed with Chandra. After correcting for incompleteness, the individual XLFs are statistically consistent with a single power-law. A break at or near L_X,Eddington , as previously reported, is not required in any individual case. The combined XLF with a reduced error, however, suggests a possible break at L_X = 5 × 10³⁸ erg s⁻¹, which may be consistent with the Eddington luminosity of neutron stars with the largest possible mass (3 M), or of He-enriched neutron star binaries. We confirm that the total X-ray luminosity of LMXBs is correlated with the the near-IR luminosities, but the scatter exceeds that expected from measurement errors. The scatter in L_X(LMXB)/L_K appears to be correlated with the specific frequency of globular clusters, as reported earlier.

We cross-correlate X-ray binaries with globular clusters determined by ground-based optical and HST observations in 6 giant elliptical galaxies. With the largest sample reported so far (300 GC LMXBs with a 5:2 ratio between red and blue GCs), we compare their X-ray properties, such as X-ray hardness, XLF and L_X/L_B and find no statistically significance difference between different groups of LMXBs. Regardless of their association with GCs, both GC and field LMXBs appear to follow the radial profile of the optical halo light, rather than that of more extended GCs. This suggests that while metallicity is a primary factor in the formation of LMXBs in GCs, there may be a secondary factor (e.g., encounter rate) playing a non-negligible role.     

Auroral signatures of transient processes in the outer magnetosphere

We present an analysis of sporadic and recurrent injections of magnetospheric ions in the midnight auroral oval during substorms and of the associated ionospheric ion outflows. The source of plasma sheet precipitating ions is determined using a simple method, based on the measured relation between the ion inverse velocity and time (l = v × t). This method is applied here to two typical passes of the Interball-Auroral (IA) satellite at distances of 3 R_E above the auroral regions. Substorm related ion injections are shown to be mainly due to time of flight effects. In contrast with particle trajectory computations (Sauvaud et al., 1999), the inverse velocity method does not require magnetic and electric field models and can thus be used systematically for the detection of time of flight dispersed ion structures (TDIS). This allowed us to build a large database of TDIS events and to perform a statistical analysis of their spatial distribution. For the cases presented here the source region of the injected ions is found at radial distances from 18 to 30 R_E near the equatorial magnetosphere. At Interball altitudes ( 3 R_E), ion injections detected at the poleward boundary of the nighside auroral oval are associated with shear Alfvén waves superimposed over large-scale quasi-static current structures. We show that the most poleward TDIS are collocated with a large outflow of ionospheric H⁺ and O⁺ displaying pitch-angle distributions peaked in the pitch-angle range 90°–120°. These ions are thus accelerated perpendicularly to the magnetic field not only in the main auroral acceleration region but also up to at least 3 R_E. The expanding auroral bulge thus constitutes a significant source of H⁺ and O⁺ ions for the mid-tail magnetosphere.     

DORIS processing at the European Space Operations Centre

This paper gives an overview of the DORIS related activities at the Navigation Support Office of the European Space Operations Centre. The DORIS activities were started in 2002 because of the launch of the Envisat satellite where ESOC is responsible for the validation of the Envisat Precise Orbits and a brief overview of the key Envisat activities at ESOC is given. Typical orbit comparison RMS values between the CNES POE (GDR-C) and the ESOC POD solution is 6.5, 18.8 and 23.1 mm in radial-, along- and cross-track direction. In the framework of the generation of the ITRF2008 ESOC participated in the reprocessing of all three space geodetic techniques; DORIS, SLR, and GPS. Here the main results of our DORIS reprocessing, in the framework of the International DORIS Service (IDS), are given. The WRMS of the weekly ESOC solution (esawd03) for the 2004–2009 period compared to the IDS-1 combined solution is of the order of 12 mm. Based on the long time series of homogeneously processed data a closer look is taken at the estimated solar radiation pressure parameters of the different satellites used in this DORIS analysis. The main aim being the stabilization of the Z-component of the geocentre estimates. We conclude that the ESOC participation to the IDS ITRF2008 contribution has been beneficial for both ESOC and the IDS. ESOC has profited significantly from the very open and direct communications and comparisons that took place within the IDS during the reprocessing campaign.     

Progress towards a new global foF2 model for the International Reference Ionosphere (IRI)

A new neural network (NN) based global empirical model for the foF2 parameter, which represents the peak electron density has been developed using extended temporal and spatial geophysical relevant inputs. The first results from this new model were presented at the International Reference Ionosphere (IRI) 2006 workshop in Buenos Aires, Argentina, and showed that this new model would be a suitable replacement for the URSI and CCIR maps currently used within the IRI model for the purpose of F2 peak electron density predictions. Measured ground based ionosonde data, from 85 global stations, spanning the period 1995–2005 and, for a few stations from 1976 to 1986, obtained from various resources of the World Data Centre (WDC) archives (Space Physics Interactive Data Resource SPIDR, the Digital Ionogram Database, DIDBase, and IPS Radio and Space Services) have been used for training a NN.     

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一,国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据,用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化,对比分析了IRI-2016与IRI-2012模型的预测结果,同时分析了IRI-2016模型输出H_mF₂的性能.结果表明,IRI模型在中高纬度地区的输出结果高于COSMIC反演结果,而赤道及低纬地区则大都偏低.与IRI-2012模型相比,IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高,但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.      

Determination of orbits with Torun Orbit Processor system

The Torun Orbit Processor system (TOP) is a multipurpose system designed to provide the research requirements of the various orbital problems for artificial satellites. The analysis of mathematical model for the TOP system began in 1991 and coding for the IBM PC and Sun SparcStation 10 computers was completed in 1993. TOP was first used for laser data processing of Lageos satellite and for simulations of satellite-to-satellite tracking technique (SST). This paper gives an overview of TOP system and highlights some problems connected with orbit recovery.