Soil density evaluated by spectral reflectance as an evidence of compaction effects

Soil compaction, reflected by high bulk density, is an environmental degradation process and new technologies are being developed for its detection. Despite the proven efficiency of remote sensing, it has not been widely used for soil density. Our objective was to evaluate the density of two soils: a Typic Quartzpisament (TQ) and a Rhodic Paleudalf (RP), using spectral reflectance obtained by a laboratory spectroradiometer between 450 and 2500 nm. Undisturbed samples were taken at two depths (0-20 and 60-80 cm), and were artificially compacted. Spectral data, obtained before and after compaction, were compared for both wet and dried compacted samples. Results demonstrated that soil density was greater in RP than in TQ at both depths due to its clayey texture. Spectral data detected high density (compacted) from low density (non-compacted) clayey soils under both wet and dry conditions. The detection of density in sandy soils by spectral reflectance was not possible. The intensity of spectral reflectance of high soil bulk density (compacted) samples was higher than for low density (non-compacted) soils due to changes in soil structure and porosity. Dry samples with high bulk density showed differences in the spectral intensity, but not in the absorption features. Wet samples in equal condition had statistically higher reflectance intensity than that of the low soil bulk density (non-compacted), and absorption differences at 1920 nm, which was due to the altered position of the water molecules. Soil line and spectral reflectance used together could detect soil bulk density variations for the clay soil. This technique could assist in the detection of high soil density in the laboratory by providing new soil information.

Characterization of nicosulfuron availability in aged soils

Sorption-desorption interactions of pesticides with soil determine their availability for transport, plant uptake, and microbial degradation. These interactions are affected by the physical-chemical properties of the pesticide and soil, and for some pesticides, their residence time in the soil. This research evaluated changes in sorption/availability of nicosulfuron (2-[[[[(4,6-dimethoxy-2-pyrimidinyl]aminolcarbonyl]amino]sulfonyl]-N,N-dimethyl-3-pyridinecarboxamide) herbicide with aging in different soils, using a radiolabeled (C-14) tracer. Aging significantly increased sorption. For instance, after the 41-day incubation, calculated K-d,K-app increased by a factor of 2 to 3 in Mollisols from the Midwestern United States and by a factor of 5 to 9 in Oxisols from Brazil and Hawaii, as compared to freshly treated soils. In view of this outcome, potential transport of nicosulfuron would be overpredicted if freshly treated soil Kd values were used to predict transport. The fact that the nicosulfuron solution concentration decreased faster than the soil concentration with time suggested that the increase in sorption was because the rate of degradation in solution and on labile sites was faster than the rate of desorption of the neutral species from the soil particles. It may have also been due to nicosulfuron anion diffusion to less accessible sites with time, leaving the more strongly bound neutral molecules for the sorption characterization. Regardless of the mechanism, these results are further evidence that increases in sorption during pesticide aging should be taken into account during the characterization of the sorption process for mathematical models of pesticide degradation and transport.

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and interannually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R(2) = 0.80 and R(2) = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages. (C) 2010 Elsevier Inc. All rights reserved.

MODIS NDVI time-series allow the monitoring of Eucalyptus plantation biomass

The use of remote sensing is necessary for monitoring forest carbon stocks at large scales. Optical remote sensing, although not the most suitable technique for the direct estimation of stand biomass, offers the advantage of providing large temporal and spatial datasets. In particular, information on canopy structure is encompassed in stand reflectance time series. This study focused on the example of Eucalyptus forest plantations, which have recently attracted much attention as a result of their high expansion rate in many tropical countries. Stand scale time-series of Normalized Difference Vegetation Index (NDVI) were obtained from MODIS satellite data after a procedure involving un-mixing and interpolation, on about 15,000 ha of plantations in southern Brazil. The comparison of the planting date of the current rotation (and therefore the age of the stands) estimated from these time series with real values provided by the company showed that the root mean square error was 35.5 days. Age alone explained more than 82% of stand wood volume variability and 87% of stand dominant height variability. Age variables were combined with other variables derived from the NDVI time series and simple bioclimatic data by means of linear (Stepwise) or nonlinear (Random Forest) regressions. The nonlinear regressions gave r-square values of 0.90 for volume and 0.92 for dominant height, and an accuracy of about 25 m(3)/ha for volume (15% of the volume average value) and about 1.6 m for dominant height (8% of the height average value). The improvement including NDVI and bioclimatic data comes from the fact that the cumulative NDVI since planting date integrates the interannual variability of leaf area index (LAI), light interception by the foliage and growth due for example to variations of seasonal water stress. The accuracy of biomass and height predictions was strongly improved by using the NDVI integrated over the two first years after planting, which are critical for stand establishment. These results open perspectives for cost-effective monitoring of biomass at large scales in intensively-managed plantation forests. (C) 2011 Elsevier Inc. All rights reserved.

Relating MODIS vegetation index time-series with structure, light absorption and stem production of fast-growing Eucalyptus plantations

By allowing the estimation of forest structural and biophysical characteristics at different temporal and spatial scales, remote sensing may contribute to our understanding and monitoring of planted forests. Here, we studied 9-year time-series of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) on a network of 16 stands in fast-growing Eucalyptus plantations in Sao Paulo State, Brazil. We aimed to examine the relationships between NDVI time-series spanning entire rotations and stand structural characteristics (volume, dominant height, mean annual increment) in these simple forest ecosystems. Our second objective was to examine spatial and temporal variations of light use efficiency for wood production, by comparing time-series of Absorbed Photosynthetically Active Radiation (APAR) with inventory data. Relationships were calibrated between the NDVI and the fractions of intercepted diffuse and direct radiation, using hemispherical photographs taken on the studied stands at two seasons. APAR was calculated from the NDVI time-series using these relationships. Stem volume and dominant height were strongly correlated with summed NDVI values between planting date and inventory date. Stand productivity was correlated with mean NDVI values. APAR during the first 2 years of growth was variable between stands and was well correlated with stem wood production (r(2) = 0.78). In contrast, APAR during the following years was less variable and not significantly correlated with stem biomass increments. Production of wood per unit of absorbed light varied with stand age and with site index. In our study, a better site index was accompanied both by increased APAR during the first 2 years of growth and by higher light use efficiency for stem wood production during the whole rotation. Implications for simple process-based modelling are discussed. (C) 2009 Elsevier B.V. All rights reserved.

Processing optimization of probiotic yogurt containing glucose oxidase using response surface methodology

Exposure to oxygen may induce a lack of functionality of probiotic dairy foods because the anaerobic metabolism of probiotic bacteria compromises during storage the maintenance of their viability to provide benefits to consumer health. Glucose oxidase can constitute a potential alternative to increase the survival of probiotic bacteria in yogurt because it consumes the oxygen permeating to the inside of the pot during storage, thus making it possible to avoid the use of chemical additives. This research aimed to optimize the processing of probiotic yogurt supplemented with glucose oxidase using response surface methodology and to determine the levels of glucose and glucose oxidase that minimize the concentration of dissolved oxygen and maximize the Bifidobacterium longum count by the desirability function. Response surface methodology mathematical models adequately described the process, with adjusted determination coefficients of 83% for the oxygen and 94% for the B. longum. Linear and quadratic effects of the glucose oxidase were reported for the oxygen model, whereas for the B. longum count model an influence of the glucose oxidase at the linear level was observed followed by the quadratic influence of glucose and quadratic effect of glucose oxidase. The desirability function indicated that 62.32 ppm of glucose oxidase and 4.35 ppm of glucose was the best combination of these components for optimization of probiotic yogurt processing. An additional validation experiment was performed and results showed acceptable error between the predicted and experimental results.

Mapping changes in coastline geomorphic features using Landsat TM and ETM+ imagery: examples in southeastern Brazil

Multitemporal Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) imagery was used to assess coastline morphological changes in southeastern Brazil. A spectral linear mixing approach (SLMA) was used to estimate fraction imagery representing amounts of vegetation, clean water (a proxy for shade) and soil. Fraction abundances were related to erosive and depositional features. Shoreline, sandy banks (including emerged and submerged banks) and sand spits were highlighted mainly by clean water and soil fraction imagery. To evaluate changes in the coastline geomorphic features, the fraction imagery generated for each data set was classified in a contextual approach using a segmentation technique and ISOSEG, an unsupervised classification. Evaluation of the classifications was performed visually and by an error matrix relating ground-truth data to classification results. Comparison of the classification results revealed an intense transformation in the coastline, and that erosive and depositional features are extremely dynamic and subject to change in short periods of time.

Pharmacodynamics of PEG-IFN-alpha-2a and HCV Response as a Function of IL28B Polymorphism in HIV/HCV-Coinfected Patients

We examined the association between IL28B single-nucleotide polymorphism rs12979860, hepatitis C virus (HCV) kinetic, and pegylated interferon alpha-2a pharmacodynamic parameters in HIV/HCV-coinfected patients from South America. Twenty-six subjects received pegylated interferon alpha-2a + ribavirin. Serum HCV-RNA and interferon concentrations were measured frequently during the first 12 weeks of therapy and analyzed using mathematical models. African Americans and whites had a similar distribution of IL28B genotypes (P = 0.5). The IL28B CC genotype was overrepresented (P = 0.015) in patients infected with HCV genotype-3 compared with genotype-1. In both genotype-1 and genotype-3, the first-phase viral decline and the average pegylated interferon-alpha-2a effectiveness during the first week of therapy were larger (trend P <= 0.12) in genotype-CC compared with genotypes-TC/TT. In genotype-1 patients, the second slower phase of viral decline (days 2-29) and infected cells loss rate, delta, were larger (P = 0.02 and 0.11, respectively) in genotype-CC than in genotypes-TC/TT. These associations were not observed in genotype-3 patients.

THE TEMPO AND MODE OF EVOLUTION OF TRANSPOSABLE ELEMENTS AS REVEALED BY MOLECULAR PHYLOGENIES RECONSTRUCTED FROM MOSQUITO GENOMES

Although many mathematical models exist predicting the dynamics of transposable elements (TEs), there is a lack of available empirical data to validate these models and inherent assumptions. Genomes can provide a snapshot of several TE families in a single organism, and these could have their demographics inferred by coalescent analysis, allowing for the testing of theories on TE amplification dynamics. Using the available genomes of the mosquitoes Aedes aegypti and Anopheles gambiae, we indicate that such an approach is feasible. Our analysis follows four steps: (1) mining the two mosquito genomes currently available in search of TE families; (2) fitting, to selected families found in (1), a phylogeny tree under the general time-reversible (GTR) nucleotide substitution model with an uncorrelated lognormal (UCLN) relaxed clock and a nonparametric demographic model; (3) fitting a nonparametric coalescent model to the tree generated in (2); and (4) fitting parametric models motivated by ecological theories to the curve generated in (3).

A decision support system to facilitate management of patients with acute gastrointestinal bleeding

Objective: To develop a model to predict the bleeding source and identify the cohort amongst patients with acute gastrointestinal bleeding (GIB) who require urgent intervention, including endoscopy. Patients with acute GIB, an unpredictable event, are most commonly evaluated and managed by non-gastroenterologists. Rapid and consistently reliable risk stratification of patients with acute GIB for urgent endoscopy may potentially improve outcomes amongst such patients by targeting scarce health-care resources to those who need it the most. Design and methods: Using ICD-9 codes for acute GIB, 189 patients with acute GIB and all. available data variables required to develop and test models were identified from a hospital medical records database. Data on 122 patients was utilized for development of the model and on 67 patients utilized to perform comparative analysis of the models. Clinical data such as presenting signs and symptoms, demographic data, presence of co-morbidities, laboratory data and corresponding endoscopic diagnosis and outcomes were collected. Clinical data and endoscopic diagnosis collected for each patient was utilized to retrospectively ascertain optimal management for each patient. Clinical presentations and corresponding treatment was utilized as training examples. Eight mathematical models including artificial neural network (ANN), support vector machine (SVM), k-nearest neighbor, linear discriminant analysis (LDA), shrunken centroid (SC), random forest (RF), logistic regression, and boosting were trained and tested. The performance of these models was compared using standard statistical analysis and ROC curves. Results: Overall the random forest model best predicted the source, need for resuscitation, and disposition with accuracies of approximately 80% or higher (accuracy for endoscopy was greater than 75%). The area under ROC curve for RF was greater than 0.85, indicating excellent performance by the random forest model Conclusion: While most mathematical models are effective as a decision support system for evaluation and management of patients with acute GIB, in our testing, the RF model consistently demonstrated the best performance. Amongst patients presenting with acute GIB, mathematical models may facilitate the identification of the source of GIB, need for intervention and allow optimization of care and healthcare resource allocation; these however require further validation. (c) 2007 Elsevier B.V. All rights reserved.

Canine population dynamics: the potential effect of sterilization campaigns

Objective. To analyze, through mathematical modeling, the potential ability of sterilization campaigns to reduce the population density of pet dogs. Methods. Mathematical models were constructed to simulate the canine population dynamics and project the results of control strategies based on several sterilization rates. Results. Even at high sterilization rates (for example, 0.80 year(-1)), it would take approximately 5 years to reduce density by 20%. Even so, other sources of population growth, such as the importing of dogs from other geographic areas, could outweigh the effects of a sterilization program. Conclusions. A program`s effectiveness is contingent upon not only on the sterilization rate, but also the rate of population growth. Sterilization campaigns may potentially reduce population density, but this reduction may not be immediately evident.

Estimation of R(0) from the initial phase of an outbreak of a vector-borne infection

The magnitude of the basic reproduction ratio R(0) of an epidemic can be estimated in several ways, namely, from the final size of the epidemic, from the average age at first infection, or from the initial growth phase of the outbreak. In this paper, we discuss this last method for estimating R(0) for vector-borne infections. Implicit in these models is the assumption that there is an exponential phase of the outbreaks, which implies that in all cases R(0) > 1. We demonstrate that an outbreak is possible, even in cases where R(0) is less than one, provided that the vector-to-human component of R(0) is greater than one and that a certain number of infected vectors are introduced into the affected population. This theory is applied to two real epidemiological dengue situations in the southeastern part of Brazil, one where R(0) is less than one, and other one where R(0) is greater than one. In both cases, the model mirrors the real situations with reasonable accuracy.