Comparison of time-lapse GPR data collected under natural and forced infiltration conditions to estimate vadose zone hydraulic parameters

Time-lapse crosshole ground-penetrating radar (GPR) data, collected while infiltration occurs, can provide valuable information regarding the hydraulic properties of the unsaturated zone. In particular, the stochastic inversion of such data provides estimates of parameter uncertainties, which are necessary for hydrological prediction and decision making. Here, we investigate the effect of different infiltration conditions on the stochastic inversion of time-lapse, zero-offset-profile, GPR data. Inversions are performed using a Bayesian Markov-chain-Monte-Carlo methodology. Our results clearly indicate that considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Official Positions for FRAX(®) clinical regarding glucocorticoids: the impact of the use of glucocorticoids on the estimate by FRAX(®) of the 10 year risk of fracture from Joint Official Positions Development Conference of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX(®).

Given the significant impact the use of glucocorticoids can have on fracture risk independent of bone density, their use has been incorporated as one of the clinical risk factors for calculating the 10-year fracture risk in the World Health Organization's Fracture Risk Assessment Tool (FRAX(®)). Like the other clinical risk factors, the use of glucocorticoids is included as a dichotomous variable with use of steroids defined as past or present exposure of 3 months or more of use of a daily dose of 5 mg or more of prednisolone or equivalent. The purpose of this report is to give clinicians guidance on adjustments which should be made to the 10-year risk based on the dose, duration of use and mode of delivery of glucocorticoids preparations. A subcommittee of the International Society for Clinical Densitometry and International Osteoporosis Foundation joint Position Development Conference presented its findings to an expert panel and the following recommendations were selected. 1) There is a dose relationship between glucocorticoid use of greater than 3 months and fracture risk. The average dose exposure captured within FRAX(®) is likely to be a prednisone dose of 2.5-7.5 mg/day or its equivalent. Fracture probability is under-estimated when prednisone dose is greater than 7.5 mg/day and is over-estimated when the prednisone dose is less than 2.5 mg/day. 2) Frequent intermittent use of higher doses of glucocorticoids increases fracture risk. Because of the variability in dose and dosing schedule, quantification of this risk is not possible. 3) High dose inhaled glucocorticoids may be a risk factor for fracture. FRAX(®) may underestimate fracture probability in users of high dose inhaled glucocorticoids. 4) Appropriate glucocorticoid replacement in individuals with adrenal insufficiency has not been found to increase fracture risk. In such patients, use of glucocorticoids should not be included in FRAX(®) calculations.

Plasma and urine profiles of Delta9-tetrahydrocannabinol and its metabolites 11-hydroxy-Delta9-tetrahydrocannabinol and 11-nor-9-carboxy-Delta9-tetrahydrocannabinol after cannabis smoking by male volunteers to estimate recent consumption by athletes.

Since 2004, cannabis has been prohibited by the World Anti-Doping Agency for all sports competitions. In the years since then, about half of all positive doping cases in Switzerland have been related to cannabis consumption. In doping urine analysis, the target analyte is 11-nor-9-carboxy-Delta(9)-tetrahydrocannabinol (THC-COOH), the cutoff being 15 ng/mL. However, the wide urinary detection window of the long-term metabolite of Delta(9)-tetrahydrocannabinol (THC) does not allow a conclusion to be drawn regarding the time of consumption or the impact on the physical performance. The purpose of the present study on light cannabis smokers was to evaluate target analytes with shorter urinary excretion times. Twelve male volunteers smoked a cannabis cigarette standardized to 70 mg THC per cigarette. Plasma and urine were collected up to 8 h and 11 days, respectively. Total THC, 11-hydroxy-Delta(9)-tetrahydrocannabinol (THC-OH), and THC-COOH were determined after hydrolysis followed by solid-phase extraction and gas chromatography/mass spectrometry. The limits of quantitation were 0.1-1.0 ng/mL. Eight puffs delivered a mean THC dose of 45 mg. Plasma levels of total THC, THC-OH, and THC-COOH were measured in the ranges 0.2-59.1, 0.1-3.9, and 0.4-16.4 ng/mL, respectively. Peak concentrations were observed at 5, 5-20, and 20-180 min. Urine levels were measured in the ranges 0.1-1.3, 0.1-14.4, and 0.5-38.2 ng/mL, peaking at 2, 2, and 6-24 h, respectively. The times of the last detectable levels were 2-8, 6-96, and 48-120 h. Besides high to very high THC-COOH levels (245 +/- 1,111 ng/mL), THC (3 +/- 8 ng/mL) and THC-OH (51 +/- 246 ng/mL) were found in 65 and 98% of cannabis-positive athletes' urine samples, respectively. In conclusion, in addition to THC-COOH, the pharmacologically active THC and THC-OH should be used as target analytes for doping urine analysis. In the case of light cannabis use, this may allow the estimation of more recent consumption, probably influencing performance during competitions. However, it is not possible to discriminate the intention of cannabis use, i.e., for recreational or doping purposes. Additionally, pharmacokinetic data of female volunteers are needed to interpret cannabis-positive doping cases of female athletes.

Use of variability in national and regional data to estimate the prevalence of lymphangioleiomyomatosis.

Background: Understanding the true prevalence of lymphangioleiomyomatosis (LAM) is important in estimating disease burden and targeting specific interventions. As with all rare diseases, obtaining reliable epidemiological data is difficult and requires innovative approaches.Aim: To determine the prevalence and incidence of LAM using data from patient organizations in seven countries, and to use the extent to which the prevalence of LAM varies regionally and nationally to determine whether prevalence estimates are related to health-care provision.Methods: Numbers of women with LAM were obtained from patient groups and national databases from seven countries (n = 1001). Prevalence was calculated for regions within countries using female population figures from census data. Incidence estimates were calculated for the USA, UK and Switzerland. Regional variation in prevalence and changes in incidence over time were analysed using Poisson regression and linear regression.Results: Prevalence of LAM in the seven countries ranged from 3.4 to 7.8/million women with significant variation, both between countries and between states in the USA. This variation did not relate to the number of pulmonary specialists in the region nor the percentage of population with health insurance, but suggests a large number of patients remain undiagnosed. The incidence of LAM from 2004 to 2008 ranged from 0.23 to 0.31/million women/per year in the USA, UK and Switzerland.Conclusions: Using this method, we have found that the prevalence of LAM is higher than that previously recorded and that many patients with LAM are undiagnosed.

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Changes in volume, surface estimate, three-dimensional shape and total number of neurons of the human primary visual cortex from midgestation until old age.

Macroscopic features such as volume, surface estimate, thickness and caudorostral length of the human primary visual cortex (Brodman's area 17) of 46 human brains between midgestation and 93 years were studied by means of camera lucida drawings from serial frontal sections. Individual values were best fitted by a logistic function from midgestation to adulthood and by a regression line between adulthood and old age. Allometric functions were calculated to study developmental relationships between all the features. The three-dimensional shape of area 17 was also reconstructed from the serial sections in 15 cases and correlated with the sequence of morphological events. The sulcal pattern of area 17 begins to develop around 21 weeks of gestation but remains rather simple until birth, while it becomes more convoluted, particularly in the caudal part, during the postnatal period. Until birth, a large increase in cortical thickness (about 83% of its mean adult value) and caudorostral length (69%) produces a moderate increase in cortical volume (31%) and surface estimate (40%) of area 17. After birth, the cortical volume and surface undergo their maximum growth rate, in spite of a rather small increase in cortical thickness and caudorostral length. This is due to the development of the pattern of gyrification within and around the calcarine fissure. All macroscopic features have reached the mean adult value by the end of the first postnatal year. With aging, the only features to undergo significant regression are the cortical surface estimate and the caudorostral length. The total number of neurons in area 17 shows great interindividual variability at all ages. No decrease in the postnatal period or in aging could be demonstrated.

A new method to estimate the infilling of alluvial sediment of glacial valleys using a sloping local base level

A new method is used to estimate the volumes of sediments of glacial valleys. This method is based on the concept of sloping local base level and requires only a digital terrain model and the limits of the alluvial valleys as input data. The bedrock surface of the glacial valley is estimated by a progressive excavation of the digital elevation model (DEM) of the filled valley area. This is performed using an iterative routine that replaces the altitude of a point of the DEM by the mean value of its neighbors minus a fixed value. The result is a curved surface, quadratic in 2D. The bedrock surface of the Rhone Valley in Switzerland was estimated by this method using the free digital terrain model Shuttle Radar Topography Mission (SRTM) (~92 m resolution). The results obtained are in good agreement with the previous estimations based on seismic profiles and gravimetric modeling, with the exceptions of some particular locations. The results from the present method and those from the seismic interpretation are slightly different from the results of the gravimetric data. This discrepancy may result from the presence of large buried landslides in the bottom of the Rhone Valley.

A new technique with high reproducibility to estimate renal oxygenation using BOLD-MRI in chronic kidney disease.

OBJECTIVES: To assess inter-observer variability of renal blood oxygenation level-dependent MRI (BOLD-MRI) using a new method of analysis, called the concentric objects (CO) technique, in comparison with the classical ROI (region of interest)-based technique. METHODS: MR imaging (3T) was performed before and after furosemide in 10 chronic kidney disease (CKD) patients (mean eGFR 43±24ml/min/1.73m(2)) and 10 healthy volunteers (eGFR 101±28ml/min1.73m(2)), and R2* maps were determined on four coronal slices. In the CO-technique, R2* values were based on a semi-automatic procedure that divided each kidney in six equal layers, whereas in the ROI-technique, all circles (ROIs) were placed manually in the cortex and medulla. The mean R2*values as assessed by two independent investigators were compared. RESULTS: With the CO-technique, inter-observer variability was 0.7%-1.9% across all layers in non-CKD, versus 1.6%-3.8% in CKD. With the ROI-technique, median variability for cortical and medullary R2* values was 3.6 and 6.8% in non-CKD, versus 4.7 and 12.5% in CKD; similar results were observed after furosemide. CONCLUSION: The CO-technique offers a new, investigator-independent, highly reproducible alternative to the ROI-based technique to estimate renal tissue oxygenation in CKD.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.