Application of compositional data analysis to geochemical data of marine sediments

In an earlier investigation (Burger et al., 2000) five sediment cores near the Rodrigues Triple Junction in the Indian Ocean were studied applying classical statistical methods (fuzzy c-means clustering, linear mixing model, principal component analysis) for the extraction of endmembers and evaluating the spatial and temporal variation of geochemical signals. Three main factors of sedimentation were expected by the marine geologists: a volcano-genetic, a hydro-hydrothermal and an ultra-basic factor. The display of fuzzy membership values and/or factor scores versus depth provided consistent results for two factors only; the ultra-basic component could not be identified. The reason for this may be that only traditional statistical methods were applied, i.e. the untransformed components were used and the cosine-theta coefficient as similarity measure. During the last decade considerable progress in compositional data analysis was made and many case studies were published using new tools for exploratory analysis of these data. Therefore it makes sense to check if the application of suitable data transformations, reduction of the D-part simplex to two or three factors and visual interpretation of the factor scores would lead to a revision of earlier results and to answers to open questions . In this paper we follow the lines of a paper of R. Tolosana- Delgado et al. (2005) starting with a problem-oriented interpretation of the biplot scattergram, extracting compositional factors, ilr-transformation of the components and visualization of the factor scores in a spatial context: The compositional factors will be plotted versus depth (time) of the core samples in order to facilitate the identification of the expected sources of the sedimentary process. Kew words: compositional data analysis, biplot, deep sea sediments

Distribució espacial i ús de l'hàbitat de la comunitat de peixos a l'estany de Banyoles

En el marc d'un projecte més ampli sobre la comunitat de peixos de la conca lacustre de Banyoles, s'ha estudiat la distribució espacial de les diferents espècies així com l'ús que aquestes fan de l'hàbitat. El poblament piscícola de l'estany de Banyoles és el resultat d'un llarg historial d'introducció d'espècies exòtiques i extinció de les poblacions autòctones locals. S'ha revisat la seva composició actual detectant un total de 18 espècies (4 autòctones i 14 introduïdes) de les quals només 13 presenten una població estable. L'estudi de l'hàbitat s'ha centrat a l'Estany, l'element principal del sistema lacustre, analitzant per separat la zona litoral i la zona limnètica. En el primer cas s'han realitzat les captures d'individus mitjançant pesca elèctrica des d'una embarcació adaptada per aquesta pràctica. La totalitat del litoral ha estat dividida en trams de característiques homogènies on s'han obtingut les captures per unitat d'esforç per cada espècie. El mostreig s'ha desenvolupat entre l'estiu de 1997 i la primavera del 2000 realitzant un total de 10 campanyes de pesca. Les espècies més abundants al litoral són la perca americana (Micropterus salmoides) i el peix sol (Lepomis gibbosus), essent també presents la perca (Perca fluviatilis), carpa (Cyprinus carpio) i el gardí (Scardinius erythrophthalmus). S'han capturat altres espècies com ara anguila (Anguilla anguilla), bagra (Squalius cephalus), sandra (Sander lucioperca), carpí (Carassius auratus) i madrilleta vera (Rutilus rutilus), però són molt menys abundants en nombre. S'ha examinat, per cadascuna de les espècies, si existeix selecció de l'hàbitat i en cas afirmatiu, quin és el preferent en base a la classificació del litoral en sis tipus de vegetació predominant. Les espècies més abundants, perca americana i peix sol, ocupen tots els hàbitats disponibles però amb una densitat diferent. La perca mostra també una clara selecció de l'hàbitat a favor de les zones molt estructurades amb abundant presència de jonca litoral. Carpa i gardí seleccionen els ambients més fondos amb major presència de matèria orgànica d'origen vegetal procedent del bogar. En general els individus ocupen les zones amb una densitat de vegetació intermèdia, majoritàriament zones de jonca a l'estiu i zones amb mansega a l'hivern, on troben refugi i els recursos tròfics necessaris. La perca americana, a més, presenta una elevada fidelitat a un mateix punt del litoral al llarg de la seva vida. La zona limnètica ha estat prospectada mensualment mitjançant ecosondació, realitzant transectes perpendiculars a l'eix principal de l'Estany, cobrint la seva totalitat. La composició d'espècies s'ha obtingut a partir de les captures fetes amb xarxes (tresmalls) amb periodicitat estacional. L'anàlisi geoestadística de la densitat de peixos ha permès descriure l'estructura espacial d'aquesta a partir dels variogrames, així com la seva variabilitat tant espacial com temporal, i obtenir els mapes de densitat. A l'hivern, la densitat de peixos a la zona limnètica assoleix els seus valors mínims i els individus es troben formant agregats dispersos, pels diferents estrats de fondària. A partir de la primavera la densitat augmenta, pel reclutament i la major freqüència d'individus que abandonen el litoral; la densitat esdevé més homogènia a les primer capes de fondària. A l'estiu la densitat és màxima i l'estrat més homogeni coincideix amb la posició de la termoclina. Aquest estructura varia en disminuir la temperatura i barrejar-se la columna d'aigua, tornant a la situació hivernal. La perca i la madrilleta vera són les espècies predominants en aquest ambient, juntament amb la carpa. La seva distribució no és homogènia i respon a les característiques limnològiques de les diferents cubetes de l'Estany. Una particularitat d'aquest, relacionada amb el seu origen càrstic, es la formació d'una ploma hidrotermal que afecta la distribució dels peixos, probablement en augmentar la terbolesa. S'ha integrat l'ús de l'hàbitat de les espècies que ocupen tant la zona limnètica com la litoral a partir del seguiment d'individus, concretament de perca i bagra. S'ha utilitzat un sistema automàtic de posicionament que estima la localització dels individus marcats amb transmissor de telemetria acústica. L'anàlisi dels desplaçaments mostra un rang superior per la bagra en comparació amb la perca. Ambdues espècies mostren una orientació en els seus desplaçaments. La perca ocupa el litoral a la nit i es desplaça a la zona limnètica de dia, amb un ritme d'activitat marcat per dos màxims coincidint amb la sortida i posta de sol; en canvi la bagra mostra una major activitat nocturna amb zones de repòs properes al litoral. S'ha estimat igualment els dominis vital de cada individu marcat.

Age-specific survival of male Golden-cheeked Warblers on the Fort Hood Military Reservation, Texas

Population models are essential components of large-scale conservation and management plans for the federally endangered Golden-cheeked Warbler (Setophaga chrysoparia; hereafter GCWA). However, existing models are based on vital rate estimates calculated using relatively small data sets that are now more than a decade old. We estimated more current, precise adult and juvenile apparent survival (Φ) probabilities and their associated variances for male GCWAs. In addition to providing estimates for use in population modeling, we tested hypotheses about spatial and temporal variation in Φ. We assessed whether a linear trend in Φ or a change in the overall mean Φ corresponded to an observed increase in GCWA abundance during 1992-2000 and if Φ varied among study plots. To accomplish these objectives, we analyzed long-term GCWA capture-resight data from 1992 through 2011, collected across seven study plots on the Fort Hood Military Reservation using a Cormack-Jolly-Seber model structure within program MARK. We also estimated Φ process and sampling variances using a variance-components approach. Our results did not provide evidence of site-specific variation in adult Φ on the installation. Because of a lack of data, we could not assess whether juvenile Φ varied spatially. We did not detect a strong temporal association between GCWA abundance and Φ. Mean estimates of Φ for adult and juvenile male GCWAs for all years analyzed were 0.47 with a process variance of 0.0120 and a sampling variance of 0.0113 and 0.28 with a process variance of 0.0076 and a sampling variance of 0.0149, respectively. Although juvenile Φ did not differ greatly from previous estimates, our adult Φ estimate suggests previous GCWA population models were overly optimistic with respect to adult survival. These updated Φ probabilities and their associated variances will be incorporated into new population models to assist with GCWA conservation decision making.

The illusion of environmental iodine deficiency

The distribution of endemic goitre in England and Wales was compared with the distribution of environmental iodine (atmospheric deposition, soil, surface water). Despite a very clear goitre belt through the west of England and Wales there was no patterning in the environmental iodine distribution. A clear seasonal variation in depositional iodine exists, with an unusually high concentration of iodine in March 1997. The temporal variation in iodine concentration is determined at the monthly and not the annual level. The presence of endemic goitre is no indicator of how iodine is distributed in the environment or vice versa!

Evaluation of ozone and water vapor fields from the ECMWF reanalysis ERA-40 during 1991-1999 in comparison with UARS satellite and MOZAIC aircraft observations

[ 1] The European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year Reanalysis (ERA-40) ozone and water vapor reanalysis fields during the 1990s have been compared with independent satellite data from the Halogen Occultation Experiment (HALOE) and Microwave Limb Sounder (MLS) instruments on board the Upper Atmosphere Research Satellite (UARS). In addition, ERA-40 has been compared with aircraft data from the Measurements of Ozone and Water Vapour by Airbus In-Service Aircraft (MOZAIC) program. Overall, in comparison with the values derived from the independent observations, the upper stratosphere in ERA-40 has about 5 - 10% more ozone and 15 - 20% less water vapor. This dry bias in the reanalysis appears to be global and extends into the middle stratosphere down to 40 hPa. Most of the discrepancies and seasonal variations between ERA-40 and the independent observations occur within the upper troposphere over the tropics and the lower stratosphere over the high latitudes. ERA-40 reproduces a weaker Antarctic ozone hole, and of less vertical extent, than the independent observations; values in the ozone maximum in the tropical stratosphere are lower for the reanalysis. ERA-40 mixing ratios of water vapor are considerably larger than those for MOZAIC, typically by 20% in the tropical upper troposphere, and they may exceed 60% in the lower stratosphere over high latitudes. The results imply that the Brewer-Dobson circulation in the ECMWF reanalysis system is too fast, as is also evidenced by deficiencies in the way ERA-40 reproduces the water vapor "tape recorder'' signal in the tropical stratosphere. Finally, the paper examines the biases and their temporal variation during the 1990s in the way ERA-40 compares to the independent observations. We also discuss how the evaluation results depend on the instrument used, as well as on the version of the data.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part II: Rockfall activity and talus deformation

During many lava dome-forming eruptions, persistent rockfalls and the concurrent development of a substantial talus apron around the foot of the dome are important aspects of the observed activity. An improved understanding of internal dome structure, including the shape and internal boundaries of the talus apron, is critical for determining when a lava dome is poised for a major collapse and how this collapse might ensue. We consider a period of lava dome growth at the Soufrière Hills Volcano, Montserrat, from August 2005 to May 2006, during which a 100 × 106 m3 lava dome developed that culminated in a major dome-collapse event on 20 May 2006. We use an axi-symmetrical Finite Element Method model to simulate the growth and evolution of the lava dome, including the development of the talus apron. We first test the generic behaviour of this continuum model, which has core lava and carapace/talus components. Our model describes the generation rate of talus, including its spatial and temporal variation, as well as its post-generation deformation, which is important for an improved understanding of the internal configuration and structure of the dome. We then use our model to simulate the 2005 to 2006 Soufrière Hills dome growth using measured dome volumes and extrusion rates to drive the model and generate the evolving configuration of the dome core and carapace/talus domains. The evolution of the model is compared with the observed rockfall seismicity using event counts and seismic energy parameters, which are used here as a measure of rockfall intensity and hence a first-order proxy for volumes. The range of model-derived volume increments of talus aggraded to the talus slope per recorded rockfall event, approximately 3 × 103–13 × 103 m3 per rockfall, is high with respect to estimates based on observed events. From this, it is inferred that some of the volumetric growth of the talus apron (perhaps up to 60–70%) might have occurred in the form of aseismic deformation of the talus, forced by an internal, laterally spreading core. Talus apron growth by this mechanism has not previously been identified, and this suggests that the core, hosting hot gas-rich lava, could have a greater lateral extent than previously considered.

Farm-level performance of genetically modified cotton - A frontier analysis of cotton production in Maharashtra

In this paper, the yield increases resulting from the cultivation of Bt cotton in Maharashtra, India, are analysed. The study relies on commercial farm, rather than trial, data and is among the first of its kind to be based on real farm and market conditions. Findings show that since its commercial release in 2002, Bt cotton has had a significant positive impact on yields and on the economic performance of cotton growers in Maharashtra. This difference remains even after controlling for different soil and insecticide inputs in the production of Bt cotton. There is also significant spatial and temporal variation in this 'benefit', and much depends upon where production is taking place and on the season.

Nectar resource diversity organises flower-visitor community structure

Communities of nectar-producing plants show high spatio-temporal variation in the patterns of volume and concentration presentation. We illustrate a novel approach for quantifying nectar reward structures in complex communities, demonstrating that nectar resource diversity (defined as the variety of nectar volume-concentration combinations available) may be a fundamental factor organising nectarivore communities. In a series of diverse bee and entomophilous flower communities in Israel, our measure of nectar resource diversity alone explains the majority of variation in bee species richness, while other nectar variables (volume, concentration, energy value, and water content) have little predictive value per se. The new measure of nectar resource diversity is highly correlated with floral species richness and particularly with the species richness of annuals, yet it is additive in its effect on bee diversity. We conclude that relying solely upon measurements of mean nectar volume and mean nectar concentration overlooks a key characteristic of community-level reward structure, nectar resource diversity, so that previous studies may have failed to identify an important determinant of flower-visitor community structure.

Pliocene seasonality across the North Atlantic inferred from cheilostome bryozoans

Previous studies have shown an inverse correlation between zooid size in cheilostome bryozoans and ambient water temperature. This relationship underlies the MART technique which uses intracolonial variation in zooid size to predict mean annual range in temperature experienced by bryozoan colonies during their life. Here we apply the MART technique to study Early and Mid Pliocene bryozoans from Central America (Panama, Costa Rica), the USA (Florida, South Carolina, North Carolina, Virginia) and the UK (Suffolk) to reconstruct palaeoseasonality across a range of latitudes for the North Atlantic during the Pliocene Epoch. Compared to the present-day, our analyses suggest greater seasonality (ca 4.5 degrees C) in the southern Caribbean at the time of Cayo Agua Formation deposition (ca 4.25 Ma), in keeping with inferred upwelling prior to the closure of the isthmian barrier at 2.7 Ma. Bryozoans also indicate seasonal upwelling on the Gulf Coast of Florida in a similar manner to the present-day. Because upwelling can be highly localised and prone to spatial and temporal variation in the Gulf of Mexico today, it contributes little to a broad understanding of Pliocene North Atlantic waters. However, MART estimates for the coastal plain region indicate a general reduction in the annual range in temperature relative to the present, suggesting that the colder surface waters that today reach south to Cape Hatteras had less influence in Early to Mid Pliocene times. These results, along with evidence from other proxies, strongly support reduced seasonality and warmer conditions along the eastern seaboard of the USA in the Early to Mid Pliocene. Finally, the MART estimates amongst Coralline Crag localities provide evidence for an increased annual range in temperature in the southern North Sea than at present. Our study shows that bryozoan MART estimates provide a powerful, independent proxy for palaeoseasonality and is the first to demonstrate that the MART technique can be applied to infer palaeoclimates across a wide range of latitudes focusing on a variety of geological formations and geographical regions. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

The effects of temperature and light integral on early vegetative growth and chlorophyll fluorescence of four contrasting genotypes of cacao (Theobroma cacao)

The effect of temperature on early vegetative growth, leaf chlorophyll fluorescence and chlorophyll content was examined on four genotypes of cacao (Amelonado, AMAZ 15–15, SCA 6 and SPEC 54/1). A controlled environment glasshouse was used to simulate the temperature conditions of three cacao-growing regions (Bahia, Brazil; Tafo, Ghana and Lower Perak, Malaysia) over the course of a year. Base temperatures calculated from increments in main stem growth varied from 18.6°C for AMAZ 15/15 to 20.8°C for SPEC 54/1. Temporal variation in Fv/Fm observed for two of the clones (SCA 6 and SPEC 54/1) in two of the compartments were correlated with temperature differences over time. Significant differences were also recorded between genotypes in leaf chlorophyll content. It was shown that variation over time in leaf chlorophyll content could be quantified accurately as a function of temperature and light integral. The results imply that genetic variability exists in cacao in response to temperature stress.

Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

Impact of a drier Early-Mid-Holocene climate upon Amazonian forests

This paper uses a palaeoecological approach to examine the impact of drier climatic conditions of the Early-Mid-Holocene (ca 8000-4000 years ago) upon Amazonia's forests and their fire regimes. Palaeovegetation (pollen data) and palaeofire (charcoal) records are synthesized from 20 sites within the present tropical forest biome, and the underlying causes of any emergent patterns or changes are explored by reference to independent palaeoclimate data and present-day patterns of precipitation, forest cover and fire activity across Amazonia. During the Early-Mid-Holocene, Andean cloud forest taxa were replaced by lowland tree taxa as the cloud base rose while lowland ecotonal areas, which are presently covered by evergreen rainforest, were instead dominated by savannahs and/or semi-deciduous dry forests. Elsewhere in the Amazon Basin there is considerable spatial and temporal variation in patterns of vegetation disturbance and fire, which probably reflects the complex heterogeneous patterns in precipitation and seasonality across the basin, and the interactions between climate change, drought- and fire susceptibility of the forests, and Palaeo-Indian land use. Our analysis shows that the forest biome in most parts of Amazonia appears to have been remarkably resilient to climatic conditions significantly drier than those of today, despite widespread evidence of forest burning. Only in ecotonal areas is there evidence of biome replacement in the Holocene. From this palaeoecological perspective, we argue against the Amazon forest 'dieback' scenario simulated for the future.

A reprocessing for climate of sea surface temperature from the Along-Track Scanning Radiometers: basis in radiative transfer

We present new radiative transfer simulations to support determination of sea surface temperature (SST) from Along Track Scanning Radiometer (ATSR) imagery. The simulations are to be used within the ATSR Reprocessing for Climate project. The simulations are based on the “Reference Forward Model” line-by-line model linked with a sea surface emissivity model that accounts for wind speed and temperature, and with a discrete ordinates scattering model (DISORT). Input to the forward model is a revised atmospheric profile dataset, based on full resolution ERA-40, with a wider range of high-latitude profiles to address known retrieval biases in those regions. Analysis of the radiative impacts of atmospheric trace gases shows that geographical and temporal variation of N2O, CH4, HNO3, and CFC-11 and CFC-12 have effects of order 0.05, 0.2, 0.1 K on the 3.7, 11, 12 μm channels respectively. In addition several trace gases, neglected in previous studies, are included using fixed profiles contributing ~ 0.04 K to top-of-atmosphere BTs. Comparison against observations for ATSR2 and AATSR indicates that forward model biases have been reduced from 0.2 to 0.5 K for previous simulations to ~ 0.1 K.

Cloud optical and microphysical properties derived from ground-based and satellite sensors over a site in the Yangtze Delta region

Comprehensive surface-based retrievals of cloud optical and microphysical properties were made at Taihu, a highly polluted site in the central Yangtze Delta region, during a research campaign from May 2008 to December 2009. Cloud optical depth (COD), effective radius (Re), and liquid water path (LWP) were retrieved from measurements made with a suite of ground-based and spaceborne instruments, including an Analytical Spectral Devices spectroradiometer, a multi␣lter rotating shadowband radiometer, a multichannel microwave radiometer profiler, and the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua satellites. Retrievals from zenith radiance measurements capture better the temporal variation of cloud properties than do retrievals from hemispherical fluxes. Annual mean LWP, COD, and Re are 115.8 ± 90.8 g/m2, 28.5 ± 19.2, and 6.9 ± 4.2 microns. Over 90% of LWP values are less than 250 g/m2. Most of the COD values (>90%) fall between 5 and 60, and ~80% of Re values are less than 10 microns. Maximum (minimum) values of LWP and Re occur in summer (winter); COD is highest in winter and spring. Raining and nonraining clouds have signi␣cant differences in LWP, COD, and Re. Rainfall frequency is best correlated with LWP, followed by COD and Re. Cloud properties retrieved from multiple ground-based instruments are also compared with those from satellite retrievals. On average, relative to surface retrievals, mean differences of satellite retrievals in cloud LWP, COD, and Re were -33.6 g/m2 (-26.4%), -5.8 (-31.4%), and 2.9 ␣m (29.3%) for 11 MODIS-Terra overpasses and -43.3 g/m2 (-22.3%), -3.0 (-10.0%), and -1.3 ␣m (-12.0%) for 8 MODIS-Aqua overpasses, respectively. These discrepancies indicate that MODIS cloud products still suffer from large uncertainties in this region.