Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Spatial and temporal analysis of land cover and landscape structure change in Zagros Forests

The Zagros oak forests in Western Iran are critically important to the sustainability of the region. These forests have undergone dramatic declines in recent decades. We evaluated the utility of the non-parametric Random Forest classification algorithm for land cover classification of Zagros landscapes, and selected the best spatial and spectral predictive variables. The algorithm resulted in high overall classification accuracies (>85%) and also equivalent classification accuracies for the datasets from the three different sensors. We evaluated the associations between trends in forest area and structure with trends in socioeconomic and climatic conditions, to identify the most likely driving forces creating deforestation and landscape structure change. We used available socioeconomic (urban and rural population, and rural income), and climatic (mean annual rainfall and mean annual temperature) data for two provinces in northern Zagros. The most correlated driving force of forest area loss was urban population, and climatic variables to a lesser extent. Landscape structure changes were more closely associated with rural population. We examined the effects of scale changes on the results from spatial pattern analysis. We assessed the impacts of eight years of protection in a protected area in northern Zagros at two different scales (both grain and extent). The effects of protection on the amount and structure of forests was scale dependent. We evaluated the nature and magnitude of changes in forest area and structure over the entire Zagros region from 1972 to 2009. We divided the Zagros region in 167 Landscape Units and developed two measures— Deforestation Sensitivity (DS) and Connectivity Sensitivity (CS) — for each landscape unit as the percent of the time steps that forest area and ECA experienced a decrease of greater than 10% in either measure. A considerable loss in forest area and connectivity was detected, but no sudden (nonlinear) changes were detected at the spatial and temporal scale of the study. Connectivity loss occurred more rapidly than forest loss due to the loss of connecting patches. More connectivity was lost in southern Zagros due to climatic differences and different forms of traditional land use.

Spatial and temporal patterns in ungulate-ecosystem interactions

Ungulates are important components of a variety of ecosystems worldwide. This dissertation integrates aspects of ungulate and forest ecology to increase our understanding of how they work together in ways that are of interest to natural resource managers, educators, and those who are simply curious about nature. Although animal ecology and ecosystem ecology are often studied separately, one of the general goals of this dissertation is to examine how they interact across spatial and temporal scales. Forest ecosystems are heterogeneous across a range of scales. Spatial and temporal habitat use patterns of forest ungulates tend to be congregated in patches where food and/or cover are readily available. Ungulates interact with ecosystem processes by selectively foraging on plants and excreting waste products in concentrated patches. Positive feedbacks may develop where these activities increase the value of habitat through soil fertilization or the alteration of plant chemistry and architecture. Heterogeneity in ecosystem processes and plant community structure, observed at both stand and local scales, may be the integrated outcome of feedbacks between ungulate behavior and abiotic resource gradients. The first chapter of this dissertation briefly discusses pertinent background information on ungulate ecology, with a focus on white-tailed deer (Odocoileus virginianus) in the Upper Great Lakes region and moose (Alces acles) in Isle Royale National Park, Michigan, USA. The second chapter demonstrates why ecological context is important for studying ungulate ecology in forest ecosystems. Excluding deer from eastern hemlock (Tsuga canadensis) stands, which deer use primarily as winter cover, resulted in less spatial complexity in soil reactive nitrogen and greater complexity in diffuse light compared to unfenced stands. The spatial patterning of herbaceous-layer cover was more similar to nitrogen where deer were present, and was a combination of nitrogen and light within deer exclosures. This relationship depends on the seasonal timing of deer habitat use because deer fertilize the soil during winter, but leave during the growing season. The third chapter draws upon an eight-year, 39-stand data set of deer fecal pellet counts in hemlock stands to estimate the amount of nitrogen that deer are depositing in hemlock stands each winter. In stands of high winter deer use, deer-excreted nitrogen inputs consistently exceeded those of atmospheric deposition at the stand scale. At the neighborhood scale, deer-excreted nitrogen was often in excess of atmospheric deposition due to the patchy distribution of deer habitat use. Spatial patterns in habitat use were consistent over the eight-year study at both stand and neighborhood scales. The fourth chapter explores how foraging selectivity by moose interacts with an abiotic resource gradient to influence forest structure and composition. Soil depth on Isle Royale varies from east to west according to glacial history. Fir saplings growing in deeper soils on the west side are generally more palatable forage for moose (lower foliar C:N) than those growing in shallower soils on the east side. Therefore, saplings growing in better conditions are less likely to reach the canopy due to moose browsing, and fir is a smaller overstory component on the west side. Lastly, chapter five focuses on issues surrounding eastern hemlock regeneration failure, which is a habitat type that is important to many wildlife species. Increasing hemlock on the landscape is complicated by several factors including disturbance regime and climate change, in addition to the influence of deer.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

Spatial Variability in Biogenic Gas Accumulations in Peat Soils Is Revealed By Ground Penetrating Radar (GPR)

We performed surface and borehole ground penetrating radar (GPR) tests, together with moisture probe measurements and direct gas sampling to detect areas of biogenic gas accumulation in a northern peatland. The main findings are: (1) shadow zones (signal scattering) observed in surface GPR correlate with areas of elevated CH4 and CO2 concentration; (2) high velocities in zero offset profiles and lower water content inferred from moisture probes correlate with surface GPR shadow zones; (3) zero offset profiles depict depth variable gas accumulation from 0-10% by volume; (4) strong reflectors may represent confining layers restricting upward gas migration. Our results have implications for defining the spatial distribution, volume and movement of biogenic gas in peatlands at multiple scales.

Spatial Patterns in Seasonal and Interannual Variability of Chlorophyll and Sea Surface Temperature in the California Current

Physical forcing and biological response within the California Current System (CCS) are highly variable over a wide range of scales. Satellite remote sensing offers the only feasible means of quantifying this variability over the full extent of the CCS. Using six years (1997-2003) of daily SST and chlorophyll imagery, we map the spatial dependence of dominant temporal variability at resolutions sufficient to identify recurrent mesoscale circulation and local pattern associated with coastal topography. Here we describe mean seasonal cycles and interannual variation; intraseasonal variability is left to a companion paper ( K. R. Legaard and A. C. Thomas, manuscript in preparation, 2006). Coastal upwelling dictates seasonality along north-central California, where weak cycles of SST fluctuate between spring minima and late summer maxima and chlorophyll peaks in early summer. Off northern California, chlorophyll maxima are bounded offshore by the seasonally recurrent upwelling jet. Seasonal cycles differ across higher latitudes and in the midlatitude Southern California Bight, where upwelling winds are less vigorous and/or persistent. Seasonality along south-central Baja is strongly affected by processes other than upwelling, despite year-round upwelling-favorable winds. Interannual variation is generally dominated by El Nino and La Nina conditions. Interannual SST variance is greatest along south-central Baja, although interannual variability constitutes a greater fraction of total variance inshore along southern Oregon and much of California. Patterns of interannual chlorophyll variance are consistent with dominant forcing through the widespread depression and elevation of the nutricline during El Nino and La Nina, respectively. Interannual variability constitutes a greater fraction of total chlorophyll variance offshore.

Spatial Patterns of Intraseasonal Variability of Chlorophyll and Sea Surface Temperature in the California Current

Six years of daily satellite data are used to quantify and map intraseasonal variability of chlorophyll and sea surface temperature (SST) in the California Current. We define intraseasonal variability as temporal variation remaining after removal of interannual variability and stationary seasonal cycles. Semivariograms are used to quantify the temporal structure of residual time series. Empirical orthogonal function (EOF) analyses of semivariograms calculated across the region isolate dominant scales and corresponding spatial patterns of intraseasonal variability. The mode 1 EOFs for both chlorophyll and SST semivariograms indicate a dominant timescale of similar to 60 days. Spatial amplitudes and patterns of intraseasonal variance derived from mode 1 suggest dominant forcing of intraseasonal variability through distortion of large scale chlorophyll and SST gradients by mesoscale circulation. Intraseasonal SST variance is greatest off southern Baja and along southern Oregon and northern California. Chlorophyll variance is greatest over the shelf and slope, with elevated values closely confined to the Baja shelf and extending farthest from shore off California and the Pacific Northwest. Intraseasonal contributions to total SST variability are strongest near upwelling centers off southern Oregon and northern California, where seasonal contributions are weak. Intraseasonal variability accounts for the majority of total chlorophyll variance in most inshore areas save for southern Baja, where seasonal cycles dominate. Contributions of higher EOF modes to semivariogram structure indicate the degree to which intraseasonal variability is shifted to shorter timescales in certain areas. Comparisons of satellite-derived SST semivariograms to those calculated from co-located and concurrent buoy SST time series show similar features.

Global Lithological Map Database v1.0 (gridded to 0.5° spatial resolution)

Lithology describes the geochemical, mineralogical, and physical properties of rocks. It plays a key role in many processes at the Earth surface, especially the fluxes of matter to soils, ecosystems, rivers, and oceans. Understanding these processes at the global scale requires a high resolution description of lithology. A new high resolution global lithological map (GLiM) was assembled from existing regional geological maps translated into lithological information with the help of regional literature. The GLiM represents the rock types of the Earth surface using 1,235,400 polygons. The lithological classification consists of three levels. The first level contains 16 lithological classes comparable to previously applied definitions in global lithological maps. The additional two levels contain 12 and 14 subclasses, respectively, which describe more specific rock attributes. According to the GLiM, the Earth is covered by 64 % sediments (a third of which is carbonates), 13 % metamorphics, 7 % plutonics, and 6 % volcanics, and 10% are covered by water or ice. The high resolution of the GLiM allows observation of regional lithological distributions which often vary from the global average. The GLiM enables regional analysis of Earth surface processes at global scales.

Tie points and log-adjusted depth scales for ODP Leg 207 Sites

Multiple copies of Cretaceous black shales extending from the early Cenomanian to the end of the Santonian were recovered at five sites on Demerara Rise during Leg 207 of the Ocean Drilling Program. These sediments are primarily composed of laminated organic-rich claystones interbedded with coarser, lightly laminated foraminferal-bearing packstones and wackestones. The black shales represent the local expression of widespread organic-rich sedimentation in the Atlantic during the mid-Cretaceous. However, incomplete recovery prevented construction of continuous composite sections, resulting in uncertainties concerning the correct stratigraphic placement of individual cores. By combining high-resolution measurements of bulk density collected shipboard on the multisensor track with continuous downhole measurements of formation resistivity using the Formation MicroScanner, an equivalent logging depth scale was constructed for black shales recovered from Sites 1258, 1260, and 1261. The integrated depths approach centimeter-scale resolution and are supported by comparisons of coarser resolution natural gamma ray emissions collected on cores and through downhole logging operations. The new depths highlight the extent of both intra- and intercore gaps and provide an opportunity to further constrain temporal and spatial paleoceanographic changes captured in proxy records from these sediments.

Benthic cover map of Heron Reef derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Providing accurate maps of coral reefs where the spatial scale and labels of the mapped features correspond to map units appropriate for examining biological and geomorphic structures and processes is a major challenge for remote sensing. The objective of this work is to assess the accuracy and relevance of the process used to derive geomorphic zone and benthic community zone maps for three western Pacific coral reefs produced from multi-scale, object-based image analysis (OBIA) of high-spatial-resolution multi-spectral images, guided by field survey data. Three Quickbird-2 multi-spectral data sets from reefs in Australia, Palau and Fiji and georeferenced field photographs were used in a multi-scale segmentation and object-based image classification to map geomorphic zones and benthic community zones. A per-pixel approach was also tested for mapping benthic community zones. Validation of the maps and comparison to past approaches indicated the multi-scale OBIA process enabled field data, operator field experience and a conceptual hierarchical model of the coral reef environment to be linked to provide output maps at geomorphic zone and benthic community scales on coral reefs. The OBIA mapping accuracies were comparable with previously published work using other methods; however, the classes mapped were matched to a predetermined set of features on the reef.

On the spatial distribution of woody plant species in a tropical montane cloud forest

Abstract The cloud forest is a special type of forest ecosystem that depends on suitable conditions of humidity and temperature to exist; hence, it is a very fragile ecosystem. The cloud forest is also one of the richest ecosystems in terms of species diversity and rate of endemism. However, today, it is one of the most threatened ecosystems in the world. Little is known about tree species distribution and coexistence among cloud forest trees. Trees are essential to understanding ecosystem functioning and maintenance because they support the ecosystem in important ways. For this dissertation, an analysis of woody plant species distribution at a small scale in a north-Peruvian Andean cloud forest was performed, and some of the factors implicated in the observed patterns were identified. Towards that end, different natural factors acting on species distribution within the forest were investigated: (i) intra-specific arrangements, (ii) heterospecific spatial relationships and (iii) relationships with external environmental factors. These analyses were conducted first on standing woody plants and then on seedlings. The woody plants were found to be clumped in the forest, either considering all the species together or each species separately. However, each species presented a specific pattern and specific spatial relationship among different-age individuals. Dispersal mode, growth form and shade tolerance played roles in the final distribution of the species. Furthermore, spatial associations among species, either positive or negative, were observed. These associations were more numerous when considering individuals of the interacting species at different developmental stages, i.e., younger individuals from one species and older individuals from another. Accordingly, competition and facilitation are asymmetric processes and vary throughout the life of an individual. Moreover, some species appear to prefer certain habitat conditions and avoid other habitats. The habitat definition that best explains species distribution is that which includes both environmental and stand characteristics; thus, a combination of these factors is necessary to understanding species' niche preferences. Seedling distribution was also associated with habitat conditions, but these conditions explained less than the 30% of the spatial variation. The position of conspecific adult individuals also affected seedling distribution; although the seedlings of many tree species avoid the vicinity of conspecifics, a few species appeared to prefer the formation of cohorts around their parent trees. The importance of habitat conditions and distance dependence with conspecifics varied among regions within the forest as well as on the developmental stage of the stand. The results from this thesis suggest that different species can coexist within a given space, forming a “puzzle” of species as a result of the intra- and interspecific spatial relationships along with niche preferences and adaptations that operate at different scales. These factors not only affect each species in a different way, but specific preferences also vary throughout species' lifespans. Resumen Resumen El bosque de niebla es uno de los ecosistemas más amenazados del mundo además de ser uno de los más frágiles. Son formaciones azonales que dependen de la existencia de unas condiciones de humedad y temperatura que permitan la formación de nubes que cubran el bosque; lo que dificulta en gran medida su conservación. También es uno de los ecosistemas con mayor riqueza de especies además de tener uno de los mayores porcentajes de endemismos. Uno de los aspectos más importantes para entender el ecosistema, es identificar y entender los elementos que lo componen y los mecanismos que regulan las relaciones entre ellos. Los árboles son el soporte del ecosistema. Sin embargo, apenas hay información sobre la distribución y coexistencia de los árboles en los bosques de niebla. Esta tesis presenta un análisis de la distribución a pequeña escala de las plantas leñosas en un bosque de niebla situado en la cordillera andina del norte de Perú; así como el análisis de algunos de los factores que pueden estar implicados en que se origine la distribución observada. Para este propósito se estudia cómo influyen factores de diferente naturaleza en la distribución de las especies (i) organización intra-específica (ii) relaciones espaciales heterospecíficas y (iii) relación con factores ambientales externos. En estos análisis se estudiaron primero las plantas jóvenes y las adultas, y después las plántulas. Los árboles aparecieron agregados en el bosque, tanto considerando todos a la vez como cuando se estudió cada especie por separado. Sin embargo, cada especie mostró un patrón distinto así como una particular relación espacial entre individuos jóvenes y adultos. El modo de dispersión, la forma de vida y la tolerancia de la especies estuvieron relacionados con el patrón general observado. Se vio también que ciertas especies aparecían relacionadas con otras, tanto de forma positiva (compartiendo zonas) como negativa (apareciendo en áreas distintas). Las asociaciones fueron mucho más numerosas cuando se consideraron los pares de especies en diferente estado de desarrollo, es decir, individuos jóvenes de una especie e individuos mayores de la otra. Eso indicaría que los procesos de competencia y facilitación son asimétricos y además varían durante la vida de la planta. Por otro lado, algunas especies aparecen preferentemente bajo ciertas condiciones de hábitat y evitan otras. La definición de hábitat a la que mejor responden las especies es cuando se incluyen tanto variables ambientales como de masa; así que ambos tipos de variables son necesarias para entender la preferencia de las especies por ciertos nichos. La distribución de las plántulas también estuvo relacionada con condiciones de hábitat, pero eso sólo llegaba a explicar hasta un 30% de la variabilidad espacial. La posición de los adultos de la misma especie también afectó a la distribución de las plántulas. En bastantes especies las plántulas evitan la cercanía de adultos de su misma especie, padres potenciales, aunque algunas especies aisladas mostraron el patrón contrario y aparecieron preferentemente en las mismas áreas que sus padres. La importancia de las condiciones de hábitat y posición de los adultos en la disposición de las plántulas varía de una zona a otra del bosque y además también varía según el estado de desarrollo de la masa.

Lacunarity of the Spatial Distributions of Soil Types in Europe.

Lacunarity as a means of quantifying textural properties of spatial distributions suggests a classification into three main classes of the most abundant soils that cover 92% of Europe. Soils with a well-defined self-similar structure of the linear class are related to widespread spatial patterns that are nondominant but ubiquitous at continental scale. Fractal techniques have been increasingly and successfully applied to identify and describe spatial patterns in natural sciences. However, objects with the same fractal dimension can show very different optical properties because of their spatial arrangement. This work focuses primary attention on the geometrical structure of the geographical patterns of soils in Europe. We made use of the European Soil Database to estimate lacunarity indexes of the most abundant soils that cover 92% of the surface of Europe and investigated textural properties of their spatial distribution. We observed three main classes corresponding to three different patterns that displayed the graphs of lacunarity functions, that is, linear, convex, and mixed. They correspond respectively to homogeneous or self-similar, heterogeneous or clustered and those in which behavior can change at different ranges of scales. Finally, we discuss the pedological implications of that classification.

AquaData and AquaGIS: Two computer utilities for temporal and spatial simulations of water-limited yield with AquaCrop.

The crop simulation model AquaCrop, recently developed by FAO can be used for a wide range of purposes. However, in its present form, its use over large areas or for applications that require a large number of simulations runs (e.g., long-term analysis), is not practical without developing software to facilitate such applications. Two tools for managing the inputs and outputs of AquaCrop, named AquaData and AquaGIS, have been developed for this purpose and are presented here. Both software utilities have been programmed in Delphi v. 5 and in addition, AquaGIS requires the Geographic Information System (GIS) programming tool MapObjects. These utilities allow the efficient management of input and output files, along with a GIS module to develop spatial analysis and effect spatial visualization of the results, facilitating knowledge dissemination. A sample of application of the utilities is given here, as an AquaCrop simulation analysis of impact of climate change on wheat yield in Southern Spain, which requires extensive input data preparation and output processing. The use of AquaCrop without the two utilities would have required approximately 1000 h of work, while the utilization of AquaData and AquaGIS reduced that time by more than 99%. Furthermore, the use of GIS, made it possible to perform a spatial analysis of the results, thus providing a new option to extend the use of the AquaCrop model to scales requiring spatial and temporal analyses.