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: Calibration of the passive cavity aerosol spectrometer probe for airborne determination of the size distribution Atmospheric Measurement Techniques, 6, 2349-2358, 2013 Author(s): Y. Cai, J. R. Snider, and P. Wechsler This work describes calibration methods for the particle sizing and particle concentration systems of the passive cavity aerosol spectrometer probe (PCASP). Laboratory calibrations conducted over six years, in support of the deployment of a PCASP on a cloud physics research aircraft, are analyzed. Instead of using the many calibration sizes recommended by the PCASP manufacturer, a relationship between particle diameter and scattered light intensity is established using three sizes of mobility-selected polystyrene latex particles, one for each amplifier gain stage. In addition, studies of two factors influencing the PCASP's determination of the particle size distribution – amplifier baseline and particle shape – are conducted. It is shown that the PCASP-derived size distribution is sensitive to adjustments of the sizing system's baseline voltage, and that for aggregates of spheres, a PCASP-derived particle size and a sphere-equivalent particle size agree within uncertainty dictated by the PCASP's sizing resolution. Robust determinations of aerosol concentration, and size distribution, also require calibration of the PCASP's aerosol flowrate sensor. Sensor calibrations, calibration drift, and the sensor's non-linear response are documented.

Description: Experimental quantification of contact freezing in an electrodynamic balance Atmospheric Measurement Techniques, 6, 2373-2382, 2013 Author(s): N. Hoffmann, A. Kiselev, D. Rzesanke, D. Duft, and T. Leisner Heterogeneous nucleation of ice in a supercooled water droplet induced by external contact with a dry aerosol particle has long been known to be more effective than freezing induced by the same nucleus immersed in the droplet. However, the experimental quantification of contact freezing is challenging. Here we report an experimental method to determine the temperature-dependent ice nucleation probability of size-selected aerosol particles. The method is based on the suspension of supercooled charged water droplets in a laminar flow of air containing aerosol particles as contact freezing nuclei. The rate of droplet–particle collisions is calculated numerically with account for Coulomb attraction, drag force and induced dipole interaction between charged droplet and aerosol particles. The calculation is verified by direct counting of aerosol particles collected by a levitated droplet. By repeating the experiment on individual droplets for a sufficient number of times, we are able to reproduce the statistical freezing behavior of a large ensemble of supercooled droplets and measure the average rate of freezing events. The freezing rate is equal to the product of the droplet–particle collision rate and the probability of freezing on a single contact, the latter being a function of temperature, size and composition of the contact ice nuclei. Based on these observations, we show that for the types of particles investigated so far, contact freezing is the dominating freezing mechanism on the timescale of our experiment.

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: Can we determine what controls the spatio-temporal distribution of d-excess and 17 O-excess in precipitation using the LMDZ general circulation model? Climate of the Past, 9, 2173-2193, 2013 Author(s): C. Risi, A. Landais, R. Winkler, and F. Vimeux Combined measurements of the H 2 18 O and HDO isotopic ratios in precipitation, leading to second-order parameter D-excess, have provided additional constraints on past climates compared to the H 2 18 O isotopic ratio alone. More recently, measurements of H 2 17 O have led to another second-order parameter: 17 O-excess. Recent studies suggest that 17 O-excess in polar ice may provide information on evaporative conditions at the moisture source. However, the processes controlling the spatio-temporal distribution of 17 O-excess are still far from being fully understood. We use the isotopic general circulation model (GCM) LMDZ to better understand what controls d-excess and 17 O-excess in precipitation at present-day (PD) and during the last glacial maximum (LGM). The simulation of D-excess and 17 O-excess is evaluated against measurements in meteoric water, water vapor and polar ice cores. A set of sensitivity tests and diagnostics are used to quantify the relative effects of evaporative conditions (sea surface temperature and relative humidity), Rayleigh distillation, mixing between vapors from different origins, precipitation re-evaporation and supersaturation during condensation at low temperature. In LMDZ, simulations suggest that in the tropics convective processes and rain re-evaporation are important controls on precipitation D-excess and 17 O-excess. In higher latitudes, the effect of distillation, mixing between vapors from different origins and supersaturation are the most important controls. For example, the lower d-excess and 17 O-excess at LGM simulated at LGM are mainly due to the supersaturation effect. The effect of supersaturation is however very sensitive to a parameter whose tuning would require more measurements and laboratory experiments. Evaporative conditions had previously been suggested to be key controlling factors of d-excess and 17 O-excess, but LMDZ underestimates their role. More generally, some shortcomings in the simulation of 17 O-excess by LMDZ suggest that general circulation models are not yet the perfect tool to quantify with confidence all processes controlling 17 O-excess.

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: Validation and empirical correction of MODIS AOT and AE over ocean Atmospheric Measurement Techniques, 6, 2455-2475, 2013 Author(s): N. A. J. Schutgens, M. Nakata, and T. Nakajima We present a validation study of Collection 5 MODIS level 2 Aqua and Terra AOT (aerosol optical thickness) and AE (Ångström exponent) over ocean by comparison to coastal and island AERONET (AErosol RObotic NETwork) sites for the years 2003–2009. We show that MODIS (MODerate-resolution Imaging Spectroradiometer) AOT exhibits significant biases due to wind speed and cloudiness of the observed scene, while MODIS AE, although overall unbiased, exhibits less spatial contrast on global scales than the AERONET observations. The same behaviour can be seen when MODIS AOT is compared against Maritime Aerosol Network (MAN) data, suggesting that the spatial coverage of our datasets does not preclude global conclusions. Thus, we develop empirical correction formulae for MODIS AOT and AE that significantly improve agreement of MODIS and AERONET observations. We show these correction formulae to be robust. Finally, we study random errors in the corrected MODIS AOT and AE and show that they mainly depend on AOT itself, although small contributions are present due to wind speed and cloud fraction in AOT random errors and due to AE and cloud fraction in AE random errors. Our analysis yields significantly higher random AOT errors than the official MODIS error estimate (0.03 + 0.05 τ), while random AE errors are smaller than might be expected. This new dataset of bias-corrected MODIS AOT and AE over ocean is intended for aerosol model validation and assimilation studies, but also has consequences as a stand-alone observational product. For instance, the corrected dataset suggests that much less fine mode aerosol is transported across the Pacific and Atlantic oceans.

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: 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: A neural network algorithm for cloud fraction estimation using NASA-Aura OMI VIS radiance measurements Atmospheric Measurement Techniques, 6, 2301-2309, 2013 Author(s): G. Saponaro, P. Kolmonen, J. Karhunen, J. Tamminen, and G. de Leeuw The discrimination of cloudy from cloud-free pixels is required in almost any estimate of a parameter retrieved from satellite data in the ultraviolet (UV), visible (VIS) or infrared (IR) parts of the electromagnetic spectrum. In this paper we report on the development of a neural network (NN) algorithm to estimate cloud fractions using radiances measured at the top of the atmosphere with the NASA-Aura Ozone Monitoring Instrument (OMI). We present and discuss the results obtained from the application of two different types of neural networks, i.e., extreme learning machine (ELM) and back propagation (BP). The NNs were trained with an OMI data sets existing of six orbits, tested with three other orbits and validated with another two orbits. The results were evaluated by comparison with cloud fractions available from the MODerate Resolution Imaging Spectrometer (MODIS) flying on Aqua in the same constellation as OMI, i.e., with minimal time difference between the OMI and MODIS observations. The results from the ELM and BP NNs are compared. They both deliver cloud fraction estimates in a fast and automated way, and they both performs generally well in the validation. However, over highly reflective surfaces, such as desert, or in the presence of dust layers in the atmosphere, the cloud fractions are not well predicted by the neural network. Over ocean the two NNs work equally well, but over land ELM performs better.

Description: Validation of stratospheric and mesospheric ozone observed by SMILES from International Space Station Atmospheric Measurement Techniques, 6, 2311-2338, 2013 Author(s): Y. Kasai, H. Sagawa, D. Kreyling, E. Dupuy, P. Baron, J. Mendrok, K. Suzuki, T. O. Sato, T. Nishibori, S. Mizobuchi, K. Kikuchi, T. Manabe, H. Ozeki, T. Sugita, M. Fujiwara, Y. Irimajiri, K. A. Walker, P. F. Bernath, C. Boone, G. Stiller, T. von Clarmann, J. Orphal, J. Urban, D. Murtagh, E. J. Llewellyn, D. Degenstein, A. E. Bourassa, N. D. Lloyd, L. Froidevaux, M. Birk, G. Wagner, F. Schreier, J. Xu, P. Vogt, T. Trautmann, and M. Yasui We observed ozone (O 3 ) in the vertical region between 250 and 0.0005 hPa (~ 12–96 km) using the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on the Japanese Experiment Module (JEM) of the International Space Station (ISS) between 12 October 2009 and 21 April 2010. The new 4 K superconducting heterodyne receiver technology of SMILES allowed us to obtain a one order of magnitude better signal-to-noise ratio for the O 3 line observation compared to past spaceborne microwave instruments. The non-sun-synchronous orbit of the ISS allowed us to observe O 3 at various local times. We assessed the quality of the vertical profiles of O 3 in the 100–0.001 hPa (~ 16–90 km) region for the SMILES NICT Level 2 product version 2.1.5. The evaluation is based on four components: error analysis; internal comparisons of observations targeting three different instrumental setups for the same O 3 625.371 GHz transition; internal comparisons of two different retrieval algorithms; and external comparisons for various local times with ozonesonde, satellite and balloon observations (ENVISAT/MIPAS, SCISAT/ACE-FTS, Odin/OSIRIS, Odin/SMR, Aura/MLS, TELIS). SMILES O 3 data have an estimated absolute accuracy of better than 0.3 ppmv (3%) with a vertical resolution of 3–4 km over the 60 to 8 hPa range. The random error for a single measurement is better than the estimated systematic error, being less than 1, 2, and 7%, in the 40–1, 80–0.1, and 100–0.004 hPa pressure regions, respectively. SMILES O 3 abundance was 10–20% lower than all other satellite measurements at 8–0.1 hPa due to an error arising from uncertainties of the tangent point information and the gain calibration for the intensity of the spectrum. SMILES O 3 from observation frequency Band-B had better accuracy than that from Band-A. A two month period is required to accumulate measurements covering 24 h in local time of O 3 profile. However such a dataset can also contain variation due to dynamical, seasonal, and latitudinal effects.

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 Atmospheric radiation measurement (ARM) program network of microwave radiometers: instrumentation, data, and retrievals Atmospheric Measurement Techniques, 6, 2359-2372, 2013 Author(s): M. P. Cadeddu, J. C. Liljegren, and D. D. Turner The Climate Research Facility of the US Department of Energy's Atmospheric Radiation Measurement (ARM) Program operates a network of ground-based microwave radiometers. Data and retrievals from these instruments have been available to the scientific community for almost 20 yr. In the past five years the network has expanded to include a total of 22 microwave radiometers deployed in various locations around the world. The new instruments cover a frequency range between 22 and 197 GHz and are consistently and automatically calibrated. The latest addition to the network is a new generation of three-channel radiometers, currently in the early stage of deployment at all ARM sites. The network has been specifically designed to achieve increased accuracy in the retrieval of precipitable water vapor (PWV) and cloud liquid water path (LWP) with the long-term goal of providing the scientific community with reliable, calibrated radiometric data and retrievals of important geophysical quantities with well-characterized uncertainties. The radiometers provide high-quality, continuous datasets that can be utilized in a wealth of applications and scientific studies. This paper presents an overview of the microwave instrumentation, calibration procedures, data, and retrievals that are available for download from the ARM data archive.

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: Airborne sun photometer PLASMA: concept, measurements, comparison of aerosol extinction vertical profile with lidar Atmospheric Measurement Techniques, 6, 2383-2389, 2013 Author(s): Y. Karol, D. Tanré, P. Goloub, C. Vervaerde, J. Y. Balois, L. Blarel, T. Podvin, A. Mortier, and A. Chaikovsky A 15-channel airborne sun-tracking photometer has been developed. The instrument provides aerosol optical depths over a wide spectral range (0.34–2.25 μm) with an accuracy (ΔAOD) of approximately 0.01. Taking measurements at different altitudes allow us to derive the aerosol extinction vertical profile. Thanks to the wide spectral range of the instrument, information on the aerosol size distribution along the vertical is also available.

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: McClear: a new model estimating downwelling solar radiation at ground level in clear-sky conditions Atmospheric Measurement Techniques, 6, 2403-2418, 2013 Author(s): M. Lefèvre, A. Oumbe, P. Blanc, B. Espinar, B. Gschwind, Z. Qu, L. Wald, M. Schroedter-Homscheidt, C. Hoyer-Klick, A. Arola, A. Benedetti, J. W. Kaiser, and J.-J. Morcrette A new fast clear-sky model called McClear was developed to estimate the downwelling shortwave direct and global irradiances received at ground level under clear skies. It is a fully physical model replacing empirical relations or simpler models used before. It exploits the recent results on aerosol properties, and total column content in water vapour and ozone produced by the MACC project (Monitoring Atmosphere Composition and Climate). It accurately reproduces the irradiance computed by the libRadtran reference radiative transfer model with a computational speed approximately 10 5 times greater by adopting the abaci, or look-up table, approach combined with interpolation functions. It is therefore suited for geostationary satellite retrievals or numerical weather prediction schemes with many pixels or grid points, respectively. McClear irradiances were compared to 1 min measurements made in clear-sky conditions at several stations within the Baseline Surface Radiation Network in various climates. The bias for global irradiance comprises between −6 and 25 W m −2 . The RMSE ranges from 20 W m −2 (3% of the mean observed irradiance) to 36 W m −2 (5%) and the correlation coefficient ranges between 0.95 and 0.99. The bias for the direct irradiance comprises between −48 and +33 W m −2 . The root mean square error (RMSE) ranges from 33 W m −2 (5%) to 64 W m −2 (10%). The correlation coefficient ranges between 0.84 and 0.98. This work demonstrates the quality of the McClear model combined with MACC products, and indirectly the quality of the aerosol properties modelled by the MACC reanalysis.

Description: Ten years of MIPAS measurements with ESA Level 2 processor V6 – Part 1: Retrieval algorithm and diagnostics of the products Atmospheric Measurement Techniques, 6, 2419-2439, 2013 Author(s): P. Raspollini, B. Carli, M. Carlotti, S. Ceccherini, A. Dehn, B. M. Dinelli, A. Dudhia, J.-M. Flaud, M. López-Puertas, F. Niro, J. J. Remedios, M. Ridolfi, H. Sembhi, L. Sgheri, and T. von Clarmann The MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) instrument on the Envisat (Environmental satellite) satellite has provided vertical profiles of the atmospheric composition on a global scale for almost ten years. The MIPAS mission is divided in two phases: the full resolution phase, from 2002 to 2004, and the optimized resolution phase, from 2005 to 2012, which is characterized by a finer vertical and horizontal sampling attained through a reduction of the spectral resolution. While the description and characterization of the products of the ESA processor for the full resolution phase has been already described in previous papers, in this paper we focus on the performances of the latest version of the ESA (European Space Agency) processor, named ML2PP V6 (MIPAS Level 2 Prototype Processor), which has been used for reprocessing the entire mission. The ESA processor had to perform the operational near real time analysis of the observations and its products needed to be available for data assimilation. Therefore, it has been designed for fast, continuous and automated analysis of observations made in quite different atmospheric conditions and for a minimum use of external constraints in order to avoid biases in the products. The dense vertical sampling of the measurements adopted in the second phase of the MIPAS mission resulted in sampling intervals finer than the instantaneous field of view of the instrument. Together with the choice of a retrieval grid aligned with the vertical sampling of the measurements, this made ill-conditioned the retrieval problem of the MIPAS operational processor. This problem has been handled with minimal changes to the original retrieval approach but with significant improvements nonetheless. The Levenberg–Marquardt method, already present in the retrieval scheme for its capability to provide fast convergence for nonlinear problems, is now also exploited for the reduction of the ill-conditioning of the inversion. An expression specifically designed for the regularizing Levenberg–Marquardt method has been implemented for the computation of the covariance matrices and averaging kernels of the retrieved products. The regularization of the Levenberg–Marquardt method is controlled by the convergence criteria and is deliberately kept weak. The resulting oscillations of the retrieved profile are a posteriori damped by an innovative self-adapting Tikhonov regularization. The convergence criteria and the weakness of the self-adapting regularization ensure that minimum constraints are used and the best vertical resolution obtainable from the measurements is achieved in all atmospheric conditions. Random and systematic errors, as well as vertical and horizontal resolution are compared in the two phases of the mission for all products, namely: temperature, H 2O , O 3 , HNO 3 , CH 4 , N 2O , NO 2 , CFC-11, CFC-12, N 2 O 5 and ClONO 2 . The use in the two phases of the mission of different optimized sets of spectral intervals ensures that, despite the different spectral resolutions, comparable performances are obtained in the whole MIPAS mission in terms of random and systematic errors, while the vertical resolution and the horizontal resolution are significantly better in the case of the optimized resolution measurements.

Description: Ground-based stratospheric O 3 and HNO 3 measurements at Thule, Greenland: an intercomparison with Aura MLS observations Atmospheric Measurement Techniques, 6, 2441-2453, 2013 Author(s): I. Fiorucci, G. Muscari, L. Froidevaux, and M. L. Santee In response to the need for improving our understanding of the evolution and the interannual variability of the winter Arctic stratosphere, in January 2009 a Ground-Based Millimeter-wave Spectrometer (GBMS) was installed at the Network for the Detection of Atmospheric Composition Change (NDACC) site in Thule (76.5° N, 68.8° W), Greenland. In this work, stratospheric GBMS O 3 and HNO 3 vertical profiles obtained from Thule during the winters 2010 (HNO 3 only), 2011 and 2012 are characterized and intercompared with co-located measurements of the Aura Microwave Limb Sounder (MLS) experiment. Using a recently developed algorithm based on Optimal Estimation, we find that the GBMS O 3 retrievals show good sensitivity (〉 80%) to atmospheric variations between ~ 17 and ~ 50 km, where their 1σ uncertainty is estimated to be the larger of ~ 11% or 0.2 ppmv. Similarly, HNO 3 profiles can be considered for scientific use between ~ 17 and ~ 45 km altitude, with a 1σ uncertainty that amounts to the larger of 15% or 0.2 ppbv. Comparisons with Aura MLS version 3.3 observations show that, on average, GBMS O 3 mixing ratios are biased negatively with respect to MLS throughout the stratosphere, with differences ranging between ~ 0.3 ppmv (8%) and 0.9 ppmv (18%) in the 17–50 km vertical range. GBMS HNO 3 values display instead a positive bias with respect to MLS up to 26 km, reaching a maximum of ~ 1 ppbv (10%) near the mixing ratio profile peak. O 3 and HNO 3 values from the two datasets prove to be well correlated at all altitudes, although their correlations worsen at the lower end of the altitude ranges considered. Column contents of GBMS and MLS O 3 (from 20 km upwards) and HNO 3 (from 17 km upwards) correlate very well and indicate that GBMS measurements can provide valuable estimates of column interannual and seasonal variations for these compounds.

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: Inter-calibration of polar imager solar channels using SEVIRI Atmospheric Measurement Techniques, 6, 2495-2508, 2013 Author(s): J. F. Meirink, R. A. Roebeling, and P. Stammes Accurate calibration of satellite imagers is a prerequisite for using their measurements in climate applications. Here we present a method for the inter-calibration of geostationary and polar-orbiting imager solar channels based on regressions of collocated near-nadir reflectances. Specific attention is paid to correcting for differences in spectral response between instruments. The method is used to calibrate the solar channels of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on the geostationary Meteosat satellite with corresponding channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the polar-orbiting Aqua satellite. The SEVIRI operational calibration is found to be stable during the years 2004 to 2009, but offset by −8, −6, and +3.5 % for channels 1 (0.6 μm), 2 (0.8 μm), and 3 (1.6 μm), respectively. These results are robust for a range of choices that can be made regarding data collocation and selection, as long as the viewing and illumination geometries of the two instruments are matched. Uncertainties in the inter-calibration method are estimated to be 1 % for channel 1 and 1.5 % for channels 2 and 3. A specific application of our method is the inter-calibration of polar imagers using SEVIRI as a transfer instrument. This offers an alternative to direct inter-calibration, which in general has to rely on high-latitude collocations. Using this method we have tied MODIS-Terra and Advanced Very High Resolution Radiometer (AVHRR) instruments on National Oceanic and Atmospheric Administration (NOAA) satellites 17 and 18 to MODIS-Aqua for the years 2007 to 2009. While reflectances of the two MODIS instruments differ less than 2 % for all channels considered, deviations of an existing AVHRR calibration from MODIS-Aqua reach −3.5 and +2.5 % for the 0.8 and 1.6 μm channels, respectively.

Description: Using ocean-glint scattered sunlight as a diagnostic tool for satellite remote sensing of greenhouse gases Atmospheric Measurement Techniques, 6, 2509-2520, 2013 Author(s): A. Butz, S. Guerlet, O. P. Hasekamp, A. Kuze, and H. Suto Spectroscopic measurements of sunlight backscattered by the Earth's surface is a technique widely used for remote sensing of atmospheric constituent concentrations from space. Thereby, remote sensing of greenhouse gases poses particularly challenging accuracy requirements for instrumentation and retrieval algorithms which, in general, suffer from various error sources. Here, we investigate a method that helps disentangle sources of error for observations of sunlight backscattered from the glint spot on the ocean surface. The method exploits the backscattering characteristics of the ocean surface, which is bright for glint geometry but dark for off-glint angles. This property allows for identifying a set of clean scenes where light scattering due to particles in the atmosphere is negligible such that uncertain knowledge of the lightpath can be excluded as a source of error. We apply the method to more than 3 yr of ocean-glint measurements by the Thermal And Near infrared Sensor for carbon Observation (TANSO) Fourier Transform Spectrometer (FTS) onboard the Greenhouse Gases Observing Satellite (GOSAT), which aims at measuring carbon dioxide (CO 2 ) and methane (CH 4 ) concentrations. The proposed method is able to clearly monitor recent improvements in the instrument calibration of the oxygen (O 2 ) A-band channel and suggests some residual uncertainty in our knowledge about the instrument. We further assess the consistency of CO 2 retrievals from several absorption bands between 6400 cm −1 (1565 nm) and 4800 cm −1 (2100 nm) and find that the absorption bands commonly used for monitoring of CO 2 dry air mole fractions from GOSAT allow for consistency better than 1.5 ppm. Usage of other bands reveals significant inconsistency among retrieved CO 2 concentrations pointing at inconsistency of spectroscopic parameters.

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: 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: Global and regional sea surface temperature trends during Marine Isotope Stage 11 Climate of the Past, 9, 2231-2252, 2013 Author(s): Y. Milker, R. Rachmayani, M. F. G. Weinkauf, M. Prange, M. Raitzsch, M. Schulz, and M. Kučera The Marine Isotope Stage (MIS) 11 (424–374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic ( N = 28) or planktonic ( N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3–6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties 〉1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

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: Parameterizing radiative transfer to convert MAX-DOAS dSCDs into near-surface box-averaged mixing ratios Atmospheric Measurement Techniques, 6, 1521-1532, 2013 Author(s): R. Sinreich, A. Merten, L. Molina, and R. Volkamer We present a novel parameterization method to convert multi-axis differential optical absorption spectroscopy (MAX-DOAS) differential slant column densities (dSCDs) into near-surface box-averaged volume mixing ratios. The approach is applicable inside the planetary boundary layer under conditions with significant aerosol load, and builds on the increased sensitivity of MAX-DOAS near the instrument altitude. It parameterizes radiative transfer model calculations and significantly reduces the computational effort, while retrieving ~ 1 degree of freedom. The biggest benefit of this method is that the retrieval of an aerosol profile, which usually is necessary for deriving a trace gas concentration from MAX-DOAS dSCDs, is not needed. The method is applied to NO 2 MAX-DOAS dSCDs recorded during the Mexico City Metropolitan Area 2006 (MCMA-2006) measurement campaign. The retrieved volume mixing ratios of two elevation angles (1° and 3°) are compared to volume mixing ratios measured by two long-path (LP)-DOAS instruments located at the same site. Measurements are found to agree well during times when vertical mixing is expected to be strong. However, inhomogeneities in the air mass above Mexico City can be detected by exploiting the different horizontal and vertical dimensions probed by the MAX-DOAS and LP-DOAS instruments. In particular, a vertical gradient in NO 2 close to the ground can be observed in the afternoon, and is attributed to reduced mixing coupled with near-surface emission inside street canyons. The existence of a vertical gradient in the lower 250 m during parts of the day shows the general challenge of sampling the boundary layer in a representative way, and emphasizes the need of vertically resolved measurements.

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: Polarization data from SCIAMACHY limb backscatter observations compared to vector radiative transfer model simulations Atmospheric Measurement Techniques, 6, 1503-1520, 2013 Author(s): P. Liebing, K. Bramstedt, S. Noël, V. Rozanov, H. Bovensmann, and J. P. Burrows SCIAMACHY is a passive imaging spectrometer onboard ENVISAT designed to obtain trace gas abundances from measured radiances and irradiances in the UV to SWIR range in nadir-, limb- and occultation-viewing modes. Its grating spectrometer introduces a substantial sensitivity to the polarization of the incoming light with nonnegligible effects on the radiometric calibration. To be able to correct for the polarization sensitivity, SCIAMACHY utilizes broadband Polarization Measurement Devices (PMDs). While for the nadir-viewing mode the measured atmospheric polarization has been validated against POLDER data (Tilstra and Stammes, 2007, 2010), a similar validation study regarding the limb-viewing mode has not yet been performed. This paper aims at an assessment of the quality of the SCIAMACHY limb polarization data. Since limb polarization measurements by other air/spaceborne instruments in the spectral range of SCIAMACHY are not available, a comparison with radiative transfer simulations by SCIATRAN V3.1 (Rozanov et al., 2013) using a wide range of atmospheric parameters is performed. SCIATRAN is a vector radiative transfer model (VRTM) capable of performing calculations of the multiply scattered radiance in a spherically symmetric atmosphere. The study shows that the limb polarization data exhibit a large time-dependent bias that decreases with wavelength. Possible reasons for this bias are a still unknown combination of insufficient accuracy or inconsistencies of the on-ground calibration data, scan mirror degradation and stress induced changes of the polarization response of components inside the optical bench of the instrument. It is shown that it should in principle be feasible to recalibrate the effective polarization sensitivity of the instrument using the in-flight data and VRTM simulations.

Description: Evaluating calibration strategies for isotope ratio infrared spectroscopy for atmospheric 13 CO 2 / 12 CO 2 measurement Atmospheric Measurement Techniques, 6, 1491-1501, 2013 Author(s): X.-F. Wen, Y. Meng, X.-Y. Zhang, X.-M. Sun, and X. Lee Isotope ratio infrared spectroscopy (IRIS) provides an in situ technique for measuring δ 13 C in atmospheric CO 2 . A number of methods have been proposed for calibrating the IRIS measurements, but few studies have systematically evaluated their accuracy for atmospheric applications. In this study, we carried out laboratory and ambient measurements with two commercial IRIS analyzers and compared the accuracy of four calibration strategies. We found that calibration based on the 12 C and 13 C mixing ratios (Bowling et al., 2003) and on linear interpolation of the measured delta using the mixing ratio of the major isotopologue (Lee et al., 2005) yielded accuracy better than 0.06‰. Over a 7-day atmospheric measurement in Beijing, the two analyzers agreed to within −0.02 ± 0.18‰ after proper calibration. However, even after calibration the difference between the two analyzers showed a slight correlation with concentration, and this concentration dependence propagated through the Keeling analysis, resulting in a much larger difference of 2.44‰ for the Keeling intercept. The high sensitivity of the Keeling analysis to the concentration dependence underscores the challenge of IRIS for atmospheric research.

Description: Cavity ring-down spectroscopy sensor development for high-time-resolution measurements of gaseous elemental mercury in ambient air Atmospheric Measurement Techniques, 6, 1477-1489, 2013 Author(s): A. Pierce, D. Obrist, H. Moosmüller, X. Faïn, and C. Moore We describe further development of a previous laboratory prototype pulsed cavity ring-down spectroscopy (CRDS) sensor into a field-deployable system for high-time-resolution, continuous, and automated measurement of gaseous elemental mercury (GEM) concentrations in ambient air. We employed an external, isotopically enriched Hg cell for automated locking and stabilization of the laser wavelength on the GEM peak absorption during measurements. Further, we describe implementation of differential absorption measurements via a piezoelectric tuning element for pulse-by-pulse tuning of the laser wavelength onto and off of the GEM absorption line. This allowed us to continuously correct (at 25 Hz) for system baseline extinction losses unrelated to GEM absorption. Extensive measurement and calibration data obtained with the system were based on spike addition in both GEM-free air and ambient air. Challenges and interferences that occurred during measurements (particularly in ambient air) are discussed including temperature and ozone (O 3 ) concentration fluctuations, and steps taken to reduce these. CRDS data were highly linear ( r 2 ≥ 0.98) with data from a commercial Tekran 2537 Hg analyzer across a wide range of GEM concentrations (0 to 127 ng m −3 ) in Hg-free and ambient air. Measurements during periods of stable background GEM concentrations provided a conservative instrument sensitivity estimate of 0.35 ng m −3 for the CRDS system when time averaged for 5 min. This sensitivity, along with concentration patterns observed in ambient air (with the CRDS system and verified with the Tekran analyzer), showed that the sensor was capable of characterizing GEM fluctuations in ambient air. The value of fast-response GEM measurements was shown by a series of GEM spike additions – highlighting that high-temporal-resolution measurement allowed for detailed characterization of fast concentration fluctuations not possible with traditional analyzers.

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: 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: 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: Skill and reliability of climate model ensembles at the Last Glacial Maximum and mid-Holocene Climate of the Past, 9, 811-823, 2013 Author(s): J. C. Hargreaves, J. D. Annan, R. Ohgaito, A. Paul, and A. Abe-Ouchi Paleoclimate simulations provide us with an opportunity to critically confront and evaluate the performance of climate models in simulating the response of the climate system to changes in radiative forcing and other boundary conditions. Hargreaves et al. (2011) analysed the reliability of the Paleoclimate Modelling Intercomparison Project, PMIP2 model ensemble with respect to the MARGO sea surface temperature data synthesis (MARGO Project Members, 2009) for the Last Glacial Maximum (LGM, 21 ka BP). Here we extend that work to include a new comprehensive collection of land surface data (Bartlein et al., 2011), and introduce a novel analysis of the predictive skill of the models. We include output from the PMIP3 experiments, from the two models for which suitable data are currently available. We also perform the same analyses for the PMIP2 mid-Holocene (6 ka BP) ensembles and available proxy data sets. Our results are predominantly positive for the LGM, suggesting that as well as the global mean change, the models can reproduce the observed pattern of change on the broadest scales, such as the overall land–sea contrast and polar amplification, although the more detailed sub-continental scale patterns of change remains elusive. In contrast, our results for the mid-Holocene are substantially negative, with the models failing to reproduce the observed changes with any degree of skill. One cause of this problem could be that the globally- and annually-averaged forcing anomaly is very weak at the mid-Holocene, and so the results are dominated by the more localised regional patterns in the parts of globe for which data are available. The root cause of the model-data mismatch at these scales is unclear. If the proxy calibration is itself reliable, then representativity error in the data-model comparison, and missing climate feedbacks in the models are other possible sources of error.

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: 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: Stable isotopic evidence of El Niño-like atmospheric circulation in the Pliocene western United States Climate of the Past, 9, 903-912, 2013 Author(s): M. J. Winnick, J. M. Welker, and C. P. Chamberlain Understanding how the hydrologic cycle has responded to warmer global temperatures in the past is especially important today as concentrations of CO 2 in the atmosphere continue to increase due to human activities. The Pliocene offers an ideal window into a climate system that has equilibrated with current atmospheric p CO 2 . During the Pliocene the western United States was wetter than modern, an observation at odds with our current understanding of future warming scenarios, which involve the expansion and poleward migration of the subtropical dry zone. Here we compare Pliocene oxygen isotope profiles of pedogenic carbonates across the western US to modern isotopic anomalies in precipitation between phases of the El Niño–Southern Oscillation (ENSO). We find that when accounting for seasonality of carbonate formation, isotopic changes through the late Pliocene match modern precipitation isotopic anomalies in El Niño years. Furthermore, isotopic shifts through the late Pliocene mirror changes through the early Pleistocene, which likely represents the southward migration of the westerly storm track caused by growth of the Laurentide ice sheet. We propose that the westerly storm track migrated northward through the late Pliocene with the development of the modern cold tongue in the east equatorial Pacific, then returned southward with widespread glaciation in the Northern Hemisphere – a scenario supported by terrestrial climate proxies across the US. Together these data support the proposed existence of background El Niño-like conditions in western North America during the warm Pliocene. If the earth behaves similarly with future warming, this observation has important implications with regard to the amount and distribution of precipitation in western North America.

Description: Mass-movement and flood-induced deposits in Lake Ledro, southern Alps, Italy: implications for Holocene palaeohydrology and natural hazards Climate of the Past, 9, 825-840, 2013 Author(s): A. Simonneau, E. Chapron, B. Vannière, S. B. Wirth, A. Gilli, C. Di Giovanni, F. S. Anselmetti, M. Desmet, and M. Magny High-resolution seismic profiles and sediment cores from Lake Ledro combined with soil and riverbed samples from the lake's catchment area are used to assess the recurrence of natural hazards (earthquakes and flood events) in the southern Italian Alps during the Holocene. Two well-developed deltas and a flat central basin are identified on seismic profiles in Lake Ledro. Lake sediments have been finely laminated in the basin since 9000 cal. yr BP and frequently interrupted by two types of sedimentary events (SEs): light-coloured massive layers and dark-coloured graded beds. Optical analysis (quantitative organic petrography) of the organic matter present in soil, riverbed and lacustrine samples together with lake sediment bulk density and grain-size analysis illustrate that light-coloured layers consist of a mixture of lacustrine sediments and mainly contain algal particles similar to the ones observed in background sediments. Light-coloured layers thicker than 1.5 cm in the main basin of Lake Ledro are synchronous to numerous coeval mass-wasting deposits remoulding the slopes of the basin. They are interpreted as subaquatic mass-movements triggered by historical and pre-historical regional earthquakes dated to AD 2005, AD 1891, AD 1045 and 1260, 2545, 2595, 3350, 3815, 4740, 7190, 9185 and 11 495 cal. yr BP. Dark-coloured SEs develop high-amplitude reflections in front of the deltas and in the deep central basin. These beds are mainly made of terrestrial organic matter (soils and lignocellulosic debris) and are interpreted as resulting from intense hyperpycnal flood event. Mapping and quantifying the amount of soil material accumulated in the Holocene hyperpycnal flood deposits of the sequence allow estimating that the equivalent soil thickness eroded over the catchment area reached up to 5 mm during the largest Holocene flood events. Such significant soil erosion is interpreted as resulting from the combination of heavy rainfall and snowmelt. The recurrence of flash flood events during the Holocene was, however, not high enough to affect pedogenesis processes and highlight several wet regional periods during the Holocene. The Holocene period is divided into four phases of environmental evolution. Over the first half of the Holocene, a progressive stabilization of the soils present through the catchment of Lake Ledro was associated with a progressive reforestation of the area and only interrupted during the wet 8.2 event when the soil destabilization was particularly important. Lower soil erosion was recorded during the mid-Holocene climatic optimum (8000–4200 cal. yr BP) and associated with higher algal production. Between 4200 and 3100 cal. yr BP, both wetter climate and human activities within the drainage basin drastically increased soil erosion rates. Finally, from 3100 cal. yr BP to the present-day, data suggest increasing and changing human land use.

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: Broadband measurements of aerosol extinction in the ultraviolet spectral region Atmospheric Measurement Techniques, 6, 861-877, 2013 Author(s): R. A. Washenfelder, J. M. Flores, C. A. Brock, S. S. Brown, and Y. Rudich Aerosols influence the Earth's radiative budget by scattering and absorbing incoming solar radiation. The optical properties of aerosols vary as a function of wavelength, but few measurements have reported the wavelength dependence of aerosol extinction cross sections and complex refractive indices. We describe a new laboratory instrument to measure aerosol optical extinction as a function of wavelength, using cavity enhanced spectroscopy with a broadband light source. The instrument consists of two broadband channels which span the 360–390 and 385–420 nm spectral regions using two light emitting diodes (LED) and a grating spectrometer with charge-coupled device (CCD) detector. We determined aerosol extinction cross sections and directly observed Mie scattering resonances for aerosols that are purely scattering (polystyrene latex spheres and ammonium sulfate), slightly absorbing (Suwannee River fulvic acid), and strongly absorbing (nigrosin dye). We describe an approach for retrieving refractive indices as a function of wavelength from the measured extinction cross sections over the 360–420 nm wavelength region. The retrieved refractive indices for PSL and ammonium sulfate agree within uncertainty with the literature values for this spectral region. The refractive index determined for nigrosin is 1.78 (± 0.03) + 0.19 (± 0.08) i at 360 nm and 1.63 (± 0.03) + 0.21 (± 0.05) i at 420 nm. The refractive index determined for Suwannee River fulvic acid is 1.71 (± 0.02) + 0.07 (± 0.06) i at 360 nm and 1.66 (± 0.02) + 0.06 (± 0.04) i at 420 nm. These laboratory results support the potential for a field instrument capable of determining ambient aerosol optical extinction, average aerosol extinction cross section, and complex refractive index as a function of wavelength.

Description: Using data assimilation to investigate the causes of Southern Hemisphere high latitude cooling from 10 to 8 ka BP Climate of the Past, 9, 887-901, 2013 Author(s): P. Mathiot, H. Goosse, X. Crosta, B. Stenni, M. Braida, H. Renssen, C. J. Van Meerbeeck, V. Masson-Delmotte, A. Mairesse, and S. Dubinkina From 10 to 8 ka BP (thousand years before present), paleoclimate records show an atmospheric and oceanic cooling in the high latitudes of the Southern Hemisphere. During this interval, temperatures estimated from proxy data decrease by 0.8 °C over Antarctica and 1.2 °C over the Southern Ocean. In order to study the causes of this cooling, simulations covering the early Holocene have been performed with the climate model of intermediate complexity LOVECLIM constrained to follow the signal recorded in climate proxies using a data assimilation method based on a particle filtering approach. The selected proxies represent oceanic and atmospheric surface temperature in the Southern Hemisphere derived from terrestrial, marine and glaciological records. Two mechanisms previously suggested to explain the 10–8 ka BP cooling pattern are investigated using the data assimilation approach in our model. The first hypothesis is a change in atmospheric circulation, and the second one is a cooling of the sea surface temperature in the Southern Ocean, driven in our experimental setup by the impact of an increased West Antarctic melting rate on ocean circulation. For the atmosphere hypothesis, the climate state obtained by data assimilation produces a modification of the meridional atmospheric circulation leading to a 0.5 °C Antarctic cooling from 10 to 8 ka BP compared to the simulation without data assimilation, without congruent cooling of the atmospheric and sea surface temperature in the Southern Ocean. For the ocean hypothesis, the increased West Antarctic freshwater flux constrainted by data assimilation (+100 mSv from 10 to 8 ka BP) leads to an oceanic cooling of 0.7 °C and a strengthening of Southern Hemisphere westerlies (+6%). Thus, according to our experiments, the observed cooling in Antarctic and the Southern Ocean proxy records can only be reconciled with the reconstructions by the combination of a modified atmospheric circulation and an enhanced freshwater flux.

Description: Climatologies from satellite measurements: the impact of orbital sampling on the standard error of the mean Atmospheric Measurement Techniques, 6, 937-948, 2013 Author(s): M. Toohey and T. von Clarmann Climatologies of atmospheric observations are often produced by binning measurements according to latitude and calculating zonal means. The uncertainty in these climatological means is characterised by the standard error of the mean (SEM). However, the usual estimator of the SEM, i.e., the sample standard deviation divided by the square root of the sample size, holds only for uncorrelated randomly sampled measurements. Measurements of the atmospheric state along a satellite orbit cannot always be considered as independent because (a) the time-space interval between two nearest observations is often smaller than the typical scale of variations in the atmospheric state, and (b) the regular time-space sampling pattern of a satellite instrument strongly deviates from random sampling. We have developed a numerical experiment where global chemical fields from a chemistry climate model are sampled according to real sampling patterns of satellite-borne instruments. As case studies, the model fields are sampled using sampling patterns of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and Atmospheric Chemistry Experiment Fourier-Transform Spectrometer (ACE-FTS) satellite instruments. Through an iterative subsampling technique, and by incorporating information on the random errors of the MIPAS and ACE-FTS measurements, we produce empirical estimates of the standard error of monthly mean zonal mean model O 3 in 5° latitude bins. We find that generally the classic SEM estimator is a conservative estimate of the SEM, i.e., the empirical SEM is often less than or approximately equal to the classic estimate. Exceptions occur only when natural variability is larger than the random measurement error, and specifically in instances where the zonal sampling distribution shows non-uniformity with a similar zonal structure as variations in the sampled field, leading to maximum sensitivity to arbitrary phase shifts between the sample distribution and sampled field. The occurrence of such instances is thus very sensitive to slight changes in the sampling distribution, and to the variations in the measured field. This study highlights the need for caution in the interpretation of the oft-used classically computed SEM, and outlines a relatively simple methodology that can be used to assess one component of the uncertainty in monthly mean zonal mean climatologies produced from measurements from satellite-borne instruments.

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: Model sensitivity to North Atlantic freshwater forcing at 8.2 ka Climate of the Past, 9, 955-968, 2013 Author(s): C. Morrill, A. N. LeGrande, H. Renssen, P. Bakker, and B. L. Otto-Bliesner We compared four simulations of the 8.2 ka event to assess climate model sensitivity and skill in responding to North Atlantic freshwater perturbations. All of the simulations used the same freshwater forcing, 2.5 Sv for one year, applied to either the Hudson Bay (northeastern Canada) or Labrador Sea (between Canada's Labrador coast and Greenland). This freshwater pulse induced a decadal-mean slowdown of 10–25% in the Atlantic Meridional Overturning Circulation (AMOC) of the models and caused a large-scale pattern of climate anomalies that matched proxy evidence for cooling in the Northern Hemisphere and a southward shift of the Intertropical Convergence Zone. The multi-model ensemble generated temperature anomalies that were just half as large as those from quantitative proxy reconstructions, however. Also, the duration of AMOC and climate anomalies in three of the simulations was only several decades, significantly shorter than the duration of ~150 yr in the paleoclimate record. Possible reasons for these discrepancies include incorrect representation of the early Holocene climate and ocean state in the North Atlantic and uncertainties in the freshwater forcing estimates.

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: A mid-Holocene climate reconstruction for eastern South America Climate of the Past, 9, 2117-2133, 2013 Author(s): L. F. Prado, I. Wainer, C. M. Chiessi, M.-P. Ledru, and B. Turcq The mid-Holocene (6000 calibrated years before present) is a key period in palaeoclimatology because incoming summer insolation was lower than during the late Holocene in the Southern Hemisphere, whereas the opposite happened in the Northern Hemisphere. However, the effects of the decreased austral summer insolation over South American climate have been poorly discussed by palaeodata syntheses. In addition, only a few of the regional studies have characterised the mid-Holocene climate in South America through a multiproxy approach. Here, we present a multiproxy compilation of mid-Holocene palaeoclimate data for eastern South America. We compiled 120 palaeoclimatological datasets, which were published in 84 different papers. The palaeodata analysed here suggest a water deficit scenario in the majority of eastern South America during the mid-Holocene if compared to the late Holocene, with the exception of northeastern Brazil. Low mid-Holocene austral summer insolation caused a reduced land–sea temperature contrast and hence a weakened South American monsoon system circulation. This scenario is represented by a decrease in precipitation over the South Atlantic Convergence Zone area, saltier conditions along the South American continental margin, and lower lake levels.

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: A chemical analyzer for charged ultrafine particles Atmospheric Measurement Techniques, 6, 2339-2348, 2013 Author(s): S. G. Gonser and A. Held New particle formation is a frequent phenomenon in the atmosphere and of major significance for the Earth's climate and human health. To date the mechanisms leading to the nucleation of particles as well as to aerosol growth are not completely understood. A lack of appropriate measurement equipment for online analysis of the chemical composition of freshly nucleated particles is one major limitation. We have developed a Chemical Analyzer for Charged Ultrafine Particles (CAChUP) capable of analyzing particles with diameters below 30 nm. A bulk of size-separated particles is collected electrostatically on a metal filament, resistively desorbed and subsequently analyzed for its molecular composition in a time of flight mass spectrometer. We report on technical details as well as characterization experiments performed with the CAChUP. Our instrument was tested in the laboratory for its detection performance as well as for its collection and desorption capabilities. The manual application of defined masses of camphene (C 10 H 16 ) to the desorption filament resulted in a detection limit between 0.5 and 5 ng, and showed a linear response of the mass spectrometer. Flow tube experiments of 25 nm diameter secondary organic aerosol from ozonolysis of alpha-pinene also showed a linear relation between collection time and the mass spectrometer's signal intensity. The resulting mass spectra from the collection experiments are in good agreement with published work on particles generated by the ozonolysis of alpha-pinene. A sensitivity study shows that the current setup of CAChUP is ready for laboratory measurements and for the observation of new particle formation events in the field.

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: Glacial fluctuations of the Indian monsoon and their relationship with North Atlantic climate: new data and modelling experiments Climate of the Past, 9, 2135-2151, 2013 Author(s): C. Marzin, N. Kallel, M. Kageyama, J.-C. Duplessy, and P. Braconnot Several paleoclimate records such as from Chinese loess, speleothems or upwelling indicators in marine sediments present large variations of the Asian monsoon system during the last glaciation. Here, we present a new record from the northern Andaman Sea (core MD77-176) which shows the variations of the hydrological cycle of the Bay of Bengal. The high-resolution record of surface water δ 18 O dominantly reflects salinity changes and displays large millennial-scale oscillations over the period 40 000 to 11 000 yr BP. Their timing and sequence suggests that events of high (resp. low) salinity in the Bay of Bengal, i.e. weak (resp. strong) Indian monsoon, correspond to cold (resp. warm) events in the North Atlantic and Arctic, as documented by the Greenland ice core record. We use the IPSL_CM4 Atmosphere-Ocean coupled General Circulation Model to study the processes that could explain the teleconnection between the Indian monsoon and the North Atlantic climate. We first analyse a numerical experiment in which such a rapid event in the North Atlantic is obtained under glacial conditions by increasing the freshwater flux in the North Atlantic, which results in a reduction of the intensity of the Atlantic meridional overturning circulation. This freshwater hosing results in a weakening of the Indian monsoon rainfall and circulation. The changes in the continental runoff and local hydrological cycle are responsible for an increase in salinity in the Bay of Bengal. This therefore compares favourably with the new sea water δ 18 O record presented here and the hypothesis of synchronous cold North Atlantic and weak Indian monsoon events. Additional sensitivity experiments are produced with the LMDZ atmospheric model to analyse the teleconnection mechanisms between the North Atlantic and the Indian monsoon. The changes over the tropical Atlantic are shown to be essential in triggering perturbations of the subtropical jet over Africa and Eurasia, that in turn affect the intensity of the Indian monsoon. These relationships are also found to be valid in additional coupled model simulations in which the Atlantic meridional overturning circulation (AMOC) is forced to resume.

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: Mid-Holocene ocean and vegetation feedbacks over East Asia Climate of the Past, 9, 2153-2171, 2013 Author(s): Z. Tian and D. Jiang Mid-Holocene ocean and vegetation feedbacks over East Asia are investigated by a set of numerical experiments performed with the version 4 of the Community Climate System Model (CCSM4). With reference to the pre-industrial period, most of the mid-Holocene annual and seasonal surface-air temperature and precipitation changes are found to result from a direct response of the atmosphere to insolation forcing, while dynamic ocean and vegetation modulate regional climate of East Asia to some extent. Because of its thermal inertia, the dynamic ocean induced an additional warming of 0.2 K for the annual mean, 0.5 K in winter (December–February), 0.0003 K in summer (June–August), and 1.0 K in autumn (September–November), but a cooling of 0.6 K in spring (March–May) averaged over China, and it counteracted (amplified) the direct effect of insolation forcing for the annual mean and in winter and autumn (spring) for that period. The dynamic vegetation had an area-average impact of no more than 0.4 K on the mid-Holocene annual and seasonal temperatures over China, with an average cooling of 0.2 K for the annual mean. On the other hand, ocean feedback induced a small increase of precipitation in winter (0.04 mm day −1 ) and autumn (0.05 mm day −1 ), but a reduction for the annual mean (0.14 mm day −1 ) and in spring (0.29 mm day −1 ) and summer (0.34 mm day −1 ) over China, while it also suppressed the East Asian summer monsoon rainfall. The effect of dynamic vegetation on the mid-Holocene annual and seasonal precipitation was comparatively small, ranging from −0.03 mm day −1 to 0.06 mm day −1 averaged over China. In comparison, the CCSM4 simulated annual and winter cooling over China agrees with simulations within the Paleoclimate Modeling Intercomparison Project (PMIP), but the results are contrary to the warming reconstructed from multiple proxy data for the mid-Holocene. Ocean feedback narrows this model–data mismatch, whereas vegetation feedback plays an opposite role but with a level of uncertainty.

Description: Comparison of AOD between CALIPSO and MODIS: significant differences over major dust and biomass burning regions Atmospheric Measurement Techniques, 6, 2391-2401, 2013 Author(s): X. Ma, K. Bartlett, K. Harmon, and F. Yu Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) provide global vertical profiles of aerosol optical properties for the first time. In this study, we employed about 6 yr (2006–2011) of CALIPSO level 3 monthly mean gridded aerosol optical depth (AOD) products (daytime and nighttime) for cloud-free conditions, to compare with the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra/Aqua level 3 monthly mean AOD dataset for the same time period. While the spatial distribution and seasonal variability of CALIPSO AOD is generally consistent with that of MODIS, CALIPSO is overall lower than MODIS as MODIS has higher frequency than CALIPSO for most bins of AOD. The correlation between MODIS and CALIPSO is better over ocean than over land. We focused on four regions that have large systematic differences: two over dust regions (the Sahara and Northwest China) and two over biomass burning regions (South Africa and South America). It is found that CALIPSO AOD is significantly lower than MODIS AOD over dust regions during the whole time period, with a maximum difference of 0.3 over the Saharan region and 0.25 over Northwest China. For biomass burning regions, CALIPSO AOD is significantly higher than MODIS AOD over South Africa, with a maximum difference of 0.25. Additionally CALIPSO AOD is slightly higher than MODIS AOD over South America for most of the time period, with a few exceptions in 2006, 2007, and 2010, when biomass burning is significantly stronger than during other years. We analyzed the impact of the satellite spatial and temporal sampling issue by using level 2 CALIPSO and MODIS products, and these systematic differences can still be found. The results of this study indicate that systematic differences of CALIPSO relative to MODIS are closely associated with aerosol types, which vary by location and season. Large differences over dust and biomass burning regions may suggest that assumptions made in satellite retrievals, such as the assumed lidar ratios for CALIPSO retrievals over dust and biomass burning regions or the surface reflectance information and/or the aerosol model utilized by the MODIS algorithm, are not appropriate.

Description: Water-soluble organic carbon in snow and ice deposited at Alpine, Greenland, and Antarctic sites: a critical review of available data and their atmospheric relevance Climate of the Past, 9, 2195-2211, 2013 Author(s): M. Legrand, S. Preunkert, B. Jourdain, J. Guilhermet, X. Fa{ï}n, I. Alekhina, and J. R. Petit While it is now recognized that organic matter dominates the present-day atmospheric aerosol load over continents, its sources remain poorly known. The studies of organic species or organic fractions trapped in ice cores may help to overcome this lack of knowledge. Available data on the dissolved (or total) organic carbon (DOC or TOC) content of snow and ice often appear largely inconsistent, and, until now, no critical review has been conducted to understand the causes of these inconsistencies. To draw a more consistent picture of the organic carbon amount present in solid precipitation that accumulates on cold glaciers, we here review available data and, when needed, complete the data set with analyses of selected samples. The different data sets are then discussed by considering the age (modern versus pre-industrial, Holocene versus Last glacial Maximum) and type (surface snow, firn, or ice) of investigated samples, the deployed method, and the applied contamination control. Finally, the OC (DOC or TOC) levels of Antarctic, Greenland, and Alpine ice cores are compared and discussed with respect to natural (biomass burning, vegetation emissions) and anthropogenic sources (fossil fuel combustion) contributing to atmospheric OC aerosol.

Description: Microwave radiometer to retrieve temperature profiles from the surface to the stratopause Atmospheric Measurement Techniques, 6, 2477-2494, 2013 Author(s): O. Stähli, A. Murk, N. Kämpfer, C. Mätzler, and P. Eriksson TEMPERA (TEMPERature RAdiometer) is a new ground-based radiometer which measures in a frequency range from 51–57 GHz radiation emitted by the atmosphere. With this instrument it is possible to measure temperature profiles from ground to about 50 km. This is the first ground-based instrument with the capability to retrieve temperature profiles simultaneously for the troposphere and stratosphere. The measurement is done with a filterbank in combination with a digital fast Fourier transform spectrometer. A hot load and a noise diode are used as stable calibration sources. The optics consist of an off-axis parabolic mirror to collect the sky radiation. Due to the Zeeman effect on the emission lines used, the maximum height for the temperature retrieval is about 50 km. The effect is apparent in the measured spectra. The performance of TEMPERA is validated by comparison with nearby radiosonde and satellite data from the Microwave Limb Sounder on the Aura satellite. In this paper we present the design and measurement method of the instrument followed by a description of the retrieval method, together with a validation of TEMPERA data over its first year, 2012.

Description: Major dust events in Europe during marine isotope stage 5 (130–74 ka): a climatic interpretation of the "markers" Climate of the Past, 9, 2213-2230, 2013 Author(s): D.-D. Rousseau, M. Ghil, G. Kukla, A. Sima, P. Antoine, M. Fuchs, C. Hatté, F. Lagroix, M. Debret, and O. Moine At present, major dust storms are occurring at mid-latitudes in the Middle East and Asia, as well as at low latitudes in Northern Africa and in Australia. Western Europe, though, does not experience such dramatic climate events, except for some African dust reaching it from the Sahara. This modern situation is of particular interest, in the context of future climate projections, since the present interglacial is usually interpreted, in this context, as an analog of the warm Eemian interval. European terrestrial records show, however, major dust events during the penultimate interglacial and early glacial. These events are easily observed in loess records by their whitish-color deposits, which lie above and below dark chernozem paleosols in Central European records of Marine Isotope Stage (MIS) 5 age. We describe here the base of the Dolni Vestonice (DV) loess sequence, Czech Republic, as the reference of such records. The dust is deposited during intervals that are characterized by poor vegetation – manifested by high δ 13 C values and low magnetic susceptibility – while the fine sand and clay in the deposits shows grain sizes that are clearly different from the overlying pleniglacial loess deposits. Some of these dust events have been previously described as "Markers" or Marker Silts (MS) by one of us (G. Kukla), and are dated at about 111–109 ka and 93–92 ka, with a third and last one slightly visible at about 75–73 ka. Other events correspond to the loess material of Kukla's cycles, and are described as eolian silts (ES); they are observed in the same DV sequence and are dated at about 106–105 ka, 88–86 ka, and 78.5–77 ka. These dates are determined by considering the OSL ages with their errors measured on the studied sequence, and the comparison with Greenland ice-core and European speleothem chronologies. The fine eolian deposits mentioned above, MS as well as ES, correspond to short events that lasted about 2 ka; they are synchronous with re-advances of the polar front over the North Atlantic, as observed in marine sediment cores. These deposits also correlate with important changes observed in European vegetation. Some ES and MS events appear to be coeval with significant dust peaks recorded in the Greenland ice cores, while others are not. This decoupling between the European eolian and Greenland dust depositions is of considerable interest, as it differs from the fully glacial situation, in which the Eurasian loess sedimentation mimics the Greenland dust record. Previous field observations supported an interpretation of MS events as caused by continental dust storms. We show here, by a comparison with speleothems of the same age found in the northern Alps, that different atmospheric-circulation modes seem to be responsible for the two categories of dust events, MS vs. ES.

Description: Retrieval of nitric oxide in the mesosphere and lower thermosphere from SCIAMACHY limb spectra Atmospheric Measurement Techniques, 6, 2521-2531, 2013 Author(s): S. Bender, M. Sinnhuber, J. P. Burrows, M. Langowski, B. Funke, and M. López-Puertas We use the ultra-violet (UV) spectra in the range 230–300 nm from the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) to retrieve the nitric oxide (NO) number densities from atmospheric emissions in the gamma-bands in the mesosphere and lower thermosphere. Using 3-D ray tracing, a 2-D retrieval grid, and regularisation with respect to altitude and latitude, we retrieve a whole semi-orbit simultaneously for the altitude range from 60 to 160 km. We present details of the retrieval algorithm, first results, and initial comparisons to data from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Our results agree on average well with MIPAS data and are in line with previously published measurements from other instruments. For the time of available measurements in 2008–2011, we achieve a vertical resolution of 5–10 km in the altitude range 70–140 km and a horizontal resolution of about 9° from 60° S–60° N. With this we have independent measurements of the NO densities in the mesosphere and lower thermosphere with approximately global coverage. This data can be further used to validate climate models or as input for them.

Description: Development and field testing of a rapid and ultra-stable atmospheric carbon dioxide spectrometer Atmospheric Measurement Techniques, 7, 4445-4453, 2014 Author(s): B. Xiang, D. D. Nelson, J. B. McManus, M. S. Zahniser, R. A. Wehr, and S. C. Wofsy We present field test results for a new spectroscopic instrument to measure atmospheric carbon dioxide (CO 2 ) with high precision (0.02 μmol mol −1 , or ppm at 1 Hz) and demonstrate high stability (within 0.1 ppm over more than 8 months), without the need for hourly, daily, or even monthly calibration against high-pressure gas cylinders. The technical novelty of this instrument (ABsolute Carbon dioxide, ABC) is the spectral null method using an internal quartz reference cell with known CO 2 column density. Compared to a previously described prototype, the field instrument has better stability and benefits from more precise thermal control of the optics and more accurate pressure measurements in the sample cell (at the mTorr level). The instrument has been deployed at a long-term ecological research site (the Harvard Forest, USA), where it has measured for 8 months without on-site calibration and with minimal maintenance, showing drift bounds of less than 0.1 ppm. Field measurements agree well with those of a commercially available cavity ring-down CO 2 instrument (Picarro G2301) run with a standard calibration protocol. This field test demonstrates that ABC is capable of performing high-accuracy, unattended, continuous field measurements with minimal use of reference gas cylinders.

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: Laminated sediments in the Bering Sea reveal atmospheric teleconnections to Greenland climate on millennial to decadal timescales during the last deglaciation Climate of the Past, 10, 2215-2236, 2014 Author(s): H. Kuehn, L. Lembke-Jene, R. Gersonde, O. Esper, F. Lamy, H. Arz, G. Kuhn, and R. Tiedemann During the last glacial termination, the upper North Pacific Ocean underwent dramatic and rapid changes in oxygenation that lead to the transient intensification of oxygen minimum zones (OMZs), recorded by the widespread occurrence of laminated sediments on circum-Pacific continental margins. We present a new laminated sediment record from the mid-depth (1100 m) northern Bering Sea margin that provides insight into these deglacial OMZ maxima with exceptional, decadal-scale detail. Combined ultrahigh-resolution micro-X-ray-fluorescence (micro-XRF) data and sediment facies analysis of laminae reveal an alternation between predominantly terrigenous and diatom-dominated opal sedimentation. The diatomaceous laminae are interpreted to represent spring/summer productivity events related to the retreating sea ice margin. We identified five laminated sections in the deglacial part of our site. Lamina counts were carried out on these sections and correlated with the Bølling–Allerød and Preboreal phases in the North Greenland Ice Core (NGRIP) oxygen isotope record, indicating an annual deposition of individual lamina couplets (varves). The observed rapid decadal intensifications of anoxia, in particular within the Bølling–Allerød, are tightly coupled to short-term warm events through increases in regional export production. This dependence of laminae formation on warmer temperatures is underlined by a correlation with published Bering Sea sea surface temperature records and δ 18 O data of planktic foraminifera from the Gulf of Alaska. The rapidity of the observed changes strongly implies a close atmospheric teleconnection between North Pacific and North Atlantic regions. We suggest that concomitant increases in export production and subsequent remineralization of organic matter in the Bering Sea, in combination with oxygen-poor waters entering the Being Sea, drove down oxygen concentrations to values below 0.1 mL L −1 and caused laminae preservation. Calculated benthic–planktic ventilation ages show no significant variations throughout the last deglaciation, indicating that changes in formation rates or differing sources of North Pacific mid-depth waters are not prime candidates for strengthening the OMZ at our site. The age models established by our correlation procedure allow for the determination of calendar age control points for the Bølling–Allerød and the Preboreal that are independent of the initial radiocarbon-based chronology. Resulting surface reservoir ages range within 730–990 yr during the Bølling–Allerød, 800–1100 yr in the Younger Dryas, and 765–775 yr for the Preboreal.

Description: Implication of methodological uncertainties for mid-Holocene sea surface temperature reconstructions Climate of the Past, 10, 2237-2252, 2014 Author(s): I. Hessler, S. P. Harrison, M. Kucera, C. Waelbroeck, M.-T. Chen, C. Anderson, A. de Vernal, B. Fréchette, A. Cloke-Hayes, G. Leduc, and L. Londeix We present and examine a multi-sensor global compilation of mid-Holocene (MH) sea surface temperatures (SST), based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst census counts. We assess the uncertainties originating from using different methodologies and evaluate the potential of MH SST reconstructions as a benchmark for climate-model simulations. The comparison between different analytical approaches (time frame, baseline climate) shows the choice of time window for the MH has a negligible effect on the reconstructed SST pattern, but the choice of baseline climate affects both the magnitude and spatial pattern of the reconstructed SSTs. Comparison of the SST reconstructions made using different sensors shows significant discrepancies at a regional scale, with uncertainties often exceeding the reconstructed SST anomaly. Apparent patterns in SST may largely be a reflection of the use of different sensors in different regions. Overall, the uncertainties associated with the SST reconstructions are generally larger than the MH anomalies. Thus, the SST data currently available cannot serve as a target for benchmarking model simulations. Further evaluations of potential subsurface and/or seasonal artifacts that may contribute to obscure the MH SST reconstructions are urgently needed to provide reliable benchmarks for model evaluations.

Description: Investigating uptake of N 2 O in agricultural soils using a high-precision dynamic chamber method Atmospheric Measurement Techniques, 7, 4455-4462, 2014 Author(s): N. J. Cowan, D. Famulari, P. E. Levy, M. Anderson, D. S. Reay, and U. M. Skiba Uptake (or negative flux) of nitrous oxide (N 2 O) in agricultural soils is a controversial issue which has proved difficult to investigate in the past due to constraints such as instrumental precision and methodological uncertainties. Using a recently developed high-precision quantum cascade laser gas analyser combined with a closed dynamic chamber, a well-defined detection limit of 4 μg N 2 O-N m −2 h −1 could be achieved for individual soil flux measurements. 1220 measurements of N 2 O flux were made from a variety of UK soils using this method, of which 115 indicated uptake by the soil (i.e. a negative flux in the micrometeorological sign convention). Only four of these apparently negative fluxes were greater than the detection limit of the method, which suggests that the vast majority of reported negative fluxes from such measurements are actually due to instrument noise. As such, we suggest that the bulk of negative N 2 O fluxes reported for agricultural fields are most likely due to limits in detection of a particular flux measurement methodology and not a result of microbiological activity consuming atmospheric N 2 O.

Description: Rapid, optical measurement of the atmospheric pressure on a fast research aircraft using open-path TDLAS Atmospheric Measurement Techniques, 7, 3653-3666, 2014 Author(s): B. Buchholz, A. Afchine, and V. Ebert Because of the high travel speed, the complex flow dynamics around an aircraft, and the complex dependency of the fluid dynamics on numerous airborne parameters, it is quite difficult to obtain accurate pressure values at a specific instrument location of an aircraft's fuselage. Complex simulations using computational fluid dynamics (CFD) models can in theory computationally "transfer" pressure values from one location to another. However, for long flight patterns, this process is inconvenient and cumbersome. Furthermore, these CFD transfer models require a local experimental validation, which is rarely available. In this paper, we describe an integrated approach for a spectroscopic, calibration-free, in-flight pressure determination in an open-path White cell on an aircraft fuselage using ambient, atmospheric water vapour as the "sensor species". The presented measurements are realised with the HAI (Hygrometer for Atmospheric Investigations) instrument, built for multiphase water detection via calibration-free TDLAS (tunable diode laser absorption spectroscopy). The pressure determination is based on raw data used for H 2 O concentration measurement, but with a different post-flight evaluation method, and can therefore be conducted at deferred time intervals on any desired flight track. The spectroscopic pressure is compared in-flight with the static ambient pressure of the aircraft avionic system and a micro-mechanical pressure sensor, located next to the open-path cell, over a pressure range from 150 to 800 hPa, and a water vapour concentration range of more than 3 orders of magnitude. The correlation between the micro-mechanical pressure sensor measurements and the spectroscopic pressure measurements shows an average deviation from linearity of only 0.14% and a small offset of 9.5 hPa. For the spectroscopic pressure evaluation we derive measurement uncertainties under laboratory conditions of 3.2 and 5.1% during in-flight operation on the HALO airplane. Under certain flight conditions we quantified, for the first time, stalling-induced, dynamic pressure deviations of up to 30% (at 200 hPa) between the avionic sensor and the optical and mechanical pressure sensors integrated in HAI. Such severe local pressure deviations from the typically used avionic pressure are important to take into account for other airborne sensors employed on such fast flying platforms as the HALO aircraft.

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: Mixing-layer height retrieval with ceilometer and Doppler lidar: from case studies to long-term assessment Atmospheric Measurement Techniques, 7, 3685-3704, 2014 Author(s): J. H. Schween, A. Hirsikko, U. Löhnert, and S. Crewell Aerosol signatures observed by ceilometers are frequently used to derive mixing-layer height (MLH) which is an essential variable for air quality modelling. However, Doppler wind lidar measurements of vertical velocity can provide a more direct estimation of MLH via simple thresholding. A case study reveals difficulties in the aerosol-based MLH retrieval during transition times when the mixing layer builds up in the morning and when turbulence decays in the afternoon. The difficulties can be explained by the fact that the aerosol distribution is related to the history of the mixing process and aerosol characteristics are modified by humidification. The results of the case study are generalized by evaluating one year of joint measurements by a Vaisala CT25K and a HALO Photonics Streamline wind lidar. On average the aerosol-based retrieval gives higher MLH than the wind lidar with an overestimation of MLH by about 300 m (600 m) in the morning (late afternoon). Also, the daily aerosol-based maximum MLH is larger and occurs later during the day and the average morning growth rates are smaller than those derived from the vertical wind. In fair weather conditions classified by less than 4 octa cloud cover the mean diurnal cycle of cloud base height corresponds well to the mixing-layer height showing potential for a simplified MLH estimation.

Description: Oxygen stable isotopes during the Last Glacial Maximum climate: perspectives from data–model ( i LOVECLIM) comparison Climate of the Past, 10, 1939-1955, 2014 Author(s): T. Caley, D. M. Roche, C. Waelbroeck, and E. Michel We use the fully coupled atmosphere–ocean three-dimensional model of intermediate complexity i LOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 years). By using a model that is able to explicitly simulate the sensor (δ 18 O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that i LOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean δ 18 O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice sheet. We observe a general pattern toward more enriched values for continental calcite δ 18 O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data–model comparison for sea surface temperature indicates that i LOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data–model comparison for calcite δ 18 O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a strong mean annual cooling in the area south of Iceland and Greenland between the LGM and present (〉 6 °C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data–model comparison also reveals that large positive calcite δ 18 O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive δ 18 Osw anomaly for the north Indian Ocean that contrasts with a large negative δ 18 Osw anomaly in the China Sea between the LGM and the present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in δ 18 Osw between the LGM and the present are not spatially homogeneous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite δ 18 O data. Nonetheless, our data–model comparison supports a heterogeneous cooling of a few degrees (2–4 °C) in the LGM Ocean.

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: Deciphering the spatio-temporal complexity of climate change of the last deglaciation: a model analysis Climate of the Past, 7, 591-602, 2011 Author(s): D. M. Roche, H. Renssen, D. Paillard, and G. Levavasseur Understanding the sequence of events occuring during the last major glacial to interglacial transition (21 ka BP to 9 ka BP) is a challenging task that has the potential to unveil the mechanisms behind large scale climate changes. Though many studies have focused on the understanding of the complex sequence of rapid climatic change that accompanied or interrupted the deglaciation, few have analysed it in a more theoretical framework with simple forcings. In the following, we address when and where the first significant temperature anomalies appeared when using slow varying forcing of the last deglaciation. We used here coupled transient simulations of the last deglaciation, including ocean, atmosphere and vegetation components to analyse the spatial timing of the deglaciation. To keep the analysis in a simple framework, we did not include freshwater forcings that potentially cause rapid climate shifts during that time period. We aimed to disentangle the direct and subsequent response of the climate system to slow forcing and moreover, the location where those changes are more clearly expressed. In a data – modelling comparison perspective, this could help understand the physically plausible phasing between known forcings and recorded climatic changes. Our analysis of climate variability could also help to distinguish deglacial warming signals from internal climate variability. We thus are able to better pinpoint the onset of local deglaciation, as defined by the first significant local warming and further show that there is a large regional variability associated with it, even with the set of slow forcings used here. In our model, the first significant hemispheric warming occurred simultaneously in the North and in the South and is a direct response to the obliquity forcing.

Description: Radio occultation bending angle anomalies during tropical cyclones Atmospheric Measurement Techniques, 4, 1053-1060, 2011 Author(s): R. Biondi, T. Neubert, S. Syndergaard, and J. K. Nielsen The tropical deep convection affects the radiation balance of the atmosphere changing the water vapor mixing ratio and the temperature of the upper troposphere lower stratosphere. The aim of this work is to better understand these processes and to investigate if severe storms leave a significant signature in radio occultation profiles in the tropical tropopause layer. Using tropical cyclone best track database and data from different GPS radio occultation missions (COSMIC, GRACE, CHAMP, SACC and GPSMET), we selected 1194 profiles in a time window of 3 h and a space window of 300 km from the eye of the cyclone. We show that the bending angle anomaly of a GPS radio occultation signal is typically larger than the climatology in the upper troposphere and lower stratosphere and that a double tropopause during deep convection can easily be detected using this technique. Comparisons with co-located radiosondes, climatology of tropopause altitudes and GOES analyses are also shown to support the hypothesis that the bending angle anomaly can be used as an indicator of convective towers. The results are discussed in connection to the GPS radio occultation receiver which will be part of the Atomic Clock Ensemble in Space (ACES) payload on the International Space Station.

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: Assimilation of GPS radio occultation data at DWD Atmospheric Measurement Techniques, 4, 1105-1113, 2011 Author(s): H. Anlauf, D. Pingel, and A. Rhodin We describe the status of the assimilation of bending angles from GPS radio occultations in the 3D-Var for DWD's operational global forecast model GME ("Global Model for Europe"). Experiments show that the assimilation of GPSRO data leads to a significant reduction of biases in the analyses of temperature, humidity and wind in the upper troposphere and the stratosphere, as well as a better r. m. s. fit in the comparison to radiosondes. The impact on forecasts is most prominent in the data sparse Southern Hemisphere, but is also quite notable in the Northern Hemisphere extra-tropics. The positive results found in the impact experiments lead to the implementation of the assimilation of GPS radio occultations from GRACE-A, FORMOSAT-3/COSMIC and GRAS/MetOp-A into the operational suite on 3 August 2010. We also show some initial results from assimilation experiments using radio occultation data from the German research satellite TerraSAR-X.

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.