Reliability of satellite-based precipitation products in capturing extreme precipitation indices over Iran

The goal of this study is to assess the performance of four widely-used satellite precipitation products in capturing extreme precipitation indices across Iran over the period 2001–2018; these products include GPM IMERG (Integrated Multi-Satellite Retrievals for Global Precipitation Measurement), TRMM 3B42 (Tropical Rainfall Measuring Mission), CHIRPS (Climate Hazards Center InfraRed Precipitation with Station data), and PERSIANN-CDR (Precipitation Estimation From Remotely Sensed Information Using Artificial Neural Networks-Climate Data Record). For this aim, a national gridded precipitation dataset was developed using a dense network of rain gauges as a reference dataset. The results suggest that the IMERG product outperforms the other three precipitation products in capturing extreme precipitation indices both temporally and spatially. TRMM 3B42 data show promising results in identifying many extreme indices, while the CHIRPS and PERSIANN-CDR products show less performance in accurately generating many of the extreme precipitation indices.     

High temporal resolution global PWV dataset of 2005–2016 by using a neural network approach to determine the mean temperature of the atmosphere

Atmospheric water vapour plays an important role in phenomena related to the global hydrologic cycle and climate change. However, the rapid temporal–spatial variation in global tropospheric water vapour has not been well investigated due to a lack of long-term, high-temporal-resolution precipitable water vapour (PWV). Accordingly, this study generates an hourly PWV dataset for 272 ground-based International Global Navigation Satellite System (GNSS) Service (IGS) stations over the period of 2005–2016 using the zenith troposphere delay (ZTD) derived from global-scale GNSS observation. The root mean square (RMS) of the hourly ZTD obtained from the IGS tropospheric product is approximately 4 mm. A fifth-generation reanalysis dataset of the European Centre for Medium-range Weather Forecasting (ECMWF ERA5) is used to obtain hourly surface temperature (T) and pressure (P), which are first validated with GNSS synoptic station data and radiosonde data, respectively. Then, T and P are used to calculate the water vapour-weighted atmospheric mean temperature (Tm) and zenith hydrostatic delay (ZHD), respectively. T and P at the GNSS stations are obtained via an interpolation in the horizontal and vertical directions using the grid-based ERA5 reanalysis dataset. Here, Tm is calculated using a neural network model, whereas ZHD is obtained using an empirical Saastamoinen model. The RMS values of T and P at the collocated 693 radiosonde stations are 1.6 K and 3.1 hPa, respectively. Therefore, the theoretical error of PWV caused by the errors in ZTD, T and P is on the order of approximately 2.1 mm. A practical comparison experiment is performed using 97 collocated radiosonde stations and 23 GNSS stations equipped with meteorological sensors. The RMS and bias of the hourly PWV dataset are 2.87/?0.16 and 2.45/0.55 mm, respectively, when compared with radiosonde and GNSS stations equipped with meteorological sensors. Additionally, preliminary analysis of the hourly PWV dataset during the EI Niño event of 2014–2016 further indicates the capability of monitoring the daily changes in atmospheric water vapour. This finding is interesting and significant for further climate research.     

Integrated drought monitoring index: A tool to monitor agricultural drought by using time-series datasets of space-based earth observation satellites

In this study, integrated drought monitoring index (IDMI) was proposed as a tool to assess and monitor the spatio-temporal dynamics of agricultural drought during the northeast monsoon season for the period from 2000 to 2016 in Tamil Nadu state, south-eastern part of Indian peninsula. The IDMI is characterized as the principal component of precipitation condition index (PCI), soil moisture condition index (SMCI), temperature condition index (TCI), and vegetation condition index (VCI) derived from time-series satellite observations of climate hazards group infra-red precipitation with stations (CHIRPS), European space agency climate change initiative (ESA-CCI) and moderate resolution imaging spectroradiometer (MODIS). The study shows that in the year 2016, about 44.4 and 17.8% of Tamil Nadu state was under extreme and severe drought conditions, respectively. Sensitivity analysis of the study shows that PCI is the most influential parameter to IDMI, followed by VCI and TCI. The validation of IDMI with 3-month standardized precipitation index (SPI) by using Pearson correlation test shows a strong positive correlation between IDMI and 3-month SPI with correlation coefficient (r) value of 0.73 and 0.77 for the wet (2005) and dry year (2016), respectively. The study clearly demonstrates the potential of IDMI derived from time-series datasets of earth observation satellites as a tool in assessment and monitoring of spatio-temporal dynamics of agricultural drought. The proposed IDMI could be effectively used as a reliable tool to monitor agricultural drought and develop its mitigation strategies to minimise the adverse effects of drought on agriculture, water resources, and livelihoods of the people.     

Evaluation of three high-resolution satellite precipitation estimates: Potential for monsoon monitoring over Pakistan

Multi-sensor precipitation datasets including two products from the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) and estimates from Climate Prediction Center Morphing Technique (CMORPH) product were quantitatively evaluated to study the monsoon variability over Pakistan. Several statistical and graphical techniques are applied to illustrate the nonconformity of the three satellite products from the gauge observations. During the monsoon season (JAS), the three satellite precipitation products captures the intense precipitation well, all showing high correlation for high rain rates (>30 mm/day). The spatial and temporal satellite rainfall error variability shows a significant geo-topography dependent distribution, as all the three products overestimate over mountain ranges in the north and coastal region in the south parts of Indus basin. The TMPA-RT product tends to overestimate light rain rates (approximately 100%) and the bias is low for high rain rates (about ±20%). In general, daily comparisons from 2005 to 2010 show the best agreement between the TMPA-V7 research product and gauge observations with correlation coefficient values ranging from moderate (0.4) to high (0.8) over the spatial domain of Pakistan. The seasonal variation of rainfall frequency has large biases (100–140%) over high latitudes (36N) with complex terrain for daily, monsoon, and pre-monsoon comparisons. Relatively low uncertainties and errors (Bias ±25% and MAE 1–10 mm) were associated with the TMPA-RT product during the monsoon-dominated region (32–35N), thus demonstrating their potential use for developing an operational hydrological application of the satellite-based near real-time products in Pakistan for flood monitoring.     

Effects observed in the equatorial and low latitude ionospheric F-region in the Brazilian sector during low solar activity geomagnetic storms and comparison with the COSMIC measurements

The main objective of the present investigation has been to compare the ionospheric parameters (NmF2 and hmF2) observed by two ground-based ionospheric sounders (one at PALMAS- located near the magnetic equator and the other at Sao Jose dos Campos-located in the low-latitude region) in the Brazilian sector with that by the satellite FORMOSAT-3/COSMIC radio occultation (RO) measurements during two geomagnetic storms which occurred in December 2006 and July 2009. It should be pointed out that in spite of increasing the latitude (to 10°) and longitude (to 20°) around the stations; we had very few common observations. It has been observed that both the peak electron density (NmF2) and peak height (hmF2) observed by two different techniques (space-borne COSMIC and ground-based ionosondes) during both the geomagnetic storm events compares fairly well (with high correlation coefficients) at the two stations in the Brazilian sector. It should be pointed out that due to equatorial spread F (ESF) in the first storm (December 2006) and no-reflections from the ionosphere during nighttime in the second storm (July 2009), we had virtually daytime data from the two ionosondes.     

A balloon-borne very long baseline interferometry experiment in the stratosphere: Systems design and developments

The balloon-borne very long baseline interferometry (VLBI) experiment is a technical feasibility study for performing radio interferometry in the stratosphere. The flight model has been developed. A balloon-borne VLBI station will be launched to establish interferometric fringes with ground-based VLBI stations distributed over the Japanese islands at an observing frequency of approximately 20?GHz as the first step. This paper describes the system design and development of a series of observing instruments and bus systems. In addition to the advantages of avoiding the atmospheric effects of absorption and fluctuation in high frequency radio observation, the mobility of a station can improve the sampling coverage (“uv-coverage”) by increasing the number of baselines by the number of ground-based counterparts for each observation day. This benefit cannot be obtained with conventional arrays that solely comprise ground-based stations. The balloon-borne VLBI can contribute to a future progress of research fields such as black holes by direct imaging.     

Dynamics of the dayside aurora

Observations of auroral particles by rockets and satellites show the existence of structured “inverted V” type precipitation events with electron characteristic energies ranging up to several hundred eV in the post-noon dayside auroral region. Two station ground-based measurements from Cape Parry, NWT and Sachs Harbour, NWT, Canada from meridian scanners and auroral all-sky TV cameras are used in conjunction with rocket-borne electron measurements to study the temporal and spatial characteristics of the associated auroral emissions. The frequent occurrence of narrow, transient auroral arcs with lifetimes of 1–2 minutes, resulting from the inverted V events, is discussed.     

Evaluation of GPM sampling error over major basins in Indian subcontinent using bootstrap technique

The study investigates the evaluation and comparison of sampling error for the Global Precipitation Measurement (GPM) mission orbital data products by implementing a bootstrap technique over the two major basins in the Indian subcontinent i.e the Ganga and the Mahanadi basin. The relative sampling error evaluated over both the Ganges and Mahanadi basins showed commendable results thus giving the confidence to adopt the bootstrap technique to evaluate the sampling error. The region over India with large seasonal rainfall seems to have less sampling uncertainty and vice versa with some regions showing exceptions which might be due to the difference in precipitation variability and space-time correlation length. The scale dependence was verified for four grid sizes along with seasonal time scale. Results indicate that the relative sampling error estimates are inversely proportional to the scale of the grid size. The comparative study of evaluation of sampling uncertainty to different precipitation types resulted to have maximum sampling error in GPM Microwave Imager (GMI) in comparison to Dual Precipitation Radar (DPR) convective and DPR total precipitation. Thus, the comparable results of sampling uncertainty between the major basins in the Indian sub-continent provides the user a decision making criteria before utilizing the GPM orbital products in any applications.     

Application of ERA5 Reanalysis to the Construction of Initial Conditions for WACCM Simulations

This study uses ECMWF fifth-generation reanalysis, ERA5, which extends to the mesopause, to construct the Initial Conditions (IC) for WACCM (Whole Atmosphere Community Climate Model) simulations. Because the biases between ERA5 and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperature data are within ±5 K below the lower mesosphere, ERA5 reanalysis is used to construct IC in the lower atmosphere. Four experiments are performed to simulate a Stratospheric Sudden Warming (SSW) event from 5 to 15 February 2016. The simulation using the WACCM default climatic IC cannot represent the sharp meteorological variation during SSW. In contrast, the 0~4 d forecast results driven by ERA5-constructed IC is consistent with ERA5 reanalysis below the middle mesosphere. Comparing with WACCM climatology ICs scheme, the ICs constructing method based on ERA5 reanalysis can obtain 67%, 40%, 22%, 4% and 6% reduction of temperature forecast RMSE at 10 hPa, 1 hPa, 0.1 hPa, 0.01 hPa and 0.001 hPa respectively. However, such improvement is not shown in the lower thermosphere.      

Recent Advances in Observation and Research of the Chinese Meridian Project

The Chinese Meridian Space Weather Monitoring Project (Meridian Project) is a ground-based geospace monitoring chain in China. It consists of 15 ground-based observation stations located roughly along 120°E longitude and 30°N latitude. In recent two years, using data from the Meridian Project, significant progress has been made in space weather and space physics research. These advances are mainly in four aspects:regional characteristics of space environment above China or along 120°E meridian line, coupling between space spheres at different heights and different physical processes, space weather disturbance and its propagation along the meridian chain, and space weather effects on ground technical facilities.      

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Conditional sea surface statistics and their impact on geophysical sea surface parameters retrieved from SAR altimetry signals

In this study we present an analytical formulation of synthetic-aperture radar (SAR) altimetry signals including narrow banded nonlinear wave fields and conditional statistics between wave elevation displacements, horizontal wave slopes and vertical wave particle velocities. Considering the wave elevation displacements coskewness with respect to horizontal slopes leads to an analytical formulation of the electromagnetic bias within a SAR-mode altimeter stack. This formulation can be either parametrized by the significant wave height (SWH) and mean wave steepness, or in terms of the variance of vertical wave velocities. The effect of conditional vertical wave particle velocity variances with respect to the observed horizontal wave slopes close to nadir incidence angles leads to an effective reduction of the azimuth blurring of SAR-mode stacks. We present here a formulation of this effect by examining JONSWAP ocean wave spectra. In most cases this effect reduces the azimuth blurring by 10% to 30%. Additionally we investigate the effect of a nonlinear wave elevation displacement probability density function (PDF) on estimated geophysical parameters. We were able to show that including an elevation displacement skewness of 0.13 improves significantly the SWH consistency between altimetry and ECMWF Reanalysis v5 ERA5 results.All of these effects are validated with respect to ERA5 model data in the North East Atlantic region and in situ data located in the German Bight and Baltic Sea.The developed model can be used in both SAR and conventional altimetry retrackers.     

Refinements of a classical technique of airglow spectroscopy

We have performed ground-based airglow observations with a tilting filter spectrometer to determine rotational temperatures in the mesopause region since 1984. Although quite simple in principle, the success of the technique depends on a number of non-trivial details, including corrections for measurements under less favorable conditions, as, for instance, full moon nights. Successive improvements have finally led to nearly continuous observations, since August 1997, with more than 200 nights of useful data per year. The technique is described in some detail and recent results are discussed.     

Determination of crustal deformations based on GPS observing-session duration in Marmara region,Turkey

This paper reports on the estimation of changes of a GPS-based crustal deformation field as function of the observing session duration. This investigation was carried out on the Marmara Continuous GPS Network (Marmara region) by using GPS data collected on a set of 10 stations during the period 2002.50–2005.47 The GPS observations were processed in the ITRF 2005 reference frame using Bernese v5.0 software.     

民航GPS地基区域完好性监视系统设计与实现

结合中国地域和空中交通流量分布的特点,提出了中国民航全球卫星定位系统GPS(Global Positioning System)地基区域完好性监视系统GRIMS(Ground-based Regional Integrity Monitoring System)的概念.一个包含7个监测站的地面监测网监测视界内的GPS卫星,判断其可用状况,并将判断结果和重要的监测数据发送到主控站.主控站计算中国范围内的接收机自主完好性监测RAIM(Receiver Autonomous Integrity Monitoring)空洞,同时推算全国民航飞机的实时位置,并向所有位于RAIM空洞中的飞机及时发送完好性告警信息.从卫星出现故障到飞机收到告警信息之间的告警延迟小于7s.该系统已经应用于中国民航,在保证飞行安全的同时提高了飞行效率.      

Large-scale traveling ionospheric disturbances of auroral origin according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered using measurements of total electron content (TEC) during the magnetic storms on October 29–31, 2003, and on November 7–11, 2004, had been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region f_oF2 were used for this purpose. The variations of TEC and f_oF2 were similar for all events mentioned above. The previous assumption that the ionospheric region with vertical extension from 150 to 200 km located near the F-layer maximum mainly contributes to the TEC variations was confirmed for the cases when the electron density disturbance at the F region maximum was not more than 50%. However, this region probably becomes vertically more extended when the electron density disturbance in the ionospheric F region is about 85%.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

夏季强对流活动对东亚低纬电离层不规则体影响的事件分析

利用海南台站和东南亚地区的多种地基和天基观测手段,对2014年7月28日夜间观测到的东亚低纬F区不规则体事件的时空变化及其物理过程进行分析。结果表明,海南台站观测到了罕见的长时间持续的F区电离层不规则体,不同手段观测到的电离层不规则体存在明显的形态差异。不同台站观测到的电离层不规则体活动存在明显的差异。海南台站经度区南北异常峰附近的TEC起伏活动在日落后至午夜附近明显增强,在午夜后明显减弱。C/NOFS卫星轨迹午夜后逐渐接近于磁赤道,且处于较低高度上,几乎总会观测到弱等离子体扰动/泡的发生,与该区域地基观测的弱电离层不规则体活动存在明显的联系。SWARM卫星在黎明海南台站附近经度区仍观测到较强的赤道异常双峰结构,且西侧异常峰区附近仍存在明显的等离子体密度耗空/泡结构。海南台站西侧磁赤道区附近（中南半岛）强对流活动（MCC）激发的重力波种子扰动对东亚低纬区等离子体泡及准周期结构的产生发挥了重要作用。      

Detection of aerosol pollution sources during sandstorms in Northwestern China using remote sensed and model simulated data

The present paper has used a comprehensive approach to study atmosphere pollution sources including the study of vertical distribution characteristics, the epicenters of occurrence and transport of atmospheric aerosol in North-West China under intensive dust storm registered in all cities of the region in April 2014. To achieve this goal, the remote sensing data using Moderate Resolution Imaging Spectroradiometer satellite (MODIS) as well as model-simulated data, were used, which facilitate tracking the sources, routes, and spatial extent of dust storms. The results of the study have shown strong territory pollution with aerosol during sandstorm. According to ground-based air quality monitoring stations data, concentrations of PM₁₀ and PM_2.5 exceeded 400?μg/m³ and 150?μg/m³, respectively, the ratio PM_2.5/PM₁₀ being within the range of 0.123–0.661. According to MODIS/Terra Collection 6 Level-2 aerosol products data and the Deep Blue algorithm data, the aerosol optical depth (AOD) at 550?nm in the pollution epicenter was within 0.75–1. The vertical distribution of aerosols indicates that the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) 532?nm total attenuates backscatter coefficient ranges from 0.01 to 0.0001?km^?1?×?sr^?1 with the distribution of the main types of aerosols in the troposphere of the region within 0–12.5?km, where the most severe aerosol contamination is observed in the lower troposphere (at 3–6?km). According to satellite sounding and model-simulated data, the sources of pollution are the deserted regions of Northern and Northwestern China.