59 results
Search Results
Now showing 1 - 10 of 59
- ItemFracture-induced softening for large-scale ice dynamics(München : European Geopyhsical Union, 2014) Albrecht, T.; Levermann, A.Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e.g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.
- ItemModeling of two different water uptake approaches for mono-and mixed-species forest stands(Basel : MDPI, 2015) Gutsch, Martin; Lasch-Born, Petra; Suckow, Felicitas; Reyer, Christopher P.O.To assess how the effects of drought could be better captured in process-based models, this study simulated and contrasted two water uptake approaches in Scots pine and Scots pine-Sessile oak stands. The first approach consisted of an empirical function for root water uptake (WU1). The second approach was based on differences of soil water potential along a soil-plant-atmosphere continuum (WU2) with total root resistance varying at low, medium and high total root resistance levels. Three data sets on different time scales relevant for tree growth were used for model evaluation: Two short-term datasets on daily transpiration and soil water content as well as a long-term dataset on annual tree ring increments. Except WU2 with high total root resistance, all transpiration outputs exceeded observed values. The strongest correlation between simulated and observed annual tree ring width occurred with WU2 and high total root resistance. The findings highlighted the importance of severe drought as a main reason for small diameter increment. However, if all three data sets were taken into account, no approach was superior to the other. We conclude that accurate projections of future forest productivity depend largely on the realistic representation of root water uptake in forest model simulations.
- ItemNCIO 1.0: A simple Fortran NetCDF interface(München : European Geopyhsical Union, 2015) Robinson, A.; Perrette, M.The NetCDF (Network Common Data Form) library has become an indispensable tool for data and model output management in geoscience. However for simple tasks, particularly in Fortran, the complexity of native NetCDF functionality can be cumbersome. The NCIO (NetCDF Input/Output) module has been designed as an interface to the NetCDF library with simplicity and ease of use in mind. While this implies that some NetCDF functionality is masked from the user, the subroutines provided here are adequate for basic serial reading and writing tasks of up to 6-D data arrays along with corresponding data attributes. The code is available online via a GitHub repository (http://www.github.com/alex-robinson/ncio), which includes an example program to illustrate the approach.
- ItemFrom cyclic ice streaming to Heinrich-like events: the grow-and-surge instability in the parallel ice sheet model(München : European Geopyhsical Union, 2017) Feldmann, Johannes; Levermann, AndersHere we report on a cyclic, physical ice-discharge instability in the Parallel Ice Sheet Model, simulating the flow of a three-dimensional, inherently buttressed ice-sheet-shelf system which periodically surges on a millennial timescale. The thermomechanically coupled model on 1 km horizontal resolution includes an enthalpy-based formulation of the thermodynamics, a nonlinear stress-balance-based sliding law and a very simple subglacial hydrology. The simulated unforced surging is characterized by rapid ice streaming through a bed trough, resulting in abrupt discharge of ice across the grounding line which is eventually calved into the ocean. We visualize the central feedbacks that dominate the subsequent phases of ice buildup, surge and stabilization which emerge from the interaction between ice dynamics, thermodynamics and the subglacial till layer. Results from the variation of surface mass balance and basal roughness suggest that ice sheets of medium thickness may be more susceptible to surging than relatively thin or thick ones for which the surge feedback loop is damped. We also investigate the influence of different basal sliding laws (ranging from purely plastic to nonlinear to linear) on possible surging. The presented mechanisms underlying our simulations of self-maintained, periodic ice growth and destabilization may play a role in large-scale ice-sheet surging, such as the surging of the Laurentide Ice Sheet, which is associated with Heinrich events, and ice-stream shutdown and reactivation, such as observed in the Siple Coast region of West Antarctica.
- ItemEnhanced economic connectivity to foster heat stress-related losses(Washington, DC : American Association for the Advancement of Science, 2016) Wenz, Leonie; Levermann, AndersAssessing global impacts of unexpected meteorological events in an increasingly connected world economy is important for estimating the costs of climate change. We show that since the beginning of the 21st century, the structural evolution of the global supply network has been such as to foster an increase of climate-related production losses. We compute first- and higher-order losses from heat stress–induced reductions in productivity under changing economic and climatic conditions between 1991 and 2011. Since 2001, the economic connectivity has augmented in such a way as to facilitate the cascading of production loss. The influence of this structural change has dominated over the effect of the comparably weak climate warming during this decade. Thus, particularly under future warming, the intensification of international trade has the potential to amplify climate losses if no adaptation measures are taken.
- ItemGlobal gridded crop model evaluation: Benchmarking, skills, deficiencies and implications(München : European Geopyhsical Union, 2017) Müller, Christoph; Elliott, Joshua; Chryssanthacopoulos, James; Arneth, Almut; Balkovic, Juraj; Ciais, Philippe; Deryng, Delphine; Folberth, Christian; Glotter, Michael; Hoek, Steven; Iizumi, Toshichika; Izaurralde, Roberto C.; Jones, Curtis; Khabarov, Nikolay; Lawrence, Peter; Liu, Wenfeng; Olin, Stefan; Pugh, Thomas A.M.; Ray, Deepak K.; Reddy, Ashwan; Rosenzweig, Cynthia; Ruane, Alex C.; Sakurai, Gen; Schmid, Erwin; Skalsky, Rastislav; Song, Carol X.; Wang, Xuhui; de Wit, Allard; Yang, HongCrop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.
- ItemEvaluating spatially resolved influence of soil and tree water status on quality of European plum grown in semi-humid climate(Lausanne : Frontiers Media, 2017) Käthner, Jana; Ben-Gal, Alon; Gebbers, Robin; Peeters, Aviva; Herppich, Werner B.; Zude-Sasse, ManuelaIn orchards, the variations of fruit quality and its determinants are crucial for resource effective measures. In the present study, a drip-irrigated plum production (Prunus domestica L. “Tophit plus”/Wavit) located in a semi-humid climate was studied. Analysis of the apparent electrical conductivity (ECa) of soil showed spatial patterns of sand lenses in the orchard. Water status of sample trees was measured instantaneously by means of leaf water potential, Ψleaf [MPa], and for all trees by thermal imaging of canopies and calculation of the crop water stress index (CWSI). Methods for determining CWSI were evaluated. A CWSI approach calculating canopy and reference temperatures from the histogram of pixels from each image itself was found to suit the experimental conditions. Soil ECa showed no correlation with specific leaf area ratio and cumulative water use efficiency (WUEc) derived from the crop load. The fruit quality, however, was influenced by physiological drought stress in trees with high crop load and, resulting (too) high WUEc, when fruit driven water demand was not met. As indicated by analysis of variance, neither ECa nor the instantaneous CWSI could be used as predictors of fruit quality, while the interaction of CWSI and WUEc did succeed in indicating significant differences. Consequently, both WUEc and CWSI should be integrated in irrigation scheduling for positive impact on fruit quality.
- ItemA new climate dataset for systematic assessments of climate change impacts as a function of global warming(München : European Geopyhsical Union, 2012) Heinke, J.; Ostberg, S.; Schaphoff, S.; Frieler, K.; Müller, C.; Gerten, D.; Meinshausen, M.; Lucht, W.Interactions between food demand, biomass energy and forest preservation are driving both food prices and land-use changes, regionally and globally. This study presents a new model called Nexus Land-Use version 1.0 which describes these interactions through a generic representation of agricultural intensification mechanisms within agricultural lands. The Nexus Land-Use model equations combine biophysics and economics into a single coherent framework to calculate crop yields, food prices, and resulting pasture and cropland areas within 12 regions inter-connected with each other by international trade. The representation of cropland and livestock production systems in each region relies on three components: (i) a biomass production function derived from the crop yield response function to inputs such as industrial fertilisers; (ii) a detailed representation of the livestock production system subdivided into an intensive and an extensive component, and (iii) a spatially explicit distribution of potential (maximal) crop yields prescribed from the Lund-Postdam-Jena global vegetation model for managed Land (LPJmL). The economic principles governing decisions about land-use and intensification are adapted from the Ricardian rent theory, assuming cost minimisation for farmers. In contrast to the other land-use models linking economy and biophysics, crops are aggregated as a representative product in calories and intensification for the representative crop is a non-linear function of chemical inputs. The model equations and parameter values are first described in details. Then, idealised scenarios exploring the impact of forest preservation policies or rising energy price on agricultural intensification are described, and their impacts on pasture and cropland areas are investigated.
- ItemPALADYN v1.0, a comprehensive land surface-vegetation-carbon cycle model of intermediate complexity(München : European Geopyhsical Union, 2016) Willeit, Matteo; Ganopolski, AndreyPALADYN is presented; it is a new comprehensive and computationally efficient land surface–vegetation–carbon cycle model designed to be used in Earth system models of intermediate complexity for long-term simulations and paleoclimate studies. The model treats in a consistent manner the interaction between atmosphere, terrestrial vegetation and soil through the fluxes of energy, water and carbon. Energy, water and carbon are conserved. PALADYN explicitly treats permafrost, both in physical processes and as an important carbon pool. It distinguishes nine surface types: five different vegetation types, bare soil, land ice, lake and ocean shelf. Including the ocean shelf allows the treatment of continuous changes in sea level and shelf area associated with glacial cycles. Over each surface type, the model solves the surface energy balance and computes the fluxes of sensible, latent and ground heat and upward shortwave and longwave radiation. The model includes a single snow layer. Vegetation and bare soil share a single soil column. The soil is vertically discretized into five layers where prognostic equations for temperature, water and carbon are consistently solved. Phase changes of water in the soil are explicitly considered. A surface hydrology module computes precipitation interception by vegetation, surface runoff and soil infiltration. The soil water equation is based on Darcy's law. Given soil water content, the wetland fraction is computed based on a topographic index. The temperature profile is also computed in the upper part of ice sheets and in the ocean shelf soil. Photosynthesis is computed using a light use efficiency model. Carbon assimilation by vegetation is coupled to the transpiration of water through stomatal conductance. PALADYN includes a dynamic vegetation module with five plant functional types competing for the grid cell share with their respective net primary productivity. PALADYN distinguishes between mineral soil carbon, peat carbon, buried carbon and shelf carbon. Each soil carbon type has its own soil carbon pools generally represented by a litter, a fast and a slow carbon pool in each soil layer. Carbon can be redistributed between the layers by vertical diffusion and advection. For the vegetated macro surface type, decomposition is a function of soil temperature and soil moisture. Carbon in permanently frozen layers is assigned a long turnover time which effectively locks carbon in permafrost. Carbon buried below ice sheets and on flooded ocean shelves is treated differently. The model also includes a dynamic peat module. PALADYN includes carbon isotopes 13C and 14C, which are tracked through all carbon pools. Isotopic discrimination is modelled only during photosynthesis. A simple methane module is implemented to represent methane emissions from anaerobic carbon decomposition in wetlands (including peatlands) and flooded ocean shelf. The model description is accompanied by a thorough model evaluation in offline mode for the present day and the historical period.
- ItemDeciphering the imprint of topology on nonlinear dynamical network stability(Bristol : IOP Publishing, 2017) Nitzbon, J.; Schultz, P.; Heitzig, J.; Hellmann, F.Coupled oscillator networks show complex interrelations between topological characteristics of the network and the nonlinear stability of single nodes with respect to large but realistic perturbations. We extend previous results on these relations by incorporating sampling-based measures of the transient behaviour of the system, its survivability, as well as its asymptotic behaviour, its basin stability. By combining basin stability and survivability we uncover novel, previously unknown asymptotic states with solitary, desynchronized oscillators which are rotating with a frequency different from their natural one. They occur almost exclusively after perturbations at nodes with specific topological properties. More generally we confirm and significantly refine the results on the distinguished role tree-shaped appendices play for nonlinear stability. We find a topological classification scheme for nodes located in such appendices, that exactly separates them according to their stability properties, thus establishing a strong link between topology and dynamics. Hence, the results can be used for the identification of vulnerable nodes in power grids or other coupled oscillator networks. From this classification we can derive general design principles for resilient power grids. We find that striving for homogeneous network topologies facilitates a better performance in terms of nonlinear dynamical network stability. While the employed second-order Kuramoto-like model is parametrised to be representative for power grids, we expect these insights to transfer to other critical infrastructure systems or complex network dynamics appearing in various other fields.