Investigation of vertical profile of rain microstructure at Ahmedabad in Indian tropical region

The microstructure of rain has been studied with observations using a vertical looking Micro Rain Radar (MRR) at Ahmedabad (23.06°N, 72.62°E), a tropical location in the Indian region. The rain height, derived from the bright band signature of melting layer of radar reflectivity profile, is found to be variable between the heights 4600 m and 5200 m. The change in the nature of rain, classified on the basis of radar reflectivity, is also observed through the MRR. It has been found that there are three types of rain, namely, convective, mixed and stratiform rain, prevailing with different vertical rain microstructures, such as, Drop Size Distribution (DSD), mean drop size, rain rate, liquid water content and average fall speed of the drops at different heights. It is observed that the vertical DSD profile is more inhomogeneous for mixed and stratiform type rain than for convective type rain. It is also found that the large number of drops of size <0.5 mm is present in convective rain whereas in stratiform rain, drops concentration is appreciable up to 1 mm. A comparison of measurements taken by ground based Disdrometer and that from the 200 m level obtained from MRR shows good agreement for rain rate and DSD at smaller rain rate values. The results may be useful for understanding rain structures over this region.     

Rainfall estimation over a Mediterranean region using a method based on various spectral parameters of SEVIRI-MSG

The ultimate objective of this paper is the estimation of rainfall over an area in Algeria using data from the SEVIRI radiometer (Spinning Enhanced Visible and Infrared Imager). To achieve this aim, we use a new Convective/Stratiform Rain Area Delineation Technique (CS-RADT). The satellite rainfall retrieval technique is based on various spectral parameters of SEVIRI that express microphysical and optical cloud properties. It uses a multispectral thresholding technique to distinguish between stratiform and convective clouds. This technique (CS-RADT) is applied to the complex situation of the Mediterranean climate of this region. The tests have been conducted during the rainy seasons of 2006/2007 and 2010/2011 where stratiform and convective precipitation is recorded. The developed scheme (CS-RADT) is calibrated by instantaneous meteorological radar data to determine thresholds, and then rain rates are assigned to each cloud type by using radar and rain gauge data. These calibration data are collocated with SEVIRI data in time and space.     

A multilayer perceptron and multiclass support vector machine based high accuracy technique for daily rainfall estimation from MSG SEVIRI data

The current paper introduces a new multilayer perceptron (MLP) and support vector machine (SVM) based approach to improve daily rainfall estimation from the Meteosat Second Generation (MSG) data. In this study, the precipitation is first detected and classified into convective and stratiform rain by two MLP models, and then four multi-class SVM algorithms were used for daily rainfall estimation. Relevant spectral and textural input features of the developed algorithms were derived from the spectral MSG SEVIRI radiometer channels. The models were trained using radar rainfall data set colected over north Algeria. Validation of the proposed daily rainfall estimation technique was performed by rain gauge network data set recorded over north Algeria. Thus, several statistical scores were calculated, such as correlation coefficient (r), root mean square error (RMSE), mean error (Bias), and mean absolute error (MAE). The findings given by: (r = 0.97, bias = 0.31 mm, RMSE = 2.20 mm and MAE = 1.07 mm), showed a quite satisfactory relationship between the estimation and the respective observed daily precipitation. Moreover, the comparison of the results with those of two advanced techniques based on random forests (RF) and weighted ‘k’ nearest neighbor (WkNN) showed higher accuracy obtained by the proposed model.     

Development of a regional rain retrieval algorithm for exclusive mesoscale convective systems over peninsular India

The present study emphasize the development of a region specific rain retrieval algorithm by taking into accounts the cloud features. Brightness temperatures (T_bs) from various TRMM Microwave Imager (TMI) channels are calibrated with near surface rain intensity as observed from the TRMM – Precipitation Radar. It shows that T_b–R relations during exclusive-Mesoscale Convective System (MCS) events have greater dynamical range compared to combined events of non-MCS and MCS. Increased dynamical range of T_b–R relations for exclusive-MCS events have led to the development of an Artificial Neural Network (ANN) based regional algorithm for rain intensity estimation. By using the exclusive MCSs algorithm, reasonably good improvement in the accuracy of rain intensity estimation is observed. A case study of a comparison of rain intensity estimation by the exclusive-MCS regional algorithm and the global TRMM 2A12 rain product with a Doppler Weather Radar shows significant improvement in rain intensity estimation by the developed regional algorithm.     

Time dependent charging of layer clouds in the global electric circuit

There is much observational data consistent with the hypothesis that the ionosphere-earth current density (J_z) in the global electric circuit, which is modulated by both solar activity and thunderstorm activity, affects atmospheric dynamics and cloud cover. One candidate mechanism involves J_z causing the accumulation of space charge on droplets and aerosol particles, that affects the rate of scavenging of the latter, notably those of Cloud Condensation Nuclei (CCN) and Ice Forming Nuclei (IFN) (  and ). Space charge is the difference, per unit volume, between total positive and total negative electrical charge that is on droplets, aerosol particles (including the CCN and IFN) and air ions. The cumulative effects of the scavenging in stratiform clouds and aerosol layers in an air mass over the lifetime of the aerosol particles of 1–10 days affects the concentration and size distribution of the CCN, so that in subsequent episodes of cloud formation (including deep convective clouds) there can be effects on droplet size distribution, coagulation, precipitation processes, and even storm dynamics.     

Variation of rain drop size distribution with rain rate at a few coastal and high altitude stations in southern peninsular India

Rain drop size distribution (DSD) was measured at four places in Southern India {Thiruvananthapuram, Kochi, Munnar and Sriharikota (SHAR)} using a Joss–Waldvogel (JW) impact type disdrometer. The data for each minute were corrected for dead time errors and rain rate was computed from the corrected data. The data for a whole month were then sorted according to rain rate (R) into several classes ranging from 0.1 to >100 mm/h. The average DSD in each class was computed, and the lognormal distribution function was fitted to the average. In all the cases, the function fitted the data very well. The fit parameters were found to have dependence on rain rate. The total number of drops (N_T), the geometric mean diameter (D_g) and the standard geometric deviation (σ) were also computed from the fit parameters. The standard geometric deviation (σ) was found to be more or less constant with rain rate at all the sites and in all months. The other two parameters (N_T and D_g) were found to vary exponentially with rain rate except in Munnar, a high altitude station. At Thiruvananthapuram, in most of the months, N_T increased exponentially with rain rate up to some value of R, which was different in different months, and then remained more or less constant or decrease slightly. In all cases, the variation of N_T and D_g was such that N_TD_g³ increased linearly with rain rate.     

An Artificial Neural Network based approach for estimation of rain intensity from spectral moments of a Doppler Weather Radar

By using a Doppler Weather Radar (DWR) at Shriharikota (13.66°N & 80.23°E), an Artificial Neural Network (ANN) based technique is proposed to improve the accuracy of rain intensity estimation. Three spectral moments of a Doppler spectra are utilized as an input data to an ANN. Rain intensity, as measured by the tipping bucket rain gauges around the DWR station, are considered as a target values for the given inputs. Rain intensity as estimated by the developed ANN model is validated by the rain gauges measurements. With the help of a developed technique, reasonable improvement in the estimation of rain intensity is observed. By using the developed technique, root mean square error and bias are reduced in the range of 34–18% and 17–3% respectively, compared to Z–R approach.     

Investigation of GPS-TEC measurements using ANN method indicating seismo-ionospheric anomalies around the time of the Chile (Mw = 8.2) earthquake of 01 April 2014

In 94 km NW of Iquique in Chile (19.610°S, 70.776°W) a powerful earthquake of M_w = 8.2 took place at 23:46:47 UTC (20:46:47 LT) on April 01, 2014. Using GPS-TEC (Total Electron Content) measurements, potential unusual variations around the time and location of the Chile earthquake have been detected based on the median and Artificial Neural Network (ANN) methods. The indices D_st, K_p, A_p and F10.7 were used to distinguish pre-earthquake anomalies from the other anomalies related to the solar-geomagnetic activities. Using the median method, striking anomalies in time series of TEC data are observed 4 days before the earthquake at 14:00 and 16:00 UTC. The ANN method detected a number of anomalies, 4 (02:00 and 16:00 UTC) and 13 (24:00 UTC) days preceding the earthquake. The results indicate that the ANN method due to its capability of non linear learning is quite promising and deserves serious attention as a robust predictor tool for seismo-ionospheric anomalies detection.     

Propagation of secondary antiprotons and cosmic rays in the Galaxy

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper, we discuss one possibility to overcome these difficulties. Using the measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local “unprocessed” component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z ⩽ 28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three “radioactive clock” isotopes in CR, ¹⁰Be, ²⁶Al, ³⁶Cl, are all consistent and indicate a halo size z_h ∼ 4 kpc based on the most accurate data taken by the ACE spacecraft.     

Seismo ionospheric anomalies in Turkey associated with Mw ≥ 6.0 earthquakes detected by GPS stations and GIM TEC

Earthquake (EQ) anomalies in the form of enhancement and depletion in ionospheric Total Electron Content (TEC) from Global Positioning System (GPS) may considerably alarm about short and long term precursors of the impending main shock. In this paper, TEC anomalies are investigated from permanent GPS ground-stations in Turkey associated to M_w ≥ 6.0 EQs occurred in 2011–2012. Temporal and spatial analyses of TEC at 2 h sampling have shown significant evidences about EQ induced ionospheric anomalies during 10–14 h of UT (Universal Time) within 5 days before M_w 6.0 Greece, and M_w 7.1, Turkish EQ. Spatial analyses have manifested arrival of TEC anomalies at UT = 10 h to epicenter of both EQs, which linger above epicenter during UT = 12–14 h and left seismogenic zone after UT = 14 h before every EQ during K_p < 3 and Dst = 0 nT. Meanwhile, a geomagnetic storm (Dst < -100 nT) induce perturbation two days after the M_w 7.1 Turkish EQ, showing no relation with epicenter during spatial analysis. It also shows that TEC can be useful to distinguish geomagnetic storm variations to successfully detect EQ precursors. These anomalies during quiet storm (K_p < 3; Dst = 0 nT) conditions may be effective to link the lithosphere and ionosphere in severe seismic zones to detect EQ precursors before future EQs. Interpretation of TEC anomalies and it enhancements over EQ epicenters during UT = 12–14 h for both EQs have shown that EQs anomalies only occurred in particular time. Whereas, geomagnetic storm effect occurred during whole abnormal day over the Earth.     

Simulation of depth–dose distributions for various ions in polyethylene medium

Study of depth–dose distributions for intermediate energy ion beams in tissue-like media such as polyethylene (CH₂)_n provides a good platform for further improvements in the fields of hadrontherapy and space radiation shielding. The depth–dose distributions for ¹²C ions at various energies and for light and intermediate ion beams (³He, ¹⁶O, ²⁰Ne and ²⁸Si) as well as for heavy ions ⁵⁶Fe in polyethylene were estimated by using simulation toolkit: Geant4. Calculations were performed mainly by considering two different combinations of standard electromagnetic (EM), binary cascade (BIC), statistical multifragmentation (SMF) and Fermi breakup (FB) models. The energies of the ion beams were selected to achieve the Bragg peaks at predefined position (∼60 mm) and as per their availability. Variations of peak-to-entrance ratio (from 7.44 ± 0.05 to 8.87 ± 0.05), entrance dose (from 2.89 ± 0.01 to 203.71 ± 0.63 MeV/mm) and entrance stopping power (from 3.608 to 208.858 MeV/mm, calculated by SRIM) with atomic number (Z) were presented in a systematic manner. The better peak-to-entrance ratio and less entrance dose in the region Z = 2 to 8 (i.e. ³He to ¹⁶O) may provide the suitability of the ion beams for hadrontherapy.     

An evaluation of the IRI-2007 storm time model at low latitude stations

This paper discusses the ability of the International Reference Ionosphere IRI-2007 storm time model to predict foF2 ionospheric parameter during geomagnetic storm periods. Experimental data (based on availability) from two low latitude stations: Vanimo (geographic coordinates, 2.7 °S, 141.3 °E, magnetic coordinates, 12.3 °S, 212.50 °E) and Darwin (geographic coordinates, 12.45 °S, 130.95 °E, magnetic coordinates, 22.9 °S, 202.7 °E) during nine storms that occurred in 2000 (Rz₁₂ = 119), 2001(Rz₁₂ = 111) and 2003 (Rz₁₂ = 64) are compared with those obtained by the IRI-2007 storm model. The results obtained show that the percentage deviation between the experimental and IRI predicted foF2 values during these storm periods is as high as 100% during the main and recovery phases. Based on the values of “relative deviation module mean” (RDMM) obtained (i.e. between 0.08 and 0.60), it is observed that there is a reasonable to poor agreement between measured foF2 values and the IRI-storm model prediction values during main and recovery phases of the storms under investigation. As a result, in addition to other studies that have been carried out from different sectors, more studies are required to be carried out. This will enable IRI community to improve on the present performance of the model. In general the IRI-storm model predictions follow normal trend of the foF2 measured values but does not reproduce well the measured values.     

Estimation of sediment settling velocity in estuarine and coastal waters using optical remote sensing data

A robust method has been developed for estimating sediment settling velocity (w_s) from high resolution optical remote sensing data in estuarine, coastal and harbor waters. This method estimates settling velocity as a function of the drag coefficient (C_d), Reynolds number (Re), grain size (D₅₀), specific gravity (Δ_SG) and grain shape (in terms of the Corey Shape Factor – CSF). These parameters were derived from the particulate inherent optical properties such as backscattering (b_bp), beam attenuation (c_p), suspended sediment concentration and turbidity using Landsat 8 OLI and HICO data. Preliminary results for the Gulf of Cambay in the eastern Arabian Sea and Yangtze river estuary in the East China Sea, showed that satellite-retrieved settling velocities (m?s^?1) varied from very low values in clear oceanic waters, intermediate values in coastal waters, to very high values in river plumes and sediment-laden coastal waters. The remote sensing retrievals of sediment properties and their settling velocities were generally consistent with the field and laboratory results, which indicate that the proposed methodology will have important implications in various coastal engineering, environmental and management studies.     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

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.     

Plasma chemical reactions in C2H2/N2, C2H4/N2, and C2H6/N2 gas mixtures of a laboratory dielectric barrier discharge

Plasma chemical reactions in C₂H₂/N₂, C₂H₄/N₂, and C₂H₆/N₂ gas mixtures have been studied by means of mass spectrometry at a medium pressure of 300 mbar in a laboratory dielectric barrier discharge. A major reaction scheme is production of larger hydrocarbons like C_nH_m with n up to 12 including formation of functional CN groups.     

Influence of the equinoctial asymmetry on the semi-annual anomaly in NmF2 near noon in the northern geographic hemisphere: A statistical study

Hourly values of NmF2 measured by 72 ionosondes near noon from 1957 to 2012 at low and middle geomagnetic latitudes of the northern geographic hemisphere are used to study the spring and autumn semi-annual anomalies in NmF2. The spring/summer, autumn/summer, spring/winter, and autumn/winter geomagnetically quiet NmF2 ratios, X_M, X_S, Y_M, and Y_S, respectively, taken near noon over each ionosonde for approximately the same winter, spring, summer, and autumn solar activity conditions are analyzed. The probabilities of occurrences of X_M, X_S, Y_M, and Y_S in intervals of change of these parameters, together with its most frequent and the mean expected values, and the probabilities to observe the F2-region spring and autumn semi-annual anomalies are calculated and studied for the first time for low, moderate, and high solar activity conditions. The geomagnetic latitude range between 10 and 55 degrees, where the ionosondes are located, is divided into 9 intervals of the same length of 5 degrees, the statistical parameters are averaged over each 5 degree interval, and the trends in these averaged statistical characteristics of the NmF2 spring and autumn semi-annual anomalies are calculated and studied for the first time. It is shown that the NmF2 equinoctial asymmetry can significantly affect the statistical parameters of the semi-annual anomaly in NmF2.     

Prediction of clear-sky biologically effective erythematic radiation (EER) From global solar radiation (250–2800 nm) at Cairo,Egypt

The hourly and daily measured clear-sky global solar radiation (G) and biologically important effective erythematic radiation (EER) incident on a horizontal surface at Cairo, Egypt (latitude 30° 05′ N & Longitude 31° 15′ E), during the period from January 1995 to December 2005 are used in this paper. The relationship between daily integrated totals of EER and the daily totals of broadband global solar radiation (250–2800 nm) is established. The temporal variability of the percentage ratio of the total daily erythema to total daily broadband solar global irradiation (EER/G) is determined. The monthly and the seasonal averages of the extraterrestrial UVB solar radiation, Mesurad and estímated UVB solar radiation and clearness index K_tUVB of UVB radiation are discussed. The average monthly mean variation of slant ozone (Z) and UVB transmission (K_tUVB) at the present work are found. The two variables show an opposite seasonal behavior, and the average monthly of slant ozone column and UVB transmission values shows the relationship between them in a clearer way than those of daily values. The estimated values of UVB solar radiation a good agreement with the measured values of the UVB solar radiation, the difference between the estimated and measured values of UVB solar radiation varies from 1.2% to 2.8%. The effect of the annual cycles of solar zenith angle (SZA) and total column ozone (TCO) on the ratios (EER/G) are presented and the correction factors are determined for removal of the ozone cycle. The seasonal variability of EER/G is also discussed. The effect of the annual cycles of solar zenith angle (SZA) and total column ozone (TCO) on the ratios (EER/G) is presented and the correction factors are determined for removal of the ozone cycle.     

Longitudinal behaviors of the IRI-B parameters of the equatorial electron density profiles retrieved from FORMOSAT-3/COSMIC radio occultation measurements

The electron density profiles in the bottomside F₂-layer ionosphere are described by the thickness parameter B0 and the shape parameter B1 in the International Reference Ionosphere (IRI) model. We collected the ionospheric electron density (N_e) profiles from the FORMOSAT-3/COSMIC (F3/C) radio occultation measurements from DoY (day number of year) 194, 2006 to DoY 293, 2008 to investigate the daytime behaviors of IRI-B parameters (B0 and B1) in the equatorial regions. Our fittings confirm that the IRI bottomside profile function can well describe the averaged profiles in the bottomside ionosphere. Analysis of the equatorial electron density profile datasets provides unprecedented detail of the behaviors of B0 and B1 parameters in equatorial regions at low solar activity. The longitudinal averaged B1 has values comparable with IRI-2007 while it shows little seasonal variation. In contrast, the observed B0 presents semiannual variation with maxima in solstice months and minima in equinox months, which is not reproduced by IRI-2007. Moreover, there are complicated longitudinal variations of B0 with patterns varying with seasons. Peaks are distinct in the wave-like longitudinal structure of B0 in equinox months. An outstanding feature is that a stable peak appears around 100°E in four seasons. The significant longitudinal variation of B0 provides challenges for further improving the presentations of the bottomside ionosphere in IRI.     

Occurrence of sporadic-E layer over the ionospheric station of Rome: Analysis of data for thirty-two years

Hourly systematic measurements of the highest frequency reflected by the sporadic-E layer (foEs) recorded at the Rome ionospheric observatory (Italy, 41.8° N, 12.5° E), were considered during the period January 1976–December 2007, to calculate the percentage of occurrence of sporadic-E layer with frequencies foEs greater than a given threshold value f_T, P(foEs > f_T).