The program system “Potsdam-4” for differential improvement of orbital elements and other parameters

The “Potsdam-4” program system can model satellite orbits with a high accuracy by means of numerical integration. It takes into account gravitational forces of the earth, the moon, the sun, and non-gravitational forces, and realizes an inertial system with a high accuracy. On the basis of this orbit calculation, geodetic and geodynamic parameters can be determined by means of photographic, laser and Doppler observations.     

Limits of photosynthesis in extrasolar planetary systems for earth-like planets.

We present a general modeling scheme for investigating the possibility of photosynthesis-based life on extrasolar planets. The scheme focuses on the identification of the habitable zone in main-sequence-star planetary systems with planets of Earth mass and size. Our definition of habitability is based on the long-term possibility of photosynthetic biomass production as a function of mean planetary surface temperature and atmospheric CO2-content. All the astrophysical, climatological, biogeochemical, and geodynamic key processes involved in the generation of photosynthesis-driven life conditions are taken into account. Implicitly, a co-genetic origin of the central star and the orbiting planet is assumed. The numerical solution of an advanced geodynamic model yields realistic look-up diagrams for determining the limits of photosynthesis in extrasolar planetary systems, assuming minimum CO2 levels set by the demand of C4 photosynthesis.     

Calibration of Envisat radar altimeter over Lake Issykkul

This study presents the results of calibration/validation (C/V) of Envisat satellite radar altimeter over Lake Issykkul located in Kyrgyzstan, which was chosen as a dedicated radar altimetry C/V site in 2004. The objectives are to estimate the absolute altimeter bias of Envisat and its orbit based on cross-over analysis with TOPEX/Poseidon (T/P), Jason-1 and Jason-2 over the ocean. We have used a new method of GPS data processing in a kinematic mode, developed at the Groupe de Recherche de Geodesie Spatiale (GRGS), which allows us to calculate the position of the GPS antenna without needing a GPS reference station. The C/V is conducted using various equipments: a local GPS network, a moving GPS antenna along the satellites tracks over Lake Issykkul, In Situ level gauges and weather stations. The absolute bias obtained for Envisat from field campaigns conducted in 2009 and 2010 is between 62.1 and 63.4 ± 3.7 cm, using the Ice-1 retracking algorithm, and between 46.9 and 51.2 cm with the ocean retracking algorithm. These results differ by about 10 cm from previous studies, principally due to improvement of the C/V procedure. Apart from the new algorithm for GPS data processing and the orbit error reduction, more attention has been paid to the GPS antenna height calculation, and we have reduced the errors induced by seiche over Lake Issykkul. This has been assured using cruise data along the Envisat satellite track at the exact date of the pass of the satellite for the two campaigns. The calculation of the Envisat radar altimeter bias with respect to the GPS levelling is essential to allow the continuity of multi-mission data on the same orbit, with the expected launch of SARAL/Altika mission in 2012. Implications for hydrology in particular, will be to produce long term homogeneous and reliable time series of lake levels worldwide.     

Comparison between African equatorial station ground-based inferred vertical E × B drift,Jicamarca direct measured drift,and IRI model

The Incoherent Scatter Radar measurement over Jicamarca, together with the IRI model-2007 measurements were compared with ground-based digisonde inferred E × B drift over Ilorin in the African region during year of solar minima (F10.7 = 81). Seasonally, Ilorin pre-reversal enhancement (PRE) had peak drift velocities of 7.2, 3.7 and 7.9 m/s for March equinox, September equinox and December solstice respectively, while Jicamarca drifts indicated 13.0, 10.5 and 5.2 m/s; as well as the IRI model with 14.3, 8.4 and 0.7 m/s in similar order. PRE value was insignificant during June solstice. The PRE magnitude of the IRI-model during the equinoxes is twice the value obtained at Ilorin. The daytime E × B drift peaked over Ilorin 1–2 h earlier than both the modeled and Jicamarca observations. This could be due to the difference in sunset time at the conjugate points corresponding to the altitude of the observation. During the evening time PRE, the respective correlation coefficients (R) for Vz–F10.7 relation over Jicamarca, Ilorin and the modeled observations are −0.5559, 0.4796 and −0.4979. Similarly, the Vz–Ap relation exhibit excellent anti-correlation coefficient (R = −0.8637) for the IRI-model, −0.4827 over Jicamarca and 0.3479 for Ilorin. Annual mean drift velocities over Jicamarca, Ilorin and IRI model measurements respectively are 10, 5.6 and 10 m/s for the peak PRE observation; 15, 16 and 21 m/s for the daytime pre-sunrise peak values; and −21, −9 and −16 m/s for the nighttime downward reversals. The root-mean square (RMS) deviation between IRI-model and the Ilorin drift between 2000 and 0500 h is 4.37, 2.03, 3.71 and 2.42 m/s for March equinox, June solstice, September equinox and December solstice respectively. For Jicamarca–Ilorin drift relation, RMS deviation is 5.48, 2.30, 3.47 and 1.27 m/s in the same order respectively. Annual hmF2 inferred drift over Ilorin during daytime is higher by a factor of ≈2 and 3 at Jicamarca and IRI model measurements respectively; and by a factor of ≈5 for both during the night-time period. The limitations in using hmF2 to infer drifts are discussed.     

Mapping the earth's magnetic and gravity fields from space: Current status and future prospects

Orbital potential field measurements are sensitive to regional variations in earth density and magnetization that occur over scales of a few hundred kilometers or greater. Global field models currently available are able to distinguish gravity variations of ±5 milligal over distances of ～1,000 km and magnetic variations of ±6 gamma over distances of ～300 km at the earth's surface. Regional variations in field strength have been detected in orbital measurements that are not apparent in higher resolution, low altitude surveys. NASA is presently studying a spacecraft mission known as GRAVSAT/MAGSAT, which would be the first satellite mission to perform a simultaneous survey of the earth's gravity and magnetic fields at low orbital altitudes. GRAVSAT/MAGSAT has been proposed for launch during the latter nineteen-eighties, and it would measure gravity field strength to an accuracy of 1 milligal and magnetic field strength to an accuracy of 2 gamma (scalar)/5 gamma (vector components) over a distance of roughly 100 km. Even greater improvements in the accuracy and spatial resolution of orbital surveys are anticipated during the nineteen-nineties with the development of potential field gradiometers and a tethered satellite system that can be deployed from the Space Shuttle to altitudes of 120 km above the earth's surface.     

Observations of ionospheric irregularities and its correspondence with sporadic e occurrence over South Korea and Japan

This study analyzed the occurrence of ionospheric irregularities over South Korea and Japan (mid-latitudes) during the years 2010–2015. The irregularities were quantified using the rate of change of total electron content (TEC) index (ROTI), which detects irregularities with scale sizes in the range of 400 m–2.5 km. The ROTI threshold for an ionospheric irregularity to have occurred was set as 0.1 TECU/min. Results showed that ionospheric irregularities mostly occur during night-time and around local noon-time in the seasons of spring and summer. In addition, the percentage of ionospheric irregularities had a high positive correlation with solar flux (F10.7) (r > 0.72). For the first time, we showed good correspondence between ionospheric irregularities measured by the ROTI index and sporadic E (Es). The median ROTI associated with ionospheric irregularities over a South Korea station (DAEJ) and a Japan station (KGNI) were 0.131 and 0.125 TECU/min, respectively. However, in severe cases of ionospheric irregularities, the ROTI values over DAEJ (KGNI) can reach 0.246 (0.217) and 0.314 (0.339) TECU/min during day and night, respectively. These critical ROTI values can be important in interpreting and monitoring ionospheric irregularity occurrence over South Korea and Japan.     

The energy spectra of protons and helium measured with the ATIC experiment

The Advanced Thin Ionization Calorimeter (ATIC) balloon experiment is designed to investigate the composition and energy spectra of cosmic rays at the highest energies currently accessible by direct measurements, i.e., the region up to 100 TeV. The instrument consists of a silicon matrix for charge measurement, a graphite target (0.75 nuclear interaction length) to induce hadronic interactions, three layers of scintillator strip hodoscopes for triggering and trajectory reconstruction, and a Bismuth Germanate (BGO) crystal calorimeter (18 radiation lengths) to measure particle energies. ATIC has had two successful Long Duration Balloon (LDB) flights from McMurdo, Antarctica: one from 12/28/00 to 01/13/01 and the other from 12/29/02 to 01/18/03. We present the energy spectra of protons and helium extracted from the first flight, over the energy range from 100 GeV to 100 TeV, and compare them with the results from other experiments at both the lower and higher energies. ATIC-1 results do not indicate significant differences in spectral shape between protons and helium over the investigated energy range.     

Ozone mini-hole observation over the Balkan Peninsula in March 2005

Areas with dimensions of 1000–3000 km in which the total ozone content (TOC) decreases fast are called ozone mini-holes. They are generated mainly dynamically in two ways, either by poor-ozone air mass transport from the tropics to higher latitudes by planetary wave activity or, they are connected with strong adiabatic uplifting of the tropopause height. An ozone mini-hole, generated by the second mechanism, was observed over the Balkan Peninsula on 19/21 March 2005. In the middle of March, the polar vortex was strongly disturbed by Rossby waves, reaching up to the lower stratosphere. Warming episodes over a geographical area, covering the Barents Sea and the Polar Sea north from Central Siberia, displaced a polar vortex fragment extremely southwards. However, the vorticity was weak and the stratospheric temperatures did not reach low values, providing conditions for ozone chemical destruction via heterogenic reactions. At the same time, a Rossby wave ridge was located below the European polar fragment. In the period from 13 to 19 March, the thermal tropopause over Sofia was uplifted almost by 3 km. Ozone distributions observed by the SCIAMACHY instrument on 18–21 March show a fast TOC decrease westwards from Ireland, which was moving eastwards during the next days, increasing the area in which the ozone content decreased. On 20/21 March low ozone content was observed above the Stara Zagora (42°N, 25°E) ground-based station by means of the GASCOD instrument, using DOAS technique. The TOMS Earth probe instrument detected 237 DU over Sofia. This is a record low March value from the beginning of the TOMS instrument measurements in 1978. In March/April the ozone distribution was characterized by its mean annual maximum of 360 DU at 42°N.     

An annual variation analysis of the ionospheric spatial gradient over a regional area for GNSS applications

An ionospheric spatial gradient represents the ionosphere delay difference between different locations, and its variation over a specific area is important for implementing differential GNSS systems. An estimation method for the ionospheric spatial gradient over a small regional area is proposed. A plate map model is implemented for the direct estimation of the gradients. Nine years of GPS data were processed to figure out the annual variation of the mean gradient at the mid-geomagnetic latitude of 30° N. Gradients along the north–south direction have a mean of 0.65 mm/km and follow solar-cycle variations.     

Study of satellite retrieved CO2 and CH4 concentration over India

This paper reports a study of spatial and temporal variations of columnar averaged concentration of CO₂ and CH₄ over India using SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) and Greenhouse gas Observing SATellite (GOSAT) data. Comparison of these data with the global view National Oceanic and Atmospheric Administration (NOAA) land data and also location specific flask data is made. The temporal variation in column averaged global CO₂ is similar to that over India and it is also similar to the NOAA surface flask data and global view. The variation in NOAA surface CH₄ is location dependent and its global view appears to vary seasonally in opposite phase with the column averaged CH₄ values from satellites, reflecting the limited comparability of surface and column averaged data. Over India the CO₂ maximum is in May and minimum in August/September while for CH₄ the maximum is in September and minimum in February/March. The seasonal variation of CH₄ over India is correlated with the eastern coastal rice cultivation.     

Rio de la Plata estuarine system: Relationship between river flow and frontal variability

This paper deals with the application of SeaWIFS images to characterize spatial and temporal variability of fronts in the Rio de la Plata estuarine system over the period 2000–2003. We aim to depict the relationship between river outflow and variability of fronts’ loci on monthly to ENSO-related timescales and the influence of the winds along Rio de la Plata (axial winds) on the abrupt changes in frontal dynamics over synoptic timescales. During the studied period both La Niña (July 1999–June 2000) and El Niño (April 2002–May 2003) events induced significant displacements of fronts. Three distinct fronts were analyzed between river, estuarine, coastal and marine waters of the Rio de la Plata: Main Turbidity Front, Main Marine Front, and Secondary Marine Front. Their number, location and separation seem to be mainly related to river outflow and second, to fresh (>8 m/s) axial winds. During low discharge periods (i.e. summer time and/or La Niña events) these winds induce abrupt changes in the location of fronts (100–200 km) and greater separation between them over synoptic timescale, whereas during high river discharge or ENSO years some of the variability of fronts location is explained by the river’s outflow fluctuations, especially by the high variability of the River Uruguay discharge.     

Comparison of Landsat OLI,ASTER, and Sentinel 2A data in lithological mapping : A Case study of Rich area (Central High Atlas,Morocco)

The eastern part of the Rich area consists of the massive Paleozoic and Meso-Cenozoic cover formations that present the geodynamic development of the study area, where is characterized by various carbonate facies of Jurassic age. The geographical characteristic of the study area leaves the zone difficult to map by conventional methods. The objective of this work focuses on the mapping of the constituent lithological units of the study area using multispectral data of Landsat OLI, ASTER, and Sentinel 2A MSI. The processing of these data is based on a precise methodology that distinguishs and highlights the limits of the different lithological units that have an approximate similarity of spectral signature. Three techniques were used to enhance the image including Principal Component Analysis (PCA), Minimum Noise Fraction (MNF), and Independent Component Analysis (ICA). Lithological mapping was performed using two types of supervised classification : Maximum likelihood classifier (MLC) and Support Vector Machine (SVM).The results of processing data show the effectiveness of Sentinel 2A data in mapping of lithological units than the ASTER and Landsat OLI data. The classification evaluation of two methods of the Sentinel 2A MSI image showed that the SVM method give a better classification with an overall accuracy of 93,93% and a Kappa coefficient of 0.93, while the MLC method present an overall accuracy of 82,86% and a Kappa coefficient of 0.80. The results of mapping obtained show a good correlation with the geological map of the study area as well as the efficiency of remote sensing in identification of different lithological units in the Central High Atlas.     

First quality assessment of the Cryosat-2 altimetric system over ocean

This paper presents the results of a Calval analysis performed for the Cryosat-2 mission over ocean. The data set used in this analysis consists of products generated by the Cryosat-2 Processor Prototype developed by CNES (Center National d’Etudes Spatiales). This data set has been analysed focusing on LRM (Low Resolution Mode) mode only. One major objective of this paper is to illustrate the potential of Cryosat-2 data over ocean, mainly for waves and Sea Level Anomaly applications. All the results indicate very good performances of the SIRAL (SAR/Interferometric Altimeter) altimeter over ocean. Crossover standard deviation is close to 6.5 cm over the analysed period (3 months) which is close to the Jason-2 and ENVISAT performance. All these results confirm that Cryosat’s altimeter can provide data almost as valuable as other flying altimetric missions, and that it has the potential to contribute to oceanography (e.g. multi-mission climate record, mesoscale monitoring in near real time) and to geodesy (e.g. mean sea surface, bathymetry).     

Correlations between neutral and ionized solar wind

We report results of a statistical study correlating ionized solar wind (ISW) fluxes observed by ACE during late 2000 and throughout 2001 with neutral solar wind (NSW) fluxes observed by IMAGE/LENA over the same period. The average correlation coefficient between the neutral and ionized solar wind is 0.66 with correlations greater than 0.80 occurring about 29% of the time. Correlations appear to be driven by high solar wind flux variability, similar to results obtained by in situ multi-spacecraft correlation studies. In this study, however, IMAGE remains inside the magnetosphere on over 95% of its orbits. As a function of day of year, or equivalently ecliptic longitude, the slope of the relationship between the neutral solar wind flux and the ionized solar wind flux shows an enhancement near the upstream direction, but the symmetry point appears shifted toward higher ecliptic longitudes than the interstellar neutral (ISN) flow direction by about 20°. The estimated peak interstellar neutral upstream density inside of 1 AU is about 7 × 10⁻³ cm⁻³.     

Characterization of the low latitude plasma density irregularities observed using C/NOFS and SCINDA data

Complex electrodynamic processes over the low latitude region often result in post sunset plasma density irregularities which degrade satellite communication and navigation. In order to forecast the density irregularities, their occurrence time, duration and location need to be quantified. Data from the Communication/Navigation Outage Forecasting System (C/NOFS) satellite was used to characterize the low latitude ion density irregularities from 2011 to 2013. This was supported by ground based data from the SCIntillation Network Decision Aid (SCINDA) receivers at Makerere (Geographic coordinate 32.6°E, 0.3°N, and dip latitude ?9.3°N) and Nairobi (Geographic coordinate 36.8°E, ?1.3°N, and dip latitude ?10.8°N). The results show that irregularities in ion density have a daily pattern with peaks from 20:00 to 24:00 Local Time (LT). Scintillation activity at L band and VHF over East Africa peaked in 2011 and 2012 from 20:00 to 24:00 LT, though in many cases scintillation at VHF persisted longer than that at L band. A longitudinal pattern in ion density irregularity occurrence was observed with peaks over 135–180°E and 270–300°E. The likelihood of ion density irregularity occurrence decreased with increasing altitude. Analysis of C/NOFS zonal ion drift velocities showed that the largest nighttime and daytime drifts were in 270–300°E and 300–330°E longitude regions respectively. Zonal irregularity drift velocities over East Africa were for the first time estimated from L-band scintillation indices. The results show that the velocity of plasma density irregularities in 2011 and 2012 varied daily, and hourly in the range of 50–150 m s^?1. The zonal drift velocity estimates from the L-band scintillation indices had good positive correlation with the zonal drift velocities derived from VHF receivers by the spaced receiver technique.     

Calibration of IRI–ITU-R peak density and height over the oceans with topside sounding data

Accuracy of IRI electron density profile depends on the F2 layer peak density and height converted by empirical formulae from the critical frequency and M3000F2 factor provided by the ITU-R (former CCIR). The CCIR/ITU-R maps generated from ground-based ionosonde measurements suffer from model assumptions, in particular, over the oceans where relatively few measurements are available due to a scarcity of ground-based ionosondes. In the present study a grid-point calibration of IRI/ITU-R maps for the foF2 and hmF2 over the oceans is proposed using modeling results based on the topside true-height profiles provided by ISIS1, ISIS2, IK-19 and Cosmos-1809 satellites for the period of 1969–1987. Topside soundings results are compared with IRI and the Russian standard model of ionosphere, SMI, and grouped to provide an empirical calibration coefficient to the peak density and height generated from ITU-R maps. The grid-point calibration coefficients maps are produced in terms of the solar activity, geodetic latitude and longitude, universal time and season allowing update of IRI–ITU-R predictions of the F2 layer peak parameters.     

Geomagnetic pulsation over conjugate locations during geomagnetic storms and substorm

Geomagnetic data collected during magnetic storm over magnetically conjugate pair (according to IGRF 2000) of high latitude stations viz., Maitri (70° 45′ S, 11° 42′ E) and Tromso (69° 40′ N, 18° 56′ E) reveal that amplitudes of Pc6 pulsation characteristically differ. The amplitude obtained from horizontal magnetic field for the Pc6 pulsation frequency range between 0.6 and 1.6 mHz significantly differs in time corresponding to peak amplitude. The relative differences in its time of occurrence found to gradually increase around initial phase of storm and remain exactly out of phase at peak amplitude of storm. Thence, it is found to be in agreement in phase gradually until storm unwinds. This indicates that simultaneous amplitude of Pc6 pulsation at conjugate pair of stations and its time of occurrence could be a key factor to infer storm arrivals somewhat prior to its peak effects. The emphasize remains on prediction of storm arrivals only by utilizing ground based magnetometer observations. However, it is necessary to understand differences on the basis of weak, moderate, strong, and super strong cases and more exactly how they behave along the line of magnetic Meridian. Nevertheless, the analysis implies that geo-effective magnetic ejecta/clouds/CIRs/sheaths/CMEs/ICMEs giving rise to geomagnetic storm can be predicted ahead of its peak effects by having magnetometer data over conjugate locations.     

Ionospheric mapping by regional spherical harmonic analysis: New developments

The method of Spherical Cap Harmonic Analysis (SCHA) /1/ has been applied /2/ to the critical frequency of the F₂ layer (f₀F₂) for mapping and modelling it over Europe. The model was based on longitudinal expansion in Fourier series, and fractional Legendre colatitudinal functions over a spherical cap including Europe. Here a new and simpler technique, previously developed for the regional modelling of the geomagnetic field /3/, is introduced and described. The basic improvement of the new method, called Adjusted Spherical Harmonic Analysis (ASHA), implies the use of conventional Spherical Harmonic (SH) functions after the colatitude interval is adjusted to that of a hemisphere. Examples are shown dealing with the application of ASHA to retrospective mapping and modelling of the monthly medians of f₀F₂ over Europe.     

Revised upper limit on the internal magnetic moment of Venus

Over four Venus years of low altitude nightside PVO magnetometer observations are used to establish a new upper limit for the magnetic moment of Venus. Improvements over previous studies include data coverage and new instrument calibration information. The upper limit on an internal dipole moment is determined to be 8.4 × 10¹⁰ T m³.