Monitoring of active layer dynamics at two transects on Svalbard using multi-channel ground-penetrating radar

Multi-channel ground-penetrating radar is used to investigate the late-summer evolution of the thaw depth and the average soil water content of the thawed active layer at a high-arctic continuous permafrost site on Svalbard, Norway. Between mid of August and mid of September 2008, five surveys have been conducted over transect lengths of 130 and 175 m each. The maximum thaw depths range from 1.6 m to 2.0 m, so that they are among the deepest thaw depths recorded for Svalbard so far. The thaw depths increase by approximately 0.2 m between mid of August and beginning of September and subsequently remain constant until mid of September. The thaw rates are approximately constant over the entire length of the transects within the measurement accuracy of about 5 to 10 cm. The average volumetric soil water content of the thawed soil varies between 0.18 and 0.27 along the investigated transects. While the measurements do not show significant changes in soil water content over the first four weeks of the study, strong precipitation causes an increase in average soil water content of up to 0.04 during the last week. These values are in good agreement with evapotranspiration and precipitation rates measured in the vicinity of the the study site. While we cannot provide conclusive reasons for the detected spatial variability of the thaw depth at the study site, our measurements show that thaw depth and average soil water content are not directly correlated. The study demonstrates the potential of multi-channel ground-penetrating radar for mapping thaw depth in permafrost areas. The novel non-invasive technique is particularly useful when the thaw depth exceeds 1.5 m, so that it is hardly accessible by manual probing. In addition, multi-channel ground-penetrating radar holds potential for mapping the latent heat content of the active layer and for estimating weekly to monthly averages of the ground heat flux during the thaw period.

Fidelity bias in mollusk assemblages from coastal lagoons of Southern Brazil

The South America southern coast exhibits many outcrops with abundant shell beds, from the Pleistocene through the Recent. How much biological information is preserved within these shell beds? Or, what is the actual probability a living community has to leave a fossil record corresponding to these shell deposits? Although ecological and biogeographical aspects might had been pointed, considering these temporal scales, up to the moment there is no taphonomically-oriented studies available. Quantitative comparisons between living (LAs), death (DAs) and fossil assemblages (FAs) are important not only in strictly taphonomic studies, but have grown a leading tool for conservation paleobiology analysis. Comparing LAs, DAs and FAs from estuaries and lagoons in the Rio Grande do Sul Coastal Plain makes possible to quantitatively understand the nature and quantity of biological information preserved in fossil associations in Holocene lagoon facies. As already noted by several authors, spatial scale parts the analysis, but we detected that the FAs refl ects live ones, rather than dead ones, as previously not realized. The results herein obtained illustrates that species present in DA are not as good preserved in recent (Holocene) fossil record as originally thought. Strictly lagoon species are most prone to leave fossil record. The authors consider that the fi delity pattern here observed for estuarine mollusks to be driven by (i) high temporal and spatial variability in the LAs, (ii) spatial mixing in the DA and (iii) differential preservation of shells, due to long residence times in the taphonomically active zone.

Spring snowmelt and midwinter thaw and refreeze north of 60°N based on Seawinds QuikScat 2000-2009

The scatterometer SeaWinds on QuikSCAT provided regular measurements at Ku-band from 1999 to 2009. Although it was designed for ocean applications, it has been frequently used for the assessment of seasonal snowmelt patterns aside from other terrestrial applications such as ice cap monitoring, phenology and urban mapping. This paper discusses general data characteristics of SeaWinds and reviews relevant change detection algorithms. Depending on the complexity of the method, parameters such as long-term noise and multiple event analyses were incorporated. Temporal averaging is a commonly accepted preprocessing step with consideration of diurnal, multi-day or seasonal averages.

Size-partitioned phytoplankton carbon concentrations retrieved from ocean color data, links to data in netCDF format

Owing to their important roles in biogeochemical cycles, phytoplankton functional types (PFTs) have been the aim of an increasing number of ocean color algorithms. Yet, none of the existing methods are based on phytoplankton carbon (C) biomass, which is a fundamental biogeochemical and ecological variable and the "unit of accounting" in Earth system models. We present a novel bio-optical algorithm to retrieve size-partitioned phytoplankton carbon from ocean color satellite data. The algorithm is based on existing methods to estimate particle volume from a power-law particle size distribution (PSD). Volume is converted to carbon concentrations using a compilation of allometric relationships. We quantify absolute and fractional biomass in three PFTs based on size - picophytoplankton (0.5-2 µm in diameter), nanophytoplankton (2-20 µm) and microphytoplankton (20-50 µm). The mean spatial distributions of total phytoplankton C biomass and individual PFTs, derived from global SeaWiFS monthly ocean color data, are consistent with current understanding of oceanic ecosystems, i.e., oligotrophic regions are characterized by low biomass and dominance of picoplankton, whereas eutrophic regions have high biomass to which nanoplankton and microplankton contribute relatively larger fractions. Global climatological, spatially integrated phytoplankton carbon biomass standing stock estimates using our PSD-based approach yield - 0.25 Gt of C, consistent with analogous estimates from two other ocean color algorithms and several state-of-the-art Earth system models. Satisfactory in situ closure observed between PSD and POC measurements lends support to the theoretical basis of the PSD-based algorithm. Uncertainty budget analyses indicate that absolute carbon concentration uncertainties are driven by the PSD parameter No which determines particle number concentration to first order, while uncertainties in PFTs' fractional contributions to total C biomass are mostly due to the allometric coefficients. The C algorithm presented here, which is not empirically constrained a priori, partitions biomass in size classes and introduces improvement over the assumptions of the other approaches. However, the range of phytoplankton C biomass spatial variability globally is larger than estimated by any other models considered here, which suggests an empirical correction to the No parameter is needed, based on PSD validation statistics. These corrected absolute carbon biomass concentrations validate well against in situ POC observations.

Results of wave modelling for Moreton Bay, Southeast Queensland, and a reef lagoon off Lizard Island, Great Barrier Reef, Australia

The distribution, abundance, behaviour, and morphology of marine species is affected by spatial variability in the wave environment. Maps of wave metrics (e.g. significant wave height Hs, peak energy wave period Tp, and benthic wave orbital velocity URMS) are therefore useful for predictive ecological models of marine species and ecosystems. A number of techniques are available to generate maps of wave metrics, with varying levels of complexity in terms of input data requirements, operator knowledge, and computation time. Relatively simple "fetch-based" models are generated using geographic information system (GIS) layers of bathymetry and dominant wind speed and direction. More complex, but computationally expensive, "process-based" models are generated using numerical models such as the Simulating Waves Nearshore (SWAN) model. We generated maps of wave metrics based on both fetch-based and process-based models and asked whether predictive performance in models of benthic marine habitats differed. Predictive models of seagrass distribution for Moreton Bay, Southeast Queensland, and Lizard Island, Great Barrier Reef, Australia, were generated using maps based on each type of wave model. For Lizard Island, performance of the process-based wave maps was significantly better for describing the presence of seagrass, based on Hs, Tp, and URMS. Conversely, for the predictive model of seagrass in Moreton Bay, based on benthic light availability and Hs, there was no difference in performance using the maps of the different wave metrics. For predictive models where wave metrics are the dominant factor determining ecological processes it is recommended that process-based models be used. Our results suggest that for models where wave metrics provide secondarily useful information, either fetch- or process-based models may be equally useful.

Grain-size analysis of surface sediments, continental margin off NW Africa

The terrigenous fraction of seabed sediments recovered along the north-west African continental margin illustrates spatial variability in grain size attributed to different transport mechanisms. Three subpopulations are determined from the grain-size analyses (n = 78) of the carbonate-free silt fraction applying an end-member modelling algorithm (G. J. Weltje, 1997). The two coarsest end-members are interpreted as representing aeolian dust, and the fine-grained end-member is related to fluvial supply. The end-member model thus allows aeolian fallout to be distinguished from fluvial-sourced mud in this area. The relative contributions of the end-members show distinct regional variations that can be related to different transport processes and pathways. Understanding present-day sediment dispersal and mixing is important for a better understanding of older sedimentary records and palaeoclimate reconstructions in the region.

(Table 3) Oxygen and carbon isotopic record of living (Rose Bengal Stained) benthic foraminifera of North Atlantic surface sediments

Variation of the d13C of living (Rose Bengal stained) deep-sea benthic foraminifera is documented from two deep-water sites (~2430 and ~3010 m) from a northwest Atlantic Ocean study area 275 km south of Nantucket Island. The carbon isotopic data of Hoeglundina elegans and Uvigerina peregrina from five sets of Multicorer and Soutar Box Core samples taken over a 10-month interval (March, May, July, and October 1996 and January 1997) are compared with an 11.5 month time series of organic carbon flux to assess the effect of organic carbon flux on the carbon isotopic composition of dominant taxa. Carbon isotopic data of Hoeglundina elegans at 3010 m show 0.3 per mil lower mean values following an organic carbon flux maximum resulting from a spring phytoplankton bloom. This d13C change following the spring bloom is suggested to be due to the presence of a phytodetritus layer on the seafloor and the subsequent depletion of d13C in the pore waters within the phytodetritus and overlying the sediment surface. Carbon isotopic data of H. elegans from the 2430 m site show an opposite pattern to that found at 3010 m with a d13C enrichment following the spring bloom. This different pattern may be due to spatial variation in phytodetritus deposition and resuspension or to a limited number of specimens recovered from the March 1996 cruise. The d13C of Uvigerina peregrina at 2430 m shows variation over the 10 month interval, but an analysis of variance shows that the variability is more consistent with core and subcore variability than with seasonal changes. The isotopic analyses are grouped into 100 µm size classes on the basis of length measurements of individual specimens to evaluate d13C ontogenetic changes of each species. The data show no consistent patterns between size classes in the d13C of either H. elegans or U. peregrina. These results suggest that variation in organic carbon flux does not preferentially affect particular size classes, nor do d13C ontogenetic changes exist within the >250 to >750 µm size range for these species at this locality. On the basis of the lack of ontogenetic changes a range of sizes of specimens from a sample can be used to reconstruct d13C in paleoceanographic studies. The prediction standard deviation, which is composed of cruise, core, subcore, and residual (replicate) variability, provides an estimate of the magnitude of variability in fossil d13C data; it is 0.27 per mil for H. elegans at 3010 m and 0.4 per mil for U. peregrina at the 2430 m site. Since these standard deviations are based on living specimens, they should be regarded as minimum estimates of variability for fossil data based on single specimen analyses. Most paleoceanographic reconstructions are based on the analysis of multiple specimens, and as a result, the standard error would be expected to be reduced for any particular sample. The reduced standard error resulting from the analysis of multiple specimens would result in the seasonal and spatial variability observed in this study having little impact on carbon isotopic records.

Paleoclimate investiations on sediment core PG1984 from Lake Sysy-Kyuele, northern Siberia

In northeastern Siberia, Russia, a 1.2 m sediment core was retrieved and radiocarbon dated from a small and shallow lake located at the western side of the lower Lena River (N 69°24', E 123°50', 81 m a.s.l.). The objective of this paper is to reconstruct the palaeoenvironmental variability and to infer major palaeoclimate trends that have occurred since ~ 13.3 cal. kyrs BP. We analysed the diatom assemblages, sedimentology (grain size, total organic carbon (TOC), total nitrogen (TN)), and the elemental and mineralogical composition using X-ray fluorescence (XRF) and X-ray diffractometry (XRD) of the sediment core. Our results show parallel changes in the diatom species composition and sediment characteristics. Enhanced minerogenic sediment input and the occurrence of pyrite is indicative of a cold period between ~ 12.7-11.6 cal. kyrs BP. The diatom data enable a qualitative inference about the local ecological conditions to be made, and reveal an oligotrophic lake system with alkaline and cold conditions during the earliest Holocene. Moderately warmer climates are inferred for the period from ~ 9.1 to 5.7 cal. kyrs BP. The major shift in the diatom assemblage, from dominance of small benthic fragilarioid taxa to a more complex diatom flora with an influx of several achnanthoid and naviculoid diatom species, occurred after a transitional period of about 1400 years (7.1 to 5.7 cal. kyrs BP) at ~ 5.7 cal. kyrs BP, indicating a circumneutral and warmer hydrological regime during the Holocene thermal maximum (HTM). Diatom valve concentrations declined starting ~ 2.8 cal. kyrs BP, but have been rising again since less than or equalt to 600 cal. years BP. This has occurred in parallel to the increased presence of acidophilous diatom taxa (e.g. Eunotia spp.) and decreased presence of small benthic fragilarioid species in the most recent sediments, which is interpreted as the result of neoglacial cooling and subsequent recent climate warming. Our findings are compared to other lake-inferred climate reconstructions along the Lena River. We conclude that the timing and spatial variability of the HTM in the lower Lena River area reveal a temporal delay from north to south.

Underway measurements of dissolved gases along POLARSTERN cruise track ANT-XXVII/2

We present measurements of pCO2, O2 concentration, biological oxygen saturation (Delta O2/Ar) and N2 saturation (Delta N2) in Southern Ocean surface waters during austral summer, 2010-2011. Phytoplankton biomass varied strongly across distinct hydrographic zones, with high chlorophyll a (Chla) concentrations in regions of frontal mixing and sea-ice melt. pCO2 and Delta O2 /Ar exhibited large spatial gradients (range 90 to 450 µatm and -10 to 60%, respectively) and co-varied strongly with Chla. However, the ratio of biological O2 accumulation to dissolved inorganic carbon (DIC) drawdown was significantly lower than expected from photosynthetic stoichiometry, reflecting the differential time-scales of O2 and CO2 air-sea equilibration. We measured significant oceanic CO2 uptake, with a mean air-sea flux (~ -20 mmol m-2 d-1) that significantly exceeded regional climatological values. N2 was mostly supersaturated in surface waters (mean Delta N2 of +2.5 %), while physical processes resulted in both supersaturation and undersaturation of mixed layer O2 (mean Delta O2phys = 2.1 %). Box model calculations were able to reproduce much of the spatial variability of Delta N2 and Delta O2phys along the cruise track, demonstrating significant effects of air-sea exchange processes (e.g. atmospheric pressure changes and bubble injection) and mixed layer entrainment on surface gas disequilibria. Net community production (NCP) derived from entrainment-corrected surface Delta O2 /Ar data, ranged from ~ -40 to > 300 mmol O2 m-2 d-1 and showed good coherence with independent NCP estimates based on seasonal mixed layer DIC deficits. Elevated NCP was observed in hydrographic frontal zones and regions of sea-ice melt with shallow mixed layer depths, reflecting the importance of mixing in controlling surface water light and nutrient availability.

(Table S1) Stable oxygen isotopic values and Mg/Ca ratios of Globigerinoides ruber w from sediment core AAS9_21

Variations in sea surface temperature (SST), d18O of sea water (?18Ow), and salinity were reconstructed for the past 68 ka using a sediment core (AAS9/21) from the eastern Arabian Sea (EAS) in order to understand the changes in evaporation and precipitation associated with the monsoon system. The Mg/Ca-derived SST record varies by ~4°C; it shows that marine isotope stage (MIS) 4 was warmer than MIS 3, that the Last Glacial Maximum was 4°C cooler than the present, and that there was a 2°C increase within the Holocene. MIS 4 records higher d18Ow and salinity values than MIS 2, suggesting variable flow of low-salinity Bay of Bengal flow into the EAS during glacial periods. The transition from MIS 4 to MIS 3 was marked with a conspicuous shift from higher to lower d18Ow values, which reflects a decrease in the evaporation-precipitation budget in the EAS, perhaps due to the strengthening of southwest monsoon. Monsoon reconstructions based on d18Ow reveal that monsoon-driven precipitation was higher during MIS 3 and MIS 1 and was lower during MIS 2 and MIS 4. This is consistent with earlier monsoon reconstructions based on upwelling indices from the western Arabian Sea. However, the amplitude of monsoon fluctuations derived through upwelling indices and d18Ow varies significantly, which may indicate spatial variability of monsoon rainfall.

Mediterranean cold-water coral measurements

This study presents aggradation rates supplemented for the first time by carbonate accumulation rates from Mediterranean cold-water coral sites considering three different regional and geomorphological settings: (i) a cold-water coral ridge (eastern Melilla coral province, Alboran Sea), (ii) a cold-water coral rubble talus deposit at the base of a submarine cliff (Urania Bank, Strait of Sicily) and (iii) a cold-water coral deposit rooted on a predefined topographic high overgrown by cold-water corals (Santa Maria di Leuca coral province, Ionian Sea). The mean aggradation rates of the respective cold-water coral deposits vary between 10 and 530 cm kyr?1 and the mean carbonate accumulation rates range between 8 and 396 g cm?2 kyr?1 with a maximum of 503 g cm?2 kyr?1 reached in the eastern Melilla coral province. Compared to other deep-water depositional environments the Mediterranean cold-water coral sites reveal significantly higher carbonate accumulation rates that were even in the range of the highest productive shallow-water Mediterranean carbonate factories (e.g. Cladocora caespitosa coral reefs). Focusing exclusively on cold-water coral occurrences, the carbonate accumulation rates of the Mediterranean cold-water coral sites are in the lower range of those obtained for the prolific Norwegian coral occurrences, but exhibit much higher rates than the cold-water coral mounds off Ireland. This study clearly indicates that cold-water corals have the potential to act as important carbonate factories and regional carbonate sinks within the Mediterranean Sea. Moreover, the data highlight the potential of cold-water corals to store carbonate with rates in the range of tropical shallow-water reefs. In order to evaluate the contribution of the cold-water coral carbonate factory to the regional or global carbonate/carbon cycle, an improved understanding of the temporal and spatial variability in aggradation and carbonate accumulation rates and areal estimates of the respective regions is needed.

Remote Sensing of Vegetation Change Across a Latitudinal Gradient in the Canadian Arctic

Global air surface temperatures and precipitation have increased over the last several decades resulting in a trend of greening across the Circumpolar Arctic. The spatial variability of warming and the inherent effects on plant communities has not proven to be uniform or homogeneous on global or local scales. We can apply remote sensing vegetation indices such as the Normalized Difference Vegetation Index (NDVI) to map and monitor vegetation change (e.g., phenology, greening, percent cover, and biomass) over time. It is important to document how Arctic vegetation is changing, as it will have large implications related to global carbon and surface energy budgets. The research reported here examined vegetation greening across different spatial and temporal scales at two disparate Arctic sites: Apex River Watershed (ARW), Baffin Island, and Cape Bounty Arctic Watershed Observatory (CBAWO), Melville Island, NU. To characterize the vegetation in the ARW, high spatial resolution WorldView-2 data were processed to create a supervised land-cover classification and model percent vegetation cover (PVC) (a similar process had been completed in a previous study for the CBAWO). Meanwhile, NDVI data spanning the past 30 years were derived from intermediate resolution Landsat data at the two Arctic sites. The land-cover classifications at both sites were used to examine the Landsat NDVI time series by vegetation class. Climate variables (i.e., temperature, precipitation and growing season length (GSL) were examined to explore the potential relationships of NDVI to climate warming. PVC was successfully modeled using high resolution data in the ARW. PVC and plant communities appear to reside along a moisture and altitudinal gradient. The NDVI time series demonstrated an overall significant increase in greening at the CBAWO (High Arctic site), specifically in the dry and mesic vegetation type. However, similar overall greening was not observed for the ARW (Low Arctic site). The overall increase in NDVI at the CBAWO was attributed to a significant increase in July temperatures, precipitation and GSL.

RELIBILITY ANALYSIS OF SIMPLE SLOPES AND SOIL-STRUCTURES WITH LINEAR LIMIT STATES

The first objective of this research was to develop closed-form and numerical probabilistic methods of analysis that can be applied to otherwise conventional methods of unreinforced and geosynthetic reinforced slopes and walls. These probabilistic methods explicitly include random variability of soil and reinforcement, spatial variability of the soil, and cross-correlation between soil input parameters on probability of failure. The quantitative impact of simultaneously considering the influence of random and/or spatial variability in soil properties in combination with cross-correlation in soil properties is investigated for the first time in the research literature. Depending on the magnitude of these statistical descriptors, margins of safety based on conventional notions of safety may be very different from margins of safety expressed in terms of probability of failure (or reliability index). The thesis work also shows that intuitive notions of margin of safety using conventional factor of safety and probability of failure can be brought into alignment when cross-correlation between soil properties is considered in a rigorous manner. The second objective of this thesis work was to develop a general closed-form solution to compute the true probability of failure (or reliability index) of a simple linear limit state function with one load term and one resistance term expressed first in general probabilistic terms and then migrated to a LRFD format for the purpose of LRFD calibration. The formulation considers contributions to probability of failure due to model type, uncertainty in bias values, bias dependencies, uncertainty in estimates of nominal values for correlated and uncorrelated load and resistance terms, and average margin of safety expressed as the operational factor of safety (OFS). Bias is defined as the ratio of measured to predicted value. Parametric analyses were carried out to show that ignoring possible correlations between random variables can lead to conservative (safe) values of resistance factor in some cases and in other cases to non-conservative (unsafe) values. Example LRFD calibrations were carried out using different load and resistance models for the pullout internal stability limit state of steel strip and geosynthetic reinforced soil walls together with matching bias data reported in the literature.