Determination of extreme floods using a distributed hydrological model

Hydrological models developed for extreme precipitation of PMP type are difficult to calibrate because of the scarcity of available data for these events. This article presents the process and results of calibration for a distributed hydrological model at fine scale developed for the estimation of probable maximal floods in the case of a PMP. This calibration is done on two Swiss catchments for two events of summer storms. The calculation done is concentrated on the estimation of the parameters of the model, divided in two parts. The first is necessary for the computation of flow speeds while the second is required for the determination of the initial and final infiltration capacities for each terrain type. The results, validated with the Nash equation show a good correlation between the simulated and observed flows. We also apply this model on two Romanian catchments, showing the river network and estimated flow.

The three horsemen of riches: Plague, war and urbanization in early modern Europe

How did Europe escape the "Iron Law of Wages?" We construct a simple Malthusian model withtwo sectors and multiple steady states, and use it to explain why European per capita incomes andurbanization rates increased during the period 1350-1700. Productivity growth can only explain a smallfraction of the rise in output per capita. Population dynamics changes of the birth and death schedules were far more important determinants of steady states. We show how a major shock to population cantrigger a transition to a new steady state with higher per-capita income. The Black Death was such ashock, raising wages substantially. Because of Engel's Law, demand for urban products increased, andurban centers grew in size. European cities were unhealthy, and rising urbanization pushed up aggregatedeath rates. This effect was reinforced by diseases spread through war, financed by higher tax revenues.In addition, rising trade also spread diseases. In this way higher wages themselves reduced populationpressure. We show in a calibration exercise that our model can account for the sustained rise in Europeanurbanization as well as permanently higher per capita incomes in 1700, without technological change.Wars contributed importantly to the "Rise of Europe", even if they had negative short-run effects. We thustrace Europe s precocious rise to economic riches to interactions of the plague shock with the belligerentpolitical environment and the nature of cities.

Parallelizing remote sensing image geometric correction

Remote sensing spatial, spectral, and temporal resolutions of images, acquired over a reasonably sized image extent, result in imagery that can be processed to represent land cover over large areas with an amount of spatial detail that is very attractive for monitoring, management, and scienti c activities. With Moore's Law alive and well, more and more parallelism is introduced into all computing platforms, at all levels of integration and programming to achieve higher performance and energy e ciency. Being the geometric calibration process one of the most time consuming processes when using remote sensing images, the aim of this work is to accelerate this process by taking advantage of new computing architectures and technologies, specially focusing in exploiting computation over shared memory multi-threading hardware. A parallel implementation of the most time consuming process in the remote sensing geometric correction has been implemented using OpenMP directives. This work compares the performance of the original serial binary versus the parallelized implementation, using several multi-threaded modern CPU architectures, discussing about the approach to nd the optimum hardware for a cost-e ective execution.

A wearable system for multi-segment foot kinetics measurement.

This study aims to design a wearable system for kinetics measurement of multi-segment foot joints in long-distance walking and to investigate its suitability for clinical evaluations. The wearable system consisted of inertial sensors (3D gyroscopes and 3D accelerometers) on toes, forefoot, hindfoot, and shank, and a plantar pressure insole. After calibration in a laboratory, 10 healthy elderly subjects and 12 patients with ankle osteoarthritis walked 50m twice wearing this system. Using inverse dynamics, 3D forces, moments, and power were calculated in the joint sections among toes, forefoot, hindfoot, and shank. Compared to those we previously estimated for a one-segment foot model, the sagittal and transverse moments and power in the ankle joint, as measured via multi-segment foot model, showed a normalized RMS difference of less than 11%, 14%, and 13%, respectively, for healthy subjects, and 13%, 15%, and 14%, for patients. Similar to our previous study, the coronal moments were not analyzed. Maxima-minima values of anterior-posterior and vertical force, sagittal moment, and power in shank-hindfoot and hindfoot-forefoot joints were significantly different between patients and healthy subjects. Except for power, the inter-subject repeatability of these parameters was CMC>0.90 for healthy subjects and CMC>0.70 for patients. Repeatability of these parameters was lower for the forefoot-toes joint. The proposed measurement system estimated multi-segment foot joints kinetics with acceptable repeatability but showed difference, compared to those previously estimated for the one-segment foot model. These parameters also could distinguish patients from healthy subjects. Thus, this system is suggested for outcome evaluations of foot treatments.

Stable isotope compositions of mammoth teeth from Niederweningen, Switzerland: Implications for the Late Pleistocene climate, environment, and diet

Oxygen and carbon isotope compositions of well-preserved mammoth teeth from the Middle Wurmian (40-70 ka) peat layer of Niederweningen, the most important mammoth site in Switzerland, were analysed to reconstruct Late Pleistocene palaeoclimatic and palaeoenvironmental conditions. Drinking water (delta(18)O values of approximately -12.3 +/- 0.9 parts per thousand were calculated front oxygen isotope compositions of mammoth tooth enamel apatite using a species-specific calibration for modern elephants. These delta(18)O(H2O) values reflect the mean oxygen isotope composition of the palaeo-precipitation and are similar to those directly measured for fate Pleistocene groundwater from aquifers in northern Switzerland and southern Germany. Using a present-day delta(18)O(H2)o-precipitation-air temperature relation for Switzerland, a mean annual air temperature (MAT) of around 4.3 +/- 2.1 degrees C can be calculated for the Middle Wurmian at this site. This MAT is in good agreement with palaeotemperature estimates on the basis of Middle Wurmian groundwater recharge temperatures and beetle assemblages. Hence, the climatic conditions in this region were around 4 degrees C cooler during the Middle Wurmian interstadial phase, around 45-50ka BP, than they are today. During this period the mammoths from Niederweningen lived in an open tundra-like, C(3) plant-dominated environment as indicated by enamel (delta(13)C values of -11.5 +/- 0.3 parts per thousand and pollen and macroplant fossils found in the embedding peat. The low variability of enamel delta(13)C and delta(18)O values from different mammoth teeth reflects similar environmental conditions and supports a relatively small time frame for the fossil assemblage. (C) 2006 Elsevier Ltd and INQUA. All rights reserved.

Chemical labelling of oyster shells used for time-calibrated high-resolution Mg/Ca ratios: A tool for estimation of past seasonal temperature variations

The geochemical compositions of biogenic carbonates are increasingly used for palaeoenvironmental reconstructions. The skeletal delta O-18 temperature relationship is dependent on water salinity, so many recent studies have focused on the Mg/Ca and Sr/Ca ratios because those ratios in water do not change significantly on short time scales. Thus, those elemental ratios are considered to be good palaeotemperature proxies in many biominerals, although their use remains ambiguous in bivalve shells. Here, we present the high-resolution Mg/Ca ratios of two modern species of juvenile and adult oyster shells, Crassostrea gigas and Ostrea edulis. These specimens were grown in controlled conditions for over one year in two different locations. In situ monthly Mn-marking of the shells has been used for day calibration. The daily Mg/Ca.ratios in the shell have been measured with an electron microprobe. The high frequency Mg/Ca variation of all specimens displays good synchronism with lunar cycles, suggesting that tides strongly influence the incorporation of Mg/Ca into the shells. Highly significant correlation coefficients (0.70<R<0.83, p<0.0001) between the Mg/Ca ratios and the seawater temperature are obtained only for juvenile C. gigas samples, while metabolic control of Mg/Ca incorporation and lower shell growth rates preclude the use of the Mg/Ca ratio in adult shells as a palaeothermometer. Data from three juvenile C. gigas shells from the two study sites are selected to establish a relationship: T = 3.77Mg/Ca + 1.88, where T is in degrees C and Mg/Ca in mmol/mol. (c) 2012 Elsevier B.V. All rights reserved.

Interpolation-free Coregistration and Phase-Correction of Airborne SAR Interferograms

This letter discusses the detection and correction ofresidual motion errors that appear in airborne synthetic apertureradar (SAR) interferograms due to the lack of precision in the navigationsystem. As it is shown, the effect of this lack of precision istwofold: azimuth registration errors and phase azimuth undulations.Up to now, the correction of the former was carried out byestimating the registration error and interpolating, while the latterwas based on the estimation of the phase azimuth undulations tocompensate the phase of the computed interferogram. In this letter,a new correction method is proposed, which avoids the interpolationstep and corrects at the same time the azimuth phase undulations.Additionally, the spectral diversity technique, used to estimateregistration errors, is critically analyzed. Airborne L-bandrepeat-pass interferometric data of the German Aerospace Center(DLR) experimental airborne SAR is used to validate the method

Estimation of percentage body fat in 6- to 13-year-old children by skinfold thickness, body mass index and waist circumference.

We evaluated the accuracy of skinfold thicknesses, BMI and waist circumference for the prediction of percentage body fat (PBF) in a representative sample of 372 Swiss children aged 6-13 years. PBF was measured using dual-energy X-ray absorptiometry. On the basis of a preliminary bootstrap selection of predictors, seven regression models were evaluated. All models included sex, age and pubertal stage plus one of the following predictors: (1) log-transformed triceps skinfold (logTSF); (2) logTSF and waist circumference; (3) log-transformed sum of triceps and subscapular skinfolds (logSF2); (4) log-transformed sum of triceps, biceps, subscapular and supra-iliac skinfolds (logSF4); (5) BMI; (6) waist circumference; (7) BMI and waist circumference. The adjusted determination coefficient (R² adj) and the root mean squared error (RMSE; kg) were calculated for each model. LogSF4 (R² adj 0.85; RMSE 2.35) and logSF2 (R² adj 0.82; RMSE 2.54) were similarly accurate at predicting PBF and superior to logTSF (R² adj 0.75; RMSE 3.02), logTSF combined with waist circumference (R² adj 0.78; RMSE 2.85), BMI (R² adj 0.62; RMSE 3.73), waist circumference (R² adj 0.58; RMSE 3.89), and BMI combined with waist circumference (R² adj 0.63; RMSE 3.66) (P < 0.001 for all values of R² adj). The finding that logSF4 was only modestly superior to logSF2 and that logTSF was better than BMI and waist circumference at predicting PBF has important implications for paediatric epidemiological studies aimed at disentangling the effect of body fat on health outcomes.

1st radiometric dating of a paleontologically dated Bathonian level from Georgia (USSR) - Use of the cathodoluminescence for selection of suitable plagioclases

Jurassic volcanic formations interlayered with (ammonite-bearing) sediments are common in the Caucasus area; this situation is of interest for the numerical calibration of the poorly documented Jurassic portion of the time scale. However, following petrographic study on thin sections no whole-rocks can be considered reliable geochronometers due to subsequent alteration; from about 20 samples, two were selected for plagioclase dating; one (V134) is probably early Kimmeridgian in age; the other (V136) is probably located in the Lower Bathonian stage according to diagnostic ammonites. Cathodoluminescence (CTL) study has shown that sample V136 was similar to usual volcanic feldspars (blue to green colour); however, the lack of CTL of the V134 plagioclase is a character common to diagenetic feldspars; consequently, in spite of a good optical preservation, this geo-chronometer cannot give an age representative of the time of emplacement of the lava flow. We have combined CTL observation with microprobe analysis in order to document the poorly known CTL behaviour of volcanic feldspars; cations Ti4+ and Fe2+ play a major role in the CTL colour of plagioclases and are able to document the growing history of these feldspars ; phenocrysts are initially rich in Fe2+ (core of the crystals, green in colour), then richer in Ti toward the exterior; microcrysts are even richer in Ti (blue to bright blue). We have also observed that natural CTL colour was modified resulting from acid ``cleaning'' of the separated feldspars : the initial blue or green colour tends to change to yellow or violet, respectively, after acid treatment probably due to oxydation of Fe2+ toward Fe3+. X-ray and microprobe analyses both indicated that plagioclases from sample V134 was near the sodic end member (albite) suggesting a diagenetic origin in this andesitic basalt; In contrast, sample V136 contains a calcic plagioclase of common composition for a doleritic basalt. The K-Ar conventional technique was applied as a preliminary tool for radiometric analysis. The Kimmeridgian Na-plagioclase sample gave a ``rejuvenated'' (85 Ma) apparent age which confirms a late genesis for the separated plagioclase phase; this interpretation is based on CTL observation, X-ray analysis, and microprobe analysis ; these techniques are able to distinguish samples which have been submitted to diagenetic alteration from those which have not. An age consistent with the stratigraphic location has been obtained from sample V136. This age of 161 +/- 3 (2-sigma) Ma, is the first one available from a sample palaeontologically located with reasonable precision within the mid Jurassic time.

Determination of plasma levels of citalopram and its demethylated and deaminated metabolites by gas chromatography and gas chromatography-mass spectrometry.

Sensitive and specific methods based on gas chromatography (GC) and gas chromatography-mass spectrometry (GC-MS) for the determination of levels of citalopram, desmethylcitalopram and didesmethylcitalopram in the plasma of patients treated with citalopram are presented, as well as a GC-MS procedure for the assay of the citalopram propionic acid derivative. After addition of a separate internal standard for each drug, liquid-solvent extraction is used to separate the basic compounds from the acid compounds. The demethylated amines are derivatized with trifluoroacetic anhydride, and the acid metabolite with methyl iodide. GC-MS is performed in the electron impact mode, as mass spectrometry by the (positive-ion) chemical ionization mode (methane and ammonia) appeared to be unsuitable. The limits of quantification were 1 ng/ml for citalopram and desmethylcitalopram and 2 ng/ml for the other metabolites. The correlation coefficients for the calibration curves (range 10-500 ng/ml) were &gt; or = 0.999 for all compounds, whether determined by GC or GC-MS.

The stable hydrogen and oxygen isotope variation of water stored in polyethylene terephthalate (PET) bottles

A set of bottled waters from a single natural spring distributed worldwide in polyethylene terephthalate (PET) bottles has been used to examine the effects of storage in plastic polymer material on the isotopic composition (delta(18)O and delta(2)H values) of the water. All samples analyzed were subjected to the same packaging procedure but experienced different conditions of temperature and humidity during storage. Water sorption and the diffusive transfer of water and water vapor through the wall of the PET bottle may cause isotopic exchange between water within the bottle and water vapor in air near the PET-water interface. Changes of about +4 parts per thousand for delta(2)H and +0.7 parts per thousand for delta(18)O have been measured for water after 253 days of storage within the PET bottle. The results of this study clearly indicate the need to use glass bottles for storing water samples for isotopic studies. It is imperative to transfer PET-bottled natural waters to glass bottles for their use as calibration material or potential international working standards. Copyright (C) 2008 John Wiley & Sons, Ltd.

A new maximum likelihood method for luminosity calibrations

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.

CP2 stars as viewed by the uvby Hbeta system

A new statistical parallax method using the Maximum Likelihood principle is presented, allowing the simultaneous determination of a luminosity calibration, kinematic characteristics and spatial distribution of a given sample. This method has been developed for the exploitation of the Hipparcos data and presents several improvements with respect to the previous ones: the effects of the selection of the sample, the observational errors, the galactic rotation and the interstellar absorption are taken into account as an intrinsic part of the formulation (as opposed to external corrections). Furthermore, the method is able to identify and characterize physically distinct groups in inhomogeneous samples, thus avoiding biases due to unidentified components. Moreover, the implementation used by the authors is based on the extensive use of numerical methods, so avoiding the need for simplification of the equations and thus the bias they could introduce. Several examples of application using simulated samples are presented, to be followed by applications to real samples in forthcoming articles.

Hydrologic Simulations of the Maquoketa River Watershed Using SWAT Working Paper 09-WP 49,June 2009

This paper describes the application of the Soil and Water Assessment Tool (SWAT) model to the Maquoketa River watershed, located in northeast Iowa. The inputs to the model were obtained from the Environmental Protection Agency’s geographic information/database system called Better Assessment Science Integrating Point and Nonpoint Sources (BASINS). Climatic data from six weather stations located in and around the watershed, and measured streamflow data from a U.S. Geological Survey gage station at the watershed outlet were used in the sensitivity analysis of SWAT model parameters as well as its calibration and validation for watershed hydrology and streamflow. A sensitivity analysis was performed using an influence coefficient method to evaluate surface runoff and base flow variations in response to changes in model input hydrologic parameters. The curve number, evaporation compensation factor, and soil available water capacity were found to be the most sensitive parameters among eight selected parameters when applying SWAT to the Maquoketa River watershed. Model calibration, facilitated by the sensitivity analysis, was performed for the period 1988 through 1993, and validation was performed for 1982 through 1987. The model performance was evaluated by well-established statistical methods and was found to explain at least 86% and 69% of the variability in the measured stream flow data for the calibration and validation periods, respectively. This initial hydrologic modeling analysis will facilitate future applications of SWAT to the Maquoketa River watershed for various watershed analysis, including water quality.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.