Novel methods improve prediction of species' distributions from occurrence data

Prediction of species' distributions is central to diverse applications in ecology, evolution and conservation science. There is increasing electronic access to vast sets of occurrence records in museums and herbaria, yet little effective guidance on how best to use this information in the context of numerous approaches for modelling distributions. To meet this need, we compared 16 modelling methods over 226 species from 6 regions of the world, creating the most comprehensive set of model comparisons to date. We used presence-only data to fit models, and independent presence-absence data to evaluate the predictions. Along with well-established modelling methods such as generalised additive models and GARP and BIOCLIM, we explored methods that either have been developed recently or have rarely been applied to modelling species' distributions. These include machine-learning methods and community models, both of which have features that may make them particularly well suited to noisy or sparse information, as is typical of species' occurrence data. Presence-only data were effective for modelling species' distributions for many species and regions. The novel methods consistently outperformed more established methods. The results of our analysis are promising for the use of data from museums and herbaria, especially as methods suited to the noise inherent in such data improve.

Time-scale and other invariants of integrative mechanical behavior in living cells.

In dealing with systems as complex as the cytoskeleton, we need organizing principles or, short of that, an empirical framework into which these systems fit. We report here unexpected invariants of cytoskeletal behavior that comprise such an empirical framework. We measured elastic and frictional moduli of a variety of cell types over a wide range of time scales and using a variety of biological interventions. In all instances elastic stresses dominated at frequencies below 300 Hz, increased only weakly with frequency, and followed a power law; no characteristic time scale was evident. Frictional stresses paralleled the elastic behavior at frequencies below 10 Hz but approached a Newtonian viscous behavior at higher frequencies. Surprisingly, all data could be collapsed onto master curves, the existence of which implies that elastic and frictional stresses share a common underlying mechanism. Taken together, these findings define an unanticipated integrative framework for studying protein interactions within the complex microenvironment of the cell body, and appear to set limits on what can be predicted about integrated mechanical behavior of the matrix based solely on cytoskeletal constituents considered in isolation. Moreover, these observations are consistent with the hypothesis that the cytoskeleton of the living cell behaves as a soft glassy material, wherein cytoskeletal proteins modulate cell mechanical properties mainly by changing an effective temperature of the cytoskeletal matrix. If so, then the effective temperature becomes an easily quantified determinant of the ability of the cytoskeleton to deform, flow, and reorganize.

Applicability of a distributed watershed pollution model in a data-poor environment in Santa Catarina State, Brazil

Intensification of agricultural production without a sound management and regulations can lead to severe environmental problems, as in Western Santa Catarina State, Brazil, where intensive swine production has caused large accumulations of manure and consequently water pollution. Natural resource scientists are asked by decision-makers for advice on management and regulatory decisions. Distributed environmental models are useful tools, since they can be used to explore consequences of various management practices. However, in many areas of the world, quantitative data for model calibration and validation are lacking. The data-intensive distributed environmental model AgNPS was applied in a data-poor environment, the upper catchment (2,520 ha) of the Ariranhazinho River, near the city of Seara, in Santa Catarina State. Steps included data preparation, cell size selection, sensitivity analysis, model calibration and application to different management scenarios. The model was calibrated based on a best guess for model parameters and on a pragmatic sensitivity analysis. The parameters were adjusted to match model outputs (runoff volume, peak runoff rate and sediment concentration) closely with the sparse observed data. A modelling grid cell resolution of 150 m adduced appropriate and computer-fit results. The rainfall runoff response of the AgNPS model was calibrated using three separate rainfall ranges (< 25, 25-60, > 60 mm). Predicted sediment concentrations were consistently six to ten times higher than observed, probably due to sediment trapping along vegetated channel banks. Predicted N and P concentrations in stream water ranged from just below to well above regulatory norms. Expert knowledge of the area, in addition to experience reported in the literature, was able to compensate in part for limited calibration data. Several scenarios (actual, recommended and excessive manure applications, and point source pollution from swine operations) could be compared by the model, using a relative ranking rather than quantitative predictions.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Perspectives for Forensic Intelligence in anti-doping: Thinking outside of the box.

Today's approach to anti-doping is mostly centered on the judicial process, despite pursuing a further goal in the detection, reduction, solving and/or prevention of doping. Similarly to decision-making in the area of law enforcement feeding on Forensic Intelligence, anti-doping might significantly benefit from a more extensive gathering of knowledge. Forensic Intelligence might bring a broader logical dimension to the interpretation of data on doping activities for a more future-oriented and comprehensive approach instead of the traditional case-based and reactive process. Information coming from a variety of sources related to doping, whether directly or potentially, would feed an organized memory to provide real time intelligence on the size, seriousness and evolution of the phenomenon. Due to the complexity of doping, integrating analytical chemical results and longitudinal monitoring of biomarkers with physiological, epidemiological, sociological or circumstantial information might provide a logical framework enabling fit for purpose decision-making. Therefore, Anti-Doping Intelligence might prove efficient at providing a more proactive response to any potential or emerging doping phenomenon or to address existing problems with innovative actions or/and policies. This approach might prove useful to detect, neutralize, disrupt and/or prevent organized doping or the trafficking of doping agents, as well as helping to refine the targeting of athletes or teams. In addition, such an intelligence-led methodology would serve to address doping offenses in the absence of adverse analytical chemical evidence.

Rainfall erosivity and rainfall return period in the experimental watershed of Aracruz, in the coastal plain of Espirito Santo, Brazil

Knowledge on the factors influencing water erosion is fundamental for the choice of the best land use practices. Rainfall, expressed by rainfall erosivity, is one of the most important factors of water erosion. The objective of this study was to determine rainfall erosivity and the return period of rainfall in the Coastal Plains region, near Aracruz, a town in the state of Espírito Santo, Brazil, based on available data. Rainfall erosivity was calculated based on historic rainfall data, collected from January 1998 to July 2004 at 5 min intervals, by automatic weather stations of the Aracruz Cellulose S.A company. A linear regression with individual rainfall and erosivity data was fit to obtain an equation that allowed data extrapolation to calculate individual erosivity for a 30-year period. Based on this data the annual average rainfall erosivity in Aracruz was 8,536 MJ mm ha-1 h-1 yr-1. Of the total annual rainfall erosivity 85 % was observed in the most critical period October to March. Annual erosive rains accounted for 38 % of the events causing erosion, although the runoff volume represented 88 % of the total. The annual average rainfall erosivity return period was estimated to be 3.4 years.

French version of the Family Attitude Scale: Psychometric properties and relation of attitudes to the respondent's psychiatric status.

The Family Attitude Scale (FAS) is a self-report measure of critical or hostile attitudes and behaviors towards another family member, and demonstrates an ability to predict relapse in psychoses. Data are not currently available on a French version of the scale. The present study developed a French version of the FAS, using a large general population sample to test its internal structure, criterion validity and relationships with the respondents' symptoms and psychiatric diagnoses, and examined the reciprocity of FAS ratings by respondents and their partners. A total of 2072 adults from an urban population undertook a diagnostic interview and completed self-report measures, including an FAS about their partner. A subset of participants had partners who also completed the FAS. Confirmatory factor analyses revealed an excellent fit by a single-factor model, and the FAS demonstrated a strong association with dyadic adjustment. FAS scores of respondents were affected by their anxiety levels and mood, alcohol and anxiety diagnoses, and moderate reciprocity of attitudes and behaviors between the partners was seen. The French version of the FAS has similarly strong psychometric properties to the original English version. Future research should assess the ability of the French FAS to predict relapse of psychiatric disorders.

Pattern recognition ecological niche models fit to presence-only and presence-absence data

1. Identifying the boundary of a species' niche from observational and environmental data is a common problem in ecology and conservation biology and a variety of techniques have been developed or applied to model niches and predict distributions. Here, we examine the performance of some pattern-recognition methods as ecological niche models (ENMs). Particularly, one-class pattern recognition is a flexible and seldom used methodology for modelling ecological niches and distributions from presence-only data. The development of one-class methods that perform comparably to two-class methods (for presence/absence data) would remove modelling decisions about sampling pseudo-absences or background data points when absence points are unavailable. 2. We studied nine methods for one-class classification and seven methods for two-class classification (five common to both), all primarily used in pattern recognition and therefore not common in species distribution and ecological niche modelling, across a set of 106 mountain plant species for which presence-absence data was available. We assessed accuracy using standard metrics and compared trade-offs in omission and commission errors between classification groups as well as effects of prevalence and spatial autocorrelation on accuracy. 3. One-class models fit to presence-only data were comparable to two-class models fit to presence-absence data when performance was evaluated with a measure weighting omission and commission errors equally. One-class models were superior for reducing omission errors (i.e. yielding higher sensitivity), and two-classes models were superior for reducing commission errors (i.e. yielding higher specificity). For these methods, spatial autocorrelation was only influential when prevalence was low. 4. These results differ from previous efforts to evaluate alternative modelling approaches to build ENM and are particularly noteworthy because data are from exhaustively sampled populations minimizing false absence records. Accurate, transferable models of species' ecological niches and distributions are needed to advance ecological research and are crucial for effective environmental planning and conservation; the pattern-recognition approaches studied here show good potential for future modelling studies. This study also provides an introduction to promising methods for ecological modelling inherited from the pattern-recognition discipline.

Chronofunctions of Heilu soil developed from Loess in Luochuan, on the chinese Loess plateau

Soil chronofunctions are an alternative for the quantification of soil-forming processes and underlie the modeling of soil genesis. To establish soil chronofunctions of a Heilu soil profile on Loess in Luochuan, selected soil properties and the 14C ages in the Holocene were studied. Linear, logarithmic, and third-order polynomial functions were selected to fit the relationships between soil properties and ages. The results indicated that third-order polynomial function fit best for the relationships between clay (< 0.002 mm), silt (0.002-0.02 mm), sand (0.02-2 mm) and soil ages, and a trend of an Ah horizon ocurrence in the profile. The logarithmic function indicated mainly variations of soil organic carbon and pH with time (soil age). The variation in CaCO3 content, Mn/Zr, Fe/Zr, K/Zr, Mg/Zr, Ca/Zr, P/Zr, and Na/Zr ratios with soil age were best described by three-order polynomial functions, in which the trend line showed migration of CaCO3 and some elements.

Electrical conductivity of soil irrigated with swine wastewater estimated by time-domain reflectometry¹

Wastewater application to soil is an alternative for fertilization and water reuse. However, particular care must be taken with this practice, since successive wastewater applications can cause soil salinization. Time-domain reflectometry (TDR) allows for the simultaneous and continuous monitoring of both soil water content and apparent electrical conductivity and thus for the indirect measurement of the electrical conductivity of the soil solution. This study aimed to evaluate the suitability of TDR for the indirect determination of the electrical conductivity (ECse) of the saturated soil extract by using an empirical equation for the apparatus TDR Trase 6050X1. Disturbed soil samples saturated with swine wastewater were used, at soil proportions of 0, 0.45, 0.90, 1.80, 2.70, and 3.60 m³ m-3. The probes were equipped with three handmade 0.20 cm long rods. The fit of the empirical model that associated the TDR measured values of electrical conductivity (EC TDR) to ECse was excellent, indicating this approach as suitable for the determination of electrical conductivity of the soil solution.

Theoretical study of elastic electron scattering off stable and exotic nuclei

Results for elastic electron scattering by nuclei, calculated with charge densities of Skyrme forces and covariant effective Lagrangians that accurately describe nuclear ground states, are compared against experiment in stable isotopes. Dirac partial-wave calculations are performed with an adapted version of the ELSEPA package. Motivated by the fact that studies of electron scattering off exotic nuclei are intended in future facilities in the commissioned GSI and RIKEN upgrades, we survey the theoretical predictions from neutron-deficient to neutron-rich isotopes in the tin and calcium isotopic chains. The charge densities of a covariant interaction that describes the low-energy electromagnetic structure of the nucleon within the Lagrangian of the theory are used to this end. The study is restricted to medium- and heavy-mass nuclei because the charge densities are computed in mean-field approach. Because the experimental analysis of scattering data commonly involves parameterized charge densities, as a surrogate exercise for the yet unexplored exotic nuclei, we fit our calculated mean-field densities with Helm model distributions. This procedure turns out to be helpful to study the neutron-number variation of the scattering observables and allows us to identify correlations of potential interest among some of these observables within the isotopic chains.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

Main-Channel Slopes of Selected Streams in Iowa for Estimation of Flood-Frequency Discharges, November 2004

This report describes a statewide study conducted to develop main-channel slope (MCS) curves for 138 selected streams in Iowa with drainage areas greater than 100 square miles. MCS values determined from the curves can be used in regression equations for estimating flood frequency discharges. Multi-variable regression equations previously developed for two of the three hydrologic regions defined for Iowa require the measurement of MCS. Main-channel slope is a difficult measurement to obtain for large streams using 1:24,000-scale topographic maps. The curves developed in this report provide a simplified method for determining MCS values for sites located along large streams in Iowa within hydrologic Regions 2 and 3. The curves were developed using MCS values quantified for 2,058 selected sites along 138 selected streams in Iowa. A geographic information system (GIS) technique and 1:24,000-scale topographic data were used to quantify MCS values for the stream sites. The sites were selected at about 5-mile intervals along the streams. River miles were quantified for each stream site using a GIS program. Data points for river-mile and MCS values were plotted and a best-fit curve was developed for each stream. An adjustment was applied to all 138 curves to compensate for differences in MCS values between manual measurements and GIS quantification. The multi-variable equations for Regions 2 and 3 were developed using manual measurements of MCS. A comparison of manual measurements and GIS quantification of MCS indicates that manual measurements typically produce greater values of MCS compared to GIS quantification. Median differences between manual measurements and GIS quantification of MCS are 14.8 and 17.7 percent for Regions 2 and 3, respectively. Comparisons of percentage differences between flood-frequency discharges calculated using MCS values of manual measurements and GIS quantification indicate that use of GIS values of MCS for Region 3 substantially underestimate flood discharges. Mean and median percentage differences for 2- to 500-year recurrence-interval flood discharges ranged from 5.0 to 5.3 and 4.3 to 4.5 percent, respectively, for Region 2 and ranged from 18.3 to 27.1 and 12.3 to 17.3 percent for Region 3. The MCS curves developed from GIS quantification were adjusted by 14.8 percent for streams located in Region 2 and by 17.7 percent for streams located in Region 3. Comparisons of percentage differences between flood discharges calculated using MCS values of manual measurements and adjusted-GIS quantification for Regions 2 and 3 indicate that the flood-discharge estimates are comparable. For Region 2, mean percentage differences for 2- to 500-year recurrence-interval flood discharges ranged between 0.6 and 0.8 percent and median differences were 0.0 percent. For Region 3, mean and median differences ranged between 5.4 to 8.4 and 0.0 to 0.3 percent, respectively. A list of selected stream sites presented with each curve provides information about the sites including river miles, drainage areas, the location of U.S. Geological Survey stream flowgage stations, and the location of streams Abstract crossing hydro logic region boundaries or the Des Moines Lobe landforms region boundary. Two examples are presented for determining river-mile and MCS values, and two techniques are presented for computing flood-frequency discharges.

Influence of the morphology of the dural sac on surgical decision making in lumbar spinal stenosis

Introduction: Surgical decision making in lumbar spinal stenosis (LSS) takes into account primarily clinical symptoms as well as concordant radiological findings. We hypothesized that a wide variation of operative threshold would be found in particular as far as judgment of severity of radiological stenosis is concerned. Patients and methods: The number of surgeons who would proceed to decompression was studied relative to the perceived severity of radiological stenosis based either on measurements of dural sac cross sectional area (DSCA) or on the recently described morphological grading as seen on axial T2 MRI images. A link to an electronic survey page with a set of ten axial T2 MRI images taken from ten patients with either low back pain or LSS were sent to members of three national or international spine societies. Those 10 images were randomly presented initially and re-shuffled on a second page including this time DSCA measurements in mm2, ranging from 14 to 226 mm2, giving a total of 20 images to appraise. Morphological grades were ranging from grade A to D. Surgeons were asked if they would consider decompression given the radiological appearance of stenosis and that symptoms of neurological claudication were severe in patients who were otherwise fit for surgery. Fisher's exact test was performed following dichotomization of data when appropriate. Results: A total of 142 spine surgeons (113 orthopedic spine surgeons, 29 neurosurgeons) responded from 25 countries. A substantial agreement was observed in operating patients with severe (grade C) or extreme (grade D) stenosis as defined by the morphological grade compared to lesser stenosis (A&B) grades (p<0.0001). Decision to operate was not dependent on number of years in practice, medical density in practicing country or specialty although more neurosurgeons would operate on grade C stenosis (p<0.005). Disclosing the DSCA measurement did not alter the decision to operate. Although 20 surgeons only had prior knowledge of the description of the morphological grading, their responses showed no statistically significant difference with those of the remaining 122 physicians. Conclusions: This study showed that surgeons across borders are less influenced by DSCA in their decision making than by the morphological appearance of the dural sac. Classifying LSS according to morphology rather than surface measurements appears to be consistent with current clinical practice.

Duration-frequency relationships of heavy rainfall in Santa Catarina, Brazil

The purpose of this study was to adjust equations that establish relationships between rainfall events with different duration and data from weather stations in the state of Santa Catarina, Brazil. In this study, the relationships between different duration heavy rainfalls from 13 weather stations of Santa Catarina were analyzed. From series of maximum annual rainfalls, and using the Gumbel-Chow distribution, the maximum rainfall for durations between 5 min and 24 h were estimated considering return periods from 2 to 100 years. The data fit to the Gumbel-Chow model was verified by the Kolmogorov-Smirnov test at 5 % significance. The coefficients of Bell's equation were adjusted to estimate the relationship between rainfall duration t (min) and the return period T (y) in relation to the maximum rainfall with a duration of 1 hour and a 10 year return period. Likewise, the coefficients of Bell's equation were adjusted based on the maximum rainfall with a duration of 1 day and a 10 year return period. The results showed that these relationships are viable to estimate short-duration rainfall events at locations where there are no rainfall records.