Multifractal analysis as a tool for validating a rainfall model

A multifractal analysis was carried out in order to validate the simulation of hourly rainfall records of a local climate model for the Iberian Peninsula. Observed and simulated hourly rainfall data from four locations in Andalusia (southern Spain) were used to carry out the study. In order to detect the influence of the length of the data series on the results, two different sizes were used for the real data: 4 years, and 20 years. The results show that algebraic tails are required to fit the probability distribution of extreme rain event sizes, and rain and dry event durations for both kinds of rainfall data. Similar results are found for the extreme rain event sizes and dry event durations fits when the real and synthetic data are considered. Nevertheless, some differences appear in the cases of rain event durations. The detection of the presence of a first‐order multifractal phase transition associated with a critical moment in the empirical moment scaling exponent function and the results of the extreme rain event sizes fits, reveal that real rainfall is a self‐organized criticality (SOC) process. That behaviour is less evident in the simulated rainfall series. The same ‘synoptic maximum’ value was found for each place with both types of rainfall data. A time clustering analysis was carried out applying the count‐based periodogram and the Fano factor methods. Some periodicities have been detected in the periodograms, especially for the longest real rainfall data series. Copyright © 2007 John Wiley & Sons, Ltd.     

A Late‐glacial chironomid record from Hawes Water,northwest England

This paper presents the results of a high‐resolution Late‐glacial chironomid stratigraphy from Hawes Water, a small carbonate lake in northern Lancashire. The samples were from a core taken from the terrestrialised margin of the present lake, which represents an intermediate depth between the true littoral and the profundal. The chironomid assemblage showed a high degree of sensitivity to both broad‐scale and short‐term temperature changes. Comparison with an existing proxy temperature record (δ¹⁸O) for the site confirmed the presence of four temperature inversions within the Late‐glacial Interstadial. A mean July air temperature inference model, derived from acid, soft‐water lakes in Norway and Svalbard, was applied to the data. Despite the absence of carbonate lakes within the Norwegian training set, there was a close similarity between trends in estimated July air temperature and the δ¹⁸O trace, with a particularly strong correspondence in the periods of clay deposition. This suggests that this model is highly robust. The inferred maximum Interstadial temperature was 13.4°C, dropping initially to 7.5°C in the Loch Lomond Stadial. Temperatures reach a maximum of nearly 10°C in this period, cool for a short period before rising rapidly to 13.2°C at the start of the Holocene. These temperatures are similar to but slightly higher than those estimated for Whitrig Bog, southeast Scotland, and lower than those inferred from coleopteran‐based models for sites in South Wales. Copyright © 2004 John Wiley & Sons, Ltd.     

A real-time software backend for the GMRT

The new era of software signal processing has a large impact on radio astronomy instrumentation. Our design and implementation of a 32 antennae, 33 MHz, dual polarization, fully real-time software backend for the GMRT, using only off-the-shelf components, is an example of this. We have built a correlator and a beamformer, using PCI-based ADC cards and a Linux cluster of 48 nodes with dual gigabit inter-node connectivity for real-time data transfer requirements. The highly optimized compute pipeline uses cache efficient, multi-threaded parallel code, with the aid of vectorized processing. This backend allows flexibility in final time and frequency resolutions, and the ability to implement algorithms for radio frequency interference rejection. Our approach has allowed relatively rapid development of a fairly sophisticated and flexible backend receiver system for the GMRT, which will greatly enhance the productivity of the telescope. In this paper we describe some of the first lights using this software processing pipeline. We believe this is the first instance of such a real-time observatory backend for an intermediate sized array like the GMRT.     

In situ U–Pb SIMS (IN-SIMS) micro-baddeleyite dating of mafic rocks: Method with examples

An in situ U–Pb SIMS (IN-SIMS) method to date micro-baddeleyite crystals as small as 3 μm is presented with results from three samples that span a variety of ages and geologic settings. The method complements ID-TIMS geochronology by extending the range of dateable crystals to sizes smaller than can be recovered by physical separation. X-ray mapping and BSE imaging are used to locate target grains in thin section, followed by SIMS analysis on a CAMECA ims 1270, using the field aperture in the transfer column to screen out ions from host phases. Internal age precisions for the method are anticipated to range from 0.1% for Precambrian rocks to 3–7% for Phanerozoic rocks. Results establish a 2689 ± 5 Ma age for mafic dikes in the Wyoming craton, USA, a 1540 ± 30 Ma age for a subaerial lava flow from the Thelon Basin of northern Canada, and a 457 ± 34 Ma age for mafic dikes in the platform sequence of southeastern Siberia. The method is ideal for relatively non-destructive dating of small samples such as extraterrestrial rocks and precious terrestrial samples.     

Groundwater recharge in irrigated semi-arid areas: quantitative hydrological modelling and sensitivity analysis

For semi-arid regions, methods of assessing aquifer recharge usually consider the potential evapotranspiration. Actual evapotranspiration rates can be below potential rates for long periods of time, even in irrigated systems. Accurate estimations of aquifer recharge in semi-arid areas under irrigated agriculture are essential for sustainable water-resources management. A method to estimate aquifer recharge from irrigated farmland has been tested. The water-balance-modelling approach was based on VisualBALAN v. 2.0, a computer code that simulates water balance in the soil, vadose zone and aquifer. The study was carried out in the Campo de Cartagena (SE Spain) in the period 1999–2008 for three different groups of crops: annual row crops (lettuce and melon), perennial vegetables (artichoke) and fruit trees (citrus). Computed mean-annual-recharge values (from irrigation+precipitation) during the study period were 397 mm for annual row crops, 201 mm for perennial vegetables and 194 mm for fruit trees: 31.4, 20.7 and 20.5% of the total applied water, respectively. The effects of rainfall events on the final recharge were clearly observed, due to the continuously high water content in soil which facilitated the infiltration process. A sensitivity analysis to assess the reliability and uncertainty of recharge estimations was carried out.     

Modelling shallow landslide susceptibility: a new approach in logistic regression by using favourability assessment

A new method for estimating shallow landslide susceptibility by combining Geographical Information System (GIS), nonparametric kernel density estimation and logistic regression is described. Specifically, a logistic regression is applied to predict the spatial distribution by estimating the probability of occurrence of a landslide in a 16 km² area. For this purpose, a GIS is employed to gather the relevant sample information connected with the landslides. The advantages of pre-processing the explanatory variables by nonparametric density estimation (for continuous variables) and a reclassification (for categorical/discrete ones) are discussed. The pre-processing leads to new explanatory variables, namely, some functions which measure the favourability of occurrence of a landslide. The resulting model correctly classifies 98.55% of the inventaried landslides and 89.80% of the landscape surface without instabilities. New data about recent shallow landslides were collected in order to validate the model, and 92.20% of them are also correctly classified. The results support the methodology and the extrapolation of the model to the whole study area (278 km²) in order to obtain susceptibility maps.     

Impacts of uncertain river flow data on rainfall‐runoff model calibration and discharge predictions

In order to quantify total error affecting hydrological models and predictions, we must explicitly recognize errors in input data, model structure, model parameters and validation data. This paper tackles the last of these: errors in discharge measurements used to calibrate a rainfall‐runoff model, caused by stage–discharge rating‐curve uncertainty. This uncertainty may be due to several combined sources, including errors in stage and velocity measurements during individual gaugings, assumptions regarding a particular form of stage–discharge relationship, extrapolation of the stage–discharge relationship beyond the maximum gauging, and cross‐section change due to vegetation growth and/or bed movement. A methodology is presented to systematically assess and quantify the uncertainty in discharge measurements due to all of these sources. For a given stage measurement, a complete PDF of true discharge is estimated. Consequently, new model calibration techniques can be introduced to explicitly account for the discharge error distribution. The method is demonstrated for a gravel‐bed river in New Zealand, where all the above uncertainty sources can be identified, including significant uncertainty in cross‐section form due to scour and re‐deposition of sediment. Results show that rigorous consideration of uncertainty in flow data results in significant improvement of the model's ability to predict the observed flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of Canadian daily ice charts with surface observations off Newfoundland,winter 1992

Abstract

Forecast ice drift rates and thicknesses displayed on daily ice charts and forecast winds for the Canadian east coast are compared to on‐ice observations made during the second Canadian Atlantic Storm Program (CASP II) of March 1992. Observed and 24‐hour forecasts of daily ice drift rates were weakly correlated even though long‐term means closely matched observations. Daily drift rates have an RMS error of 13 cm s^‐1 relative to a 15 cm s^‐1 mean in addition to an RMS direction error of 50 degrees. Contributions towards daily drift uncertainties were: the estimation of winds, unmodelled physics of ocean and ice cover processes; and the inconsistency in the methods used by the ice forecaster. Correlation coefficients between forecast winds and on‐ice observed winds decreased from 0.8 at 0‐hour to 0.7 for the 30‐hour forecast. Similar results were found between ice drift rates from forecast winds. Histograms of ice thicknesses observed along narrow swaths using a helicopter‐towed electromagnetic sensor compared well with undeformed ice thicknesses representing large areas on ice charts, with differences mainly caused by difference in ice type representation and by co‐registering the two data sets.