Global morphology of ionospheric sporadic E layer from the FormoSat-3/COSMIC GPS radio occultation experiment

Ionospheric sporadic-E (Es) activity and global morphology were studied using the 50 Hz signal-to-noise ratio amplitude and excess phase measurements from the FormoSat-3/Constellation Observing System for Meteorology, Ionosphere and Climate (FS3/COSMIC) GPS radio occultation (RO) observations. The results are presented for data collected during the last sunspot cycle from mid-2006 to the end of 2017. The FS3/COSMIC generally performed more than 1000 complete E-region GPS RO observations per day, which were used to retrieve normalized L1-band amplitude standard deviation (SDL1) and relative electron density (N_e) profiles successfully. More or less 31% of those observations were identified as Es events based on SDL1 and peak SDL1 altitude criteria. We found that the peak Es-event i values are approximately proportional to the logarithms of the corresponding peak N_e differences. Five major geographical zones were identified, in which the seasonal and diurnal Es occurrence patterns are markedly different. These five zones include the geomagnetic equatorial zone (??5°?<?magnetic latitude (ML)?<?5°), two extended geomagnetic mid-latitude zones (15°?<?ML?<?55°, and ??55°?<?ML < ??15°), and two auroral zones (70°?<?ML, and ML < ??70°). The Es climatology, namely its variations with each identified zone, altitude, season, and local time has been documented.     

Characterization of forests and deforestation in Cambodia using ALOS/PALSAR observation

In this study, we have demonstrated the capability of full polarimetric ALOS/Phased Array L-band Synthetic Aperture Radar data for the characterization of the forests and deforestation in Cambodia, to support climate change mitigation policies of Reducing Emission from Deforestation and Forest Degradation (REDD). We have observed mean backscattering coefficient (σ°), entropy (H), alpha angle (α), anisotropy (A), pedestal height (PH), Radar Vegetation Index (RVI) and Freeman–Durden three-component decomposition parameters. The observations show that the forest types and deforested area are showing variable polarimetric and backscattering properties because of the structural difference. Evergreen forest is characterized by a high value of σ° HV (?12.96 dB) as compared with the deforested area (σ° HV=?22.2 dB). The value of polarimetric parameters such as entropy (0.93), RVI (0.91), PH (0.41) and Freeman–Durden volume scattering (0.43) is high for evergreen forest, whereas deforested area is characterized by the low values of entropy (0.36) and RVI (0.17). Based on these parameters, it is found that σ° HV, entropy, RVI and PH provide best results among other parameters.     

FORMOSAT-3/COSMIC observations of the ionospheric auroral oval development

The ionospheric radiance and electron density observed by the tiny ionospheric photometer (TIP) and GPS occultation experiment (GOX) payloads on FORMOSAT-3/COSMIC satellites are applied to determine the boundaries of the auroral oval and its width in the winter nighttime ionosphere for both hemispheres. The TIP collects ionospheric emission at 135.6 nm due to electron impact excitation, while the GOX offers ionospheric electron density profiles with radio occultation (RO) technique. Comparison between them shows similar patterns of the plasma structure in the polar caps. The mean width of the auroral bands ranges between about 2 and 11° latitude in the winter nighttime and it varies with longitudes. The comparison by month suggests that the mean radius of the auroral ovals varies with the intensity of the auroral radiance.     

Artificial ionospheric wave number 4 structure below the F2 region due to the Abel retrieval of radio occultation measurements

We analyzed the effect of the Abel inversion on the wave number 4 (WN4) structure from the GPS radio occultation (RO)–measured electron densities by using the FORMOSAT-3/COSMIC (F-3/C) observations under the equinox condition. The Abel-retrieved electron density from both the F-3/C observations and the simulated results by an empirical model with an imposed WN4 structure in the F layer are investigated. It is found that the Abel inversion can reproduce the real WN4 structure well in the F2 layer. However, it will result in pseudo and reversed-phase WN4 structure in the lower altitude (F1 and E layers). Quantitatively, relative ±15% WN4 signature in the F2 layer can produce ±40% artificial WN4 in the E and F1 layers. Analysis on the F-3/C data shows about ±15% WN4 signature in the F2 layer and ±50% WN4 with reversed-phase in the E and F1 layers. The F-3/C-observed WN4 structure in the E and F1 layers might be the combinations of the real WN4 signature and the artificial effects of Abel retrieval.     

Hyperspectral radiometric observation of the northeast Arabians Sea during April 2006

Parameters were retrieved from the hyperspectral radiometer like upwelling and downwelling radiance (Lu and Ed) upwelling and downwelling attenuation (K-Lu and K-Ed) for 9 stations in the northeast Arabian Sea between 16–26 April 2006. Data was analyzed for 5 offshore and 4 coastal stations of the cruise SS-244, on board FORV “Sagar Sampada” between latitude 9-22oN and longitude 68–74°E. The peak for all parameters was observed to be different respectively for depths 1, 5, 10, 20, 30, 40, 50 meters in coastal and offshore stations. Each peak in the respective wavelength is due to a particular composition; phytoplankton pigments have spectral peaks at 443, 490, 555, 670 nm, suspended matter, sediments have peaks at 630 and 670 nm. Detailed analysis with High Performance Liquid Chromatography (HPLC) data and comparison with the water composition of our hyperspectral radiometer results show that the marine cyanophyte, Trichodesmium bloom produces high pigment concentrations of chlorophyll-a, zeaxanthin, β-carotene and pheophytin and their absorptions are interpreted at wavelengths 443, 490, 515 and 536 nanometers, respectively. A dip around 515 nm was seen in the Ed and Lu profiles in our study.     

Global ionosphere map constructed by using total electron content from ground-based GNSS receiver and FORMOSAT-3/COSMIC GPS occultation experiment

Effects of rapidly changing ionospheric weather are critical in high accuracy positioning, navigation, and communication applications. A system used to construct the global total electron content (TEC) distribution for monitoring the ionospheric weather in near-real time is needed in the modern society. Here we build the TEC map named Taiwan Ionosphere Group for Education and Research (TIGER) Global Ionospheric Map (GIM) from observations of ground-based GNSS receivers and space-based FORMOSAT-3/COSMIC (F3/C) GPS radio occultation observations using the spherical harmonic expansion and Kalman filter update formula. The TIGER GIM (TGIM) will be published in near-real time of 4-h delay with a spatial resolution of 2.5° in latitude and 5° in longitude and a high temporal resolution of every 5 min. The F3/C TEC results in an improvement on the GIM of about 15.5%, especially over the ocean areas. The TGIM highly correlates with the GIMs published by other international organizations. Therefore, the routinely published TGIM in near-real time is not only for communication, positioning, and navigation applications but also for monitoring and scientific study of ionospheric weathers, such as magnetic storms and seismo-ionospheric anomalies.     

TRIGONOMETRICAL LEVELLING ACROSS THE INLAND ICE IN NORTH GREENLAND

Abstract

As part of the scientific work of the British North Greenland Expedition (1952–1954), a programme of trigonometrical levelling was carried out from the east to the west coast of Greenland, along a line across the inland ice between latitudes 76° 40′ N., and 78° 10′ N. The primary purpose of the work was to determine accurately the heights above sea level of a series of gravity stations, the gravity measurements being made in connection with determinations of ice thickness. For meteorological purposes it was necessary to know also the altitude of the Expedition's central station, situated in latitude 78° 04′ N., longitude 38° 29′ W. The accuracy necessary for the purpose of the gravity survey was a few metres for the altitudes, while the latitude of each gravity station had to be determined with an accuracy of ± 0.1 minute.     

GPS radio occultation measurements on ionospheric electron density from low Earth orbit

Since the proof-of-concept GPS/Meteorology (GPS/MET) experiment successfully demonstrated active limb sounding of the Earth’s neutral atmosphere and ionosphere via GPS radio occultation (RO) from low Earth orbit, the developments of electron density (n _e) retrieval techniques and powerful processing systems have made a significant progress in recent years. In this study, the researches of n _e profiling from space-based GPS RO observations are briefly reviewed. Applying to the Formosat-3/Constellation Observing System for Meteorology, Ionosphere and Climate (FS3/COSMIC) data, we also present a compensatory Abel inversion technique including the effects of large-scale horizontal gradients and/or inhomogeneous ionospheric n _e obtained from an improved near real-time phenomenological model of the TaiWan Ionospheric Model. The results were evaluated by the ionosonde foF2 and foE data and showed improvements of rms foF2 difference from 29.2 to 16.5% in relative percentage and rms foE difference from 54.2 to 32.7% over the standard Abel inversion.     

Intra-Seasonal Oscillations Associated with Indian Ocean Warm Pool and Summer Monsoon Rainfall and Their Inter-Annual Variability

Sea surface temperature (SST) and wind speed (WS) derived from Multichannel Scanning Microwave Radiometer (MSMR), onboard IRS-P4 (Oceansat-1) satellite were used to generate spatially averaged (80° to 100° E & 0° to 10° S) daily data during June to September from 1999 to 2001 along with collocated outgoing long wave radiation (OLR) data from NCEP/NCAR reanalysis product. Daily rainfall data over the peninsular India during the above period was taken from the weekly weather report published by the Indian Meteorological Department. The data were then subjected to power spectral analysis. Dominant 30 to 60 days oscillations were observed over both ocean and land during all the 3 years investigated. It is quite interesting to note that the intra-seasonal oscillation with 60 days periodicity dominates in both ocean and land during 1999 and then transforms to 30 and 40 days periodicities in 2000 and 2001 respectively with a phase lag of about 25 to 30 days. It was also observed that the time-latitude section of zonally averaged OLR and TMI derived cloud liquid water data clearly depict the propagation of convection and cloud from the equator to the north at the rate of 0.75° to 1° latitude per day which corroborates well with the rate of propagation derived from the phase lag obtained between the processes occurring at the equatorial Indian Ocean and rainfall recorded at the designated land segments. The results are in good agreement with the observations made by earlier investigators.     

Automatic detection of near surface geological and hydrological features and investigating their influence on groundwater accumulation and salinity in southwest Egypt using remote sensing and GIS

This study employs two automated algorithms, the topographic fabric and the deterministic eight-node (D8), to reveal near surface geological fractures and their associated paleodrainages network from a 90?m DEM of the Shuttle Topographic Radar Mission sensor. The topographic fabric algorithm, which calculates the slope and aspect which define a vector normal to the earth’s surface and then compute direction cosines of normal vector at each point, was used to auto-detect fault zones. The deterministic eight-node (D8) algorithm, which determines in which neighbouring pixel any water in a central pixel will flow naturally, was used to delineate paleodrainages concealed beneath sand sheets. Seven sets of geological structures were recognized. Their major trends were found to be in the S46°W, S58°W, S81°W, N107°W, N–S, E–W and N152°W and share similar trends of the revealed paleodrainages. The results suggest a strong spatial relationship between the features extracted from DEM and groundwater potential.     

Effect of seasonal spectral variations on land cover classification

The supervised classification (Maximum likelihood) on three dates of IRS (LISS III) satellite data was performed to study the effect of seasonal spectral variation on land cover classification for the study area falling in the Solan district of Himachal Pradesh at latitude 30° 50’ N to 31° ’N and longitude 77° 00’ E to 77° 15’ E. It was found that the summer dataset was better with overall classification accuracy of 76% as compared to winter and spring dataset with classification accuracy of 49 and 46%, respectively.     

利用FY3C与FY3D GPS掩星资料研究中低纬地区偶发E层形态的比较

偶发E层(sporadic E, Es)是主要发生在90~120 km高度的电子密度显著增强的电离层薄层,Es层的存在会导致掩星观测中全球导航卫星系统信号强度和相位的强烈波动。利用2019-01—2021-12风云三号C（Fengyun-3C,FY3C）和风云三号D(Fengyun-3D,FY3D)卫星GPS（global positioning system）掩星观测的50 Hz信噪比数据提取Es层信息,进而对两颗卫星数据分别反演得到的60°S~60°N中低纬地区Es层发生率的时空分布及季节变化进行比较。结果发现,虽然两颗卫星掩星资料得到的Es层发生率分布形态基本一致,均反映了Es层的发生率与地磁场和中性大气背景风场的相关性,但在大部分季节和地区,由FY3D得到的Es层发生率低于由FY3C得到的结果,北半球夏季中纬地区尤为明显,而FY3C反演结果与基于电离层与气候星座观测系统掩星数据的反演结果更为接近。导致差异的可能原因包括两颗卫星信噪比廓线的上边界高度分布和地方时覆盖上的差异、两颗卫星掩星接收机噪声水平的差异等。上述结果表明,后续融合两颗卫星的掩星数据进行Es层相关研究时,可能需要顾及两颗卫星接收机的不同噪声水平,在Es层发生的判定策略上进行针对性调整。     

GIS application to determine the effect of topography on landuse in Ashwani Khad watershed

The landuse status of Ashwani Khad watershed has been obtained using 1RS-ID satellite data for 1999 and further topographic analysis has been carried out using GIS software-ARC/ INFO and ARCVIEW. It has been found that of the total geographical area (85.30 sq. km) of the Ashwani Khad watershed which lies between 30°50′ to 31°N latitude and 77°05′ to 77° 15′E longitude in Himachal Pradesh, 54.53 % constituted wasteland, 33.55% agriculture and least 11.92 % forest. The altitude, aspect and slope have exhibited marked effect on land utilization. Agriculture and wasteland have been found maximum in mid altitude (1300-1500 m) and moderate slopes (13.2-26.4 degree), whereas, agriculture and forest have been maximum in flat and north aspect.     

Assessing antenna field of view and receiver clocks of COSMIC and GRACE satellites: lessons for COSMIC-2

We provide suggestions for the approved COSMIC-2 satellite mission regarding the field of view (FOV) and the clock stability of its future GNSS receiver based on numerical analyses using COSMIC GPS data. While the GRACE GPS receiver is mounted on the zenith direction, the precise orbit determination (POD) antennas of COSMIC are not. The COSMIC antenna design results in a narrow FOV and a reduction in the number of GPS observations. To strengthen the GPS geometry, GPS data from two POD antennas of COSMIC are used to estimate its orbits. The phase residuals of COSMIC are at the centimeter level, compared to the millimeter level of GRACE. The receiver clock corrections of COSMIC and GRACE are at the microsecond and nanosecond levels, respectively. The clock spectra of COSMIC at the frequencies of 0–0.005 Hz contain significant powers, indicating potential systematic errors in its clock corrections. The clock stability, expressed by the Allan deviation, of COSMIC ranges from 10^?9 to 10^?11 over 1 to 10⁴ s, compared to 10^?12 to 10^?14 for GRACE. Compared to USO-based clock of GRACE, the clock of COSMIC is degraded in its stability and is linked to the reduction of GPS data quality. Lessons for improvement of COSMIC-2 over COSMIC in FOV and receiver clock stability are given.     

Quality assessment of COSMIC/FORMOSAT-3 GPS radio occultation data derived from single- and double-difference atmospheric excess phase processing

This study evaluates the quality of GPS radio occultation (RO) atmospheric excess phase data derived with single- and double-difference processing algorithms. A spectral analysis of 1 s GPS clock estimates indicates that a sampling interval of 1 s is necessary to adequately remove the GPS clock error with single-difference processing. One week (May 2–8, 2009) of COSMIC/FORMOSAT-3 data are analyzed in a post-processed mode with four different processing strategies: (1) double-differencing with 1 s GPS ground data, (2) single-differencing with 30 s GPS clock estimates (standard COSMIC Data Analysis and Archival Center product), (3) single-differencing with 5 s GPS clocks, and (4) single-differencing with 1 s GPS clocks. Analyses of a common set of 5,596 RO profiles show that the neutral atmospheric bending angles and refractivities derived from single-difference processing with 1 s GPS clocks are the highest quality. The random noise of neutral atmospheric bending angles between 60 and 80 km heights is about 1.50e−6 rad for the single-difference cases and 1.74e−6 rad for double-differencing. An analysis of pairs of collocated soundings also shows that bending angles derived from single-differencing with 1 s GPS clocks are more consistent than with the other processing strategies. Additionally, the standard deviation of the differences between RO and high-resolution European Center for Medium range Weather Forecasting (ECMWF) refractivity profiles at 30 km height is 0.60% for single-differencing with 1 and 5 s GPS clocks, 0.68% for single-differencing with 30 s clocks, and 0.66% for double-differencing. A GPS clock-sampling interval of 1 s or less is required for single- and zero-difference processing to achieve the highest quality excess atmospheric phase data for RO applications.     

Electron density profiles in the equatorial ionosphere observed by the FORMOSAT-3/COSMIC and a digisonde at Jicamarca

We examine for the first time the ionospheric electron density profiles concurrently observed by the GPS occultation experiment (GOX) onboard the FORMOSAT-3/COSMIC (F3/C) and the ground-based digisonde portable sounder DPS-4 at Jicamarca (12°S, 283°W, 1°N geomagnetic) in 2007. Our results show that the F3/C generally underestimates the F2-peak electron density NmF2 and the F2-peak height hmF2. On the other hand, when the equatorial ionization anomaly (EIA) pronouncedly appears during daytime, the total electron content (TEC) derived from the radio occultation of the GPS signal recorded by the F3/C GOX is significantly enhanced. This results in the NmF2 at Jicamarca being overestimated by the Abel inversion on the enhanced TEC during the afternoon period.     

NOTES ON THE BASE LINES OF THE CEYLON TRIANGULATION

Abstract

The triangulation of Ceylon depends for its scale upon two bases, each about 5½ miles long, situated at Negombo on the West Coast (latitude 7° 10′) and at Batticaloa on the East Coast (latitude 7° 40′). Both bases are in low, flat country; brick towers up to 70 feet high had to be built over the terminals to enable observations to be taken to surrounding points. These lines have recently been re-measured.     

Algorithm for Harmonic Tidal Analysis along T/P Tracks： Taking Differences of Observed Sea Surface Heights at Adjacent Points as Observations

An algorithm （differential mode） is presented for the improvement of harmonic tidal analysis along T/P tracks, in which the differences between the observed sea surface heights at adjacent points are taken as observations. Also, the observation equations are constrained with the results of the crossover analysis; the parameter estimations are performed at 0.1° latitude intervals by the least squares. Cycle 10 to 330 T/P altimeter data covering the China Sea and the Northwest Pacific Ocean （2°-50° N,105°-150° E） are adopted for a refined along-track harmonic tidal analysis, and harmonic constants of 12 constituents in 8 474 points are obtained, which indicates that the algorithm can efficiently remove non-tidal effects in the altimeter observations, and improve the precision of tide parameters. Moreover, parameters along altimetry tracks represent a smoother distribution than those obtained by traditional algorithms. The root mean squares of the fitting errors between the tidal height model and the observations reduce from 11 cm to 1.3 cm.     

A Case Study: Heavy Rainfall Event Comparison Between Daily Satellite Rainfall Estimation Products with IMD Gridded Rainfall Over Peninsular India During 2015 Winter Monsoon

India Meteorological department (IMD) used INSAT-3D Metrological Satellite Imager data to drive two type rainfall estimation products viz-Hydro Estimate (HE) and INSAT Multi-Spectral Rainfall Algorithm (IMSRA) on half hourly rainfall rate and daily accumulated rainfall in millimeter (mm). Integrated Multi-Satellite Retrieval for GPM (IMERG) product is being derived by NASA and JAXA by using Global Precipitation Mission (GPM) satellites data. IMSRA and GPM (IMERG) are gridded data at 10 km spatial resolution and HE is available at pixel level (4 km at Nadir). IMD provides gridded rainfall data at 0.25° × 0.25° resolution which is based on wide coverage of 6955 actual observation. In present study, validation of INSAT-3D based Hydro Estimator (HE), INSAT Multi-Spectral Rainfall Algorithm (IMSRA) and Integrated Multi-Satellite Retrieval for GPM (IMERG) of Global Precipitation Mission (GPM) satellites are carried out with IMD gridded data set for heavy rainfall event during winter monsoon, over peninsular India (November–December 2015). In validation, Nash–Sutcliffe efficiencies (NSE), RMSE, Correlation, Skilled scores are calculated at grid level for heavy and very heavy rain categories and the values of NSE of HE (? 32.36, ? 3.12), GPM (? 68.67, ? 2.39) and IMSRA (? 0.02, 0.28) on 16th November 2015 and HE (? 13.65, ? 1.69), GPM (? 43.79, ? 2.94) and IMSRA (? 1.08, ? 1.60) on 2nd Dec 2015, for heavy and very heavy rainfall. On both days, HE is showing better rainfall estimate compare to GPM for Heavy rainfall and GPM showing better estimation for very heavy rainfall events. In all the cases IMSRA is underestimating, if daily rain fall exceeded 75 mm.     

Impact of GPS radio occultation data assimilation on regional weather predictions

The impact of GPS radio occultation (RO) data assimilation on severe weather predictions in East Asia is introduced and reviewed. Both the local observation operator that assimilates the retrieved refractivity as local point measurement, and the nonlocal observation operator that assimilates the integrated retrieved refractivity along a straight raypath have been utilized in WRF 3DVAR to improve the initial analysis of the model. A general evaluation of the impact of these approaches on Asian regional analysis and daily prediction is provided in this paper. In general, the GPS RO data assimilation may improve prediction of severe weather such as typhoons and Mei-yu systems when COSMIC data were available, ranging from several points in 2006 to a maximum of about 60 in 2007 and 2008 in this region. Based on a number of experiments, regional model predictions at 5 km resolution were not significantly influenced by different observation operators, although the nonlocal observation operator sometimes results in slightly better track forecast. These positive impacts are seen not only in typhoon track prediction but also in prediction of local heavy rainfall associated with severe weather over Taiwan. The impact of 56 GPS RO soundings on track prediction of Cyclone Gonu (2007) over the Indian Ocean is also appealing when compared to other tracks assimilated with different observations. From a successive evaluation of skill scores for real-time forecasts on Mei-yu frontal systems operationally conducted over a longer period and predictions of six typhoons in 2008, assimilation of GPS RO data appears to have some positive impact on regional weather predictions, on top of existent assimilation with all other observations.