Low Earth orbit constellation design using a multi-objective genetic algorithm for GNSS reflectometry missions

Spaceborne global navigation satellite system reflectometry (GNSS-R) is an innovative bistatic radar remote sensing technique utilizing low Earth orbit (LEO) based GNSS-R instruments to acquire GNSS L-band opportunistic signals for measuring geophysical parameters. A GNSS-R LEO constellation with an optimization design for its specialized missions is very significant and necessary. However, the constellation design involves multi-parameter and multi-objective optimization, and the classical analytic solution is not capable of such a complicated issue. This study proposes a multi-objective LEO constellation design method with a genetic algorithm (GA) and presents a framework for designing two GNSS-R LEO constellations, termed “lower-latitude constellation” for typhoons and hurricanes observation in the tropics and “global constellation” for global geophysical parameter measurements. Then, the observation capability of both designed constellations is evaluated in terms of the number of reflection points, spatial coverage density, and revisit time to verify the GA efficiency in LEO constellation design. Results show that the two designed LEO constellations with high fitness function values possess optimal orbit parameter set configuration and outperform the existing CyGNSS constellations in observation performance. Compared with CyGNSS, the number of reflection points observed by the lower-latitude constellation and the global constellation increases by 38% and 45%, as well as the spatial coverage density increases by 28% and 36%. The revisit time for the lower-latitude constellation is reduced by 0.29 h, whereas the revisit time for the global constellation increases by one hour.     

LEO enhanced Global Navigation Satellite System (LeGNSS) for real-time precise positioning services

Global Navigation Satellite System (GNSS) has been widely used in many geosciences areas with its Positioning, Navigation and Timing (PNT) service. However, GNSS still has its own bottleneck, such as the long initialization period of Precise Point Positioning (PPP) without dense reference network. Recently, the concept of PNTRC (Positioning, Navigation, Timing, Remote sensing and Communication) has been put forward, where Low Earth Orbit (LEO) satellite constellations are recruited to fulfill diverse missions. In navigation aspect, a number of selected LEO satellites can be equipped with a transmitter to transmit similar navigation signals to ground users, so that they can serve as GNSS satellites but with much faster geometric change to enhance GNSS capability, which is named as LEO constellation enhanced GNSS (LeGNSS). As a result, the initialization time of PPP is expected to be shortened to the level of a few minutes or even seconds depending on the number of the LEO satellites involved. In this article, we simulate all the relevant data from June 8th to 14th, 2014 and investigate the feasibility of LeGNSS with the concentration on the key issues in the whole data processing for providing real-time PPP service based on a system configuration with fourteen satellites of BeiDou Navigation Satellite System (BDS), twenty-four satellites of the Global Positioning System (GPS), and sixty-six satellites of the Iridium satellite constellations. At the server-end, Precise Orbit Determination (POD) and Precise Clock Estimation (PCE) with various operational modes are investigated using simulated observations. It is found out that GNSS POD with partial LEO satellites is the most practical mode of LeGNSS operation. At the user-end, the Geometry Dilution Of Precision (GDOP) and Signal-In-Space Ranging Error (SISRE) are calculated and assessed for different positioning schemes in order to demonstrate the performance of LeGNSS. Centimeter level SISRE can be achieved for LeGNSS.     

Minimum of PDOP and its applications in inter-satellite links (ISL) establishment of Walker-δ constellation

Within the next decade, there will be a number of GNSS (Global Navigation Satellite System) available, i.e. modernized GPS, Galileo, restored GLONASS, BeiDou and many other regional GNSS augmentation systems. Thus, measurement redundancies and geometry of the satellites can be improved. GDOP (Geometric Dilution of Precision) and PDOP (Position Dilution of Precision) are associated with the constellation geometry of satellites, and they are the geometrically determined factors that describe the effect of geometry on the relationship between measurement error and position error. GDOP and PDOP are often used as standards for selecting good satellites to meet the desired positioning precision. In this paper, the related conclusions of minimum of GDOP which was discussed are given, and it is used to study the minimum of PDOP for two cases that the receiver is on the earth’s surface and the receiver is on satellite. The corresponding theorem and constructive solutions of minimum of PDOP are given. Then, the rationality of the ISL (inter-satellite link) establishment criteria in Walker-δ constellation is discussed by using the theory of minimum of PDOP. Finally, the minimum of PDOP is calculated when the number of satellites is 4–10, and these results are verified by using Monte Carlo method.     

Combining GPS and GLONASS for time and frequency transfer

Geodetic time and frequency transfer (TFT) consists in a comprehensive modeling of code and carrier phase observations from Global Navigation Satellite System (GNSS) in order to determine the synchronization errors between two remote clocks connected to GNSS receivers. Using either common view (CV), or Precise Point Positioning (PPP), current GNSS time transfer uses only GPS measurements. This study combines GPS and GLONASS observations in geodetic TFT in order to determine the added value of the GLONASS data in the results. Using the software Atomium, we demonstrate on one hand that using both constellations improves the solution for both CV and PPP results when analysing short data batches. In that case, there are not enough GPS code data to calibrate the solution, and additional GLONASS code data allows us to retrieve a correct absolute value for the solution. On the other hand, the CV results of frequency transfer are not significantly affected by adding GLONASS data, while in PPP the combination with GLONASS modifies the frequency transfer results, and in particular the daily frequency offset, with maximum differences of 150 ps between the TFT solutions obtained with GPS-only or GPS + GLONASS.     

Iridium/ORBCOMM机会信号融合定位技术

在全球导航卫星系统(GNSS)不可用情况下，低地球轨道(LEO)卫星机会信号(SOP)定位技术是一种有效的导航定位解决方案。单LEO星座机会信号定位技术面临星座构型不足或可见卫星偏少等问题，多LEO星座机会信号融合定位技术可有效解决该问题。通过分析瞬时多普勒定位原理，建立了Iridium/ORBCOMM机会信号融合定位模型，引入基于Helmert方差估计的加权最小二乘算法进一步提高定位精度。实测数据表明:基于Helmert方差估计的Iridium/ORBCOMM机会信号融合定位精度优于70 m，验证了多LEO星座机会信号融合定位的可行性和有效性。      

Ambiguity resolution performance with GPS and BeiDou for LEO formation flying

The evolving BeiDou Navigation Satellite System constellation brings new opportunities for high-precision applications. In this contribution the focus will be on one such application, namely precise and instantaneous relative navigation of a formation of LEO satellites. The aim is to assess the ambiguity resolution performance with the future GPS and BeiDou constellations depending on system choice (GPS, BeiDou, or GPS+BeiDou), single- or dual-frequency observations, receiver noise, and uncertainties in ionosphere modelling. In addition, for the GPS+BeiDou constellation it will be shown how the growing BeiDou constellation in the years to come can already bring an important performance improvement compared to the GPS-only case. The performance will be assessed based on the percentage of time that the required precision can be obtained with a partial ambiguity resolution strategy.     

In-orbit assessment of laser retro-reflector efficiency onboard high orbiting satellites

The navigation and geodetic satellites that orbit the Earth at altitudes of approximately 20,000 km are tracked routinely by many of the Satellite Laser Ranging (SLR) stations of the International Laser Ranging Service (ILRS). In order to meet increasing demands on SLR stations for daytime and nighttime observations, any new mission needs to ensure a strong return signal so that the target is easily acquirable. The ILRS has therefore set a minimum effective cross-section of 100 million square metres for the on-board laser retro-reflector arrays (LRAs) and further recommends the use of ‘uncoated’ cubes in the arrays. Given the large number of GNSS satellites that are currently supported by SLR, it is informative to make an assessment of the relative efficiencies of the various LRAs employed. This paper uses the laser ranging observations themselves to deduce and then compare the efficiencies of the LRAs on the COMPASS-M1 navigation satellite, two satellites from the GPS and three from the GLONASS constellations, the two GIOVE test satellites from the upcoming Galileo constellation, the two Etalon geodetic spheres and the geosynchronous communications test satellite, ETS-8. All the LRAs on this set of satellites employ back-coated retro-reflector cubes, except those on the COMPASS-M1 and ETS-8 vehicles which are uncoated. A measure of return signal strength, and thus of LRA-efficiency, is calculated using the laser-range full-rate data archive from 2007 to 2010, scaled to remove the effects of variations in satellite range, atmospheric attenuation and retro-reflector target total surface area. Observations from five SLR stations are used in this study; they are Herstmonceux (UK), Yarragadee (Australia), Monument Peak and McDonald (USA) and Wettzell (Germany). Careful consideration is given to the treatment of the observations from each station in order to take account of local working practices and system upgrades. The results show that the uncoated retro-reflector cubes offer significant improvements in efficiency.     

精密进近阶段的多系统GNSS组合RAIM可用性算法及分析

给出了多系统全球卫星导航系统（GNSS）组合接收机自主完好性监测（ReceiverAutonomousIntegrityMonitoring,RAIM）可用性计算方法,在此基础上利用GPS、GLONASS实测数据与BDS、Galileo全星座仿真数据,分析了BDS、GPS、GLONASS和Galileo不同组合在精密进近阶段的RAIM可用性。通过试验分析发现,BDS的5颗地球同步轨道卫星和3颗倾斜地球同步轨道卫星对亚洲、非洲和欧洲大部分地区的RAIM可用性有很大的贡献。这些地区站星间几何观测结构得到改善,使得RAIM可用性相对于其他地区有很大幅度的提升。在亚太地区APV-I阶段单系统导航情况下,北斗导航系统RAIM可用性达到99.5%,高于其他三个导航系统。在精密进近阶段（APV-I、APV-II和CAT-I）,BDS与其他导航系统（GPS、GLONASS和Galileo）的组合导航可以满足全球大部分区域的RAIM可用性需求,大多可达到100%。     

低轨卫星的北斗反射信号海洋覆盖分析

北斗系统作为我国自行研制的导航系统，具有独特的混合星座特性，其反射事件在海洋上的空间分布和覆盖性能的研究成果较少。针对上述问题，模拟了北斗三号（Beidou System 3，BDS3）以及全球定位系统（Global Positioning System，GPS）空间星座，模拟不同性能参数的低轨接收卫星，在此仿真系统的基础上，计算反射事件在各仿真场景下的海洋覆盖率以及同时发生反射事件的卫星数量。结果表明：BDS3与GPS相比具有更大的反射事件海洋覆盖率，覆盖性能更优；仿真周期7天与1天相比,前者卫星产生的反射事件海洋覆盖率约为后者5倍；低轨道卫星（Low Earth Orbit,LEO）轨道高度越高，天线波束角度越大，反射事件海洋覆盖率越大；3颗LEO卫星下同时发生反射事件卫星数的平均值约为单颗LEO卫星的3倍。可通过设计多星组网来提高反射事件海洋覆盖率。研究结果对星载反射技术海洋遥感应用方面有一定的参考价值。     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations

Differential Code Bias (DCB) is an essential correction that must be provided to the Global Navigation Satellite System (GNSS) users for precise position determination. With the continuous deployment of Low Earth Orbit (LEO) satellites, DCB estimation using observations from GNSS receivers onboard the LEO satellites is drawing increasing interests in order to meet the growing demands on high-quality DCB products from LEO-based applications, such as LEO-based GNSS signal augmentation and space weather research. Previous studies on LEO-based DCB estimation are usually using the geometry-free combination of GNSS observations, and it may suffer from significant leveling errors due to non-zero mean of multipath errors and short-term variations of receiver code and phase biases. In this study, we utilize the uncombined Precise Point Positioning (PPP) model for LEO DCB estimation. The models for uncombined PPP-based LEO DCB estimation are presented and GPS observations acquired from receivers onboard three identical Swarm satellites from February 1 to 28, 2019 are used for the validation. The results show that the average Root Mean Square errors (RMS) of the GPS satellite DCBs estimated with onboard data from each of the three Swarm satellites using the uncombined PPP model are less than 0.18 ns when compared to the GPS satellite DCBs obtained from IGS final daily Global Ionospheric Map (GIM) products. Meanwhile, the corresponding average RMS of GPS satellite DCBs estimated with the conventional geometry-free model are 0.290, 0.210, 0.281 ns, respectively, which are significantly larger than those obtained with the uncombined PPP model. It is also noted that the estimated GPS satellite DCBs by Swarm A and C satellites are highly correlated, likely attributed to their similar orbit type and space environment. On the other hand, the Swarm receiver DCBs estimated with uncombined PPP model, with Standard Deviation (STD) of 0.065, 0.037 and 0.071 ns, are more stable than those obtained from the official Swarm Level 2 products with corresponding STD values of 0.115, 0.101, and 0.109 ns, respectively. The above indicates that high-quality DCB products can be estimated based on uncombined PPP with LEO onboard observations.     

GNSS-R open-loop difference phase altimetry: Results from a bridge experiment

The paper explores a method to obtain accurate lake surface heights using measurements of the Global Navigation Satellite System (GNSS) carrier phase reflected from the lake surface. The method is referred to as Global Navigation Satellite System-Reflection (GNSS-R) open-loop difference phase altimetry method. It consists of two key technologies: one is the open-loop tracking method to track the GNSS-R signals, where the direct GNSS signal’s frequency is used as a reference frequency to obtain the carrier phases of the GNSS-R signals; the other key technology is time difference phase altimetry method to invert the lake surface heights using two or more carrier phases of GNSS-R signals received simultaneously. A validation experiment is carried out on the SANYING bridge over GUANTING lake using a GNSS-R receiver developed by the Center for Space Science and Applied Research (CSSAR), processing the data with GNSS-R open-loop difference phase altimetry method. The lake surface height results are consistent with the height results of GPS dual-frequency differential positioning altimetry. The results show that we can achieve centimeter level height in one minute average, by using 11 minutes carrier phase data of three GNSS-R signals received simultaneously.     

NeQuick G model based scale factor determination for using SBAS ionosphere corrections at low earth orbit

A space-based augmentation system (SBAS) provides real-time correction data for global navigation satellite system (GNSS) users near ground. In order to use the SBAS ionosphere correction for low Earth orbit (LEO) satellites, the correction should be scaled down for the LEO altitude. This scale factor varies with ionosphere distribution and it is hard to determine the value at LEO in real time. We propose a real-time scale factor determination method by using Galileo GNSS’s NeQuick G model. A LEO satellite GPS data and SBAS data received on ground were used to evaluate the performance of the NeQuick G derived variable scale factor. The NeQuick G derived scale factor shows a significant accuracy improvement over NeQuick G model or pre-determined constant scale factor. It improves a vertical positioning accuracy of the LEO satellite. The error mean reductions of the vertical positioning over NeQuick G and the constant scale factor are 31.5% and 11.7%, respectively.     

A real-time GNSS-R system based on software-defined radio and graphics processing units

Reflected signals of the Global Navigation Satellite System (GNSS) from the sea or land surface can be utilized to deduce and monitor physical and geophysical parameters of the reflecting area. Unlike most other remote sensing techniques, GNSS-Reflectometry (GNSS-R) operates as a passive radar that takes advantage from the increasing number of navigation satellites that broadcast their L-band signals. Thereby, most of the GNSS-R receiver architectures are based on dedicated hardware solutions. Software-defined radio (SDR) technology has advanced in the recent years and enabled signal processing in real-time, which makes it an ideal candidate for the realization of a flexible GNSS-R system. Additionally, modern commodity graphic cards, which offer massive parallel computing performances, allow to handle the whole signal processing chain without interfering with the PC’s CPU. Thus, this paper describes a GNSS-R system which has been developed on the principles of software-defined radio supported by General Purpose Graphics Processing Units (GPGPUs), and presents results from initial field tests which confirm the anticipated capability of the system.     

GRAS radio occultation on-board of Metop

The GRAS radio occultation instrument is flying on Metop-A and belongs to the EPS (EUMETSAT Polar System). GRAS observes GPS satellites in occultation. Within this work, validation of GRAS closed-loop bending angle data against co-located ECMWF profiles extracted from model fields and occultations from the COSMIC constellation of radio occultation instruments is shown. Results confirm the high data quality and robustness, where GRAS shows lower bending angle noise against ECMWF than COSMIC and in terms of occultations per day, one GRAS ≈ two COSMIC satellites. This is partly due to the operational setup of EPS. For the investigation we focus on two observation periods where updates in the ECMWF (March 2009) and COSMIC processing (October 2009) have improved the statistics further. Bending angles biases agree to within 0.5% against ECMWF and to within 0.1% against COSMIC after the updates for altitudes between 8 and 40 km. In addition, we also analyze the impact of the Metop orbit processing on the derived GRAS bending angle data, where different GPS and Metop orbit solutions are analyzed. Results show that a batch based orbit processing would improve in particular the bending angle bias behavior at higher altitudes. Requirements for the operational processing of GRAS data are briefly outlined, options to ease the use of other positioning system satellites in the near future are discussed. A simplified analysis on the observation of several of these systems, e.g. GPS and Galileo, from one platform shows that about 16% of occultations are found within 300 km, ±3 h, thus providing similar information. A constellation of 2 GRAS like instruments would have only about 10% close-by.     

Assessing the performance of BDS-3 for multi-GNSS static and kinematic PPP-AR

As of 2021, a total of four different GNSS constellations – namely, GPS, GLONASS, Galileo, and BDS-3 – can be used with Full Operational Capability (FOC). In this work, the contribution of BDS-3 FOC to GPS + GLONASS + Galileo (GRE) PPP-AR is investigated, considering the three different cut-off angles (7°, 30°, and 45°) and different lengths of static observation sessions (24-, 12-, 6-, 3-, 1-, 0.5-, 0.25-hour). The data of 31 IGS-MGEX stations is processed with GRE PPP-AR and GREC3 (GPS + GLONASS (using float mode) + Galileo + BDS-3) PPP-AR modes. The results showed that BDS-3 degraded the horizontal (except for 24-h sessions) and vertical accuracy of static GRE PPP-AR solutions regardless of the elevation cutoff angle and observation time. The kinematic results showed that BDS-3 significantly contributed to the accuracy of GRE kinematic PPP-AR for 30° and 45° cutoff angles. The convergence time analysis showed that BDS-3 only contributes to GRE kinematic PPP-AR for the vertical component.     

基于TLE的Starlink星座第一阶段部署情况分析

近年来,随着卫星技术的快速发展和低轨（low earth orbit,LEO）卫星宽带互联网建设需求的不断增加,低轨大规模星座发展日新月异。针对Starlink星座初始化部署问题,首先论述了“星链”（Starlink）星座现状,分析在轨卫星高度变化。然后利用公开的两行轨道根数（two-line element,TLE）,从卫星发射入轨、轨道面分布两个方面,简要分析了Starlink星座的部署情况,给出升交点的变化规律;同时仿真分析了Starlink星座对地面的覆盖性能。最后,给出星座轨道面和相位分布、故障卫星处置以及可见卫星数量。所分析的结果以期为中国未来部署大规模LEO星座的建设提供借鉴。     

Compatibility of Terrestrial Reference Frames used in GNSS broadcast messages during an 8 week period of 2019

The operational Terrestrial Reference Frames (TRFs) realized through the evaluation of broadcast ephemerides for GPS, GLONASS, Galileo, BeiDou-2 and BeiDou-3 have been compared to IGS14, the TRF realized by the International GNSS Service (IGS). The TRFs realized by the GPS, GLONASS, Galileo, and BeiDou-2 and BeiDou-3 broadcast ephemerides are the orbital realizations of WGS 84 (G1762′), PZ90.11, GTRF19v01, and BDCS respectively. These TRFs are compared using up to 56 days of data (21 July-14 Sept 2019) at a 5 or 15-min rate. The operational TRFs are compared to IGS14 in a 7-parameter similarity (Helmert) transformation. Numerical results show that the operational GNSS TRFs differ from IGS14 at a level no greater than 4 cm for Galileo, 6 cm for GPS and BeiDou-3, 13 cm for GLONASS, and 48 cm for a limited set of BeiDou-2 Medium Earth Orbit (MEO) vehicles.     

Coastal sea level measurements using a single geodetic GPS receiver

This paper presents a method to derive local sea level variations using data from a single geodetic-quality Global Navigation Satellite System (GNSS) receiver using GPS (Global Positioning System) signals. This method is based on multipath theory for specular reflections and the use of Signal-to-Noise Ratio (SNR) data. The technique could be valuable for altimeter calibration and validation. Data from two test sites, a dedicated GPS tide gauge at the Onsala Space Observatory (OSO) in Sweden and the Friday Harbor GPS site of the EarthScope Plate Boundary Observatory (PBO) in USA, are analyzed. The sea level results are compared to independently observed sea level data from nearby and in situ tide gauges. For OSO, the Root-Mean-Square (RMS) agreement is better than 5 cm, while it is in the order of 10 cm for Friday Harbor. The correlation coefficients are better than 0.97 for both sites. For OSO, the SNR-based results are also compared with results from a geodetic analysis of GPS data of a two receivers/antennae tide gauge installation. The SNR-based analysis results in a slightly worse RMS agreement with respect to the independent tide gauge data than the geodetic analysis (4.8 cm and 4.0 cm, respectively). However, it provides results even for rough sea surface conditions when the two receivers/antennae installation no longer records the necessary data for a geodetic analysis.     

GNSS water vapour tomography – Expected improvements by combining GPS,GLONASS and Galileo observations

The German Research Centre for Geosciences (GFZ) operates a GNSS water vapour tomography system using about 350 German GNSS stations. The GNSS data processing at the GFZ works in near real-time and provides zenith total delays, integrated water vapour and slant delay data operationally. This large data set of more than 50,000 slant delays per hour is used to reconstruct spatially resolved humidity fields by means of tomographic techniques. It can be expected that additional observations from the future Galileo system provide more information with improved quality. A simulation study covering 12 h at 14 July 2009 was therefore started to estimate the impact of GPS, Galileo and GLONASS data on the GNSS tomography. It is shown that the spatial coverage of the atmosphere with slant paths is highly improved by combining observations from two or three satellite systems. Equally important for a reliable tomographic reconstruction is the distribution of slant path intersections as they are required to locate the integrated delay information. The number of intersection points can be increased by a factor of 4 or 8 if two or three systems are combined and their distribution will cover larger regions of the atmosphere. The combined data sets can be used to increase the spatiotemporal resolution of the reconstructed humidity fields up to 30 km horizontally, 300 m vertically and 15 min. The reconstruction quality could not be improved considerably using the currently available techniques.