Search Results

Now showing 1 - 10 of 24
  • Item
    Towards global empirical upscaling of FLUXNET eddy covariance observations: Validation of a model tree ensemble approach using a biosphere model
    (München : European Geopyhsical Union, 2009) Jung, M.; Reichstein, M.; Bondeau, A.
    Global, spatially and temporally explicit estimates of carbon and water fluxes derived from empirical up-scaling eddy covariance measurements would constitute a new and possibly powerful data stream to study the variability of the global terrestrial carbon and water cycle. This paper introduces and validates a machine learning approach dedicated to the upscaling of observations from the current global network of eddy covariance towers (FLUXNET). We present a new model TRee Induction ALgorithm (TRIAL) that performs hierarchical stratification of the data set into units where particular multiple regressions for a target variable hold. We propose an ensemble approach (Evolving tRees with RandOm gRowth, ERROR) where the base learning algorithm is perturbed in order to gain a diverse sequence of different model trees which evolves over time. We evaluate the efficiency of the model tree ensemble (MTE) approach using an artificial data set derived from the Lund-Potsdam-Jena managed Land (LPJmL) biosphere model. We aim at reproducing global monthly gross primary production as simulated by LPJmL from 1998–2005 using only locations and months where high quality FLUXNET data exist for the training of the model trees. The model trees are trained with the LPJmL land cover and meteorological input data, climate data, and the fraction of absorbed photosynthetic active radiation simulated by LPJmL. Given that we know the "true result" in the form of global LPJmL simulations we can effectively study the performance of the MTE upscaling and associated problems of extrapolation capacity. We show that MTE is able to explain 92% of the variability of the global LPJmL GPP simulations. The mean spatial pattern and the seasonal variability of GPP that constitute the largest sources of variance are very well reproduced (96% and 94% of variance explained respectively) while the monthly interannual anomalies which occupy much less variance are less well matched (41% of variance explained). We demonstrate the substantially improved accuracy of MTE over individual model trees in particular for the monthly anomalies and for situations of extrapolation. We estimate that roughly one fifth of the domain is subject to extrapolation while MTE is still able to reproduce 73% of the LPJmL GPP variability here. This paper presents for the first time a benchmark for a global FLUXNET upscaling approach that will be employed in future studies. Although the real world FLUXNET upscaling is more complicated than for a noise free and reduced complexity biosphere model as presented here, our results show that an empirical upscaling from the current FLUXNET network with MTE is feasible and able to extract global patterns of carbon flux variability.
  • Item
    Heinrich event 1: An example of dynamical ice-sheet reaction to oceanic changes
    (München : European Geopyhsical Union, 2011) Álvarez-Solas, J.; Montoya, M.; Ritz, C.; Ramstein, G.; Charbit, S.; Dumas, C.; Nisancioglu, K.; Dokken, T.; Ganopolski, A.
    Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.
  • Item
    Similarity estimators for irregular and age-uncertain time series
    (München : European Geopyhsical Union, 2014) Rehfeld, K.; Kurths, J.
    Paleoclimate time series are often irregularly sampled and age uncertain, which is an important technical challenge to overcome for successful reconstruction of past climate variability and dynamics. Visual comparison and interpolation-based linear correlation approaches have been used to infer dependencies from such proxy time series. While the first is subjective, not measurable and not suitable for the comparison of many data sets at a time, the latter introduces interpolation bias, and both face difficulties if the underlying dependencies are nonlinear. In this paper we investigate similarity estimators that could be suitable for the quantitative investigation of dependencies in irregular and age-uncertain time series. We compare the Gaussian-kernel-based cross-correlation (gXCF, Rehfeld et al., 2011) and mutual information (gMI, Rehfeld et al., 2013) against their interpolation-based counterparts and the new event synchronization function (ESF). We test the efficiency of the methods in estimating coupling strength and coupling lag numerically, using ensembles of synthetic stalagmites with short, autocorrelated, linear and nonlinearly coupled proxy time series, and in the application to real stalagmite time series. In the linear test case, coupling strength increases are identified consistently for all estimators, while in the nonlinear test case the correlation-based approaches fail. The lag at which the time series are coupled is identified correctly as the maximum of the similarity functions in around 60–55% (in the linear case) to 53–42% (for the nonlinear processes) of the cases when the dating of the synthetic stalagmite is perfectly precise. If the age uncertainty increases beyond 5% of the time series length, however, the true coupling lag is not identified more often than the others for which the similarity function was estimated. Age uncertainty contributes up to half of the uncertainty in the similarity estimation process. Time series irregularity contributes less, particularly for the adapted Gaussian-kernel-based estimators and the event synchronization function. The introduced link strength concept summarizes the hypothesis test results and balances the individual strengths of the estimators: while gXCF is particularly suitable for short and irregular time series, gMI and the ESF can identify nonlinear dependencies. ESF could, in particular, be suitable to study extreme event dynamics in paleoclimate records. Programs to analyze paleoclimatic time series for significant dependencies are included in a freely available software toolbox.
  • Item
    High-resolution numerical modeling of mesoscale island wakes and sensitivity to static topographic relief data
    (München : European Geopyhsical Union, 2015) Nunalee, C.G.; Horváth, Á.; Basu, S.
    Recent decades have witnessed a drastic increase in the fidelity of numerical weather prediction (NWP) modeling. Currently, both research-grade and operational NWP models regularly perform simulations with horizontal grid spacings as fine as 1 km. This migration towards higher resolution potentially improves NWP model solutions by increasing the resolvability of mesoscale processes and reducing dependency on empirical physics parameterizations. However, at the same time, the accuracy of high-resolution simulations, particularly in the atmospheric boundary layer (ABL), is also sensitive to orographic forcing which can have significant variability on the same spatial scale as, or smaller than, NWP model grids. Despite this sensitivity, many high-resolution atmospheric simulations do not consider uncertainty with respect to selection of static terrain height data set. In this paper, we use the Weather Research and Forecasting (WRF) model to simulate realistic cases of lower tropospheric flow over and downstream of mountainous islands using the default global 30 s United States Geographic Survey terrain height data set (GTOPO30), the Shuttle Radar Topography Mission (SRTM), and the Global Multi-resolution Terrain Elevation Data set (GMTED2010) terrain height data sets. While the differences between the SRTM-based and GMTED2010-based simulations are extremely small, the GTOPO30-based simulations differ significantly. Our results demonstrate cases where the differences between the source terrain data sets are significant enough to produce entirely different orographic wake mechanics, such as vortex shedding vs. no vortex shedding. These results are also compared to MODIS visible satellite imagery and ASCAT near-surface wind retrievals. Collectively, these results highlight the importance of utilizing accurate static orographic boundary conditions when running high-resolution mesoscale models.
  • Item
    Identifying environmental controls on vegetation greenness phenology through model-data integration
    (München : European Geopyhsical Union, 2014) Forkel, M.; Carvalhais, N.; Schaphoff, S.; v. Bloh, W.; Migliavacca, M.; Thurner, M.; Thonicke, K.
    Existing dynamic global vegetation models (DGVMs) have a limited ability in reproducing phenology and decadal dynamics of vegetation greenness as observed by satellites. These limitations in reproducing observations reflect a poor understanding and description of the environmental controls on phenology, which strongly influence the ability to simulate longer-term vegetation dynamics, e.g. carbon allocation. Combining DGVMs with observational data sets can potentially help to revise current modelling approaches and thus enhance the understanding of processes that control seasonal to long-term vegetation greenness dynamics. Here we implemented a new phenology model within the LPJmL (Lund Potsdam Jena managed lands) DGVM and integrated several observational data sets to improve the ability of the model in reproducing satellite-derived time series of vegetation greenness. Specifically, we optimized LPJmL parameters against observational time series of the fraction of absorbed photosynthetic active radiation (FAPAR), albedo and gross primary production to identify the main environmental controls for seasonal vegetation greenness dynamics. We demonstrated that LPJmL with new phenology and optimized parameters better reproduces seasonality, inter-annual variability and trends of vegetation greenness. Our results indicate that soil water availability is an important control on vegetation phenology not only in water-limited biomes but also in boreal forests and the Arctic tundra. Whereas water availability controls phenology in water-limited ecosystems during the entire growing season, water availability co-modulates jointly with temperature the beginning of the growing season in boreal and Arctic regions. Additionally, water availability contributes to better explain decadal greening trends in the Sahel and browning trends in boreal forests. These results emphasize the importance of considering water availability in a new generation of phenology modules in DGVMs in order to correctly reproduce observed seasonal-to-decadal dynamics of vegetation greenness.
  • Item
    Modeling the global emission, transport and deposition of trace elements associated with mineral dust
    (München : European Geopyhsical Union, 2015) Zhang, Y.; Mahowald, N.; Scanza, R.A.; Journet, E.; Desboeufs, K.; Albani, S.; Kok, J.F.; Zhuang, G.; Chen, Y.; Cohen, D.D.; Paytan, A.; Patey, M.D.; Achterberg, E.P.; Engelbrecht, J.P.; Fomba, K.W.
    Trace element deposition from desert dust has important impacts on ocean primary productivity, the quantification of which could be useful in determining the magnitude and sign of the biogeochemical feedback on radiative forcing. However, the impact of elemental deposition to remote ocean regions is not well understood and is not currently included in global climate models. In this study, emission inventories for eight elements primarily of soil origin, Mg, P, Ca, Mn, Fe, K, Al, and Si are determined based on a global mineral data set and a soil data set. The resulting elemental fractions are used to drive the desert dust model in the Community Earth System Model (CESM) in order to simulate the elemental concentrations of atmospheric dust. Spatial variability of mineral dust elemental fractions is evident on a global scale, particularly for Ca. Simulations of global variations in the Ca / Al ratio, which typically range from around 0.1 to 5.0 in soils, are consistent with observations, suggesting that this ratio is a good signature for dust source regions. The simulated variable fractions of chemical elements are sufficiently different; estimates of deposition should include elemental variations, especially for Ca, Al and Fe. The model results have been evaluated with observations of elemental aerosol concentrations from desert regions and dust events in non-dust regions, providing insights into uncertainties in the modeling approach. The ratios between modeled and observed elemental fractions range from 0.7 to 1.6, except for Mg and Mn (3.4 and 3.5, respectively). Using the soil database improves the correspondence of the spatial heterogeneity in the modeling of several elements (Ca, Al and Fe) compared to observations. Total and soluble dust element fluxes to different ocean basins and ice sheet regions have been estimated, based on the model results. The annual inputs of soluble Mg, P, Ca, Mn, Fe and K associated with dust using the mineral data set are 0.30 Tg, 16.89 Gg, 1.32 Tg, 22.84 Gg, 0.068 Tg, and 0.15 Tg to global oceans and ice sheets.
  • Item
    Inertia gravity waves in the upper troposphere during the MaCWAVE winter campaign - Part II: Radar investigations and modelling studies
    (München : European Geopyhsical Union, 2006) Serafimovich, A.; Zülicke, Ch.; Hoffmann, P.; Peters, D.; Dalin, P.; Singer, W.
    We present an experimental and modelling study of a strong gravity wave event in the upper troposphere/lower stratosphere near the Scandinavian mountain ridge. Continuous VHF radar measurements during the MaCWAVE rocket and ground-based measurement campaign were performed at the Norwegian Andoya Rocket Range (ARR) near Andenes (69.3° N, 16° E) in January 2003. Detailed gravity wave investigations based on PSU/NCAR Fifth-Generation Mesoscale Model (MM5) data have been used for comparison with experimentally obtained results. The model data show the presence of a mountain wave and of an inertia gravity wave generated by a jet streak near the tropopause region. Temporal and spatial dependencies of jet induced inertia gravity waves with dominant observed periods of about 13 h and vertical wavelengths of ~4.5–5 km are investigated with wavelet transform applied on radar measurements and model data. The jet induced wave packet is observed to move upstream and downward in the upper troposphere. The model data agree with the experimentally obtained results fairly well. Possible reasons for the observed differences, e.g. in the time of maximum of the wave activity, are discussed. Finally, the vertical fluxes of horizontal momentum are estimated with different methods and provide similar amplitudes. We found indications that the derived positive vertical flux of the horizontal momentum corresponds to the obtained parameters of the jet-induced inertia gravity wave, but only at the periods and heights of the strongest wave activity.
  • Item
    Historical and idealized climate model experiments: An intercomparison of Earth system models of intermediate complexity
    (München : European Geopyhsical Union, 2013) Eby, M.; Weaver, A.J.; Alexander, K.; Zickfeld, K.; Abe-Ouchi, A.; Cimatoribus, A.A.; Crespin, E.; Drijfhout, S.S.; Edwards, N.R.; Eliseev, A.V.; Feulner, G.; Fichefet, T.; Forest, C.E.; Goosse, H.; Holden, P.B.; Joos, F.; Kawamiya, M.; Kicklighter, D.; Kienert, H.; Matsumoto, K.; Mokhov, I.I.; Monier, E.; Olsen, S.M.; Pedersen, J.O.P.; Perrette, M.; Philippon-Berthier, G.; Ridgwell, A.; Schlosser, A.; Schneider von Deimling, T.; Shaffer, G.; Smith, R.S.; Spahni, R.; Sokolov, A.P.; Steinacher, M.; Tachiiri, K.; Tokos, K.; Yoshimori, M.; Zeng, N.; Zhao, F.
    Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.
  • Item
    Extreme events in gross primary production: A characterization across continents
    (München : European Geopyhsical Union, 2014) Zscheischler, J.; Reichstein, M.; Harmeling, S.; Rammig, A.; Tomelleri, E.; Mahecha, M.D.
    Climate extremes can affect the functioning of terrestrial ecosystems, for instance via a reduction of the photosynthetic capacity or alterations of respiratory processes. Yet the dominant regional and seasonal effects of hydrometeorological extremes are still not well documented and in the focus of this paper. Specifically, we quantify and characterize the role of large spatiotemporal extreme events in gross primary production (GPP) as triggers of continental anomalies. We also investigate seasonal dynamics of extreme impacts on continental GPP anomalies. We find that the 50 largest positive extremes (i.e., statistically unusual increases in carbon uptake rates) and negative extremes (i.e., statistically unusual decreases in carbon uptake rates) on each continent can explain most of the continental variation in GPP, which is in line with previous results obtained at the global scale. We show that negative extremes are larger than positive ones and demonstrate that this asymmetry is particularly strong in South America and Europe. Our analysis indicates that the overall impacts and the spatial extents of GPP extremes are power-law distributed with exponents that vary little across continents. Moreover, we show that on all continents and for all data sets the spatial extents play a more important role for the overall impact of GPP extremes compared to the durations or maximal GPP. An analysis of possible causes across continents indicates that most negative extremes in GPP can be attributed clearly to water scarcity, whereas extreme temperatures play a secondary role. However, for Europe, South America and Oceania we also identify fire as an important driver. Our findings are consistent with remote sensing products. An independent validation against a literature survey on specific extreme events supports our results to a large extent.
  • Item
    The Arctic Summer Cloud Ocean Study (ASCOS): Overview and experimental design
    (München : European Geopyhsical Union, 2014) Tjernström, M.; Leck, C.; Birch, C.E.; Bottenheim, J.W.; Brooks, B.J.; Brooks, I.M.; Bäcklin, L.; Chang, R.Y.-W.; de Leeuw, G.; Di Liberto, L.; de la Rosa, S.; Granath, E.; Graus, M.; Hansel, A.; Heintzenberg, J.; Held, A.; Hind, A.; Johnston, P.; Knulst, J.; Martin, M.; Matrai, P.A.; Mauritsen, T.; Müller, M.; Norris, S.J.; Orellana, M.V.; Orsini, D.A.; Paatero, J.; Persson, P.O.G.; Gao, Q.; Rauschenberg, C.; Ristovski, Z.; Sedlar, J.; Shupe, M.D.; Sierau, B.; Sirevaag, A.; Sjogren, S.; Stetzer, O.; Swietlicki, E.; Szczodrak, M.; Vaattovaara, P.; Wahlberg, N.; Westberg, M.; Wheeler, C.R.
    The climate in the Arctic is changing faster than anywhere else on earth. Poorly understood feedback processes relating to Arctic clouds and aerosol–cloud interactions contribute to a poor understanding of the present changes in the Arctic climate system, and also to a large spread in projections of future climate in the Arctic. The problem is exacerbated by the paucity of research-quality observations in the central Arctic. Improved formulations in climate models require such observations, which can only come from measurements in situ in this difficult-to-reach region with logistically demanding environmental conditions. The Arctic Summer Cloud Ocean Study (ASCOS) was the most extensive central Arctic Ocean expedition with an atmospheric focus during the International Polar Year (IPY) 2007–2008. ASCOS focused on the study of the formation and life cycle of low-level Arctic clouds. ASCOS departed from Longyearbyen on Svalbard on 2 August and returned on 9 September 2008. In transit into and out of the pack ice, four short research stations were undertaken in the Fram Strait: two in open water and two in the marginal ice zone. After traversing the pack ice northward, an ice camp was set up on 12 August at 87°21' N, 01°29' W and remained in operation through 1 September, drifting with the ice. During this time, extensive measurements were taken of atmospheric gas and particle chemistry and physics, mesoscale and boundary-layer meteorology, marine biology and chemistry, and upper ocean physics. ASCOS provides a unique interdisciplinary data set for development and testing of new hypotheses on cloud processes, their interactions with the sea ice and ocean and associated physical, chemical, and biological processes and interactions. For example, the first-ever quantitative observation of bubbles in Arctic leads, combined with the unique discovery of marine organic material, polymer gels with an origin in the ocean, inside cloud droplets suggests the possibility of primary marine organically derived cloud condensation nuclei in Arctic stratocumulus clouds. Direct observations of surface fluxes of aerosols could, however, not explain observed variability in aerosol concentrations, and the balance between local and remote aerosols sources remains open. Lack of cloud condensation nuclei (CCN) was at times a controlling factor in low-level cloud formation, and hence for the impact of clouds on the surface energy budget. ASCOS provided detailed measurements of the surface energy balance from late summer melt into the initial autumn freeze-up, and documented the effects of clouds and storms on the surface energy balance during this transition. In addition to such process-level studies, the unique, independent ASCOS data set can and is being used for validation of satellite retrievals, operational models, and reanalysis data sets.