Spatial and temporal patterns of dry spell lengths in the Iberian Peninsula for the second half of the twentieth century

Summary Daily pluviometric records of 43 meteorological stations across the Iberian Peninsula have permitted a detailed analysis of dry spell patterns for the period 1951–2000 by distinguishing daily amount thresholds of 0.1, 1.0, 5.0 and 10.0 mm/day. The analyses are based on three annual series, namely the number of dry spells, N, the average dry spell length, L, and the extreme dry spell length, L _max. First, the statistical significance of local trends for the annual series of N, L and L _max has been investigated by means of the Mann-Kendall test and significant field trends have been established by means of Monte Carlo simulations. Clear signs of negative field trends are detected for N (1.0 and 10.0 mm/day) and L (0.1 mm/day). Second, the Weibull model fits well the empirical distributions of dry spell lengths for all the rain gauges, whatever the daily amount threshold, with a well ranged spatial distribution of their parameters u and k. On the basis of the Weibull distribution, return period maps for 2, 5, 10, 25 and 50 years have been obtained for dry spell lengths with respect to the four daily threshold levels. While for 0.1 and 1.0 mm/day the longest dry spells are expected at the south of the Iberian Peninsula, for 5.0 and 10.0 mm/day they are mostly detected at the southeast. Finally, the elapsed time between consecutive dry spells has been analysed by considering the same rain amount thresholds and different dry spell lengths at increasing intervals of 10 days. This analysis makes evident a significant negative field trend of the elapsed time between consecutive dry spells of lengths ranging from 10 to 20 days for daily amount thresholds of 1.0, 5.0 and 10.0 mm/day. Authors’ addresses: X. Lana, C. Serra, Departament de Física i Enginyeria Nuclear, ETSEIB, Universitat Politècnica de Catalunya, Av. Diagonal 647 planta 11, 08028 Barcelona, Spain; M. D. Marínez, Departament de Física Aplicada, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain; A. Burgue?o, Departament de Meteorologia i Astronomia, Universitat de Barcelona, 08028 Barcelona, Spain; J. Martín-Vide, L. Gómez, Grup de Climatologia, Universitat de Barcelona, 08028 Barcelona, Spain.     

Comparison of Vegetation Indices from Two Ground Based Sensors

The best and commonly used ground-based sensor to monitor crop growth, ASD FieldSpecPro Spectroradiometer (Analytical Spectral Devices, Boulder, CO, USA) is a passive sensor, which can be used under adequate light condition. However, now-a-days active sensors such as GreenSeeker? (GS) handheld crop response (Trimble Agriculture division, USA) are used for monitoring crop growth and are flexible in terms of timeliness and illumination conditions besides being cheaper than the ASD. Before its wide use, the suitability and accuracy of GS should be assessed by comparing the NDVI measured by this instrument with that by ASD, under diverse wheat growing conditions of India. Keeping this in view, the present experiment was undertaken with the following objectives: (1) to find out the temporal variation of NDVI measured both by ASD and GS treatments, (2) to find out relationship between the NDVI measured through ASD and GS and, (3) to evaluate the suitability of GS for NDVI measurements. It was observed that the numerical value of NDVI as measured by GS was always significantly (P < 0.05) lower than that measured by ASD for all the experiments under study. The NDVI-ASD and NDVI-GS were significantly positively correlated (P < 0.01) with the correlation coefficients being +0.94, +0.88 and +0.87 for irrigation and nitrogen experiment, irrigation and cultivars experiment, and tillage, residue and nitrogen experiments, respectively. Further, the regression equation developed between the NDVI-ASD and NDVI-GS: [NDVI-GS = 1.070 × (NDVI-ASD ? 0.292] can be successfully used to compute the NDVI of ASD from that computed by GS.     

Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

Statistically optimal estimation of Greenland Ice Sheet mass variations from GRACE monthly solutions using an improved mascon approach

We present an improved mascon approach to transform monthly spherical harmonic solutions based on GRACE satellite data into mass anomaly estimates in Greenland. The GRACE-based spherical harmonic coefficients are used to synthesize gravity anomalies at satellite altitude, which are then inverted into mass anomalies per mascon. The limited spectral content of the gravity anomalies is properly accounted for by applying a low-pass filter as part of the inversion procedure to make the functional model spectrally consistent with the data. The full error covariance matrices of the monthly GRACE solutions are properly propagated using the law of covariance propagation. Using numerical experiments, we demonstrate the importance of a proper data weighting and of the spectral consistency between functional model and data. The developed methodology is applied to process real GRACE level-2 data (CSR RL05). The obtained mass anomaly estimates are integrated over five drainage systems, as well as over entire Greenland. We find that the statistically optimal data weighting reduces random noise by 35–69%, depending on the drainage system. The obtained mass anomaly time-series are de-trended to eliminate the contribution of ice discharge and are compared with de-trended surface mass balance (SMB) time-series computed with the Regional Atmospheric Climate Model (RACMO 2.3). We show that when using a statistically optimal data weighting in GRACE data processing, the discrepancies between GRACE-based estimates of SMB and modelled SMB are reduced by 24–47%.     

Estimation of Partial Pressure of Carbon Dioxide and Air-Sea Fluxes in Hooghly Estuary Based on In Situ and Satellite Observations

An empirical model is developed and used with remotely sensed predictors: sea surface temperature (SST) and chlorophyll-a concentration (Chl-a), to compute surface water partial pressure of carbon dioxide (pCO_2w) and air-sea fluxes of CO₂ in the Hooghly estuary and its adjacent coastal oceans. In situ observations used here were based on measurements carried out in this region during winter and summer periods in 2008. The estimated pCO_2w compares well with the in situ observations at root mean square error ±18 μatm. In winter, estimated pCO_2w ranges between 320 and 500 μatm with large values (>400 μatm) on the south-western and south-eastern flanks of the coastal domain and lower values (340–375 μatm) on the main-channel. In summer, it remained spatially uniform at 450 μatm. Extrapolation of the results over the study region based on the Moderate Imaging Specroradiometer (MODIS) measured SST and Chl-a suggests that the region is a strong source of atmospheric CO₂ during the summer with net release of 0.095 Tg C year^?1 (equivalent to mean flux of 90 molC m^?2 year^?1) and is a weak source during the winter with net release of 0.006 Tg C yr^?1 (0.5 molC m^?2 year^?1) from the geographical extent of 6000 Km² area.     

Assessment of Kalpana-1 Rainfall Product over Indian Meteorological Sub-Divisions During the Summer Monsoon Season

In this paper, Kalpana-1 derived INSAT Multispectral Rainfall Algorithm (IMSRA) rainfall estimates are compared with two multisatellite rainfall products namely, TRMM Multisatellite Precipitation Analysis (TMPA)-3B42 and Global Satellite Mapping of Precipitation (GSMaP), and India Meteorological Department (IMD) surface rain gauge (SRG)-based rainfall at meteorological sub-divisional scale over India. The performance of the summer monsoon rainfall of 2013 over Indian meteorological sub-divisions is assessed at different temporal scales. Comparison of daily accumulated rainfall over India from IMSRA shows a linear correlation of 0.72 with TMPA-3B42 and 0.70 with GSMaP estimates. IMSRA is capable to pick up daily rainfall variability over the monsoon trough region as compared to TMPA-3B42 and GSMaP products, but underestimates moderate to heavy rainfall events. Satellite-derived rainfall maps at meteorological sub-divisional scales are in reasonably good agreement with IMD-SRG based rainfall maps with some exceptions. However, IMSRA performs better than GSMaP product at meteorological sub-divisional scale and comparable with TMPA data. All the satellite-derived rainfall products underestimate orographic rainfall along the west coast, the Himalayan foothills and over the northeast India and overestimate rainfall over the southeast peninsular India. Overall results suggest that IMSRA estimates have potential for monsoon rainfall monitoring over the Indian meteorological sub-divisions and can be used for various hydro-meteorological applications.