Soil moisture mapping over the Chinese Loess Plateau using ENVISAT/ASAR data

In this paper, we used the available algorithm for soil moisture estimation over LOPEX05 (the Loess Plateau land surface process Experiment (2005)) area. The available algorithm used ENVISAT/ASAR AP mode VV polarization observational data at a low incidence angle and ground measured soil moistures. The ground measurements were performed in the summer of the 2005 during the LOPEX05 field campaign. The validated results indicate that an average difference between the soil moistures estimated from the microwave remote sensing and ground measurements is less than 0.02 cm³/cm³, with a RMS error of 2.0%, and a maximum less than 0.04 cm³/cm³. The algorithm was applied to the surface soil moisture mapping later. The results show that this algorithm is suitable for monitoring soil moisture information of the agricultural fields over the Chinese Loess Plateau, when ground land cover situation and the resolution of imagery data are taken into account. However, we also find that there are large differences over the steep slope region, the edge of mesa. The results are not acceptable for surface soil moisture estimation in these regions. Thus, the surface soil moisture retrieval in the steep slope region of the Loess Plateau need to be further investigated in the future.     

Satellite ocean colour algorithm for Prochlorococcus, Synechococcus, and picoeukaryotes concentration retrieval in the South China Sea

An algorithm for retrieval of surface waters cell concentrations (in cell/ml) for three picophytoplankton components, Prochlorococcus (Pro), Synechococcus (Syn), and picoeukaryotes (Peuk) in the South China Sea (SCS), from ocean colour satellite data was developed and tested. Level 3 merged multisensor Ocean Colour Climate Change Initiative satellite data is used. Training is performed using in situ data on abundances of the three phytoplankton components. Several predictors derived from satellite reflectance data were tested. The regression form that assures the highest accuracy of the algorithm was chosen based on cross-validation (CV). According to the CV on test data subset, the algorithm performance is characterized by the r value 0.89, 0.72, and 0.73 and MAPD 38, 71 and 51% for Peuk, Pro, and Syn respectively. This is one of the few studies aimed at the Peuk, Pro, and Syn distribution research in the northern SCS using ocean colour satellite data. This is the only research providing algorithm with accuracy estimates of the Peuk, Pro, and Syn concentrations retrieval from the ocean colour data. Analysis of the developed algorithm allows us to conclude that both mechanisms (specific spectral features caused by pigments composition and spectrum features sensitive to general primary productivity, e.g. band ratios in 443–510?nm range and spectrum absolute values) are important for getting accurate information on the picophytoplankton composition.     

Study of the NGC 2770 interstellar medium through Hα, millimetric and optical polarimetric data of SN 2008D and SN 2007uy

We present the preliminary results of the recent H_α narrow-band imaging carried out for NGC 2770 with the Gran Telescopio Canarias (GTC) equipped with OSIRIS. We put the polarization measurements reported in Gorosabel et al. (2010) for SN 2007uy and SN 2008D in the context of the morphological information inferred from the H_α imaging. We estimated the orientation of the interstellar polarization (ISP) at the position of SN 2007uy and, most interestingly, at the site of SN 2008D which has been subject of an intensive debate due to its possible connection with Gamma-Ray Bursts (GRBs). The H_α imaging reveals a clumpy interstellar medium (ISM) composed of hundreds of compact emitting regions, for which we determined their sizes. The derived size for the H_α emitting region coincident with SN 2008D is consistent with the ISM cell size limits imposed by Gorosabel et al. (2010) based on millimetric data. A deeper data analysis is under way and will be published elsewhere. This article represents the first attempt to study the site of a possible GRB-like event combining millimetric, polarimetric and narrow-band data.     

Investigation of the relationship between ionospheric foF2 and earthquakes

Variations of the ionospheric F2 region critical frequency (foF2) have been investigated statistically before earthquakes during 1980–2008 periods in Japan area. Ionosonde data was taken from Kokubunji station which is in the earthquake preparation zone for all earthquakes. Standard Deviations and Inter-Quartile Range methods are applied to the foF2 data. It is observed that there are anomalous variations in foF2 before earthquakes. These variations can be regarded as ionospheric precursors and may be used for earthquake prediction.     

Evolution of periodic orbits near the Lagrangian point L2

A study of the evolution of the periodic and the quasi-periodic orbits near the Lagrangian point L₂, which is located to the right of the smaller primary on the line joining the primaries and whose distance from the more massive primary is greater than the distance between the primaries, in the framework of restricted three-body problem for the Sun–Jupiter, Earth–Moon (relatively large mass ratio) and Saturn–Titan (relatively small mass ratio) systems is made. Two families of periodic orbits around the smaller primary are identified using the Poincaré surface of section method – family I (initially elliptical, gradually becomes egg-shaped with the increase in the Jacobi constant C and elongated towards the more massive primary) and family II (initially egg-shaped orbits elongated towards L₂ and gradually becomes elliptical with the increase in C). The family I in the Sun–Jupiter and Saturn–Titan systems contains two separatrix caused by third-order and fourth-order resonances, while the Earth–Moon system has only one separatrix which is caused by third-order resonances. Also in the Sun–Jupiter and the Saturn–Titan systems, family I merge with family II, around Jacobian constant 3.0393 and 3.0163, respectively, while in the Earth–Moon system, family II evolves separately from two different branches. The two branches merge at C = 3.184515. In the Earth–Moon system, the family II contains a separatrix due to third-order resonances which is absent in the other two systems.     

Proposal of new models of the bottom-side B0 and B1 parameters for IRI

The time series of hourly electron density profiles N(h) obtained from 27 ionosonde stations distributed world-wide have been used to obtain N(h) average profiles on a monthly basis and to extract the expected bottom-side parameters that define the IRI profile under quiet conditions. The time series embrace the time interval from 1998 to 2006, which practically contains the entire solar cycle 23. The Spherical Harmonic Analysis (SHA) has been used as an analytical technique for modeling globally the B0 and B1 parameters as general functions on a spherical surface. Due to the irregular longitudinal distribution of the stations over the globe, it has been assumed that the ionosphere remains approximately constant in form for a given day under quiet conditions for a particular coordinate system. Since the Earth rotates under a Sun-fixed system, the time differences have been considered to be equivalent to longitude differences. The time dependence has been represented by a two-degree Fourier expansion to model the annual and semiannual variations and the year-by-year analyses of the B0 and B1 have furnished nine sets of spherical harmonic coefficients for each parameter. The spatial–temporal yearly coefficients have been further expressed as linear functions of Rz12 to model the solar cycle dependence. The resultant analytical model provides a tool to predict B0 and B1 at any location distributed among the used range of latitudes (70°N–50°S) and at any time that improves the fit to the observed data with respect to IRI prediction.     

Nighttime thermospheric meridional neutral winds inferred from ionospheric h′F and hpF2 data

Nighttime thermospheric meridional winds aligned to the magnetic meridian have been inferred using h′F and hpF2 ionosonde data taken from two equatorial stations, Manaus (2.9°S, 60.0°W, dip latitude 6.0°N) and Palmas (10.17°S, 48.2°W, dip latitude 6.2°S), and one low-latitude station, Sao Jose dos Campos (23.21°S, 45.86°W, dip latitude 17.26°S), during geomagnetic quiet days of August and September, 2002. Using an extension of the ionospheric servo model and a simple formulation of the diffusive vertical drift velocity, the magnetic meridional component of the thermospheric neutral winds is inferred, respectively, at the peak (hpF2) and at the base (h′F) heights of the F region over Sao Jose dos Campos. An approach has been included in the models to derive the effects of the electrodynamic drift over Sao Jose dos Campos from the time derivative of hpF2 and h′F observed at the equatorial stations. The magnetic meridional winds inferred from the two methods, for the months of August and September, are compared with winds calculated using the HWM-90 model and with measurements from Fabry–Perot technique. The results show varying agreements and disagreements. Meridional winds calculated from hpF2 ionospheric data (servo model) may produce errors of about 59 m/s, whereas the method calculated from the F-region base height (h′F) ionospheric data gives errors of about 69 m/s during the occurrence of equatorial spread-F.     

Seasonal-to-interannual variability of chlorophyll-a bloom timing associated with physical forcing in the Gulf of Cádiz

Seasonal-to-interannual variability of the winter-spring bloom in the Gulf of Cádiz, eastern North Atlantic, has been investigated using chlorophyll-a remote sensing (CHL). These data have been obtained from the GlobColour project; the temporal coverage extends from September 1997 to December 2010. In this study we develop a generic quantitative approach for describing the temporal variability in the shape of the winter-spring bloom within a region. Variability in both the timing and magnitude of the bloom in the basin has been evaluated as a function of physical properties in the water column such as Mixed Layer Depth (MLD, GODAS model), sea surface temperature (SST, from AVHRR radiometers), photosynthetically-active radiation (PAR, from ocean color data) and euphotic depth (Z_eu, from ocean color data). The analysis indicated that the timing, size and duration of the phytoplankton bloom in this area are largely controlled by both meteorological and oceanographic conditions at different scales; this means that it is likely to vary widely from one year to another.     

A new spectral index to detect Poaceae grass abundance in Mongolian grasslands

The objectives of the present study were to develop a new index based on remotely-sensed data for detecting the abundance of grasses in the family Poaceae, which has a high palatability for livestock in Mongolia, and to map the distribution of these grasses in the semi-arid Mongolian steppes. We measured ground-based spectral reflectance of pure plant leaves – including Poaceae grasses – and soils, as well as in-situ in the Mongolian grasslands. The hyper-spectral data, taken by a spectroradiometer, were converted into four multi-spectral bands (i.e., blue, green, red, and NIR) to simulate satellite-based imagery data. In order to magnify the characteristics of the spectral signal of Poaceae, NGBDI (Normalized Green-Blue Difference Index), NGRDI (Normalized Green-Red Difference Index), NDVI (Normalized Difference Vegetation Index), NNBDI (Normalized NIR-Blue Difference Index) were calculated from the four multi-spectral reflectance values. Poaceae Abundance Index (PAI) was derived by combining these four normalized difference indices. PAI was found out to be a good indicator to discriminate Poaceae grass from the other plant spectral data.     

Seismic deformation of the Mw 8.0 Wenchuan earthquake from high-rate GPS observations

High-rate GPS positioning has been recognized as a powerful tool in estimating epoch-wise station displacement which is particularly useful for seismology. In this study, station displacements during the 12 May 2008 M_w 8.0 Wenchuan earthquake are derived from the 1-Hz GPS data collected at a set of stations in China. The impacts of integer ambiguity resolution and station environment-dependent effects are investigated in order to yield more accurate results. The position accuracy of horizontal components of better than 1 cm suggests that GPS can sense the rapid position oscillation of about 2 cm in amplitude. Temporal and spatial analysis is applied to the surface displacement at station XANY and the characteristics of the movements due to Rayleigh and Love waves are detected and discussed. The comparison of GPS-derived displacement with relevant synthetic data computed based on a recently published rapture model shows a reasonable agreement in waveform. The various differences in amplitude need further investigation and also imply that rapture inversion might be improved if GPS-derived displacement is assimilated.     

Estimation of the elastic Earth parameters (h2, l2) using SLR data

We present results for the global elastic parameters h₂ and l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two satellites LAGEOS 1 and LAGEOS 2 observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analysed using different approaches. The analysis was done separately for the two satellites and approaches to estimate the two elastic parameters independently and together were performed. We do a sequential analysis and study the stability of the estimates as a function of length of the data set used. The adjusted final values for h₂ equal to 0.6151 ± 0.0008 and 0.6152 ± 0.0008, and those for l₂ equal to 0.0886 ± 0.0003 and 0.0881 ± 0.0003 for LAGEOS 1 and LAGEOS 2 tracking data are compared to other independently derived estimates. These parameters and their errors achieve stability at about the 24 and 27 month time interval for h₂ and l₂, respectively.     

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.     

Comparison of observed hmF2 and IRI 2007 model with M(3000)F2 estimation of hmF2 at low solar activity for an equatorial station

This is to investigate ways of improving the Equatorial F2-layer peak heights estimated from M(3000)F2 ionosonde data measured using the Ionospheric Prediction Service (IPS-42) sounder at Ouagadougou, Burkina Faso (Latitude +12.4°N, Longitude +1.5°W, Dip latitude +5.9°N) during a low solar activity year (1995). For this purpose, we have compared the observed h_mF2 (h_mF2_obs) deduced using an algorithm from scaled virtual heights of quiet day ionograms and the predicted h_mF2 values which is given by the IRI 2007 model (h_mF2_{IRI 2007}) with the ionosonde measured M(3000)F2 estimation of the h_mF2 values (h_mF2_est) respectively. The correlation coefficients R² for all the seasons were found to range from 0.259 to 0.692 for h_mF2_obs values, while it ranges from 0.551 to 0.875 for the h_mF2_{IRI 2007} values. During the nighttime, estimated h_mF2 (h_mF2_est) was found to be positively correlated with the h_mF2_obs values by the post-sunset peak representation which is also represented by the h_mF2_{IRI 2007} values. We also investigated the validity of the h_mF2_est values by finding the percentage deviations when compared with the h_mF2_obs and h_mF2_{IRI 2007}.     

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

GDP spatialization in Ningbo City based on NPP/VIIRS night-time light and auxiliary data using random forest regression

Accurate spatial distribution information on gross domestic product (GDP) is of great importance for the analysis of economic development, industrial distribution and urbanization processes. Traditional administrative unit-based GDP statistics cannot depict the detailed spatial differences in GDP within each administrative unit. This paper presents a study of GDP spatialization in Ningbo City, China based on National Polar-orbiting Partnership (NPP)/Visible Infrared Imaging Radiometer Suite (VIIRS) night-time light (NTL) data and town-level GDP statistical data. The Landsat image, land cover, road network and topographic data were also employed as auxiliary data to derive independent variables for GDP modelling. Multivariate linear regression (MLR) and random forest (RF) regression were used to estimate GDP at the town scale and were assessed by cross-validation. The results show that the RF model achieved significantly higher accuracy, with a mean absolute error (MAE) of 109.46 million China Yuan (CNY)·km⁻² and a determinate coefficient (R² = 0.77) than the MLR model (MAE = 161.8 million CNY·km⁻², R² = 0.59). Meanwhile, by comparing with the estimated GDP data at the county level, the town-level estimated data showed a better performance in mapping GDP distribution (MAE decreased from 115.1 million CNY·km⁻² to 74.8 million CNY·km⁻²). Among all of the independent variables, NTL, land surface temperature (T_s) and plot ratio (PR) showed higher impacts on the GDP estimation accuracy than the other variables. The GDP density map generated by the RF model depicted the detailed spatial distribution of the economy in Ningbo City. By interpreting the spatial distribution of the GDP, we found that the GDP of Ningbo was high in the northeast and low in the southwest and formed continuous clusters in the north. In addition, the GDP of Ningbo also gradually decreased from the urban centre to its surrounding areas. The produced GDP map provides a good reference for the future urban planning and socio-economic development strategies.     

Equatorial M(3000)F2 estimation of F2-layer models peak heights during low solar activity

M(3000)F2 estimation of h_mF2 based on four different formulated models viz: (1) Shimazaki (1955) (2) Bradley and Dudeney (1973), (3) Dudeney (1974) and (4) Bilitza et al. (1979) at an equatorial station in West Africa during low solar activity period (1995) are used to validate its conformity with observed and International Reference Ionosphere (IRI) model. Local time analyses of data from fifteen (15) selected days during the January and July solstices and April and October equinoxes are used. The results obtained show that the M(3000)F2 estimation of h_mF2 from the ionosonde-measured values using the Ionospheric Prediction Service (IPS-42) sounder compared to the observed values which were deduced using an algorithm from scaled virtual heights of quiet day ionograms are highly correlated with Bilitza model. International Reference Ionosphere (IRI 2007) model for the equatorial region also agrees with the formulation developed by Bilitza et al. (1979) for the four different seasons of the year. h_mF2 is highest (425 km) in summer (June solstice) season and lowest (386 km) in autumn (September equinox) season with daytimes peaks occurring at 1100–1200 LT during the solstices and at 1000 LT during the equinoxes respectively. Also, the post-sunset peaks are highest (362 km) at the spring (March equinox) and lowest (308 km) at the summer (June solstice) both occurring between 1800 and 2000 LT.     

Using S-band dual polarized radar for convective/stratiform rain indexing and the correspondence with AMSR-E GSFC profiling algorithm

The separation of rain types in convective and stratiform regimes has long been a goal in microwave remote sensing of precipitation research. In this essence, a dual polarized radar based indexing scheme that provides information on convective and stratiform (C/S) rain regimes has been presented in correspondence with advanced microwave scanning radiometer – earth observing system (AMSR-E) GSFC profiling algorithm estimate of convective rain percentage. The dual polarized radar based C/S indexing scheme first retrieves the normalized gamma drop size distribution parameters, median volume drop diameter (D₀) and concentration parameter (N_w), from dual polarized radar measurements Z_H and Z_DR, representing reflectivity and differential reflectivity respectively, by means of the genetic programming approach. Next, the C/S rain index is calculated based on the formulation of an empirical relation in N_w–D₀ domain. The scheme has been inspected and applied on measurements from the S-band Chilbolton dual polarized radar. A considerable number of “coincident” cases from the radar and the AMSR-E observations are investigated. It has been revealed that the dual polarized radar based C/S rain indexing is in a similar pattern with the AMSR-E GSFC profiling algorithm estimate of convective rain percentage. Generally, as C/S rain index value increases, which signifies a stratiform to convective trend, the AMSR-E convective rain percentage also increases.     

Possibility to detect the cluster emission up to R200 with XMM-Newton data

We present the results of analysis XMM-Newton data of galaxy cluster CL0016+16, which enables us to trace X-ray emission and temperature profile up to the virial radius. We obtained similar results using three different backgrounds. We checked the possibility of detection of cluster emission up to the virial radius with XMM-Newton data with hydrodynamical cosmology simulation from the Adaptive Mesh Refinement technique, code RAMSES by Teyssier [Teyssier, R. Cosmological hydrodynamics with adaptive mesh refinement. A new high resolution code called RAMSES. A&A 385, 337, 2002], convolution with XMM-Newton and the data base of the spectra by Sauvageot et al. [Sauvageot, J.-L., Belsole, E., Pratt, G.W. The late merging phase of a galaxy cluster: XMM EPIC observations of A 3266. A&A, 444, 673, 2005]. For the first time we were able to compute the mass of CL0016 up to R₂₀₀, we found, assuming hydrostatic equilibrium framework, M₂₀₀ = (1.15 ± 0.11) × 10¹⁵M_⊙.