ALBERT

Description: Not all calcite ballast is created equal: differing effects of foraminiferan and coccolith calcite on the formation and sinking of aggregates Biogeosciences Discussions, 10, 14861-14885, 2013 Author(s): K. Schmidt, C. L. De La Rocha, M. Gallinari, and G. Cortese Correlation between particulate organic carbon (POC) and calcium carbonate sinking through the deep ocean has led to the idea that ballast provided by calcium carbonate is important for the export of POC from the surface ocean. While this idea is certainly to some extent true, it is worth considering in more nuance, for example, examining the different effects on the aggregation and sinking of POC of small, non-sinking calcite particles like coccoliths and large, rapidly sinking calcite like planktonic foraminiferan tests. We have done that here in a simple experiment carried out in roller tanks that allow particles to sink continuously without being impeded by container walls. Coccoliths were efficiently incorporated into aggregates that formed during the experiment, increasing their sinking speed compared to similarly sized aggregates lacking added calcite ballast. The foraminiferan tests, which sank as fast as 700 m d −1 , became associated with only very minor amounts of POC. In addition, when they collided with other, larger, foraminferan-less aggregates, they fragmented them into two smaller, more slowly sinking aggregates. While these effects were certainly exaggerated within the confines of the roller tanks, they clearly demonstrate that calcium carbonate ballast is not just calcium carbonate ballast- different forms of calcium carbonate ballast have notably different effects on POC aggregation, sinking, and export.

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

Description: Beyond multifractional Brownian motion: new stochastic models for geophysical modelling Nonlinear Processes in Geophysics, 20, 643-655, 2013 Author(s): J. Lévy Véhel Multifractional Brownian motion (mBm) has proved to be a useful tool in various areas of geophysical modelling. Although a versatile model, mBm is of course not always an adequate one. We present in this work several other stochastic processes which could potentially be useful in geophysics. The first alternative type is that of self-regulating processes : these are models where the local regularity is a function of the amplitude, in contrast to mBm where it is tuned exogenously. We demonstrate the relevance of such models for digital elevation maps and for temperature records. We also briefly describe two other types of alternative processes, which are the counterparts of mBm and of self-regulating processes when the intensity of local jumps is considered in lieu of local regularity: multistable processes allow one to prescribe the local intensity of jumps in space/time, while this intensity is governed by the amplitude for self-stabilizing processes .

Description: Phenology as a strategy for carbon optimality: a global model Biogeosciences Discussions, 10, 15107-15152, 2013 Author(s): S. Caldararu, D. W. Purves, and P. I. Palmer Phenology is essential to our understanding of biogeochemical cycles and the climate system. We develop a global mechanistic model of leaf phenology based on the hypothesis that phenology is a strategy for optimal carbon gain at the canopy level so that trees adjust leaf gains and losses in response to environmental factors such as light, temperature and soil moisture, to achieve maximum carbon assimilation. We fit this model to five years of satellite observations of leaf area index (LAI) using a Bayesian fitting algorithm. We show that our model is able to reproduce phenological patterns for all vegetation types and use it to explore variations in growing season length and the climate factors that limit leaf growth for different biomes. Phenology in wet tropical areas is limited by leaf age physiological constraints while at higher latitude leaf seasonality is limited by low temperature and light availability. Leaf growth in grassland regions is limited by water availability but often in combination with other factors. This model will advance the current understanding of phenology for ecosystem carbon models and our ability to predict future phenological behaviour.

Description: Timing of sea ice retreat can alter phytoplankton community structure in the western Arctic Ocean Biogeosciences Discussions, 10, 15153-15180, 2013 Author(s): given_name prefix surname suffix, A. Fujiwara, T. Hirawake, K. Suzuki, I. Imai, and S.-I. Saitoh This study assesses the response of phytoplankton assemblages to recent climate change, especially with regard to the shrinking of sea ice in the northern Chukchi Sea of the western Arctic Ocean. Distribution patterns of phytoplankton groups in the late summers of 2008–2010 were analyzed based on HPLC pigment signatures and, the following four major algal groups were inferred via multiple regression and cluster analyses: prasinophytes, diatoms, haptophytes and dinoflagellates. A remarkable interannual difference in the distribution pattern of the groups was found in the northern basin area. Haptophytes dominated and dispersed widely in warm surface waters in 2008, whereas prasinophytes dominated in cold water in 2009 and 2010. A difference in the onset date of sea ice retreat was evident among years – the sea ice retreat in 2008 was 1–2 months earlier than in 2009 and 2010. The spatial distribution of early sea ice retreat matched the areas in which a shift in algal community composition was observed. Steel-Dwass's multiple comparison tests were used to assess the physical, chemical and biological parameters of the four clusters. We found a statistically significant difference in temperature between the haptophyte-dominated cluster and the other clusters, suggesting that the change in the phytoplankton communities was related to the earlier sea ice retreat in 2008 and the corollary increase in sea surface temperatures. Longer periods of open water during the summer, which are expected in the future, may affect food webs and biogeochemical cycles in the western Arctic due to shifts in phytoplankton community structure.

Description: Snow on the Ross Ice Shelf: comparison of reanalyses and observations from automatic weather stations The Cryosphere, 7, 1399-1410, 2013 Author(s): L. Cohen and S. Dean Snow accumulation measurements from automatic weather stations (AWS) around the Ross Ice Shelf (RIS), Antarctica, are used to provide a new set of ground-based observations which are compared to precipitation from the ECMWF ERA-Interim and NCEP/NCAR Reanalysis-2 datasets. The high temporal resolution of the AWS snow accumulation measurements allow for an event-based comparison of reanalyses precipitation to the in situ observations. Snow accumulation records from nine AWS provide multiple years of accumulation data between 2008 and 2012 over a relatively large, homogeneous region of Antarctica, and also provide the basis for a statistical evaluation of accumulation and precipitation events. The complex effects of wind on snow accumulation (which can both limit and enhance accumulation) complicate the use of the accumulation measurements, but this analysis shows that they can provide a valuable source of ground-based observations for comparisons to modelled precipitation on synoptic timescales. The analysis shows that ERA-Interim reproduces more precipitation events than NCEP-2, and these events correspond to an average 8.2% more precipitation. Significant correlations between reanalyses and AWS event sizes are seen at several stations and show that ERA-Interim consistently produces larger precipitation events than NCEP-2.

Description: Uptake of phytodetritus by benthic foraminifera under oxygen depletion at the Indian Margin (Arabian Sea) Biogeosciences Discussions, 10, 15305-15335, 2013 Author(s): A. J. Enge, U. Witte, M. Kucera, and P. Heinz Benthic foraminifera in sediments on the Indian margin of the Arabian Sea where the oxygen minimum zone (OMZ) impinges on the continental slope are exposed to particularly severe levels of oxygen depletion. Food supply for the benthic community is high but delivered in distinct pulses during upwelling and water mixing events associated with summer and winter monsoon periods. In order to investigate the response by benthic foraminifera to such pulsed food delivery under oxygen concentrations of less than 0.1 mL L −1 (4.5 μmol L −1 ), an in situ isotope labeling experiment ( 13 C, 15 N) was performed at the western continental slope of India at 540 m water depth (OMZ core region). The assemblage of living foraminifera (〉125 μm) in the uppermost centimeter at this depth is characterized by an unexpectedly high population density of 3982 ind. 10 cm −2 and a strong dominance by few calcareous species. For the experiment, we concentrated on the nine most abundant taxa, which constitute 93% of the entire foraminifera population at 0–1 cm sediment depth. Increased concentrations of 13 C and 15 N in the cytoplasm indicate that all investigated taxa took up the labeled phytodetritus during the 4 day experimental phase. In total, these nine species had assimilated 113.8 mg C m −2 (17.5% of the total added carbon). The uptake of nitrogen by the three most abundant taxa ( Bolivina aff. B. dilatata , Cassidulina sp., Bulimina gibba ) was 2.7 mg N m −2 (2% of the total added nitrogen) and showed the successful application of 15 N as tracer in foraminiferal studies. The short-term response to the offered phytodetritus varied largely among foraminiferal species with Uvigerina schwageri being by far the most important species in short-term processing whereas the most abundant taxa Bolivina aff. B. dilatata and Cassidulina sp. showed comparably low uptake of the offered food. We suggest that the observed species-specific differences are related to individual biomass of species and to specific feeding preferences. The high numbers of living foraminifera and their rapid response to deposited fresh phytodetritus demonstrate the importance of foraminifera in short-term carbon cycling under oxygen-depleted conditions. We propose that foraminifera at the studied site benefit from unique adaptations in their metabolisms to nearly anoxic conditions as well as from the exclusion of macrofauna and the resulting relaxation of competition for food and low predation pressure.

Description: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity The Cryosphere, 7, 1473-1480, 2013 Author(s): H. Löwe, F. Riche, and M. Schneebeli Finding relevant microstructural parameters beyond density is a longstanding problem which hinders the formulation of accurate parameterizations of physical properties of snow. Towards a remedy, we address the effective thermal conductivity tensor of snow via anisotropic, second-order bounds. The bound provides an explicit expression for the thermal conductivity and predicts the relevance of a microstructural anisotropy parameter Q , which is given by an integral over the two-point correlation function and unambiguously defined for arbitrary snow structures. For validation we compiled a comprehensive data set of 167 snow samples. The set comprises individual samples of various snow types and entire time series of metamorphism experiments under isothermal and temperature gradient conditions. All samples were digitally reconstructed by micro-computed tomography to perform microstructure-based simulations of heat transport. The incorporation of anisotropy via Q considerably reduces the root mean square error over the usual density-based parameterization. The systematic quantification of anisotropy via the two-point correlation function suggests a generalizable route to incorporate microstructure into snowpack models. We indicate the inter-relation of the conductivity to other properties and outline a potential impact of Q on dielectric constant, permeability and adsorption rate of diffusing species in the pore space.

Description: Summertime canopy albedo is sensitive to forest thinning Biogeosciences Discussions, 10, 15373-15414, 2013 Author(s): J. Otto, D. Berveiller, F.-M. Bréon, N. Delpierre, G. Geppert, A. Granier, W. Jans, A. Knohl, A. Kuusk, B. Longdoz, E. Moors, M. Mund, B. Pinty, M.-J. Schelhaas, and S. Luyssaert Despite an emerging body of literature linking canopy albedo to forest management, understanding of the process is still fragmented. We combined a stand-level forest gap model with a canopy radiation transfer model and satellite-derived model parameters to quantify the effects of forest thinning, that is removing trees at a certain time during the forest rotation, on summertime canopy albedo. The effects of different forest species (pine, beech, oak) and four thinning strategies (light to intense thinning regimes) were examined. During stand establishment, summertime canopy albedo is driven by tree species. In the later stages of stand development, the effect of tree species on summertime canopy albedo decreases in favour of an increasing influence of forest thinning on summertime canopy albedo. These trends continue until the end of the rotation where thinning explains up to 50% of the variance in near-infrared canopy albedo and up to 70% of the variance in visible canopy albedo. More intense thinning lowers the summertime shortwave albedo in the canopy by as much as 0.02 compared to unthinned forest. The structural changes associated with forest thinning can be described by the change in LAI in combination with crown volume. However, forests with identical canopy structure can have different summertime albedo values due to their location: the further north a forest is situated, the more the solar zenith angle increases and thus the higher is the summertime canopy albedo, independent of the wavelength. Despite the increase of absolute summertime canopy albedo values with latitude, the difference in canopy albedo between managed and unmanaged forest decreases with increasing latitude. Forest management thus strongly altered summertime forest albedo.

Description: Influence of high-order mechanics on simulation of glacier response to climate change: insights from Haig Glacier, Canadian Rocky Mountains The Cryosphere, 7, 1527-1541, 2013 Author(s): S. Adhikari and S. J. Marshall Evolution of glaciers in response to climate change has mostly been simulated using simplified dynamical models. Because these models do not account for the influence of high-order physics, corresponding results may exhibit some biases. For Haig Glacier in the Canadian Rocky Mountains, we test this hypothesis by comparing simulation results obtained from 3-D numerical models that deal with different assumptions concerning physics, ranging from simple shear deformation to comprehensive Stokes flow. In glacier retreat scenarios, we find a minimal role of high-order mechanics in glacier evolution, as geometric effects at our site (the presence of an overdeepened bed) result in limited horizontal movement of ice (flow speed on the order of a few meters per year). Consequently, high-order and reduced models all predict that Haig Glacier ceases to exist by ca. 2080 under ongoing climate warming. The influence of high-order mechanics is evident, however, in glacier advance scenarios, where ice speeds are greater and ice dynamical effects become more important. Although similar studies on other glaciers are essential to generalize such findings, we advise that high-order mechanics are important and therefore should be considered while modeling the evolution of active glaciers. Reduced model predictions may be adequate for other glaciologic and topographic settings, particularly where flow speeds are low and where mass balance changes dominate over ice dynamics in determining glacier geometry.

Description: Meteorological drivers of ablation processes on a cold glacier in the semi-arid Andes of Chile The Cryosphere, 7, 1513-1526, 2013 Author(s): S. MacDonell, C. Kinnard, T. Mölg, L. Nicholson, and J. Abermann Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s −1 ), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m −2 ; summer hourly maximum 1354 W m −2 ) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Description: Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates The Cryosphere, 7, 1499-1512, 2013 Author(s): I. Sasgen, H. Konrad, E. R. Ivins, M. R. Van den Broeke, J. L. Bamber, Z. Martinec, and V. Klemann We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 ± 18 Gt yr −1 . Therefore, our AIS mass balance of −114 ± 23 Gt yr −1 is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (−26 ± 3 Gt yr −1 and −127 ± 7 Gt yr −1 , respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 ± 13 Gt yr −1 ), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 ± 7 Gt yr −1 (ca. 0.4 mm yr −1 ) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.

Description: Non-extensivity and long-range correlations in the earthquake activity at the West Corinth rift (Greece) Nonlinear Processes in Geophysics, 20, 713-724, 2013 Author(s): G. Michas, F. Vallianatos, and P. Sammonds In the present work the statistical properties of the earthquake activity in a highly seismic region, the West Corinth rift (Central Greece), are being studied by means of generalized statistical physics. By using a dataset that covers the period 2001–2008, we investigate the earthquake energy distribution and the distribution of the time intervals (interevent times) between the successive events. As has been reported previously, these distributions exhibit complex statistical properties and fractality. By using detrended fluctuation analysis (DFA), a well-established method for detection of long-range correlations in non-stationary signals, it is shown that long-range correlations are also present in the earthquake activity. The existence of these properties motivates us to use non-extensive statistical physics (NESP) to investigate the statistical properties of the frequency-magnitude and the interevent time distributions, along with other well-known relations in seismology, such as the gamma distribution for interevent times. The results of the analysis indicate that the statistical properties of the earthquake activity can be successfully reproduced by means of NESP and that the earthquake activity at the West Corinth rift is correlated at all-time scales.

Description: Oxygen minimum zone of the open Arabian Sea: variability of oxygen and nitrite from daily to decadal time scales Biogeosciences Discussions, 10, 15455-15517, 2013 Author(s): K. Banse, S. W. A. Naqvi, P. V. Narvekar, J. R. Postel, and D. A. Jayakumar The oxygen minimum zone (OMZ) of the Arabian Sea is the thickest of the three oceanic OMZs, which is of global biogeochemical significance because of denitrification in the upper part leading to N 2 and N 2 O production. The residence time of the OMZ water is believed to be less than a decade. The upper few hundred meters of this zone are nearly anoxic but non-sulfidic and still support animal (metazoan) pelagic life, possibly as a result of episodic injections of O 2 by physical processes. The very low O 2 values obtained with the new STOX sensor in the eastern tropical South Pacific probably also characterize the Arabian Sea OMZ, but there is no apparent reason as to why the temporal trends of the historic data should not hold. We report on discrete measurements of dissolved O 2 and NO 2 - , besides temperature and salinity, made between 1959 and 2004 well below the tops of the sharp pycno- and oxyclines near 150, 200, 300, 400, and 500 m depth. We assemble nearly all O 2 determinations (originally, 849 values, 695 in the OMZ) by the visual endpoint detection of the iodometric Winkler procedure, which in our data base yields about 0.04 mL L −1 (∼2 μM) O 2 above the endpoint from modern automated titration methods. We find 632 values acceptable (480 from 150 stations in the OMZ). The data are grouped in zonally-paired boxes of 1° lat. and 2° long. centered at 8°, 10°, 12°, 15°, 18°, 20°, and 21° N along 65° E and 67° E. The latitudes of 8–12° N, outside the OMZ, are only treated in passing. The principal results are as follows: (1) an O 2 climatology for the upper OMZ reveals a marked seasonality at 200 to 500 m depth with O 2 levels during the northeast monsoon and spring intermonsoon season elevated over those during the southwest monsoon season (median difference, 0.08 mL L −1 [3.5 μM]). The medians of the slopes of the seasonal regressions of O 2 on year for the NE and SW monsoon seasons are −0.0043 and −0.0019 mL L −1 a −1 , respectively (−0.19 and −0.08 μM a −1 ; n = 10 and 12, differing at p = 0.01); (2) four decades of statistically significant decreases of O 2 between 15° and 20° N but a trend to a similar increase near 21° N are observed. The balance of the mechanisms that more or less annually maintain the O 2 levels are still uncertain. At least between 300 and 500 m the annual reconstitution of the decrease is inferred to be due to lateral, isopycnal re-supply of O 2 , while at 200 (250?) m it is diapycnal, most likely by eddies. Similarly, recent models show large vertical advection of O 2 well below the pycno- cum -oxycline. The spatial (within drift stations) and temporal (daily) variability in hydrography and chemistry is large also below the principal pycnocline. The seasonal change of hydrography is considerable even at 500 m. There is no trend in the redox environment for a quarter of a century at a GEOSECS station near 20° N. In the entire OMZ the slopes on year within seasons for the quite variable NO 2 - (taken as an indicator of active denitrification) do not show a clear pattern. Also, future O 2 or nutrient budgets for the OMZ should not be based on single cruises or sections obtained during one season only. Steady state cannot be assumed any longer for the intermediate layers of the central Arabian Sea.

Description: Stable isotopes dissect food webs from top to the bottom Biogeosciences Discussions, 10, 14923-14952, 2013 Author(s): J. J. Middelburg Stable isotopes have been used extensively to study food web functioning, i.e. the flow of energy and matter among organisms. Traditional food-web studies are based on the natural variability of carbon and nitrogen isotopes and are limited to larger organisms that can be physically separated from their environment. Recent developments allow isotope ratio measurements of microbes and this in turn allows then measurement of entire food webs, i.e. from small producers at the bottom to large consumers at the top. Here, I provide a concise review on the use and potential of stable isotope to reconstruct end-to-end food webs. I will first discuss food web reconstruction based on natural abundances isotope data and will then show that the use of stable isotopes as deliberately added tracers provides complementary information. Finally, challenges and opportunities for end-to-end food web reconstructions in a changing world are discussed.

Description: The impact of initial spread calibration on the RELO ensemble and its application to Lagrangian dynamics Nonlinear Processes in Geophysics, 20, 621-641, 2013 Author(s): M. Wei, G. Jacobs, C. Rowley, C. N. Barron, P. Hogan, P. Spence, O. M. Smedstad, P. Martin, P. Muscarella, and E. Coelho A number of real-time ocean model forecasts were carried out successfully at Naval Research Laboratory (NRL) to provide modeling support and numerical guidance to the CARTHE GLAD at-sea experiment during summer 2012. Two RELO ensembles and three single models using NCOM and HYCOM with different resolutions were carried out. A calibrated ensemble system with enhanced spread and reliability was developed to better support this experiment. The calibrated ensemble is found to outperform the un-calibrated ensemble in forecasting accuracy, skill, and reliability for all the variables and observation spaces evaluated. The metrics used in this paper include RMS error, anomaly correlation, PECA, Brier score, spread reliability, and Talagrand rank histogram. It is also found that even the un-calibrated ensemble outperforms the single forecast from the model with the same resolution. The advantages of the ensembles are further extended to the Lagrangian framework. In contrast to a single model forecast, the RELO ensemble provides not only the most likely Lagrangian trajectory for a particle in the ocean, but also an uncertainty estimate that directly reflects the complicated ocean dynamics, which is valuable for decision makers. The examples show that the calibrated ensemble with more reliability can capture trajectories in different, even opposite, directions, which would be missed by the un-calibrated ensemble. The ensembles are applied to compute the repelling and attracting Lagrangian coherent structures (LCSs), and the uncertainties of the LCSs, which are hard to obtain from a single model forecast, are estimated. It is found that the spatial scales of the LCSs depend on the model resolution. The model with the highest resolution produces the finest, small-scale, LCS structures, while the model with lowest resolution generates only large-scale LCSs. The repelling and attracting LCSs are found to intersect at many locations and create complex mesoscale eddies. The fluid particles and drifters in the middle of these tangles are subject to attraction and repulsion simultaneously from these two kinds of LCSs. As a result, the movements of particles near the Deepwater Horizon (DWH) location are severely limited. This is also confirmed by the Lagrangian trajectories predicted by the ensembles.

Description: Multifractal properties of embedded convective structures in orographic precipitation: toward subgrid-scale predictability Nonlinear Processes in Geophysics, 20, 605-620, 2013 Author(s): M. Nogueira, A. P. Barros, and P. M. A. Miranda Rain and cloud fields produced by fully nonlinear idealized cloud resolving numerical simulations of orographic convective precipitation display statistical multiscaling behavior, implying that multifractal diagnostics should provide a physically robust basis for the downscaling and sub-grid scale parameterizations of moist processes. Our results show that the horizontal scaling exponent function (and respective multiscaling parameters) of the simulated rainfall and cloud fields varies with atmospheric and terrain properties, particularly small-scale terrain spectra, atmospheric stability, and advective timescale. This implies that multifractal diagnostics of moist processes for these simulations are fundamentally transient, exhibiting complex nonlinear behavior depending on atmospheric conditions and terrain forcing at each location. A particularly robust behavior found here is the transition of the multifractal parameters between stable and unstable cases, which has a clear physical correspondence to the transition from stratiform to organized (banded and cellular) convective regime. This result is reinforced by a similar behavior in the horizontal spectral exponent. Finally, our results indicate that although nonlinearly coupled fields (such as rain and clouds) have different scaling exponent functions, there are robust relationships with physical underpinnings between the scaling parameters that can be explored for hybrid dynamical-statistical downscaling.

Description: Responses of carbon dioxide flux and plant biomass to drought in a treed peatland in northern Alberta: a climate change perspective Biogeosciences Discussions, 10, 14999-15031, 2013 Author(s): T. M. Munir, B. Xu, M. Perkins, and M. Strack Northern peatland ecosystems represent large carbon stocks that are susceptible to changes such as accelerated mineralization due to water table lowering expected under a climate change scenario. During the growing seasons of 2011 and 2012 we monitored CO 2 fluxes and plant biomass along a microtopographic gradient (hummocks-hollows) in an undisturbed dry continental boreal treed bog (control) and a nearby site that was drained (drained) in 2001. Ten years of drainage in the bog significantly increased coverage of shrubs at hummocks and lichens at hollows. Considering measured hummock coverage and including tree incremental growth, we estimate that the control site was a larger sink in 2011 of −40 than that of −13 g C m −2 in 2012 while the drained site was a source of 144 and 140 g C m −2 over the same years. We infer that, drainage induced changes in vegetation growth led to increased biomass to counteract a portion of soil carbon losses. These results suggest that spatial variability (microtopography) and changes in vegetation community in boreal peatlands will affect how these ecosystems respond to lowered water table potentially induced by climate change.

Description: The coccolithophores Emiliania huxleyi and Coccolithus pelagicus : extant populations from the Norwegian-Iceland Sea and Fram Strait Biogeosciences Discussions, 10, 15077-15106, 2013 Author(s): C. V. Dylmer, J. Giraudeau, V. Hanquiez, and K. Husum Extant coccolithophores and their relation to the governing oceanographic features in the northern North Atlantic were investigated along two zonal transects of surface water sampling, both conducted during summer 2011 and fall 2007. The northern transects crossed Fram Strait and its two opposing boundary currents (West Spitsbergen Current and East Greenland Current), while the southern transects sampled the Norwegian and Iceland Seas (passing the island Jan Mayen) from the Lofoten Islands to the continental margin off Eastern Greenland. The distribution of the dominant coccolithophore species Emiliania huxleyi and Coccolithus pelagicus is discussed in view of both the surface hydrology at the time of sampling and the structure of the surface mixed layer. Remote-sensing images as well as CTD and ARGO profiles are used to constrain the physico-chemical state of the surface water at the time of sampling. Both transects were characterized by strong seasonal differences in bulk coccolithophore standing stocks with maximum values of 53 × 10 3 cells L −1 for the northern transect and 72 × 10 3 cells L −1 for the southern transect in fall and summer, respectively. The highest recorded bulk cell densities are essentially explained by E. huxleyi . This species shows a zonal shift in peak abundance in the Norwegian-Iceland Seas from a summer maximum in the Lofoten gyre to peak cell densities around the island Jan Mayen in fall. Vertical mixing of Atlantic waters west of Lofoten Island, a phenomenom related to pervasive summer large scale atmospheric changes in the eastern Nordic Seas, on one hand, and strengthened influence of melt-water and related surface water stratification around the island Jan Mayen during fall, on the other hand, explains the observed seasonal migration of the E. huxleyi peak production area, as well as the seasonal change in dominating species within the Iceland Sea. In addition our datasets are indicative of a well-defined maximum boundary temperature of 6 °C for the production of C. pelagicus in the northern North Atlantic. The Fram Strait transects provides, to our knowledge, a first view of the zonal distribution of extant coccolithophores in this remote setting during summer and fall. Our datasets are indicative of a seasonal change in the species community from an E. huxleyi -dominated assemblage during summer to a C. pelagicus -rich population during fall. Here, higher irradiance and increased Atlantic water influence during summer favored the production of the opportunistic species E. huxleyi close to the Arctic Front, whereas the peak production area during fall, with high concentrations of C. pelagicus , lays in true Arctic/Polar waters.

Description: Anthropogenic and natural methane fluxes in Switzerland synthesized within a spatially-explicit inventory Biogeosciences Discussions, 10, 15181-15224, 2013 Author(s): R. V. Hiller, D. Bretscher, T. DelSontro, T. Diem, W. Eugster, R. Henneberger, S. Hobi, E. Hodson, D. Imer, M. Kreuzer, T. Künzle, L. Merbold, P. A. Niklaus, B. Rihm, A. Schellenberger, M. H. Schroth, C. J. Schubert, H. Siegrist, J. Stieger, N. Buchmann, and D. Brunner We present the first high-resolution (500 m × 500 m) gridded methane (CH 4 ) emission inventory for Switzerland, which integrates the national emission totals reported to the United Nations Framework Convention on Climate Change (UNFCCC) and recent CH 4 flux studies conducted by research groups across Switzerland. In addition to anthropogenic emissions, we also include natural and semi-natural CH 4 fluxes, i.e., emissions from lakes and reservoirs, wetlands, wild animals as well as uptake by forest soils. National CH 4 emissions were disaggregated using detailed geostatistical information on source locations and their spatial extent and process- or area-specific emission factors. In Switzerland, the highest CH 4 emissions in 2011 originated from the agricultural sector (150 Gg CH 4 yr −1 ), mainly produced by ruminants and manure management, followed by emissions from waste management (15 Gg CH 4 yr −1 ) mainly from landfills and the energy sector (12 Gg CH 4 yr −1 ), which was dominated by emissions from natural gas distribution. Compared to the anthropogenic sources, emissions from natural and semi-natural sources were relatively small (6 Gg CH 4 yr −1 ), making up only 3 % of the total emissions in Switzerland. CH 4 fluxes from agricultural soils were estimated to be not significantly different from zero (between −1.5 and 0 Gg CH 4 yr −1 ), while forest soils are a CH 4 sink (approx. −2.8 Gg CH 4 yr −1 ), partially offsetting other natural emissions. Estimates of uncertainties are provided for the different sources, including an estimate of spatial disaggregation errors deduced from a comparison with a global (EDGAR v4.2) and a European CH 4 inventory (TNO/MACC). This new spatially-explicit emission inventory for Switzerland will provide valuable input for regional scale atmospheric modeling and inverse source estimation.

Description: Seasonal trends of dry and bulk concentration of nitrogen compounds over a rain forest in Ghana Biogeosciences Discussions, 10, 15225-15255, 2013 Author(s): F. Fattore, T. Bertolini, S. Materia, S. Gualdi, A. Thongo M'Bou, G. Nicolini, R. Valentini, A. De Grandcourt, D. Tedesco, and S. Castaldi African tropical forests of the equatorial belt might receive significant input of extra nitrogen derived from biomass burning occurring in the north savanna belt and transported equator wards by NE winds. In order to test this hypothesis an experiment was set up in a tropical rain forest in the National park of Ankasa (Ghana) aiming at: quantifying magnitude and seasonal variability of concentrations of N compounds, present as gas and aerosol (dry nitrogen) or in the rainfall (bulk nitrogen), over the studied forest; relating their seasonal variability to trends of local and regional winds and rainfall and to variations of fire events in the region. Three Delta systems, implemented for monthly measurements of NO 2 , were mounted over a tower at 45 m height, 20 m above forest canopy to sample gas (NH 3 , NO 2 , HNO 3 , HCl, SO 2 ) and aerosol (NH 4 + , NO 3 − , and several ions), together with three tanks for bulk rainfall collection (to analyze NH 4 + , NO 3 − and ion concentration). The tower was provided with a sonic anemometer to estimate local wind data. The experiment started in October 2011 and data up to October 2012 are presented. To interpret the observed seasonal trends of measured compounds, local and regional meteo data and regional satellite fire data were analyzed. The concentration of N compounds significantly increased from December to April, during the drier period, peaking in December-February when North Eastern winds (Harmattan) were moving dry air masses over the West central African region and the inter tropical convergence zone (ITCZ) was at its minimum latitude over the equator. This period also coincided with peaks of fire in the whole region. On the contrary, N concentration in gas, aerosol and rain decreased from May to October when prevalent winds arrived from the sea (South-East), during the Monsoon period. Both ionic compositions of rain and analysis of local wind direction showed a significant and continuous presence of see-breeze at site. The ionic composition of rain water resulted much closer to see water and poorer in N compounds from May to October.

Description: Inferences from CO 2 and CH 4 concentration profiles at the Zotino Tall Tower Observatory (ZOTTO) on local summer-time ecosystem fluxes Biogeosciences Discussions, 10, 15337-15372, 2013 Author(s): J. Winderlich, C. Gerbig, O. Kolle, and M. Heimann The Siberian region is still sparsely covered by ecosystem observatories, which motivates to exploit existing datasets to gain spatially and temporally better-resolved carbon fluxes. The Zotino Tall Tower Observatory (ZOTTO, 60°48' N, 89°21' E) observations of CO 2 and CH 4 mole fractions as well as meteorological parameters from six different heights up to 301 m allow for an additional estimate of surface-atmosphere fluxes of CO 2 and CH 4 for the Middle-Siberian region since 2009. The total carbon flux is calculated from the storage and the turbulent flux component. The gradients between the different tower levels determine the storage flux component, which dominates the local fluxes, especially during night. As a correction term, the turbulent flux component was estimated by the modified Bowen ratio method based on the sensible heat flux measurements at the top of the tower. The gained average night time fluxes (23:00 to 04:00 local time) are 2.7 ± 1.1 μmol (m 2 s) −1 for CO 2 and 5.6 ± 4.5 nmol (m 2 s) −1 for CH 4 during the summer months June-September in 2009 and 2011. During day, the method is limited due to numeric instabilities from vanishing vertical gradients; however, the derived CO 2 fluxes exhibit reasonable diurnal shape and magnitude compared to the eddy covariance technique, which become available at the site in 2012. Therefore, the tall tower data facilitates the extension of the new eddy covariance flux dataset back in time. The diurnal signal of the CH 4 flux is predominantly characterized by a strong morning transition, which is explained by local topographic effects.

Description: Live foraminiferal faunas (Rose Bengal stained) from the northern Arabian Sea: links with bottom-water oxygenation Biogeosciences Discussions, 10, 15257-15304, 2013 Author(s): C. Caulle, K. A. Koho, M. Mojtahid, G. J. Reichart, and F. J. Jorissen Live (Rose Bengal stained) benthic foraminifera from the Murray Ridge, within and below the northern Arabian Sea Oxygen Minimum Zone (OMZ), were studied in order to determine the relationship between faunal composition, bottom-water oxygenation (BWO), pore-water chemistry and organic matter (organic carbon and phytopigment) distribution. A series of multicores were recovered from a ten-station oxygen (BWO: 2–78 μM) and bathymetric (885–3010 m depth) transect during the winter monsoon in January 2009. Foraminifera were investigated from three different size fractions (63–125 μm, 125–150 μm and 〉 150 μm). The larger foraminifera (〉 125 μm) were strongly dominated by agglutinated species (e.g. Reophax spp.). In contrast, in the 63–125 μm fraction, calcareous taxa were more abundant, especially in the core of the OMZ, suggesting an opportunistic behaviour. On the basis of a Principal Component Analysis, three foraminiferal groups were identified, reflecting the environmental parameters along the study transect. The faunas from the shallowest stations, in the core of the OMZ (BWO: 2 μM), were composed of "low oxygen" species, typical of the Arabian Sea OMZ (e.g., Rotaliatinopsis semiinvoluta , Praeglobobulimina spp. , Bulimina exilis, Uvigerina peregrina type parva ). These taxa are adapted to the very low BWO conditions and to high phytodetritus supplies. The transitional group, typical for the lower part of the OMZ (BWO: 5–16 μM), is composed of more cosmopolitan taxa tolerant to low-oxygen concentrations ( Globocassidulina subglobosa , Ehrenbergina trigona ). Below the OMZ (BWO: 26–78 μM), where food availability is more limited and becomes increasingly restricted to surficial sediments, more cosmopolitan calcareous taxa were present, such as Bulimina aculeata, Melonis barleeanus, Uvigerina peregrina and Epistominella exigua . Miliolids were uniquely observed in this last group, reflecting the higher BWO. At these deeper sites, the faunas exhibit a clear depth succession of superficial, intermediate and deep-infaunal microhabitats, because of the deeper oxygen and nitrate penetration into the sediment.

Description: Evidence of meltwater retention within the Greenland ice sheet The Cryosphere, 7, 1433-1445, 2013 Author(s): A. K. Rennermalm, L. C. Smith, V. W. Chu, J. E. Box, R. R. Forster, M. R. Van den Broeke, D. Van As, and S. E. Moustafa Greenland ice sheet mass losses have increased in recent decades with more than half of these attributed to surface meltwater runoff. However, the magnitudes of englacial storage, firn retention, internal refreezing and other hydrologic processes that delay or reduce true water export to the global ocean remain less understood, partly due to a scarcity of in situ measurements. Here, ice sheet surface meltwater runoff and proglacial river discharge between 2008 and 2010 near Kangerlussuaq, southwestern Greenland were used to establish sub- and englacial meltwater storage for a small ice sheet watershed (36–64 km 2 ). This watershed lacks significant potential meltwater storage in firn, surface lakes on the ice sheet and in the proglacial area, and receives limited proglacial precipitation. Thus, ice sheet surface runoff not accounted for by river discharge can reasonably be attributed to retention in sub- and englacial storage. Evidence for meltwater storage within the ice sheet includes (1) characteristic dampened daily river discharge amplitudes relative to ice sheet runoff; (2) three cold-season river discharge anomalies at times with limited ice sheet surface melt, demonstrating that meltwater may be retained up to 1–6 months; (3) annual ice sheet watershed runoff is not balanced by river discharge, and while near water budget closure is possible as much as 54% of melting season ice sheet runoff may not escape to downstream rivers; (4) even the large meltwater retention estimate (54%) is equivalent to less than 1% of the ice sheet volume, which suggests that storage in en- and subglacial cavities and till is plausible. While this study is the first to provide evidence for meltwater retention and delayed release within the Greenland ice sheet, more information is needed to establish how widespread this is along the Greenland ice sheet perimeter.

Description: High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records The Cryosphere, 7, 1481-1498, 2013 Author(s): S. Kutuzov, M. Shahgedanova, V. Mikhalenko, P. Ginot, I. Lavrentiev, and S. Kemp The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

Description: Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls The Cryosphere, 7, 1455-1471, 2013 Author(s): A. Rabatel, A. Letréguilly, J.-P. Dedieu, and N. Eckert We present time series of equilibrium-line altitude (ELA) measured from the end-of-summer snow line altitude computed using satellite images, for 43 glaciers in the western Alps over the 1984–2010 period. More than 120 satellite images acquired by Landsat, SPOT and ASTER were used. In parallel, changes in climate variables, summer cumulative positive degree days (CPDD) and winter precipitation, were analyzed over the same time period using 22 weather stations located inside and around the study area. Assuming a continuous linear trend over the study period: (1) the average ELA of the 43 glaciers increased by about 170 m; (2) summer CPDD increased by about 150 PDD at 3000 m a.s.l.; and (3) winter precipitation remained rather stationary. Summer CPDD showed homogeneous spatial and temporal variability; winter precipitation showed homogeneous temporal variability, but some stations showed a slightly different spatial pattern. Regarding ELAs, temporal variability between the 43 glaciers was also homogeneous, but spatially, glaciers in the southern part of the study area differed from glaciers in the northern part, mainly due to a different precipitation pattern. A sensitivity analysis of the ELAs to climate and morpho-topographic variables (elevation, aspect, latitude) highlighted the following: (1) the average ELA over the study period of each glacier is strongly controlled by morpho-topographic variables; and (2) the interannual variability of the ELA is strongly controlled by climate variables, with the observed increasing trend mainly driven by increasing temperatures, even if significant nonlinear, low-frequency fluctuations appear to be driven by winter precipitation anomalies. Finally, we used an expansion of Lliboutry's approach to reconstruct fluctuations in the ELA of any glacier of the study area with respect to morpho-topographic and climate variables, by quantifying their respective weight and the related uncertainties in a consistent manner within a hierarchical Bayesian framework. This method was tested and validated using the ELA measured on the satellite images.

Description: A top-down model to generate ensembles of runoff from a large number of hillslopes Nonlinear Processes in Geophysics, 20, 683-704, 2013 Author(s): P. R. Furey, V. K. Gupta, and B. M. Troutman We hypothesize that total hillslope water loss for a rainfall–runoff event is inversely related to a function of a lognormal random variable, based on basin- and point-scale observations taken from the 21 km 2 Goodwin Creek Experimental Watershed (GCEW) in Mississippi, USA. A top-down approach is used to develop a new runoff generation model both to test our physical-statistical hypothesis and to provide a method of generating ensembles of runoff from a large number of hillslopes in a basin. The model is based on the assumption that the probability distributions of a runoff/loss ratio have a space–time rescaling property. We test this assumption using streamflow and rainfall data from GCEW. For over 100 rainfall–runoff events, we find that the spatial probability distributions of a runoff/loss ratio can be rescaled to a new distribution that is common to all events. We interpret random within-event differences in runoff/loss ratios in the model to arise from soil moisture spatial variability. Observations of water loss during events in GCEW support this interpretation. Our model preserves water balance in a mean statistical sense and supports our hypothesis. As an example, we use the model to generate ensembles of runoff at a large number of hillslopes for a rainfall–runoff event in GCEW.

Description: A mechanism for catastrophic filter divergence in data assimilation for sparse observation networks Nonlinear Processes in Geophysics, 20, 705-712, 2013 Author(s): G. A. Gottwald and A. J. Majda We study catastrophic filter divergence in data assimilation procedures whereby the forecast model develops severe numerical instabilities leading to a blow-up of the solution. Catastrophic filter divergence can occur in sparse observational grids with small observational noise for intermediate observation intervals and finite ensemble sizes. Using a minimal five-dimensional model, we establish that catastrophic filter divergence is a numerical instability of the underlying forecast model caused by the filtering procedure producing analyses which are not consistent with the true dynamics, and stiffness caused by the fast attraction of the inconsistent analyses towards the attractor during the forecast step.

Description: Landscape-scale changes in forest structure and functional traits along an Andes-to-Amazon elevation gradient Biogeosciences Discussions, 10, 15415-15454, 2013 Author(s): G. P. Asner, C. Anderson, R. E. Martin, D. E. Knapp, R. Tupayachi, T. Kennedy-Bowdoin, F. Sinca, and Y. Malhi Elevation gradients provide opportunities to explore environmental controls on forest structure and functioning, but plot-based studies have proven highly variable due to limited geographic scope. We used airborne imaging spectroscopy and LiDAR (light detection and ranging) to quantify changes in three-dimensional forest structure and canopy functional traits in a series of 25 ha landscapes distributed along a 3300 m elevation gradient from lowland Amazonia to treeline in the Peruvian Andes. Canopy greenness, photosynthetic fractional cover and exposed non-photosynthetic vegetation varied as much across lowland forests (100–200 m) as they did from the lowlands to the Andean treeline (3400 m). Elevation was positively correlated with canopy gap density and understory vegetation cover, and negatively related to canopy height and vertical profile. Increases in gap density were tightly linked to increases in understory plant cover, and larger gaps (20–200 m 2 produced 25–30 times the response in understory cover than did smaller gaps ( 〈 5 m 2 . Scaling of gap size to gap frequency was, however, relatively constant along the elevation gradient, which when combined with other canopy structural information, indicates equilibrium turnover patterns from the lowlands to treeline. Our results provide a first landscape-scale quantification of forest structure and canopy functional traits with changing elevation, thereby improving our understanding of disturbance, demography and ecosystem processes in the Andes-to-Amazon corridor.

Description: A note on the water budget of temperate glaciers The Cryosphere, 7, 1557-1564, 2013 Author(s): J. Oerlemans In this note, the total dissipative melting in temperate glaciers is studied. The analysis is based on the notion that the dissipation is determined by the loss of potential energy due to the downward motion of mass (ice, snow, meltwater and rain). A mathematical formulation of the dissipation is developed and applied to a simple glacier geometry. In the next step, meltwater production resulting from enhanced ice motion during a glacier surge is calculated. The amount of melt energy available follows directly from the lowering of the centre of gravity of the glacier. To illustrate the concept, schematic calculations are presented for a number of glaciers with different geometric characteristics. Typical dissipative melt rates, expressed as water-layer depth averaged over the glacier, range from a few centimetres per year for smaller glaciers to half a metre per year for Franz Josef Glacier, one of the most active glaciers in the world (in terms of mass turnover). The total generation of meltwater during a surge is typically half a metre. For Variegated Glacier a value of 70 cm is found, for Kongsvegen 20 cm. These values refer to water layer depth averaged over the entire glacier. The melt \textit{rate} depends on the duration of the surge. It is generally an order of magnitude greater than water production by `normal' dissipation. On the other hand, the additional basal melt rate during a surge is comparable in magnitude with the water input from meltwater and precipitation. This suggests that enhanced melting during a surge does not grossly change the total water budget of a glacier. Basal water generated by enhanced sliding is an important ingredient in many theories of glacier surges. It provides a positive feedback mechanism that actually makes the surge happen. The results found here suggest that this can only work if water generated by enhanced sliding accumulates in a part of the glacier base where surface meltwater and rain have no or very limited access. This finding seems compatible with the fact that, on many glaciers, surges are initiated in the lower accumulation zone.

Description: Springtime phytoplankton dynamics in the Arctic Krossfjorden and Kongsfjorden (Spitsbergen) as a function of glacier proximity Biogeosciences Discussions, 10, 15519-15557, 2013 Author(s): A. M.-T. Piquet, W. H. van de Poll, R. J. W. Visser, C. Wiencke, H. Bolhuis, and A. G. J. Buma The hydrographic properties of the Kongsfjorden – Krossfjorden system (79° N, Spitsbergen) are affected by Atlantic water incursions as well as glacier meltwater runoff. This results in strong physical gradients (temperature, salinity and irradiance) within the fjords. Here, we tested the hypothesis that glaciers affect phytoplankton dynamics as early as the productive spring bloom period. During two campaigns in 2007 (late spring) and 2008 (early spring) we studied hydrographic characteristics and phytoplankton variability along 2 transects in both fjords, using HPLC-CHEMTAX pigment fingerprinting, molecular fingerprinting (DGGE) and sequencing of 18S rRNA genes. The sheltered inner fjord locations remained colder during spring as opposed to the outer locations. Vertical light attenuation coefficients increased from early spring onwards, at all locations, but in particular at the inner locations. During the end of spring, meltwater input had stratified surface waters throughout the fjords. The inner fjord locations were characterized by overall lower phytoplankton biomass. Furthermore HPLC-CHEMTAX data revealed that diatoms and Phaeocystis sp. were replaced by small nano- and picophytoplankton during late spring, coinciding with low nutrient availability. The innermost stations showed higher relative abundances of nano- and picophytoplankton throughout, notably of cyanophytes and cryptophytes. Molecular fingerprinting revealed a high similarity between inner fjord samples from early spring and late spring samples from all locations, while outer samples from early spring clustered separately. We conclude that glacier influence, mediated by early meltwater input, modifies phytoplankton biomass and composition already during the spring bloom period, in favor of low biomass and small cell size communities. This may affect higher trophic levels especially when regional warming further increases the period and volume of meltwater.

Description: Simulating microbial degradation of organic matter in a simple porous system using the 3-D diffusion based model MOSAIC Biogeosciences Discussions, 10, 15613-15640, 2013 Author(s): O. Monga, P. Garnier, V. Pot, E. Coucheney, N. Nunan, W. Otten, and C. Chenu This paper deals with the simulation of microbial degradation in soil within pore space at microscopic scale. Pore space was described using sphere network coming from a geometrical modeling algorithm. The biological model was improved regarding previous work in order to include transformation of dissolved organic compounds and diffusion processes. Our model was tested using experimental results of a simple substrate decomposition (Fructose) within a simple media (the sand). Diverse microbial communities were inoculated. Separated incubations in microcosms were carried out using 5 different bacterial communities at 2 different water potentials of −10 cm and −100 cm of water. We calibrated the biological parameters by means of experimental data obtained at high water content and we tested the model without any parameters change at low water content. Same as for experimental data, our simulation results showed the decrease in water content involved the decrease of mineralisation. The model was able to simulate the decrease of connectivity between substrate and microorganism due the decrease of water content.

Description: Influence of temperature and CO 2 on the strontium and magnesium composition of coccolithophore calcite Biogeosciences Discussions, 10, 15559-15586, 2013 Author(s): M. N. Müller, M. Lebrato, U. Riebesell, J. Barcelos e Ramos, K. G. Schulz, S. Blanco-Ameijeiras, S. Sett, A. Eisenhauer, and H. M. Stoll Marine calcareous sediments provide a fundamental basis for paleoceanographic studies aiming to reconstruct past oceanic conditions and understand key biogeochemical element cycles. Calcifying unicellular phytoplankton (coccolithophores) are a major contributor to both carbon and calcium cycling by photosynthesis and the production of calcite (coccoliths) in the euphotic zone and the subsequent long-term deposition and burial into marine sediments. Here we present data from controlled laboratory experiments on four coccolithophore species and elucidate the relation between the divalent cation (Sr, Mg and Ca) partitioning in coccoliths and cellular physiology (growth, calcification and photosynthesis). Coccolithophores were cultured under different seawater temperature and carbonate chemistry conditions. The partition coefficient of strontium ( D Sr ) was positively correlated with both carbon dioxide ( p CO 2 ) and temperature but displayed no coherent relation to particulate organic and inorganic carbon production rates. Furthermore, D Sr correlated positively with cellular growth rates when driven by temperature but no correlation was present when changes in growth rates were p CO 2 -induced. The results demonstrate the complex interaction between environmental forcing and physiological control on the strontium partitioning in coccolithophore calcite. The partition coefficient of magnesium ( D Mg ) displayed species-specific differences and elevated values under nutrient limitation. No conclusive correlation between coccolith D Mg and temperature was observed but p CO 2 induced a rising trend in coccolith D Mg . Interestingly, the best correlation was found between coccolith D Mg and chlorophyll a production suggesting that chlorophyll a and calcite associated Mg originate from the same intracellular pool. These results give an extended insight into the driving factors that lead to variations in the coccolith Mg / Ca ratio and can be used for Sr / Ca and Mg / Ca paleoproxy calibration.

Description: Forest NEP is significantly driven by previous year's weather Biogeosciences Discussions, 10, 15587-15611, 2013 Author(s): S. Zielis, S. Etzold, R. Zweifel, W. Eugster, M. Haeni, and N. Buchmann Understanding the response of forest net ecosystem productivity (NEP) to environmental drivers under climate change is highly relevant for predictions of annual forest carbon (C) flux budgets. Modeling annual forest NEP with soil–vegetation–atmosphere transfer models (SVATs), however, remains challenging due to unknown responses of forests to weather of the previous year. In this study, we addressed the influence of previous year's weather on the inter-annual variability of NEP for a subalpine spruce forest in Switzerland. Analysis of long-term (1997–2011) eddy covariance measurements showed that the Norway spruce forest Davos Seehornwald was a consistent sink for atmospheric CO 2 , sequestering 210 ± 88 g C m −2 per year on average. Previous year's weather strongly affected inter-annual variability of NEP, increasing the explained variance in linear models to 53% compared to 20% without previous year's weather. Thus, our results highlight the need to consider previous year's weather in modeling annual C budgets of forests. Furthermore, soil temperature in the current year's spring played a major role controlling annual NEP, mainly by influencing gross primary productivity early in the year, with spring NEP accounting for 56% of annual NEP. Consequently, we expect an increase in net CO 2 uptake with future climate warming, as long as no other resources become limiting.

Description: Synoptic evaluation of carbon cycling in Beaufort Sea during summer: contrasting river inputs, ecosystem metabolism and air–sea CO 2 fluxes Biogeosciences Discussions, 10, 15641-15710, 2013 Author(s): A. Forest, P. Coupel, B. Else, S. Nahavandian, B. Lansard, P. Raimbault, T. Papakyriakou, Y. Gratton, L. Fortier, J.-É. Tremblay, and M. Babin The accelerated decline in Arctic sea ice combined with an ongoing trend toward a more dynamic atmosphere is modifying carbon cycling in the Arctic Ocean. A critical issue is to understand how net community production (NCP; the balance between gross primary production and community respiration) responds to changes and modulates air–sea CO 2 fluxes. Using data collected as part of the ArcticNet-Malina 2009 expedition in southeastern Beaufort Sea (Arctic Ocean), we synthesize information on sea ice, wind, river, water column properties, metabolism of the planktonic food web, organic carbon fluxes and pools, as well as air–sea CO 2 exchange, with the aim of identifying indices of ecosystem response to environmental changes. Data were analyzed to develop a non-steady-state carbon budget and an assessment of NCP against air–sea CO 2 fluxes. The mean atmospheric forcing was a mild upwelling-favorable wind (~5 km h −1 ) blowing from the N-E and a decaying ice cover ( 600 mg C m −2 d −1 ) over the shelf prior to our survey, (2) freshwater dilution by river runoff and ice melt, and (3) the presence of cold surface waters offshore. Only the Mackenzie River delta and localized shelf areas directly affected by upwelling were identified as substantial sources of CO 2 to the atmosphere (〉10mmol C m −2 d −1 ). Although generally

Description: Global glacier changes: a revised assessment of committed mass losses and sampling uncertainties The Cryosphere, 7, 1565-1577, 2013 Author(s): S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp Most glaciers and ice caps (GIC) are out of balance with the current climate. To return to equilibrium, GIC must thin and retreat, losing additional mass and raising sea level. Because glacier observations are sparse and geographically biased, there is an undersampling problem common to all global assessments. Here, we further develop an assessment approach based on accumulation-area ratios (AAR) to estimate committed mass losses and analyze the undersampling problem. We compiled all available AAR observations for 144 GIC from 1971 to 2010, and found that most glaciers and ice caps are farther from balance than previously believed. Accounting for regional and global undersampling errors, our model suggests that GIC are committed to additional losses of 32 ± 12% of their area and 38 ± 16% of their volume if the future climate resembles the climate of the past decade. These losses imply global mean sea-level rise of 163 ± 69 mm, assuming total glacier volume of 430 mm sea-level equivalent. To reduce the large uncertainties in these projections, more long-term glacier measurements are needed in poorly sampled regions.

Description: Sub-inertial modulation of nonlinear Kelvin waves in the coastal zone Nonlinear Processes in Geophysics, 20, 357-364, 2013 Author(s): D. V. Stepanov and V. V. Novotryasov Observational evidence is presented for interaction between nonlinear internal Kelvin waves at the ω t,i (where the ω t is the semidiurnal frequency and the ω i is the inertial frequency) and random oscillations of the background coastal current at the sub-inertial Ω frequency in the Japan/East Sea. Enhanced coastal currents at the sum ω + and difference ω-frequencies ω ± =ω t,i ± Ω have properties of propagating Kelvin waves, which suggests permanent energy exchange from the sub-inertial band to the mesoscale ω ± band. This interaction may be responsible for a greater-than-predicted intensification, steepening and breaking of boundary-trapped Kelvin waves. The problem of interaction between the nonlinear Kelvin wave at the frequency ω and the low-frequency narrowband noise with representative frequency Ω≪ω is investigated using the theory of nonlinear weak dispersion waves.

Description: Foraminiferal survival after long term experimentally induced anoxia Biogeosciences Discussions, 10, 9243-9284, 2013 Author(s): D. Langlet, E. Geslin, C. Baal, E. Metzger, F. Lejzerowicz, B. Riedel, M. Zuschin, J. Pawlowski, M. Stachowitsch, and F. J. Jorissen Anoxia has been successfully induced in four benthic chambers installed on the Northern Adriatic seafloor from 1 week to 10 months. To accurately determine whether benthic foraminifera can survive experimentally induced prolonged anoxia, the CellTrackerGreen method has been applied. Numerous individuals have been found living at all sampling times and at all sampling depths, showing that benthic foraminifera can survive up to 10 months of anoxia with co-occurring hydrogen sulphides. However, foraminiferal standing stocks decrease with sampling time in an irregular way. A large difference in standing stock between two cores samples in initial conditions indicates the presence of a large spatial heterogeneity of the foraminiferal faunas. An unexpected increase in standing stocks after 1 month is tentatively interpreted as a reaction to increased food availability due to the massive mortality of infaunal macrofaunal organisms. After this, standing stocks decrease again in a core sampled after 2 months of anoxia, to attain a minimum in the cores sampled after 10 months. We speculate that the trend of overall decrease of standing stocks is not due to the adverse effects of anoxia and hydrogen sulphides, but rather due to a continuous diminution of labile organic matter.

Description: Environmental controls on the Emiliania huxleyi calcite mass Biogeosciences Discussions, 10, 9285-9313, 2013 Author(s): M. T. Horigome, P. Ziveri, M. Grelaud, K.-H. Baumann, G. Marino, and P. G. Mortyn Although ocean acidification is expected to impact (bio)calcification by decreasing the seawater carbonate ion concentration, [CO 3 2− ], there exists evidence of non-uniform response of marine calcifying plankton to low seawater [CO 3 2− ]. This raises questions on the role of environmental factors other than acidification and on the complex physiological responses behind calcification. Here we investigate the synergistic effect of multiple environmental parameters, including temperature, nutrient (nitrate and phosphate) availability, and seawater carbonate chemistry on the coccolith calcite mass of the cosmopolitan coccolithophore Emiliania huxleyi , the most abundant species in the world ocean. We use a suite of surface (late Holocene) sediment samples from the South Atlantic and southwestern Indian Ocean taken from depths lying well above the modern lysocline. The coccolith calcite mass in our results presents a latitudinal distribution pattern that mimics the main oceanographic features, thereby pointing to the potential importance of phosphorus and temperature in determining coccolith mass by affecting primary calcification and possibly driving the E. huxleyi morphotype distribution. This evidence does not necessarily argue against the potentially important role of the rapidly changing seawater carbonate chemistry in the future, when unabated fossil fuel burning will likely perturb ocean chemistry beyond a critical point. Rather our study highlights the importance of evaluating the combined effect of several environmental stressors on calcifying organisms to project their physiological response(s) in a high CO 2 world and improve interpretation of paleorecords.

Description: UV-induced carbon monoxide emission from living vegetation Biogeosciences Discussions, 10, 9373-9388, 2013 Author(s): D. Bruhn, K. R. Albert, T. N. Mikkelsen, and P. Ambus The global burden of carbon monoxide (CO) is rather uncertain. In this paper we address the potential for UV-induced CO emission by living terrestrial vegetation surfaces. Real-time measurements of CO concentrations were made with a cavity enhanced laser spectrometer connected in closed loop to either an ecosystem chamber or a plant-leaf scale chamber. Leaves of all examined plant species exhibited emission of CO in response to artificial UV-radiation as well as the UV-component of natural solar radiation. The UV-induced rate of CO emission exhibited a rather low dependence on temperature, indicating an abiotic process. The emission of CO in response to the UV-component of natural solar radiation was also evident at the ecosystem scale.

Description: Nutrient dynamics in tropical rivers, estuarine-lagoons, and coastal ecosystems along the eastern Hainan Island Biogeosciences Discussions, 10, 9091-9147, 2013 Author(s): R. H. Li, S. M. Liu, Y. W. Li, G. L. Zhang, J. L. Ren, and J. Zhang Nutrient dynamics were studied along the eastern Hainan Island based on field observations during 2006–2009, to understand nutrient biogeochemical processes and to have an overview of human perturbations on coastal ecosystems in this tropical region. The concentrations of nutrients in the rivers had seasonal variations enriched with dissolved inorganic nitrogen (DIN). High riverine concentrations of nitrate were mainly originated from agricultural fertilizer input. The ratios of DIN : PO 4 3− ranged from 37 to 1063, suggesting preferential PO 4 3− relative to nitrogen in the rivers. The areal yields of dissolved silicate (DSi) varied from 76 to 448 × 10 3 mol km −2 yr −1 due to erosion over the drainage area, inducing high levels of DSi among worldwide tropical systems. Aquaculture ponds contained high concentrations of NH 4 + (up to 157 μM) and DON (up to 130 μM). Particulate phosphorus concentrations (0.5 ∼1.4 μM) were in lower level comparied with estuaries around the world. Particulate silicate levels in rivers and lagoons were lower than global average level. Nutrient biogeochemistry in coastal areas were affected by human activities (e.g. aquaculture, agriculture), as well as natural events such as typhoon. Nutrient concentrations were low because open sea water dispersed land-derived nutrients. Nutrient budgets were built based on a steady-state box model, which showed that riverine fluxes would be magnified by estuarine processes (e.g. regeneration, desorption) in the Wenchanghe/Wenjiaohe Estuary, Wanquan River estuary, and the Laoyehai Lagoon except in the Xiaohai Lagoon. Riverine and groundwater input were the major sources of nutrients to the Xiaohai Lagoon and the Laiyehai Lagoon, respectively. Riverine input and aquaculture effluent were the major sources of nutrients to the eastern coastal of Hainan Island. Nutrient inputs to the coastal ecosystem can be increased by typhoon-induced runoff of rainwater, and phytoplankton bloom in the sea would be caused.

Description: Recent changes in spring snowmelt timing in the Yukon River basin detected by passive microwave satellite data The Cryosphere, 7, 905-916, 2013 Author(s): K. A. Semmens and J. M. Ramage Spring melt is a significant feature of high latitude snowmelt dominated drainage basins influencing hydrological and ecological processes such as snowmelt runoff and green-up. Melt duration, defined as the transition period from snowmelt onset until the end of the melt refreeze, is characterized by high diurnal amplitude variations (DAV) where the snowpack is melting during the day and refreezing at night, after which the snowpack melts constantly until depletion. Determining trends for this critical period is necessary for understanding how the Arctic is changing with rising temperatures and provides a baseline from which to assess future change. To study this dynamic period, brightness temperature ( T b ) data from the Special Sensor Microwave Imager (SSM/I) 37 V-GHz frequency from 1988 to 2010 were used to assess snowmelt timing trends for the Yukon River basin, Alaska/Canada. Annual T b and DAV for 1434 Equal-Area Scalable Earth (EASE)-Grid pixels (25 km resolution) were processed to determine melt onset and melt refreeze dates from T b and DAV thresholds previously established in the region. Temporal and spatial trends in the timing of melt onset and melt refreeze, and the duration of melt were analyzed for the 13 sub-basins of the Yukon River basin with three different time interval approaches. Results show a lengthening of the melt period for the majority of the sub-basins with a significant trend toward later end of melt refreeze after which the snowpack melts day and night leading to snow clearance, peak discharge, and green-up. Earlier melt onset trends were also found in the higher elevations and northernmost sub-basins (Porcupine, Chandalar, and Koyukuk rivers). Latitude and elevation displayed the dominant controls on melt timing variability and spring solar flux was highly correlated with melt timing in middle (∼600–1600 m) elevations.

Description: Southern Hemisphere imprint for Indo–Asian summer monsoons during the last glacial period as revealed by Arabian Sea productivity records Biogeosciences Discussions, 10, 9315-9343, 2013 Author(s): T. Caley, S. Zaragosi, J. Bourget, P. Martinez, B. Malaizé, F. Eynaud, L. Rossignol, T. Garlan, and N. Ellouz-Zimmermann The monsoon is one of the most important climatic phenomena: it promotes inter-hemispheric exchange of energy and affects the economical prosperity of several countries exposed to its seasonal seesaw. Previous studies in both the Indian and Asian monsoon systems have suggested a dominant north hemispheric (NH) control on summer monsoon dynamics at the scale of suborbital-millennial climatic changes, while the forcing/response of Indian and Asian monsoons at the orbital scale remains a matter of debate. Here nine marine sediment cores distributed across the whole Arabian Sea are used to build a regional surface marine productivity signal. The productivity signal is driven by the intensity of Indian summer monsoon winds. Results demonstrate the existence of an imprint of suborbital Southern Hemisphere (SH) temperature changes (i.e., Antarctica) on the Indian summer monsoon during the last glacial period, challenging the traditional and exclusive NH forcing hypothesis. Meanwhile, during the last deglaciation, the NH plays a more significant role. The δ 18 O signal recorded in the Asian monsoon speleothem records could be exported by winds from the Indian summer monsoon region, as recently proposed in modelling exercise, explaining the SH signature observed in Asian cave speleothems. Contrary to the view of a passive response of Indian and Asian monsoons to NH anomalies, the present results strongly suggest that the Indo–Asian summer monsoon plays an active role in amplifying millennial inter-hemispheric asymmetric patterns. Additionally, this study helps to decipher the observed differences between Indian and Asian-speleothem monsoonal records at the orbital-precession scale.

Description: Seasonal controls on snow distribution and aerial ablation at the snow-patch and landscape scales, McMurdo Dry Valleys, Antarctica The Cryosphere, 7, 917-931, 2013 Author(s): J. W. Eveland, M. N. Gooseff, D. J. Lampkin, J. E. Barrett, and C. D. Takacs-Vesbach Accumulated snow in the McMurdo Dry Valleys, while limited, has great ecological significance to subnivian soil environments. Though sublimation dominates the ablation process in this region, measurable increases in soil moisture and insulation from temperature extremes provide more favorable conditions with respect to subnivian soil communities. While precipitation is not substantial, significant amounts of snow can accumulate, via wind transport, in topographic lees along the valley bottoms, forming thousands of discontinuous snow patches. These patches have the potential to act as significant sources of local meltwater, controlling biogeochemical cycling and the landscape distribution of microbial communities. Therefore, determining the spatial and temporal dynamics of snow at multiple scales is imperative to understanding the broader ecological role of snow in this region. High-resolution satellite imagery acquired during the 2009–2010 and 2010–2011 austral summers was used to quantify the distribution of snow across Taylor and Wright valleys. Extracted snow-covered area from the imagery was used as the basis for assessing inter-annual variability and seasonal controls on accumulation and ablation of snow at multiple scales. In addition to landscape analyses, fifteen 1 km 2 plots (3 in each of 5 study regions) were selected to assess the prevalence of snow cover at finer spatial scales, referred to herein as the snow-patch scale. Results confirm that snow patches tend to form in the same locations each year with some minor deviations observed. At the snow-patch scale, neighboring patches often exhibit considerable differences in aerial ablation rates, and particular snow patches do not reflect trends for snow-covered area observed at the landscape scale. These differences are presumably related to microtopographic influences acting on individual snow patches, such as wind sheltering and differences in snow depth due to the underlying topography. This highlights the importance of both the landscape and snow-patch scales in assessing the effects of snow cover on biogeochemical cycling and microbial communities.

Description: Heterogeneity of impacts of high CO 2 on the North Western European Shelf Biogeosciences Discussions, 10, 9389-9413, 2013 Author(s): Y. Artioli, J. C. Blackford, G. Nondal, R. G. J. Bellerby, S. L. Wakelin, J.T. Holt, M. Butenschön, and J. I. Allen The increase in atmospheric CO 2 is a dual threat to the marine environment: from one side it drives climate change leading to changes in water temperature, circulation patterns and stratification intensity; on the other side it causes a decrease in pH (Ocean Acidification or OA) due to the increase in dissolved CO 2 . Assessing the combined impact of climate change and OA on marine ecosystems is a challenging task: the response of the ecosystem to a single driver is highly variable and still uncertain, as well as the interaction between these that could be either synergistic or antagonistic. In this work we use the coupled oceanographic-ecosystem model POLCOMS-ERSEM driven by climate forcing to study the interaction between climate change and OA. We focus in particular on primary production and nitrogen speciation. The model has been run in three different configurations in order to separate the impacts of ocean acidification from those due to climate change. The model shows significant interaction among the drivers and high variability in the spatial response of the ecosystem. Impacts of climate change and of OA on primary production have similar magnitude, compensating in some area and exacerbating in others. On the contrary, the direct impact of OA on nitrification is much lower than the one imposed by climate change.

Description: On the role of mesoscale eddies for the biological productivity and biogeochemistry in the eastern tropical Pacific Ocean off Peru Biogeosciences Discussions, 10, 9179-9211, 2013 Author(s): L. Stramma, H. W. Bange, R. Czeschel, A. Lorenzo, and M. Frank Mesoscale eddies seem to play an important role for both the hydrography and biogeochemistry of the eastern tropical Pacific Ocean (ETSP) off Peru. However, detailed surveys of these eddies are not available, which has so far hampered an in depth understanding of their implications for nutrient distribution and biological productivity. In this study three eddies along a section at 16°45' S have been surveyed intensively during R/V Meteor cruise M90 in November 2012. A coastal mode water eddy, an open ocean mode water eddy and an open ocean cyclonic eddy have been identified and sampled in order to determine both their hydrographic properties and their influence on the biogeochemical setting of the ETSP. In the thermocline the temperature of the coastal anticyclonic eddy was up to 2 °C warmer, 0.2 more saline and the swirl velocity was up to 35 cm s –1 . The observed temperature and salinity anomalies, as well as swirl velocities of both types of eddies were about twice as large as had been described for the mean eddies in the ETSP and the observed heat and salt anomalies (AHA, ASA) show a much larger variability than the mean AHA and ASA. We found that the eddies contributed significantly to productivity by maintaining pronounced subsurface maxima of chlorophyll. Based on a comparison of the coastal (young) mode water eddy and the open ocean (old) mode water eddy we conclude that the aging of eddies when they detach from the coast and move westward to the open ocean considerably influences the eddies' properties: chlorophyll maxima are weaker and nutrients are subducted. The coastal mode water eddy was found to be a hotspot of nitrogen loss in the OMZ, whereas, the open ocean cyclonic eddy was of negligible importance for nitrogen loss. Our results show that the important role the eddies play in the ETSP can only be fully deciphered and understood through dedicated high spatial and temporal resolution oceanographic/biogeochemical surveys.

Description: Soil moisture modifies the response of soil respiration to temperature in a desert shrub ecosystem Biogeosciences Discussions, 10, 9213-9242, 2013 Author(s): B. Wang, T. S. Zha, X. Jia, B. Wu, Y. Q. Zhang, and S. G. Qin The response of soil respiration (Rs) to soil temperature and moisture have been well documented in forests, but data and information from desert shrub ecosystems are limited. Soil CO 2 efflux from a desert shrub ecosystem was measured continuously with automated chambers in Ningxia, northwest China, from June to October 2012. The responses of Rs to Ts was strongly affected diurnally by soil moisture, with the diel variation in Rs being strongly related to 10 cm soil temperature (Ts) at moderate and high soil volumetric water content (VWC), but less related to Ts at low VWC. Ts typically lagged Rs by 3–4 h, however, the lag time varied in relation to VWC, with increased lag times at low VWC. Over the seasonal cycle, daily mean Rs was positively correlated with Ts when VWC exceeded 0.08 m 3 m −3 , but became decoupled from Ts when VWC dropped below this threshold. The annual temperature sensitivity of Rs ( Q 10 ) was 1.5. The short-term sensitivity of Rs to Ts, computed using three-day windows, varied significantly over the seasonal cycle; the short-term Q 10 was negatively correlated with Ts and positively correlated with VWC. These results suggest the potential for a negative feedback to climate warming in desert ecosystems, related to the impact of low soil moisture on Rs. The results highlight the biological causes of diel hysteresis between Rs and Ts and the need for carbon cycle models to account for the interacting effects of Ts and VWC as joint determinants of Rs in desert ecosystem.

Description: Do successive climate extremes weaken the resistance of plant communities? An experimental study using plant assemblages Biogeosciences Discussions, 10, 9149-9177, 2013 Author(s): F. E. Dreesen, H. J. De Boeck, I. A. Janssens, and I. Nijs The probability that plant communities undergo successive climate extremes increases under climate change. Exposure to an extreme event might elicit acclimatory responses and thereby greater resistance to a subsequent event, but might also reduce resistance if the recovery period is too short or resilience too low. Using experimental plant assemblages, we compared the effects of two successive extremes (either two drought extremes, two heat extremes or two drought + heat extremes) to those of assemblages being exposed only to the second extreme. Additionally, the recovery period between the successive extremes was varied (2, 3.5 or 6 weeks). Among the different types of climate extremes, combined drought + heat extremes induced substantial leaf and plant mortality, while the effects of drought and heat extremes were smaller. Preceding drought + heat extremes lowered the resistance in terms of leaf survival to a subsequent drought + heat extreme if the recovery period was 2 weeks, even though the leaves had completely recovered during that interval. No reduced resistance to subsequent extremes was recorded with longer recovery times or with drought or heat extremes. Despite mortality on the short term, the drought + heat and the heat extremes increased the end-of-season aboveground biomass, independent of the number of events or the recovery period. These results show that the effect of a preceding extreme event disappears quite quickly, but that recurrent climate extremes with short time intervals can weaken the resistance of herbaceous plant assemblages. This can however be compensated afterwards through rapid recovery and secondary, positive effects in the longer term.

Description: Climate and land use change impacts on global terrestrial ecosystems, fire, and river flows in the HadGEM2-ES Earth System Model using the Representative Concentration Pathways Biogeosciences Discussions, 10, 6171-6223, 2013 Author(s): R. A. Betts, N. Golding, P. Gonzalez, J. Gornall, R. Kahana, G. Kay, L. Mitchell, and A. Wiltshire A new generation of an Earth System Model now includes a number of land surface processes directly relevant to analyzing potential impacts of climate change. This model, HadGEM2-ES, allows us to assess the impacts of climate change, multiple interactions, and feedbacks as the model is run. This paper discusses the results of century-scale HadGEM2-ES simulations from an impacts perspective–specifically, terrestrial ecosystems and water resources–for four different scenarios following the Representative Concentration Pathways (RCPs), being used for next assessment report of the Intergovernmental Panel on Climate Change (IPCC). Over the 21st Century, simulated changes in global and continential-scale terrestrial ecosystems due to climate change appear to be very similar in all 4 RCPs, even though the level of global warming by the end of the 21st Century ranges from 2 °C in the lowest scenario to 5.5° in the highest. A warming climate generally favours broadleaf trees over needleleaf, needleleaf trees over shrubs, and shrubs over herbaceous vegetation, resulting in a poleward shift of temperate and boreal forests and woody tundra in all scenarios. Although climate related changes are slightly larger in scenarios of greater warming, the largest differences between scenarios arise at regional scales as a consequence of different patterns of anthropogenic land cover change. In the model, the scenario with the lowest global warming results in the most extensive decline in tropical forest cover due to a large expansion of agriculture. Under all four RCPs, fire potential could increase across extensive land areas, particularly tropical and sub-tropical latitudes. River outflows are simulated to increase with higher levels of CO 2 and global warming in all projections, with outflow increasing with mean temperature at the end of the 21st Century at the global scale and in North America, Asia, and Africa. In South America, Europe, and Australia, the relationship with climate warming and CO 2 rise is less clear, probably as a result of land cover change exerting a dominant effect in those regions.

Description: Horizontal distribution of Fukushima-derived radiocesium in zooplankton in the northwestern Pacific Ocean Biogeosciences Discussions, 10, 6143-6170, 2013 Author(s): M. Kitamura, Y. Kumamoto, H. Kawakami, E. C. Cruz, and K. Fujikura The magnitude of the 9.0 Tohoku earthquake and the ensuing tsunami on 11 March 2011, inflicted heavy damage on the Fukushima Dai-ichi nuclear power plant (FNPP1). Fission products were emitted, falling over a broad range in the northern hemisphere, and water contaminated with radionuclides leaked into the ocean. In this study, we described the horizontal distribution of the Fukushima-derived radiocesium in zooplankton and in seawater in the western North Pacific Ocean (500–2100 km from the FNPP1) 10 months after the accident. 134 Cs and 137 Cs were detected in zooplankton and seawater from all the stations. Because of its short half-lives, 134 Cs detected in our samples could only be derived from the FNPP1 accident. The highest 137 Cs activity in zooplankton was same order of magnitude as that one month after the accident, and average activity was one or two orders of magnitude higher than 137 Cs activities observed before the accident around Japan. Horizontally, the radiocesium activity concentrations in zooplankton were high at around 25° N while those in surface seawater were high at around the transition area between the Kuroshio and the Oyashio Currents (36–40° N). We observed subsurface radiocesium maxima in density range of the North Pacific Subtropical Mode Water and occurrence of many diel vertical migratory zooplanktons. These suggested that the high activity concentrations in the subtropical zooplankton at around 25° N were connected to the subsurface radiocesium and active vertical migration of zooplankton. However, the high activity concentrations of radiocesium in subsurface seawater did not necessarily follow the higher radiocesium activity in zooplankton. Biological characteristics of zooplankton community possibly influenced how large was contamination of radiocesium in the community but it is still unknown what kind of biological factors were important.

Description: Estimating carbonate parameters from hydrographic data for the intermediate and deep waters of the Southern Hemisphere Oceans Biogeosciences Discussions, 10, 6225-6257, 2013 Author(s): H. C. Bostock, S. E. Mikaloff Fletcher, and M. J. M. Williams Using GLODAP and CLIVAR ocean carbon data, we have developed several multiple linear regression (MLR) algorithms to estimate alkalinity and dissolved inorganic carbon (DIC) in the intermediate and deep waters of the Southern Hemisphere (south of 25° S) from only hydrographic data (temperature, salinity and dissolved oxygen). A Monte Carlo experiment was used to identify a potential density (σ θ ) of 27.5 as an optimal break point between the two regimes with different MLR algorithms. The algorithms provide a good estimate of DIC ( R 2 =0.98) and alkalinity ( R 2 =0.91), and excellent agreement for aragonite and calcite saturation states ( R 2 =0.99). Combining the algorithms with the CSIRO Atlas of Regional Seas (CARS), we have been able to map the calcite saturation horizon (CSH) and aragonite saturation horizon (ASH) for the Southern Ocean at a spatial resolution of 0.5°. These maps are more detailed and more consistent with oceanography than the gridded GLODAP data. The high resolution ASH map reveals a dramatic circumpolar shoaling at the Polar Front. North of 40° S the CSH is deepest in the Atlantic (~ 4000 m) and shallower in the Pacific Ocean (~ 2750 m), while the CSH sits between 3200 and 3400 m in the Indian Ocean.

Description: Ice-shelf buttressing and the stability of marine ice sheets The Cryosphere, 7, 647-655, 2013 Author(s): G. H. Gudmundsson Ice-shelf buttressing and the stability of marine-type ice sheets are investigated numerically. Buttressing effects are analysed for a situation where a stable grounding line is located on a bed sloping upwards in the direction of flow. Such grounding-line positions are known to be unconditionally unstable in the absence of transverse flow variations. It is shown that ice-shelf buttressing can restore stability under these conditions. Ice flux at the grounding line is, in general, not a monotonically increasing function of ice thickness. This, possibly at first somewhat counterintuitive result, is found to be fully consistent with recent theoretical work. Grounding lines on retrograde slopes are conditionally stable, and the stability regime is a non-trivial function of bed and ice-shelf geometry. The stability of grounding lines cannot be assessed from considerations of local bed slope only.

Description: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data The Cryosphere, 7, 615-630, 2013 Author(s): M. Tedesco, X. Fettweis, T. Mote, J. Wahr, P. Alexander, J. E. Box, and B. Wouters A combined analysis of remote sensing observations, regional climate model (RCM) outputs and reanalysis data over the Greenland ice sheet provides evidence that multiple records were set during summer 2012. Melt extent was the largest in the satellite era (extending up to ∼97% of the ice sheet) and melting lasted up to ∼2 months longer than the 1979–2011 mean. Model results indicate that near surface temperature was ∼3 standard deviations (σ) above the 1958–2011 mean, while surface mass balance (SMB) was ∼3σ below the mean and runoff was 3.9σ above the mean over the same period. Albedo, exposure of bare ice and surface mass balance also set new records, as did the total mass balance with summer and annual mass changes of, respectively, −627 Gt and −574 Gt, 2σ below the 2003–2012 mean. We identify persistent anticyclonic conditions over Greenland associated with anomalies in the North Atlantic Oscillation (NAO), changes in surface conditions (e.g., albedo, surface temperature) and preconditioning of surface properties from recent extreme melting as major driving mechanisms for the 2012 records. Less positive if not increasingly negative SMB will likely occur should these characteristics persist.

Description: Productivity of aboveground coarse wood biomass and stand age related to soil hydrology of Amazonian forests in the Purus-Madeira interfluvial area Biogeosciences Discussions, 10, 6417-6459, 2013 Author(s): B. B. L. Cintra, J. Schietti, T. Emillio, D. Martins, G. Moulatlet, P. Souza, C. Levis, C. A. Quesada, and J. Schöngart The ongoing demand for information on forest productivity has increased the number of permanent monitoring plots across the Amazon. Those plots, however, do not comprise the whole diversity of forest types in the Amazon. The complex effects of soil, climate and hydrology on the productivity of seasonally waterlogged interfluvial wetland forests are still poorly understood. The presented study is the first field-based estimate for tree ages and wood biomass productivity in the vast interfluvial region between the Purus and Madeira rivers. We estimate stand age and wood biomass productivity by a combination of tree-ring data and allometric equations for biomass stocks of eight plots distributed along 600 km in the Purus-Madeira interfluvial area that is crossed by the BR-319 highway. We relate stand age and wood biomass productivity to hydrological and edaphic conditions. Mean productivity and stand age were 5.6 ± 1.1 Mg ha −1 yr −1 and 102 ± 18 yr, respectively. There is a strong relationship between tree age and diameter, as well as between mean diameter increment and mean wood density within a plot. Regarding the soil hydromorphic properties we find a positive correlation with wood biomass productivity and a negative relationship with stand age. Productivity also shows a positive correlation with the superficial phosphorus concentration. In addition, superficial phosphorus concentration increases with enhanced soil hydromorphic condition. We raise three hypotheses to explain these results: (1) the reduction of iron molecules on the saturated soils with plinthite layers close to the surface releases available phosphorous for the plants; (2) the poor structure of the saturated soils creates an environmental filter selecting tree species of faster growth rates and shorter life spans and (3) plant growth on saturated soil is favored during the dry season, since there should be low restrictions for soil water availability.

Description: Application of wavelet transform for evaluation of hydrocarbon reservoirs: example from Iranian oil fields in the north of the Persian Gulf Nonlinear Processes in Geophysics, 20, 231-238, 2013 Author(s): M. R. Saadatinejad and H. Hassani The Persian Gulf and its surrounding area are some of the biggest basins and have a very important role in producing huge amounts of hydrocarbon, and this potential was evaluated in order to explore the target for geoscientists and petroleum engineers. Wavelet transform is a useful and applicable technique to reveal frequency contents of various signals in different branches of science and especially in petroleum studies. We applied two major capacities of continuous mode of wavelet transform in seismic investigations. These investigations were operated to detect reservoir geological structures and some anomalies related to hydrocarbon to develop and explore new petroleum reservoirs in at least 4 oilfields in the southwest of Iran. It had been observed that continuous wavelet transform results show some discontinuities in the location of faults and are able to display them more clearly than other seismic methods. Moreover, continuous wavelet transform, utilizing Morlet wavelet, displays low-frequency shadows on 4 different iso-frequency vertical sections to identify reservoirs containing gas. By comparing these different figures, the presence of low-frequency shadows under the reservoir could be seen and we can relate these variations from anomalies at different frequencies as an indicator of the presence of hydrocarbons in the target reservoir.

Description: Sources and fate of terrestrial dissolved organic carbon in lakes of a Boreal Plains region recently affected by wildfire Biogeosciences Discussions, 10, 6093-6141, 2013 Author(s): D. Olefeldt, K. J. Devito, and M. R. Turetsky Downstream mineralization and sedimentation of terrestrial dissolved organic carbon (DOC) render lakes important for landscape carbon cycling in the boreal region, with regulating processes potentially sensitive to perturbations associated with climate change including increased occurrence of wildfire. In this study we assessed chemical composition and reactivity (during both dark and UV incubations) of DOC from lakes and terrestrial sources within a peatland-rich western boreal plains region partially affected by a recent wildfire. While wildfire was found to increase aromaticity of DOC in peat pore-water above the water table, it had no effect on concentrations or composition of DOC from peatland wells and neither affected mineral well or lake DOC characteristics. Lake DOC composition reflected a mixing of peatland and mineral groundwater, with a greater influence of mineral sources to lakes in coarse- than fine-textured settings. Peatland DOC was less biodegradable than mineral DOC, but both mineralization and sedimentation of peatland DOC increased substantially during UV incubations through selective removal of aromatic humic and fulvic acids. DOC composition in lakes with longer residence times had characteristics consistent with increased UV-mediated processing. We estimate that about half of terrestrial DOC inputs had been lost within lakes, mostly due to UV-mediated processes. The importance of within-lake losses of aromatic DOC from peatland sources through UV-mediated processes indicate that terrestrial-aquatic C linkages in the study region are largely disconnected from recent terrestrial primary productivity. Together, our results suggest that characteristics of the study region (climate, surface geology and lake morphometry) render linkages between terrestrial and aquatic C cycling insensitive to the effects of wildfire by determining dominant terrestrial sources and within-lake processes of DOC removal.

Description: Balanced conditions or slight mass gain of glaciers in the Lahaul and Spiti region (northern India, Himalaya) during the nineties preceded recent mass loss The Cryosphere, 7, 569-582, 2013 Author(s): C. Vincent, Al. Ramanathan, P. Wagnon, D. P. Dobhal, A. Linda, E. Berthier, P. Sharma, Y. Arnaud, M. F. Azam, P. G. Jose, and J. Gardelle The volume change of the Chhota Shigri Glacier (India, 32° 20 N, 77° 30' E) between 1988 and 2010 has been determined using in situ geodetic measurements. This glacier has experienced only a slight mass loss between 1988 and 2010 (−3.8 ± 2.0 m w.e. (water equivalent) corresponding to −0.17 ± 0.09 m w.e. yr −1 ). Using satellite digital elevation models (DEM) differencing and field measurements, we measure a negative mass balance (MB) between 1999 and 2010 (−4.8 ± 1.8 m w.e. corresponding to −0.44 ± 0.16 m w.e. yr −1 ). Thus, we deduce a slightly positive or near-zero MB between 1988 and 1999 (+1.0 ± 2.7 m w.e. corresponding to +0.09 ± 0.24 m w.e. yr −1 ). Furthermore, satellite DEM differencing reveals that the MB of the Chhota Shigri Glacier (−0.39 ± 0.15 m w.e. yr −1 ) has been only slightly less negative than the MB of a 2110 km 2 glaciarized area in the Lahaul and Spiti region (−0.44 ± 0.09 m w.e. yr −1 ) during 1999−2011. Hence, we conclude that the ice wastage is probably moderate in this region over the last 22 yr, with near equilibrium conditions during the nineties, and an ice mass loss after. The turning point from balanced to negative mass budget is not known but lies probably in the late nineties and at the latest in 1999. This positive or near-zero MB for Chhota Shigri Glacier (and probably for the surrounding glaciers of the Lahaul and Spiti region) during at least part of the 1990s contrasts with a recent compilation of MB data in the Himalayan range that indicated ice wastage since 1975. However, in agreement with this compilation, we confirm more negative balances since the beginning of the 21st century.

Description: Solute specific scaling of inorganic nitrogen and phosphorus uptake in streams Biogeosciences Discussions, 10, 6671-6693, 2013 Author(s): R. O. Hall Jr., M. A. Baker, E. J. Rosi-Marshall, and J. L. Tank Stream ecosystem processes such as nutrient cycling may vary with stream position in the watershed. Using a scaling approach, we examined the relationship between stream size and nutrient uptake length, which represents the mean distance that a dissolved solute travels prior to removal from the water column. Ammonium uptake length increased proportionally with stream size measured as specific discharge (discharge/stream width) with a scaling exponent = 1.01. In contrast, the scaling exponent for nitrate (NO 3 − ) was 1.19 and for soluble reactive phosphorus (SRP) was 1.35, suggesting that uptake lengths for these nutrients increased more rapidly than increases in specific discharge. Additionally, the ratio of nitrogen (N) uptake length to SRP uptake length declined with stream size; there was lower demand for SRP relative to N as stream size increased. Ammonium and NO 3 − uptake velocity positively related with stream metabolism, while SRP did not. Finally, we related the scaling of uptake length and specific discharge to that of stream length using Hack's law and downstream hydraulic geometry. Ammonium uptake length increased less than proportionally with distance from the headwaters, suggesting a strong role for larger streams and rivers in regulating nutrient transport.

Description: Responses of lower trophic-level organisms to typhoon passage on the outer shelf of the East China Sea: an incubation experiment Biogeosciences Discussions, 10, 6605-6635, 2013 Author(s): N. Yasuki, K. Suzuki, and A. Tsuda Typhoons can induce vertical mixing, upwelling, or both in the water column due to strong wind stress. These events can induce phytoplankton blooms in the oligotrophic ocean after typhoon passage. However, little is known about the responses of lower trophic-level organisms or changes in the community structure following the passage of typhoons, particularly in offshore regions. Therefore, we evaluated community succession on the outer shelf of the East China Sea through on-deck bottle incubation experiments simulating hydrographic conditions after the passage of a typhoon. Under all of the experimental conditions we tested, chlorophyll a concentrations increased more than 9-fold within 6 days, and these algal cells were mainly composed of large diatoms (〉10 μm). Ciliates also increased along with the diatom bloom. These results suggest that increases in diatom and ciliate populations may enhance biogenic carbon export in the water column. Typhoons can affect not only phytoplankton productivity, but also the composition of lower trophic-level organisms and biogeochemical processes in oligotrophic offshore regions.

Description: The geomorphological effect of cornice fall avalanches in the Longyeardalen valley, Svalbard The Cryosphere, 7, 1361-1374, 2013 Author(s): M. Eckerstorfer, H. H. Christiansen, L. Rubensdotter, and S. Vogel The study of snow avalanches and their geomorphological effect in the periglacial parts of the cryosphere is important for enhanced geomorphological process understanding as well as hazard-related studies. Only a few field studies, and particularly few in the High Arctic, have quantified avalanche sedimentation. Snow avalanches are traditionally ranked behind rockfall in terms of their significance for mass-wasting processes of rockslopes. Cornice fall avalanches are at present the most dominant snow avalanche type at two slope systems, called Nybyen and Larsbreen, in the valley Longyeardalen in central Svalbard. Both slope systems are on northwest-facing lee slopes underneath a large summit plateau, with annual cornices forming on the top. High-frequency and magnitude cornice fall avalanching is observed by daily automatic time-lapse photography. In addition, rock debris sedimentation by cornice fall avalanches was measured directly in permanent sediment traps or by snow inventories. The results from a maximum of seven years of measurements in a total of 13 catchments show maximum mean rock debris sedimentation rates ranging from 8.2 to 38.7 kg m −2 at Nybyen, and from 0.8 to 55.4 kg m −2 at Larsbreen. Correspondingly, avalanche fan surfaces accreted from 2.6 to 8.8 mm yr −1 at Nybyen, and from 0.2 to 13.9 mm yr −1 at Larsbreen. This comparably efficient rockslope mass wasting is due to collapsing cornices producing cornice fall avalanches containing large amounts of rock debris throughout the entire winter. The rock debris of different origin stems from the plateau crests, the adjacent free rock face and the transport pathway, accumulating distinct avalanche fans at both slope systems. Cornice fall avalanche sedimentation also contributed to the development of a rock glacier at the Larsbreen site during the Holocene. We have recorded present maximum rockwall retreat rates of 0.9 mm yr −1 at Nybyen, but as much as 6.7 mm yr −1 at Larsbreen, while average Holocene rockwall retreat rates of 1.1 mm yr −1 at Nybyen have been determined earlier. As cornice fall avalanches are the dominant type of avalanche in central Svalbard, the related geomorphological effect is assumed to be of significance at periglacial landscape scale. A climate-induced shift in prevailing winter wind direction could change the rockslope sedimentation effectively by changing the snow avalanche activity.

Description: Icehouse-greenhouse variations in marine denitrification Biogeosciences Discussions, 10, 14769-14813, 2013 Author(s): T. J. Algeo, P. A. Meyers, R. S. Robinson, H. Rowe, and G. Q. Jiang Long-term secular variation in the isotopic composition of seawater fixed nitrogen (N) is poorly known. Here, we document variation in the N-isotopic composition of marine sediments (δ 15 N sed ) since 660 Ma (million years ago) in order to understand major changes in the marine N cycle through time and their relationship to first-order climate variation. During the Phanerozoic, greenhouse climate modes were characterized by low δ 15 N sed (∼ −2 to +2‰) and icehouse climate modes by high δ 15 N sed (∼ +4 to +8‰). Shifts toward higher δ 15 N sed occurred rapidly during the early stages of icehouse modes, prior to the development of major continental glaciation, suggesting a potentially important role for the marine N cycle in long-term climate change. Reservoir box modeling of the marine N cycle demonstrates that secular variation in δ 15 N sed was likely due to changes in the dominant locus of denitrification, with a shift in favor of sedimentary denitrification during greenhouse modes owing to higher eustatic (global sea-level) elevations and greater on-shelf burial of organic matter, and a shift in favor of water-column denitrification during icehouse modes owing to lower eustatic elevations, enhanced organic carbon sinking fluxes, and expanded oceanic oxygen-minimum zones. The results of this study provide new insights into operation of the marine N cycle, its relationship to the global carbon cycle, and its potential role in modulating climate change at multimillion-year timescales.

Description: Technical Note: Disturbance of soil structure can lead to release of methane entrapped in glacier forefield soils Biogeosciences Discussions, 10, 14815-14834, 2013 Author(s): P. A. Nauer, E. Chiri, J. Zeyer, and M. H. Schroth Investigations of sources and sinks of atmospheric CH 4 are needed to understand the global CH 4 cycle and climate-change mitigation options. Glaciated environments might play a critical role due to potential feedbacks with global glacial meltdown. In an emerging glacier forefield, an ecological shift occurs from an anoxic, potentially methanogenic subglacial sediment to an oxic proglacial soil, in which soil-microbial consumption of atmospheric CH 4 is initiated. The development of this change in CH 4 turnover can be quantified by soil-gas profile analysis. We found evidence for CH 4 entrapped in glacier forefield soils when comparing two methods for the collection of soil-gas samples: a modified steel rod (SR) designed for one-time sampling and rapid screening (samples collected ~ 1 min after hammering the SR into the soil), and a novel multi-level sampler (MLS) for repetitive sampling through a previously installed access tube (samples collected weeks after access-tube installation). In glacier forefields on siliceous bedrock, sub-atmospheric CH 4 concentrations were observed with both methods. Conversely, elevated soil-CH 4 concentrations were observed in calcareous glacier forefields, but only in samples collected with the SR, while MLS samples all showed sub-atmospheric CH 4 concentrations. Time-series SR soil-gas sampling (additional samples collected 2, 3, 5, and 7 min after hammering) confirmed the transient nature of the elevated soil-CH 4 concentrations, which were decreasing from ~ 100 μL L −1 towards background levels within minutes. This hints towards the existence of entrapped CH 4 in calcareous glacier forefield soil that can be released when sampling soil-gas with the SR. Laboratory experiments with miniature soil cores collected from two glacier forefields confirmed CH 4 entrapment in these soils. Treatment by sonication and acidification resulted in a massive release of CH 4 from calcareous cores (on average 0.3–1.8 μg CH 4 (g d.w.) −1 ); release from siliceous cores was 1–2 orders of magnitude lower (0.02–0.03 μg CH 4 (g d.w.) −1 ). Clearly, some form of CH 4 entrapment exists in calcareous glacier forefield soils, and to a much lesser extent in siliceous glacier forefield soils. Its nature and origin remain unclear and will be subject of future investigations.

Description: Dynamics of particulate organic carbon flux in a global ocean model Biogeosciences Discussions, 10, 14715-14767, 2013 Author(s): I. D. Lima, P. J. Lam, and S. C. Doney The sinking of particulate organic carbon (POC) is a key component of the ocean carbon cycle and plays an important role in the global climate system. However, the processes controlling the fraction of primary production that is exported from the euphotic zone (export ratio) and how much of it survives respiration in the mesopelagic to be sequestered in the deep ocean (transfer efficiency) are not well understood. In this study, we use a three-dimensional, coupled physical-biogeochemical model (CCSM-BEC) to investigate the processes controlling the export of particulate organic matter from the euphotic zone and its flux to depth. We also compare model results with sediment trap data and other parameterizations of POC flux to depth to evaluate model skill and gain further insight into the causes of error and uncertainty in POC flux estimates. In the model, export ratios are mainly a function of diatom relative abundance and temperature while absolute fluxes and transfer efficiency are driven by mineral ballast composition of sinking material. The temperature dependence of the POC remineralization length scale is modulated by denitrification under low O 2 concentrations and lithogenic (dust) fluxes. Lithogenic material is an important control of transfer efficiency in the model, but its effect is restricted to regions of strong atmospheric dust deposition. In the remaining regions, CaCO 3 content of exported material is the main factor affecting transfer efficiency. The fact that mineral ballast composition is inextricably linked to plankton community structure results in correlations between export ratios and ballast minerals fluxes (opal and CaCO 3 ), and transfer efficiency and diatom relative abundance that do not necessarily reflect ballast or direct ecosystem effects, respectively. This suggests that it might be difficult to differentiate between ecosystem and ballast effects in observations. The model's skill at reproducing sediment trap observations is equal to or better than that of other parameterizations. However, the sparseness and relatively large uncertainties of sediment trap data makes it difficult to accurately evaluate the skill of the model and other parameterizations. More POC flux observations, over a wider range of ecological regimes, are necessary to thoroughly evaluate and test model results and better understand the processes controlling POC flux to depth in the ocean.

Description: Thin terrestrial sediment deposits on intertidal sandflats: effects on pore water solutes and juvenile bivalve burial behaviour Biogeosciences Discussions, 10, 14835-14860, 2013 Author(s): A. Hohaia, K. Vopel, and C. A. Pilditch Changes in land use and climate increase the supply of terrestrial sediment (hereafter, TS) to coastal waters worldwide but the effects of these sediments on benthic ecosystem functioning are not well known. Past experiments with defaunated, intertidal sediment suggested a link between the de-oxygenation of sediments underlying a thin (mm) layer of TS and reduced burial rates of juvenile macrofaunal recruits. We examined this link predicting that surficial TS deposits will still negatively affect burial when applied to sediments that are initially well oxygenated due to bioturbation (C) or depleted of organic matter (D). We observed the behaviour of post-settlement juveniles of the tellinid bivalve \textit{Macomona liliana} on the surface of four treatments; C, D, and the same sediments to which we added a thin layer of TS (CTS, DTS). Pore water analyses confirmed that the diffusive impedance of the 1.7–1.9 mm TS deposit decreased the oxygenation of the underlying intertidal sediment (CTS) but not that of the depleted sediment (DTS). Unexpectedly, (1) the application of a TS deposit significantly increased but not decreased the probability of burial, irrespectively of treatment, and (2) juveniles more likely buried into C than into D. We attribute the failure to document a negative effect of TS on the recruits' burial to the activity of the resident macroinfauna (CTS) or the absence of organic matter (DTS). Our results underline the important role of the resident macrofauna in mediating the stress response of benthic ecosystems.

Description: Calcium carbonate corrosivity in an Alaskan inland sea Biogeosciences Discussions, 10, 14887-14922, 2013 Author(s): W. Evans, J. T. Mathis, and J. N. Cross Ocean acidification is the hydrogen ion increase caused by the oceanic uptake of anthropogenic CO 2 , and is a focal point in marine biogeochemistry, in part, because this chemical reaction reduces calcium carbonate (CaCO 3 ) saturation states (Ω) to levels that are corrosive (i.e. Ω ≤ 1) to shell-forming marine organisms. However, other processes can drive CaCO 3 corrosivity; specifically, the addition of tidewater glacial melt. Carbonate system data collected in May and September from 2009 through 2012 in Prince William Sound (PWS), a semi-enclosed inland sea located on the south-central coast of Alaska that is ringed with fjords containing tidewater glaciers, reveal the unique impact of glacial melt on CaCO 3 corrosivity. Initial limited sampling was expanded in September 2011 to span large portions of the western and central sound, and included two fjords proximal to tidewater glaciers: Icy Bay and Columbia Bay. The observed conditions in these fjords affected CaCO 3 corrosivity in the upper water column (

Description: Surface circulation and upwelling patterns around Sri Lanka Biogeosciences Discussions, 10, 14953-14998, 2013 Author(s): A. de Vos, C. B. Pattiaratchi, and E. M. S. Wijeratne Sri Lanka occupies a unique location within the equatorial belt in the northern Indian Ocean with the Arabian Sea on its western side and the Bay of Bengal on its eastern side. The region is characterised by bi-annually reversing monsoon winds resulting from seasonal differential heating and cooling of the continental land mass and the ocean. This study explored elements of the dynamics of the surface circulation and coastal upwelling in the waters around Sri Lanka using satellite imagery and the Regional Ocean Modelling System (ROMS) configured to the study region and forced with ECMWF interim data. The model was run for 2 yr to examine the seasonal and shorter term (∼10 days) variability. The results confirmed the presence of the reversing current system in response to the changing wind field: the eastward flowing Southwest Monsoon Current (SMC) during the Southwest (SW) monsoon transporting 11.5 Sv and the westward flowing Northeast Monsoon Current (NMC) transporting 9.5 Sv during the Northeast (NE) monsoon, respectively. A recirculation feature located to the east of Sri Lanka during the SW monsoon, the Sri Lanka Dome, is shown to result from the interaction between the SMC and the Island of Sri Lanka. Along the eastern and western coasts, during both monsoon periods, flow is southward converging along the south coast. During the SW monsoon the Island deflects the eastward flowing SMC southward whilst along the east coast the southward flow results from the Sri Lanka Dome recirculation. The major upwelling region, during both monsoon periods, is located along the south coast and is shown to be due to flow convergence and divergence associated with offshore transport of water. Higher surface chlorophyll concentrations were observed during the SW monsoon. The location of the flow convergence and hence the upwelling centre was dependent on the relative strengths of wind driven flow along the east and west coasts: during the SW (NE) monsoon the flow along the western (eastern) coast was stronger and hence the upwelling centre was shifted to the east (west). The presence of upwelling along the south coast during both monsoon periods may explain the blue whale ( Balaenoptera musculus ) aggregations in this region.

Description: Strong sensitivity of Southern Ocean carbon uptake and nutrient cycling to wind stirring Biogeosciences Discussions, 10, 15033-15076, 2013 Author(s): K. B. Rodgers, O. Aumont, S. E. Mikaloff Fletcher, Y. Plancherel, L. Bopp, C. de Boyer Montégut, D. Iudicone, R. F. Keeling, G. Madec, and R. Wanninkhof Here we test the hypothesis that winds have an important role in determining the rate of exchange of CO 2 between the atmosphere and ocean through wind stirring over the Southern Ocean. This is tested with a sensitivity study using an ad hoc parameterization of wind stirring in an ocean carbon cycle model. The objective is to identify the way in which perturbations to the vertical density structure of the planetary boundary in the ocean impacts the carbon cycle and ocean biogeochemistry. Wind stirring leads to reduced uptake of CO 2 by the Southern Ocean over the period 2000–2006, with differences of order 0.9 Pg C yr −1 over the region south of 45° S. Wind stirring impacts not only the mean carbon uptake, but also the phasing of the seasonal cycle of carbon and other species associated with ocean biogeochemistry. Enhanced wind stirring delays the seasonal onset of stratification, and this has large impacts on both entrainment and the biological pump. It is also found that there is a strong sensitivity of nutrient concentrations exported in Subantarctic Mode Water (SAMW) to wind stirring. This finds expression not only locally over the Southern Ocean, but also over larger scales through the impact on advected nutrients. In summary, the large sensitivity identified with the ad hoc wind stirring parameterization offers support for the importance of wind stirring for global ocean biogeochemistry, through its impact over the Southern Ocean.

Description: Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram The Cryosphere, 7, 1385-1398, 2013 Author(s): R. Bhambri, T. Bolch, P. Kawishwar, D. P. Dobhal, D. Srivastava, and B. Pratap Glaciers in the Karakoram show long-term irregular behaviour with comparatively frequent and sudden advances. A glacier inventory of the upper Shyok valley situated in northeast Karakoram has been generated for the year 2002 using Landsat ETM+ and SRTM3 DEM as baseline data for the investigations and subsequent change analysis. The upper Shyok valley contained 2123 glaciers (larger than 0.02 km 2 in size) with an area of 2977.9 ± 95.3 km 2 in 2002. Out of these, 18 glaciers with an area of 1004.1 ± 32.1 km 2 showed surge-type behaviour. Change analysis based on Hexagon KH-9 (years 1973 and 1974) and Landsat TM/ETM+ (years 1989, 2002 and 2011) images had to be restricted to a subset of 136 glaciers (covering an area of 1609.7 ± 51.5 km 2 in 2002) due to adverse snow conditions. The area of the investigated glaciers, including the 18 surge-type glaciers identified, showed no significant changes during all studied periods. However, the analysis provides a hint that the overall glacier area slightly decreased until about 1989 (area 1973: 1613.6 ± 43.6 km 2 ; area 1989: 1602.0 ± 33.6 km 2 ) followed by an increase (area 2002: 1609.7 ± 51.5; area 2011: 1615.8 ± 35.5 km 2 ). Although the overall change in area is insignificant, advances in glacier tongues since the end of the 1980s are clearly visible. Detailed estimations of length changes for individual glaciers since the 1970s and for Central Rimo Glacier since the 1930s confirm the irregular retreat and advance.

Description: Distributed allocation of mobile sensing swarms in gyre flows Nonlinear Processes in Geophysics, 20, 657-668, 2013 Author(s): K. Mallory, M. A. Hsieh, E. Forgoston, and I. B. Schwartz We address the synthesis of distributed control policies to enable a swarm of homogeneous mobile sensors to maintain a desired spatial distribution in a geophysical flow environment, or workspace. In this article, we assume the mobile sensors (or robots) have a "map" of the environment denoting the locations of the Lagrangian coherent structures or LCS boundaries. Using this information, we design agent-level hybrid control policies that leverage the surrounding fluid dynamics and inherent environmental noise to enable the team to maintain a desired distribution in the workspace. We discuss the stability properties of the ensemble dynamics of the distributed control policies. Since realistic quasi-geostrophic ocean models predict double-gyre flow solutions, we use a wind-driven multi-gyre flow model to verify the feasibility of the proposed distributed control strategy and compare the proposed control strategy with a baseline deterministic allocation strategy. Lastly, we validate the control strategy using actual flow data obtained by our coherent structure experimental testbed.

Description: Weak layer fracture: facets and depth hoar The Cryosphere, 7, 1447-1453, 2013 Author(s): I. Reiweger and J. Schweizer Understanding failure initiation within weak snow layers is essential for modeling and predicting dry-snow slab avalanches. We therefore performed laboratory experiments with snow samples containing a weak layer consisting of either faceted crystals or depth hoar. During these experiments the samples were loaded with different loading rates and at various tilt angles until fracture. The strength of the samples decreased with increasing loading rate and increasing tilt angle. Additionally, we took pictures of the side of four samples with a high-speed video camera and calculated the displacement using a particle image velocimetry (PIV) algorithm. The fracture process within the weak layer could thus be observed in detail. Catastrophic failure started due to a shear fracture just above the interface between the depth hoar layer and the underlying crust.

Description: Scatter of mass changes estimates at basin scale for Greenland and Antarctica The Cryosphere, 7, 1411-1432, 2013 Author(s): V. R. Barletta, L. S. Sørensen, and R. Forsberg During the last decade, the GRACE mission has provided valuable data for determining the mass changes of the Greenland and Antarctic ice sheets. Yet, discrepancies still exist in the published mass balance results, and comprehensive analyses on the sources of errors and discrepancies are lacking. Here, we present monthly mass changes together with trends derived from GRACE data at basin scale for both the Greenland and Antarctic ice sheets, and we assess the variability and errors for each of the possible sources of discrepancies, and we do this in an unprecedented systematic way, taking into account mass inference methods, data sets and background models. We find a very good agreement between the monthly mass change results derived from two independent methods, which represents a cross validation. For the monthly solutions, we find that most of the scatter is caused by the use of the two different data sets rather than the two different methods applied. Besides the well-known GIA trend uncertainty, we find that the geocenter motion and the recent de-aliasing corrections significantly impact the trends, with contributions of +13.2 Gt yr −1 and −20 Gt yr −1 , respectively, for Antarctica, which is more affected by these than Greenland. We show differences between the use of release RL04 and the new RL05 and confirm a lower noise content in the new release. The overall scatter of the solutions well exceeds the uncertainties propagated from the data errors and the leakage (as done in the past); hence we calculate new sound total errors for the monthly solutions and the trends. We find that the scatter in the monthly solutions caused by applying different estimates of geocenter motion time series (degree-1 corrections) is significant – contributing with up to 40% of the total error. For the whole GRACE period (2003–2011) our trend estimate for Greenland is −234 ± 20 Gt yr −1 and −83 ± 36 Gt yr −1 for Antarctica (−111 ± 15 Gt yr −1 in the western part). We also find a clear (with respect to our errors) increase of mass loss in the last four years.

Description: Four-dimensional ensemble-variational data assimilation for global deterministic weather prediction Nonlinear Processes in Geophysics, 20, 669-682, 2013 Author(s): M. Buehner, J. Morneau, and C. Charette The goal of this study is to evaluate a version of the ensemble-variational data assimilation approach (EnVar) for possible replacement of 4D-Var at Environment Canada for global deterministic weather prediction. This implementation of EnVar relies on 4-D ensemble covariances, obtained from an ensemble Kalman filter, that are combined in a vertically dependent weighted average with simple static covariances. Verification results are presented from a set of data assimilation experiments over two separate 6-week periods that used assimilated observations and model configuration very similar to the currently operational system. To help interpret the comparison of EnVar versus 4D-Var, additional experiments using 3D-Var and a version of EnVar with only 3-D ensemble covariances are also evaluated. To improve the rate of convergence for all approaches evaluated (including EnVar), an estimate of the cost function Hessian generated by the quasi-Newton minimization algorithm is cycled from one analysis to the next. Analyses from EnVar (with 4-D ensemble covariances) nearly always produce improved, and never degraded, forecasts when compared with 3D-Var. Comparisons with 4D-Var show that forecasts from EnVar analyses have either similar or better scores in the troposphere of the tropics and the winter extra-tropical region. However, in the summer extra-tropical region the medium-range forecasts from EnVar have either similar or worse scores than 4D-Var in the troposphere. In contrast, the 6 h forecasts from EnVar are significantly better than 4D-Var relative to radiosonde observations for both periods and in all regions. The use of 4-D versus 3-D ensemble covariances only results in small improvements in forecast quality. By contrast, the improvements from using 4D-Var versus 3D-Var are much larger. Measurement of the fit of the background and analyzed states to the observations suggests that EnVar and 4D-Var can both make better use of observations distributed over time than 3D-Var. In summary, the results from this study suggest that the EnVar approach is a viable alternative to 4D-Var, especially when the simplicity and computational efficiency of EnVar are considered. Additional research is required to understand the seasonal dependence of the difference in forecast quality between EnVar and 4D-Var in the extra-tropics.

Description: Impact of seawater Ca 2+ on the calcification and calcite Mg/Ca of Amphistegina lessonii Biogeosciences Discussions, 11, 17463-17489, 2014 Author(s): A. Mewes, G. Langer, S. Thoms, G. Nehrke, G.-J. Reichart, L. J. de Nooijer, and J. Bijma Mg/Ca ratios in foraminiferal tests are routinely used as paleo temperature proxy, but on long timescales, also hold the potential to reconstruct past seawater Mg/Ca. Impact of both temperature and seawater Mg/Ca on Mg incorporation in foraminifera have been quantified by a number of studies. The underlying mechanism responsible for Mg incorporation in foraminiferal calcite and its sensitivity to environmental conditions, however, is not fully identified. A recently published biomineralization model (Nehrke et al., 2013) proposes a combination of transmembrane transport and seawater leakage or vacuolization to link calcite Mg/Ca to seawater Mg/Ca and explains inter-species variability in Mg/Ca ratios. To test the assumptions of this model, we conducted a culture study in which seawater Mg/Ca was manipulated by varying [Ca 2+ ] and keeping [Mg 2+ ] constant. Foraminiferal growth rates, test thickness and calcite Mg/Ca of newly formed chambers were analyzed. Results showed optimum growth rates and test thickness at Mg/Ca closest to that of ambient seawater. Calcite Mg/Ca is positively correlated to seawater Mg/Ca, indicating that not absolute seawater [Ca 2+ ] and [Mg 2+ ], but the telative ratio controls Mg/Ca in tests. These results demonstrate that the calcification process cannot be based only on seawater vacuolization, supporting the mixing model proposed by Nehrke et al. (2013). Here we, however, suggest a transmembrane transport fractionation that is not as strong as suggested by Nehrke et al. (2013).

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Description: Fluxes of carbon and nutrients to the Iceland Sea surface layer and inferred primary productivity and stoichiometry Biogeosciences Discussions, 11, 15399-15433, 2014 Author(s): E. Jeansson, R. G. J. Bellerby, I. Skjelvan, H. Frigstad, S. R. Ólafsdóttir, and J. Ólafsson Fluxes of carbon and nutrients to the upper 100 m of the Iceland Sea are evaluated. The study utilises hydro-chemical data from the Iceland Sea time-series station (68.00° N, 12.67° W), for the years between 1993 and 2006. By comparing data of dissolved inorganic carbon (DIC) and nutrients in the surface layer (upper 100 m), and a sub-surface layer (100–200 m), we calculate monthly deficits in the surface, and use these to deduce the surface layer fluxes that affect the deficits: vertical mixing, horizontal advection, air–sea exchange, and biological activity. The deficits show a clear seasonality with a minimum in winter, when the mixed layer is at the deepest, and a maximum in early autumn, when biological uptake has removed much of the nutrients. The annual vertical fluxes of DIC and nitrate amounts to 1.7 ± 0.3 and 0.23 ± 0.07 mol m −2 yr −1 , respectively, and the annual air–sea uptake of atmospheric CO 2 is 4.4 ± 1.1 mol m −2 yr −1 . The biologically driven changes in DIC during the year relates to net community production (NCP), and the net annual NCP corresponds to export production, and is here calculated to 6.1 ± 0.9 mol C m −2 yr −1 . The typical, median C : N ratio during the period of net community uptake is 11, and thus clearly higher than Redfield, but is varying during the season.

Description: Technical Note: Artificial coral reef mesocosms for ocean acidification investigations Biogeosciences Discussions, 11, 15463-15505, 2014 Author(s): J. Leblud, L. Moulin, A. Batigny, P. Dubois, and P. Grosjean The design and evaluation of replicated artificial mesocosms are presented in the context of a thirteen month experiment on the effects of ocean acidification on tropical coral reefs. They are defined here as (semi)-closed (i.e. with or without water change from the reef) mesocosms in the laboratory with a more realistic physico-chemical environment than microcosms. Important physico-chemical parameters (i.e. pH, p O 2 , p CO 2 , total alkalinity, temperature, salinity, total alkaline earth metals and nutrients availability) were successfully monitored and controlled. Daily variations of irradiance and pH were applied to approach field conditions. Results highlighted that it was possible to maintain realistic physico-chemical parameters, including daily changes, into artificial mesocosms. On the other hand, the two identical artificial mesocosms evolved differently in terms of global community oxygen budgets although the initial biological communities and physico-chemical parameters were comparable. Artificial reef mesocosms seem to leave enough degrees of freedom to the enclosed community of living organisms to organize and change along possibly diverging pathways.

Description: Phytoplankton community structure in the North Sea: coupling between remote sensing and automated in situ analysis at the single cell level Biogeosciences Discussions, 11, 15621-15662, 2014 Author(s): M. Thyssen, S. Alvain, A. Lefèbvre, D. Dessailly, M. Rijkeboer, N. Guiselin, V. Creach, and L.-F. Artigas Phytoplankton observation in the ocean can be a challenge in oceanography. Accurate estimations of their biomass and dynamics will help to understand ocean ecosystems and refine global climate models. This requires relevant datasets of phytoplankton at a functional level and on a daily and sub meso scale. In order to achieve this, an automated, high frequency, dedicated scanning flow cytometer (SFC, Cytobuoy, NL), has been developed to cover the entire size range of phytoplankton cells whilst simultaneously taking pictures of the largest of them. This cytometer was directly connected to the water inlet of a~pocket Ferry Box during a cruise in the North Sea, 8–12 May 2011 (DYMAPHY project, INTERREG IV A "2 Seas"), in order to identify the phytoplankton community structure of near surface waters (6 m) with a high resolution spacial basis (2.2 ± 1.8 km). Ten groups of cells, distinguished on the basis of their optical pulse shapes, were described (abundance, size estimate, red fluorescence per unit volume). Abundances varied depending on the hydrological status of the traversed waters, reflecting different stages of the North Sea blooming period. Comparisons between several techniques analyzing chlorophyll a and the scanning flow cytometer, using the integrated red fluorescence emitted by each counted cell, showed significant correlations. The community structure observed from the automated flow cytometry was compared with the PHYSAT reflectance anomalies over a daily scale. The number of matchups observed between the SFC automated high frequency in situ sampling and the remote sensing was found to be two to three times better than when using traditional water sampling strategies. Significant differences in the phytoplankton community structure within the two days for which matchups were available, suggest that it is possible to label PHYSAT anomalies not only with dominant groups, but at the level of the community structure.

Description: Long term effects on regional European boreal climate due to structural vegetation changes Biogeosciences Discussions, 11, 15507-15547, 2014 Author(s): J. H. Rydsaa, F. Stordal, and L. M. Tallaksen Amplified warming at high latitudes over the past decades has led to changes in the boreal and arctic climate system, such as structural changes in high latitude ecosystems and soil moisture properties. These changes trigger land-atmosphere feedbacks, through altered energy partitioning in response to changes in albedo and surface water fluxes. Local scale changes in the arctic and boreal zone may propagate to affect large scale climatic features. In this study, MODIS land surface data are used with the Weather Research and Forecasting model (WRF V3.5.1) and Noah LSM, in a series of experiments to simulate the influence of structural vegetation changes over a Northern European boreal ecosystem. Emphasis is placed on surface energy partitioning and near surface atmospheric variables, in order to investigate changes in atmospheric response due to observed and anticipated structural vegetation changes. We find that a northward migration of evergreen needle leaf forest into tundra regions causes an increase in latent rather than sensible heat fluxes, increased near surface temperatures and boundary layer height. Shrub expansion in tundra areas has only small effects on surface fluxes. However, it influences near surface wind speeds and boundary layer height. Northward migration of mixed forest across the present southern border of the boreal forest has largely opposite effects on surface fluxes and the near surface atmosphere, and acts to moderate the overall mean regional effects of boreal forest migration on the near surface atmosphere.

Description: Spatial variability and hotspots of soil N 2 O fluxes from intensively grazed grassland Biogeosciences Discussions, 11, 15327-15360, 2014 Author(s): N. J. Cowan, P. Norman, D. Famulari, P. E. Levy, D. S. Reay, and U. M. Skiba One hundred N 2 O flux measurements were made from an area of intensively managed grazed grassland in central Scotland using a high resolution dynamic chamber method. The field contained a variety of features from which N 2 O fluxes were measured including a manure heap, patches of decaying grass silage, and areas of increased sheep activity. Individual fluxes varied significantly across the field varying from 2 to 79 000 μg N 2 O-N m −2 h −1 . Soil samples were collected at 55 locations to investigate relationships between soil properties and N 2 O flux. Fluxes of N 2 O correlated strongly with soil NO 3 − concentrations. Distribution of NO 3 − and the high spatial variability of N 2 O flux across the field are shown to be linked to the distribution of waste from grazing animals and the resultant reactive nitrogen compounds in the soil which are made available for microbiological processes. Features within the field such as shaded areas and manure heaps contained significantly higher available nitrogen than the rest of the field. Although these features only represented 1.1% of the area of the field, they contributed to over 55% of the total estimated daily N 2 O flux.

Description: Size-fractionated dissolved primary production and carbohydrate composition of the coccolithophore Emiliania huxleyi Biogeosciences Discussions, 11, 15289-15325, 2014 Author(s): C. Borchard and A. Engel Extracellular release (ER) by phytoplankton is the major source of fresh dissolved organic carbon (DOC) in marine ecosystems and accompanies primary production during all growth phases. Little is known, so far, on size and composition of released molecules, and to which extent ER occurs passively, by leakage, or actively, by exudation. Here, we report on ER by the widespread and bloom-forming coccolithophore Emiliania huxleyi grown under steady state conditions in phosphorus controlled chemostats (N : P = 29, growth rate of μ = 0.2 d −1 ). 14 C incubations were accomplished to determine primary production (PP), comprised by particulate (PO 14 C) and dissolved organic carbon (DO 14 C), and the concentration and composition of particulate combined carbohydrates (pCCHO), and of high molecular weight (〉1 kDa, HMW) dissolved combined carbohydrates (dCCHO) as major components of ER. Information on size distribution of ER products was obtained by investigating distinct size classes (

Description: Autonomous profiling float observations of the high biomass plume downstream of the Kerguelen plateau in the Southern Ocean Biogeosciences Discussions, 11, 17413-17462, 2014 Author(s): M. Grenier, A. Della Penna, and T. W. Trull Natural iron fertilisation from Southern Ocean islands results in high primary production and phytoplankton biomass accumulations readily visible in satellite ocean colour observations. These images reveal great spatial complexity with highly varying concentrations of chlorophyll, presumably reflecting both variations in iron supply and conditions favouring phytoplankton accumulation. To examine the second aspect, in particular the influences of variations in temperature and stratification, we deployed four autonomous profiling floats in the Antarctic Circumpolar Current near the Kerguelen plateau in the Indian sector of the Southern Ocean. Each "bio-profiler" measured more than 250 profiles of temperature ( T ), salinity ( S ), dissolved oxygen, chlorophyll fluorescence (Chl a ), and particle backscatter in the top 300 m of the water column, sampling up to 5 profiles per day along meandering trajectories extending up to 1000 km. Comparison of surface Chl a estimates (top 50 m depth; analogous to values from satellite images) with total water column inventories revealed largely linear relationships, suggesting that dilution of chlorophyll by mixed layer depth variations plays only a minor role in the spatial distributions observed by satellite, and correspondingly that these images provide credible information on total and not just surface biomass accumulations. Regions of very high Chl a accumulation (1.5–10 μg L -1 ) were associated predominantly with a narrow T – S class of surface waters, which appears to derive from the northern Kerguelen plateau. In contrast, waters with only moderate Chl a enrichments (0.5–1.5 μg L -1 ) displayed no clear correlation with water properties, including no dependence on mixed layer depth, suggesting a diversity of sources of iron and/or its efficient dispersion across filaments of the plume. The lack of dependence on mixed layer depth also indicates a limited influence on production by light limitation. One float became trapped in a cyclonic eddy, allowing temporal evaluation of the water column in early autumn. During this period, decreasing surface Chl a inventories corresponded with decreases in oxygen inventories on sub-mixed layer density surfaces, consistent with significant export of organic matter and its respiration and storage as dissolved inorganic carbon in the ocean interior. These results are encouraging for the expanded use of autonomous observing platforms to study biogeochemical, carbon cycle, and ecological problems, although the complex blend of Lagrangian and Eulerian sampling achieved by the floats suggests that arrays rather than single floats will often be required.

Description: On the use of the post-closure method uncertainty band to evaluate the performance of land surface models against eddy covariance flux data Biogeosciences Discussions, 11, 16911-16951, 2014 Author(s): J. Ingwersen, K. Imukova, P. Högy, and T. Streck The energy balance of eddy covariance (EC) flux data is normally not closed. Therefore, at least if used for modeling, EC flux data are usually post-closed, i.e. the measured turbulent fluxes are adjusted so as to close the energy balance. At the current state of knowledge, however, it is not clear how to partition the missing energy in the right way. Eddy flux data therefore contain some uncertainty due to the unknown nature of the energy balance gap, which should be considered in model evaluation and the interpretation of simulation results. We propose to construct the post-closure method uncertainty band (PUB), which essentially designates the differences between non-adjusted flux data and flux data adjusted with the three post-closure methods (Bowen ratio, latent heat flux (LE) and sensible heat flux ( H ) method). To demonstrate this approach, simulations with the NOAH-MP land surface model were evaluated based on EC measurements conducted at a winter wheat stand in Southwest Germany in 2011, and the performance of the Jarvis and Ball–Berry stomatal resistance scheme was compared. The width of the PUB of the LE was up to 110 W m –2 (21% of net radiation). Our study shows that it is crucial to account for the uncertainty of EC flux data originating from lacking energy balance closure. Working with only a single post-closing method might result in severe misinterpretations in model-data comparisons.

Description: Carbon, oxygen and biological productivity in the Southern Ocean in and out the Kerguelen plume: CARIOCA drifter results Biogeosciences Discussions, 11, 16877-16909, 2014 Author(s): L. Merlivat, J. Boutin, and F. d'Ovidio The Kerguelen Plateau region in the Indian sector of the Southern Ocean supports annually a large-scale phytoplankton bloom which is naturally fertilized with iron. As part of the second Kerguelen Ocean and Plateau compared Study expedition (KEOPS2) in austral spring (October–November 2011), one Carioca buoy was deployed east of the Kerguelen plateau. It drifted eastward downstream in the Kerguelen plume. Hourly surface measurements of p CO 2 , O 2 and ancillary observations were collected between 1 November 2011 to 12 February 2012 with the aim of characterizing the spatial and temporal variability of the biological Net Community Production (NCP) downstream the Kerguelen plateau, assess the impact of iron-induced productivity on the biological carbon consumption and consequently on the CO 2 flux exchanged at the air–sea interface. The trajectory of the buoy until mid-December was within the longitude range, 72–83° E, close to the polar front and then in the polar frontal zone, PFZ, until 97° E. From 17 November to 16 December, the buoy drifted within the Kerguelen plume following a filament carrying dissolved iron, DFe, for a total distance of 700 km. In the first part of the trajectory, the ocean surface waters are a sink for CO 2 and a source for CO 2 , with fluxes of respective mean values equal to −8 and +38 mmol CO 2 m −2 d −1 . Eastward, as the buoy escapes the iron enriched filament, the fluxes are in opposite direction, with respective mean values of +5 and −48 mmol O 2 m −2 d −1 . These numbers clearly indicate the strong impact of biological processes on the biogeochemistry in the surface waters within the Kerguelen plume in November-mid-December, while it is undetectable eastward in the PFZ from mid-December to mid-February. While the buoy follows the Fe enriched filament, simultaneous observations of dissolved inorganic carbon, DIC, and dissolved oxygen, O 2 , highlight biological events lasting from 2 to 4 days. Stoichiometric ratios, O 2 /C, between 1.1 and 1.4 are observed indicating new and regenerated production regimes. NCP estimates range from 60 to 140 mmol C m −2 d −1 . Based on the relationship between the time a water parcel has left the plateau and its iron content, we have highlighted that the main control on the value of NCP is the availability of iron in the upper water column, with the largest NCP occurring in waters that have recently left the plateau and presented the highest iron concentrations.

Description: Comparison of UV/Vis and FDOM sensors for in situ monitoring of stream DOC concentrations Biogeosciences Discussions, 11, 16855-16876, 2014 Author(s): G.-Y. Yoo, Y. Jeong, E.-J. Lee, J.-H. Park, and N.-H. Oh Optical measurements using ultra-violet/visible (UV/Vis) spectrophotometric sensors and fluorescent dissolved organic matter (FDOM) sensors have recently been used as proxies of dissolved organic carbon (DOC) concentrations of streams and rivers at high temporal resolution. Despite of the merits of the sensors, temperature changes and particulate matter in water can interfere the sensor readings, over- or under-estimating DOC concentrations. However, little efforts have been made to compare responses of the two types of the sensors in natural conditions. We conducted both laboratory experiments and in situ monitoring with a UV/Vis sensor and a FDOM sensor during the three storm events in the fall of 2012 and the spring of 2013 in a forest stream in Korea in order to compare their performance. Laboratory experiments using the Suwannee River natural organic matter, humic acid, and fulvic acid demonstrated strong linear relationships between both the sensor signals and measured DOC concentrations with R 2 ≥ 0.98. Although temperature compensation might not be needed for the UV/Vis sensor, it was sensitive to relativley small changes in turbidity. In contrast, the FDOM sensor was insenstive to relatively low turbidity while the FDOM sensor outputs decreased significantly as temperature increased, requiring temperature compensated FDOM (e.g. FDOM 20 for 20 °C) for in situ monitoring of DOC. The results suggest that both sensors can be employed as a~proxy for stream DOC concentrations after temperature and turbidity compensation in a forest stream where terrestrially derived humic-like materials are dominant components.

Description: Modeling the impact of riverine DON removal by marine bacterioplankton on primary production in the Arctic Ocean Biogeosciences Discussions, 11, 16953-16992, 2014 Author(s): V. Le Fouest, M. Manizza, B. Tremblay, and M. Babin The planktonic and biogeochemical dynamics of the Arctic shelves exhibit a strong variability in response to Arctic warming. In this study, in order to elucidate on the processes regulating the production of phytoplankton (PP) and bacterioplankton (BP) and their interactions, we employ a biogeochemical model coupled to a pan-Arctic ocean-sea ice model (MITgcm) to explicitly simulate and quantify the contribution of usable dissolved organic nitrogen (DON) drained by the major circum-Arctic rivers on PP and BP in a scenario of melting sea ice (1998–2011). Model simulations suggest that on average between 1998 and 2011, the removal of usable RDON by bacterioplankton is responsible of a ~26% increase of the annual BP for the whole Arctic Ocean. With respect to total PP, the model simulates an increase of ~8% on an annual basis and of ~18% in summer. Recycled ammonium is responsible for the PP increase. The recycling of RDON by bacterioplankton promotes higher BP and PP but there is no significant temporal trend in the BP : PP ratio within the ice-free shelves over the 1998–2011 period. This suggests no significant evolution in the balance between autotrophy and heterotrophy in the last decade with a constant annual flux of RDON into the coastal ocean although changes in RDON supply and further reduction in sea ice cover could potentially alter this delicate balance.

Description: An algorithm for detecting Trichodesmium surface blooms in the South Western Tropical Pacific Biogeosciences Discussions, 8, 5653-5689, 2011 Author(s): C. Dupouy, D. Benielli-Gary, J. Neveux, Y. Dandonneau, and T. K. Westberry Trichodesmium , a major colonial cyanobacterial nitrogen fixer, forms large blooms in NO 3 -depleted tropical oceans and enhances CO 2 sequestration by the ocean due to its ability to fix dissolved dinitrogen. Thus, its importance in C and N cycles requires better estimates of its distribution at basin to global scales. However, existing algorithms to detect them from satellite have not yet been successful in the South Western Tropical Pacific (SWTP). Here, a novel approach based on radiance anomaly spectra (RAS) observed in SeaWiFS imagery is used to detect Trichodesmium during the austral summertime in the SWTP. Selected pixels are characterized by a restricted range of parameters quantifying RAS spectra quantitative parameters (e.g. slope, intercept, curvature). The fraction of valid pixels identified as Trichodesmium surface blooms in the region 5° S–25° S 160° E–190° E is low (between 0.01 and 0.2 %), but is about 100 times higher than suggested by previous algorithms. This represents a total surface area which varies from 1500 to 20 000 km 2 . A monthly distribution of Trichodesmium surface accumulations in the SWTP is presented which demonstrates that the number of selected pixels peaks in November–February each year, consistent with field observations. This approach was validated with in situ observations of Trichodesmium surface accumulations for the period 1998–2010.

Description: Corrigendum to "Breeding and predictability in the baroclinic rotating annulus using a perfect model" published in Nonlin. Processes Geophys., 15, 469–487, 2008 Nonlinear Processes in Geophysics, 18, 359-359, 2011 Author(s): R. M. B. Young and P. L. Read No abstract available.

Description: Rainfall patterns after fire differentially affect the recruitment of three Mediterranean shrubs Biogeosciences Discussions, 8, 5761-5786, 2011 Author(s): J. M. Moreno, E. Zuazua, B. Pérez, B. Luna, A. Velasco, and V. Resco de Dios In fire-prone environments, the "event-dependent hypothesis" states that plant population changes are driven by the unique set of conditions of a fire (e.g., fire season, climate). Climate variability, in particular changes in rainfall patterns, can be most important for seeder species, since they must regenerate after fire from seeds, and for Mediterranean shrublands, given the high yearly variability of rainfall in these ecosystems. Yet, the role of rainfall variability and its interaction with fire characteristics (e.g., fire season) on plant populations has received little attention. Here we investigated the changes in seedling emergence and recruitment of three seeder species ( Cistus ladanifer , Erica umbellata and Rosmarinus officinalis ) after fires lit during three different years and at two times during the fire season (early and late in the fire season) to account for potential changes in the soil seed-bank during the year. Three plots were burned at each season, for a total of 18 plots burned during the three years. After fire, emerged seedlings were tallied, tagged and monitored during three years (two the last burning year). Rainfall during the study period was rather variable, and in some years was well below average. Seedling emergence after fire varied by a factor of 3 to 10, depending on the species and on the burning year. The bulk of seedling emergence occurred in the first year after fire, and seedling recruitment at the end of the study period was tightly correlated with this early emergence. Seedling emergence in E umbellata and R officinalis , but not in C ladanifer , were correlated with precipitation in the fall and winter immediately after the fire, being E umbellata most sensitive to low rainfall. Fire season was generally not an important factor in controlling emergence and recruitment. We discuss how projected changes in rainfall patterns with global warming can alter the balance of species in this shrubland, and can drive some species to near local extinction.

Description: Strategies for coupling global and limited-area ensemble Kalman filter assimilation Nonlinear Processes in Geophysics, 18, 415-430, 2011 Author(s): D. Merkova, I. Szunyogh, and E. Ott This paper compares the forecast performance of four strategies for coupling global and limited area data assimilation: three strategies propagate information from the global to the limited area process, while the fourth strategy feeds back information from the limited area to the global process. All four strategies are formulated in the Local Ensemble Transform Kalman Filter (LETKF) framework. Numerical experiments are carried out with the model component of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) and the NCEP Regional Spectral Model (RSM). The limited area domain is an extended North-America region that includes part of the north-east Pacific. The GFS is integrated at horizontal resolution T62 (about 150 km in the mid-latitudes), while the RSM is integrated at horizontal resolution 48 km. Experiments are carried out both under the perfect model hypothesis and in a realistic setting. The coupling strategies are evaluated by comparing their deterministic forecast performance at 12-h and 48-h lead times. The results suggest that the limited area data assimilation system has the potential to enhance the forecasts at 12-h lead time in the limited area domain at the synoptic and sub-synoptic scales (in the global wave number range of about 10 to 40). There is a clear indication that between the forecast performance of the different coupling strategies those that cycle the limited area assimilation process produce the most accurate forecasts. In the realistic setting, at 12-h forecast time the limited area systems produce more modest improvements compared to the global system than under the perfect model hypothesis, and at 48-h forecast time the global forecasts are more accurate than the limited area forecasts.

Description: Phosphorus transformations as a function of pedogenesis: a synthesis of soil phosphorus data using Hedley fractionation method Biogeosciences Discussions, 8, 5907-5934, 2011 Author(s): X. Yang and W. M. Post In spite of the importance of phosphorus (P) as a limiting nutrient in terrestrial ecosystems, our understanding of terrestrial P dynamics and our ability to model P cycling are hampered by the lack of consistent measurements of soil P. The Hedley fractionation method provides a comprehensive assessment of soil P and has been widely used in recent decades. Here we expand an earlier study that summarized Hedley P data from the literature to create a larger Hedley P database and further investigate the relationships between distributions of different forms of P and the stages of soil development. Our expanded Hedley P database generally supports what the Walker and Syers (1976) conceptual model predicts: the gradual decrease and eventual depletion of primary mineral P (mainly apatite P); the continual increase and eventual dominance of occluded P; and the overall decrease of total P during soil development. However the analysis disagrees with Walker and Syers (1976) in that we found labile inorganic P(Pi) and secondary Pi (non-occluded P in Walker and Syers' model) to be a significant fraction of total P throughout all soil orders with different weathering stages. By analyzing the Hedley-labile P and vegetation P demand, we found that the amount of labile P is much greater than vegetation demand, even in highly weathered soils commonly considered P limited. We conclude that measured Hedley-labile P should not be defined as plant available P since most of this labile P likely ends up as immobilized by microbes. Our analysis of the database also shows that carbon (C) and nitrogen (N) in soil organic matter are closely linked in all soil orders, but P is decoupled from C and N in highly weathered soils with larger variations of nitrogen:organic P (N:Po) ratio and higher mean values of N:Po ratio, compared to slightly and intermediately weathered soils.

Description: Spectral methods for internal waves: indistinguishable density profiles and double-humped solitary waves Nonlinear Processes in Geophysics, 18, 351-358, 2011 Author(s): M. Dunphy, C. Subich, and M. Stastna Internal solitary waves are widely observed in both the oceans and large lakes. They can be described by a variety of mathematical theories, covering the full spectrum from first order asymptotic theory (i.e. Korteweg-de Vries, or KdV, theory), through higher order extensions of weakly nonlinear-weakly nonhydrostatic theory, to fully nonlinear-weakly nonhydrostatic theories and finally exact theory based on the Dubreil-Jacotin-Long (DJL) equation that is formally equivalent to the full set of Euler equations. We discuss how spectral and pseudospectral methods allow for the computation of novel phenomena in both approximate and exact theories. In particular we construct markedly different density profiles for which the coefficients in the KdV theory are very nearly identical. These two density profiles yield qualitatively different behaviour for both exact, or fully nonlinear, waves computed using the DJL equation and in dynamic simulations of the time dependent Euler equations. For exact, DJL, theory we compute exact solitary waves with two-scales, or so-called double-humped waves.

Description: Greenhouse gas flux measurements in a forestry-drained peatland indicate a large carbon sink Biogeosciences Discussions, 8, 5787-5825, 2011 Author(s): A. Lohila, K. Minkkinen, M. Aurela, J.-P. Tuovinen, T. Penttilä, and T. Laurila Drainage for forestry purposes changes the conditions in the peat and leads to increased growth of shrubs and trees. Concurrently, the production and uptake of the greenhouse gases carbon dioxide (CO 2 ), methane (CH 4 ) and nitrous oxide (N 2 O) are likely to change: due to the accelerated decomposition of oxic peat, drained peatlands are generally considered to loose peat carbon (C). We measured CO 2 exchange with the eddy covariance (EC) method above a drained nutrient-poor peatland forest in Southern Finland for 16 months in 2004–2005. The site, classified as a dwarf-shrub pine bog, had been ditched about 35 years earlier. CH 4 and N 2 O fluxes were measured at 2–5 week intervals with the chamber technique. Drainage had resulted in a relatively little change in the water table level, being on average 40 cm below the ground in 2005. The annual net ecosystem exchange was −870 g CO 2 m −2 yr −1 in the calendar year 2005, varying from −810 to −900 g CO 2 m −2 yr −1 during the 16 month period under investigation. The site was a small sink of CH 4 (−0.12 g CH 4 m −2 yr −1 ) and a small source of N 2 O (0.10 g N 2 O m −2 yr −1 ). Photosynthesis was detected throughout the year when the air temperature exceeded −3 °C. As the annual accumulation of C in the above and below ground tree biomass (550 g CO 2 m −2 ) was significantly less than the net exchange of CO 2 , about 300 g CO 2 m −2 yr −1 (~80 g C m −2 ) was likely to have accumulated as organic matter into the peat soil. This is a higher average accumulation rate than previously reported for natural northern peatlands, and the first time C accumulation has been shown, by EC measurements, to occur in a drained peatland. Our results suggest that forestry-drainage may significantly increase the CO 2 uptake rate of nutrient-poor peatland ecosystems.

Description: A nonlinear method of removing harmonic noise in geophysical data Nonlinear Processes in Geophysics, 18, 367-379, 2011 Author(s): Y. Jeng and C.-S. Chen A nonlinear, adaptive method to remove the harmonic noise that commonly resides in geophysical data is proposed in this study. This filtering method is based on the ensemble empirical mode decomposition algorithm in conjunction with the logarithmic transform. We present a synthetic model study to investigate the capability of signal reconstruction from the decomposed data, and compare the results with those derived from other 2-D adaptive filters. Applications to the real seismic data acquired by using an ocean bottom seismograph and to a shot gather of the ground penetrating radar demonstrate the robustness of this method. Our work proposes a concept that instead of Fourier-based approaches, the harmonic noise removal in geophysical data can be achieved effectively by using an alternative nonlinear adaptive data analysis method, which has been applied extensively in other scientific studies.

Description: Effect of mosaic representation of vegetation in land surface schemes on simulated energy and carbon balances Biogeosciences Discussions, 8, 5849-5879, 2011 Author(s): R. Li and V. K. Arora Energy and carbon balance implications of representing vegetation using a composite or mosaic approach in a land surface scheme are investigated. In the composite approach the attributes of different plant functional types (PFTs) present in a grid cell are aggregated in some fashion for energy and water balance calculations. The resulting physical environmental conditions (including net radiation, soil moisture and soil temperature) are common to all PFTs and affect their ecosystem processes. In the mosaic approach energy and water balance calculations are performed separately for each PFT tile using its own vegetation attributes, so each PFT "sees" different physical environmental conditions and its carbon balance evolves somewhat differently from that in the composite approach. Simulations are performed at selected boreal, temperate and tropical locations to illustrate the differences caused by using the composite versus the mosaic approaches of representing vegetation. Differences in grid averaged primary energy fluxes are generally less than 5 % between the two approaches. Grid-averaged carbon fluxes and pool sizes can, however, differ by as much as 46 %. Simulation results suggest that differences in carbon balance between the two approaches arise primarily through differences in net radiation which directly affects net primary productivity, and thus leaf area index and vegetation biomass.

Description: Comparison of correlation analysis techniques for irregularly sampled time series Nonlinear Processes in Geophysics, 18, 389-404, 2011 Author(s): K. Rehfeld, N. Marwan, J. Heitzig, and J. Kurths Geoscientific measurements often provide time series with irregular time sampling, requiring either data reconstruction (interpolation) or sophisticated methods to handle irregular sampling. We compare the linear interpolation technique and different approaches for analyzing the correlation functions and persistence of irregularly sampled time series, as Lomb-Scargle Fourier transformation and kernel-based methods. In a thorough benchmark test we investigate the performance of these techniques. All methods have comparable root mean square errors (RMSEs) for low skewness of the inter-observation time distribution. For high skewness, very irregular data, interpolation bias and RMSE increase strongly. We find a 40 % lower RMSE for the lag-1 autocorrelation function (ACF) for the Gaussian kernel method vs. the linear interpolation scheme,in the analysis of highly irregular time series. For the cross correlation function (CCF) the RMSE is then lower by 60 %. The application of the Lomb-Scargle technique gave results comparable to the kernel methods for the univariate, but poorer results in the bivariate case. Especially the high-frequency components of the signal, where classical methods show a strong bias in ACF and CCF magnitude, are preserved when using the kernel methods. We illustrate the performances of interpolation vs. Gaussian kernel method by applying both to paleo-data from four locations, reflecting late Holocene Asian monsoon variability as derived from speleothem δ 18 O measurements. Cross correlation results are similar for both methods, which we attribute to the long time scales of the common variability. The persistence time (memory) is strongly overestimated when using the standard, interpolation-based, approach. Hence, the Gaussian kernel is a reliable and more robust estimator with significant advantages compared to other techniques and suitable for large scale application to paleo-data.

Description: Contributions of agricultural plants and soils to N 2 O emission in a farmland Biogeosciences Discussions, 8, 5505-5535, 2011 Author(s): J. Li, X. Lee, Q. Yu, X. Tong, Z. Qin, and B. Macdonald The goal of this study was to quantify the roles of plants and soil in the N 2 O budget of a cropland in North China. Plant and soil N 2 O fluxes were measured with transparent and dark plant chambers and soil chambers, respectively, in three adjacent fields of fertilized cotton, fertilized maize and unfertilized soybean. During the observation period, the soil flux was 448 ± 89, 230 ± 74 and 90 ± 14 μg N 2 O m −2 h −1 in cotton, maize and soybean fields, respectively. The plant flux was 54 ± 43 and 16 ± 41 μg N 2 O m −2 h −1 , about 10 % and 26 % to the total ecosystem flux, for the cotton and the soybean field, respectively. Ignoring the contribution of plants would cause an obvious underestimation on the ecosystem N 2 O flux. The influence of sunlight on plant N 2 O flux was insignificant. However, in the cotton field, the responses of the plant N 2 O flux to air temperature and soil ammonium content were significant under sunlight but insignificant under darkness, suggesting that stomatal activity might influence the release process. In the cotton field, temperature sensitivity of plant N 2 O emission was 1.13, much lower than the value of soil flux (5.74). No relationship was found between plant N 2 O flux and soil nitrate content. It was implied that nitrate reduction in plants might not be the main source of plant N 2 O emission under field conditions. The seasonal patterns of the soil and plant N 2 O emissions were similarly affected by fertilization, indicating that plants might serve as a passive conduit transporting N 2 O produced in the soil.