75 results
Search Results
Now showing 1 - 10 of 75
- ItemImplementation of multirate time integration methods for air pollution modelling(München : European Geopyhsical Union, 2012) Schlegel, M.; Knoth, O.; Arnold, M.; Wolke, R.Explicit time integration methods are characterised by a small numerical effort per time step. In the application to multiscale problems in atmospheric modelling, this benefit is often more than compensated by stability problems and step size restrictions resulting from stiff chemical reaction terms and from a locally varying Courant-Friedrichs-Lewy (CFL) condition for the advection terms. Splitting methods may be applied to efficiently combine implicit and explicit methods (IMEX splitting). Complementarily multirate time integration schemes allow for a local adaptation of the time step size to the grid size. In combination, these approaches lead to schemes which are efficient in terms of evaluations of the right-hand side. Special challenges arise when these methods are to be implemented. For an efficient implementation, it is crucial to locate and exploit redundancies. Furthermore, the more complex programme flow may lead to computational overhead which, in the worst case, more than compensates the theoretical gain in efficiency. We present a general splitting approach which allows both for IMEX splittings and for local time step adaptation. The main focus is on an efficient implementation of this approach for parallel computation on computer clusters.
- 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.
- ItemA simple equation for the melt elevation feedback of ice sheets(Katlenburg-Lindau : European Geosciences Union, 2016) Levermann, Anders; Winkelmann, RicardaIn recent decades, the Greenland Ice Sheet has been losing mass and has thereby contributed to global sea-level rise. The rate of ice loss is highly relevant for coastal protection worldwide. The ice loss is likely to increase under future warming. Beyond a critical temperature threshold, a meltdown of the Greenland Ice Sheet is induced by the self-enforcing feedback between its lowering surface elevation and its increasing surface mass loss: the more ice that is lost, the lower the ice surface and the warmer the surface air temperature, which fosters further melting and ice loss. The computation of this rate so far relies on complex numerical models which are the appropriate tools for capturing the complexity of the problem. By contrast we aim here at gaining a conceptual understanding by deriving a purposefully simple equation for the self-enforcing feedback which is then used to estimate the melt time for different levels of warming using three observable characteristics of the ice sheet itself and its surroundings. The analysis is purely conceptual in nature. It is missing important processes like ice dynamics for it to be useful for applications to sea-level rise on centennial timescales, but if the volume loss is dominated by the feedback, the resulting logarithmic equation unifies existing numerical simulations and shows that the melt time depends strongly on the level of warming with a critical slowdown near the threshold: the median time to lose 10 % of the present-day ice volume varies between about 3500 years for a temperature level of 0.5 °C above the threshold and 500 years for 5 °C. Unless future observations show a significantly higher melting sensitivity than currently observed, a complete meltdown is unlikely within the next 2000 years without significant ice-dynamical contributions.
- ItemTobac 1.2: Towards a flexible framework for tracking and analysis of clouds in diverse datasets(Katlenburg-Lindau : Copernicus, 2019) Heikenfeld, Max; Marinescu, Peter J.; Christensen, Matthew; Watson-Parris, Duncan; Senf, Fabian; van den Heever, Susan C.; Stier, PhilipWe introduce tobac (Tracking and Object-Based Analysis of Clouds), a newly developed framework for tracking and analysing individual clouds in different types of datasets, such as cloud-resolving model simulations and geostationary satellite retrievals. The software has been designed to be used flexibly with any two-or three-dimensional timevarying input. The application of high-level data formats, such as Iris cubes or xarray arrays, for input and output allows for convenient use of metadata in the tracking analysis and visualisation. Comprehensive analysis routines are provided to derive properties like cloud lifetimes or statistics of cloud properties along with tools to visualise the results in a convenient way. The application of tobac is presented in two examples. We first track and analyse scattered deep convective cells based on maximum vertical velocity and the threedimensional condensate mixing ratio field in cloud-resolving model simulations. We also investigate the performance of the tracking algorithm for different choices of time resolution of the model output. In the second application, we show how the framework can be used to effectively combine information from two different types of datasets by simultaneously tracking convective clouds in model simulations and in geostationary satellite images based on outgoing longwave radiation. The tobac framework provides a flexible new way to include the evolution of the characteristics of individual clouds in a range of important analyses like model intercomparison studies or model assessment based on observational data. © 2019 Author(s).
- ItemRegional and inter-regional effects in evolving climate networks(Göttingen : Copernicus GmbH, 2014) Hlinka, J.; Hartman, D.; Jajcay, N.; Vejmelka, M.; Donner, R.; Marwan, N.; Kurths, J.; Paluš, M.Complicated systems composed of many interacting subsystems are frequently studied as complex networks. In the simplest approach, a given real-world system is represented by an undirected graph composed of nodes standing for the subsystems and non-oriented unweighted edges for interactions present among the nodes; the characteristic properties of the graph are subsequently studied and related to the system's behaviour. More detailed graph models may include edge weights, orientations or multiple types of links; potential time-dependency of edges is conveniently captured in so-called evolving networks. Recently, it has been shown that an evolving climate network can be used to disentangle different types of El Niño episodes described in the literature. The time evolution of several graph characteristics has been compared with the intervals of El Niño and La Niña episodes. In this study we identify the sources of the evolving network characteristics by considering a reduced-dimensionality description of the climate system using network nodes given by rotated principal component analysis. The time evolution of structures in local intra-component networks is studied and compared to evolving inter-component connectivity.
- ItemThe Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 2: Dynamic equilibrium simulation of the Antarctic ice sheet(München : European Geopyhsical Union, 2011) Martin, M.A.; Winkelmann, R.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.We present a dynamic equilibrium simulation of the ice sheet-shelf system on Antarctica with the Potsdam Parallel Ice Sheet Model (PISM-PIK). The simulation is initialized with present-day conditions for bed topography and ice thickness and then run to steady state with constant present-day surface mass balance. Surface temperature and sub-shelf basal melt distribution are parameterized. Grounding lines and calving fronts are free to evolve, and their modeled equilibrium state is compared to observational data. A physically-motivated calving law based on horizontal spreading rates allows for realistic calving fronts for various types of shelves. Steady-state dynamics including surface velocity and ice flux are analyzed for whole Antarctica and the Ronne-Filchner and Ross ice shelf areas in particular. The results show that the different flow regimes in sheet and shelves, and the transition zone between them, are captured reasonably well, supporting the approach of superposition of SIA and SSA for the representation of fast motion of grounded ice. This approach also leads to a natural emergence of sliding-dominated flow in stream-like features in this new 3-D marine ice sheet model.
- ItemApplication of a model-based rainfall-runoff database as efficient tool for flood risk management(Chichester : John Wiley and Sons Ltd, 2013) Brocca, L.; Liersch, S.; Melone, F.; Moramarco, T.; Volk, M.A framework for a comprehensive synthetic rainfall-runoff database was developed to study catchment response to a variety of rainfall events. The framework supports effective flood risk assessment and management and implements simple approaches. It consists of three flexible components, a rainfall generator, a continuous rainfallrunoff model, and a database management system. The system was developed and tested at two gauged river sections along the upper Tiber River (central Italy). One of the main questions was to investigate how simple such approaches can be applied without impairing the quality of the results. The rainfall-runoff model was used to simulate runoff on the basis of a large number of rainfall events. The resulting rainfallrunoff database stores pre-simulated events classified on the basis of the rainfall amount, initial wetness conditions and initial discharge. The real-time operational forecasts follow an analogue method that does not need new model simulations. However, the forecasts are based on the simulation results available in the rainfall-runoff database (for the specific class to which the forecast belongs). Therefore, the database can be used as an effective tool to assess possible streamflow scenarios assuming different rainfall volumes for the following days. The application to the study site shows that magnitudes of real flood events were appropriately captured by the database. Further work should be dedicated to introduce a component for taking account of the actual temporal distribution of rainfall events into the stochastic rainfall generator and to the use of different rainfall-runoff models to enhance the usability of the proposed procedure.
- ItemTreatment of non-ideality in the SPACCIM multiphase model - Part 1: Model development(München : European Geopyhsical Union, 2016) Rusumdar, A.J.; Wolke, R.; Tilgner, A.; Herrmann, H.Ambient tropospheric deliquesced particles generally comprise a complex mixture of electrolytes, organic compounds, and water. Dynamic modeling of physical and chemical processes in this complex matrix is challenging. Thus, up-to-date multiphase chemistry models generally do not consider non-ideal solution effects. Therefore, the present study was aimed at presenting further development of the SPACCIM (Spectral Aerosol Cloud Chemistry Interaction Model) through treatment of solution non-ideality, which has not been considered before. The present paper firstly describes the model developments including (i) the implementation of solution non-ideality in aqueous-phase reaction kinetics in the SPACCIM framework, (ii) the advancements in the coupling scheme of microphysics and multiphase chemistry and (iii) the required adjustments of the numerical schemes, especially in the sparse linear solver and the calculation of the Jacobian. Secondly, results of sensitivity investigations are outlined, aiming at the evaluation of different activity coefficient modules and the examination of the contributions of different intermolecular forces to the overall activity coefficients. Finally, first results obtained with the new model framework are presented. The SPACCIM parcel model was developed and, so far, applied for the description of aerosol–cloud interactions. To advance SPACCIM also for modeling physical and chemical processes in deliquesced particles, the solution non-ideality has to be taken into account by utilizing activities in reaction terms instead of aqueous concentrations. The main goal of the extended approach was to provide appropriate activity coefficients for solved species. Therefore, an activity coefficient module was incorporated into the kinetic model framework of SPACCIM. Based on an intercomparison of different activity coefficient models and the comparison with experimental data, the AIOMFAC approach was implemented and extended by additional interaction parameters from the literature for mixed organic–inorganic systems. Moreover, the performance and the capability of the applied activity coefficient module were evaluated by means of water activity measurements, literature data and results of other activity coefficient models. Comprehensive comparison studies showed that the SpactMod (SPACCIM activity coefficient module) is valuable for predicting the thermodynamic behavior of complex mixtures of multicomponent atmospheric aerosol particles. First simulations with a detailed chemical mechanism have demonstrated the applicability of SPACCIM-SpactMod. The simulations indicate that the treatment of solution non-ideality might be needed for modeling multiphase chemistry processes in deliquesced particles. The modeled activity coefficients imply that chemical reaction fluxes of chemical processes in deliquesced particles can be both decreased and increased depending on the particular species involved in the reactions. For key ions, activity coefficients on the order of 0.1–0.8 and a strong dependency on the charge state as well as the RH conditions are modeled, implying a lowered chemical processing of ions in concentrated solutions. In contrast, modeled activity coefficients of organic compounds are in some cases larger than 1 under deliquesced particle conditions and suggest the possibility of an increased chemical processing of organic compounds. Moreover, the model runs have shown noticeable differences in the pH values calculated with and without consideration of solution non-ideality. On average, the predicted pH values of the simulations considering solution non-ideality are −0.27 and −0.44 pH units lower under 90 and 70 % RH conditions, respectively. More comprehensive results of detailed SPACCIM-SpactMod studies on the multiphase processing in organic–inorganic mixtures of deliquesced particles are described in a companion paper.
- ItemScheme for calculation of multi-layer cloudiness and precipitation for climate models of intermediate complexity(München : European Geopyhsical Union, 2013) Eliseev, A.V.; Coumou, D.; Chernokulsky, A.V.; Petoukhov, V.; Petri, S.In this study we present a scheme for calculating the characteristics of multi-layer cloudiness and precipitation for Earth system models of intermediate complexity (EMICs). This scheme considers three-layer stratiform cloudiness and single-column convective clouds. It distinguishes between ice and droplet clouds as well. Precipitation is calculated by using cloud lifetime, which depends on cloud type and phase as well as on statistics of synoptic and convective disturbances. The scheme is tuned to observations by using an ensemble simulation forced by the ERA-40-derived climatology for 1979–2001. Upon calibration, the scheme realistically reproduces basic features of fields of cloud fractions, cloud water path, and precipitation. The simulated globally and annually averaged total cloud fraction is 0.59, and the simulated globally averaged annual precipitation is 100 cm yr−1. Both values agree with empirically derived values. The simulated cloud water path is too small, probably because the simulated vertical extent of stratiform clouds is too small. Geographical distribution and seasonal changes of calculated cloud fraction and precipitation are broadly realistic as well. However, some important regional biases still remain in the scheme, e.g. too little precipitation in the tropics. We discuss possibilities for future improvements in the scheme.