47 results
Search Results
Now showing 1 - 10 of 47
- 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.
- 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.
- 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.
- ItemEvaluation of long-range transport and deposition of desert dust with the CTM MOCAGE(Milton Park : Taylor & Francis, 2017) Martet, M.; Peuch, V-H.; Laurent, B.; Marticorena, B.; Bergametti, G.Desert dust modelling and forecasting attract growing interest, due to the numerous impacts of dusts on climate, numerical weather prediction, health, ecosystems, transportation, as well as on many industrial activities. The validation of numerical tools is a very important activity in this context, and we present here an example of such an effort, combining in situ (horizontal visibility in SYNOP messages, IMPROVE database) and remote-sensing data (satellite imagery, AERONET aerosol optical thickness data). Interestingly, these measurements are available routinely, and not only in the context of dedicated measurements campaign; thus, they can be used in an operational context to monitor the performances of operational forecasting systems. MOCAGE is the chemistry-transport model of Météo-France, used operationally to forecast the three-dimensional transport of dusts and their deposition. Two very long-range transport episodes of dust have been studied: one case of Saharan dust transported to East America through Asia and Pacific observed in November 2004 and one case of Saharan dust transported from West Africa to Caribbean Islands in May 2007. Episodes of geographical extension had seldom been studied, and they provide a very selective reference to compare the modelled desert dusts with. The representation of dusts in MOCAGE appears to be realistic in these two very different cases. In turn, the model simulations are used to make the link between the complementary information provided by the different measurements tools, providing a fully consistent picture of the entire episodes. The evolution of the aerosol size distribution during the episodes has also been studied. With no surprise, our study underlines that deposition processes are very sensitive to the size of dust particles. If the atmospheric cycle, in terms of mass, is very much under the influence of larger particles (some micrometres and above), only the finer particles actually travel over thousands of kilometres. This illustrates the need for an accurate representation of size distributions for this aerosol component in numerical models and advocates for using a size-resolved (bin) approach as sinks, and particularly, deposition do not affect the emitted log-normal distributions symmetrically on both sides of the median diameter. Overall, the results presented in this study provide an evaluation of Météo-France operational dust forecasting system MOCAGE.
- ItemMeteorological and trace gas factors affecting the number concentration of atmospheric Aitken (DP Combining double low line 50 nm) particles in the continental boundary layer: Parameterization using a multivariate mixed effects model(München : European Geopyhsical Union, 2011) Mikkonen, S.; Korhonen, H.; Romakkaniemi, S.; Smith, J.N.; Joutsensaari, J.; Lehtinen, K.E.J.; Hamed, A.; Breider, T.J.; Birmili, W.; Spindler, G.; Plass-Duelmer, C.; Facchini, M.C.; Laaksonen, A.Measurements of aerosol size distribution and different gas and meteorological parameters, made in three polluted sites in Central and Southern Europe: Po Valley, Italy, Melpitz and Hohenpeissenberg in Germany, were analysed for this study to examine which of the meteorological and trace gas variables affect the number concentration of Aitken (Dp= 50 nm) particles. The aim of our study was to predict the number concentration of 50 nm particles by a combination of in-situ meteorological and gas phase parameters. The statistical model needs to describe, amongst others, the factors affecting the growth of newly formed aerosol particles (below 10 nm) to 50 nm size, but also sources of direct particle emissions in that size range. As the analysis method we used multivariate nonlinear mixed effects model. Hourly averages of gas and meteorological parameters measured at the stations were used as predictor variables; the best predictive model was attained with a combination of relative humidity, new particle formation event probability, temperature, condensation sink and concentrations of SO2, NO2 and ozone. The seasonal variation was also taken into account in the mixed model structure. Model simulations with the Global Model of Aerosol Processes (GLOMAP) indicate that the parameterization can be used as a part of a larger atmospheric model to predict the concentration of climatically active particles. As an additional benefit, the introduced model framework is, in theory, applicable for any kind of measured aerosol parameter.
- ItemThe Potsdam Parallel Ice Sheet Model (PISM-PIK) - Part 1: Model description(München : European Geopyhsical Union, 2011) Winkelmann, R.; Martin, M.A.; Haseloff, M.; Albrecht, T.; Bueler, E.; Khroulev, C.; Levermann, A.We present the Potsdam Parallel Ice Sheet Model (PISM-PIK), developed at the Potsdam Institute for Climate Impact Research to be used for simulations of large-scale ice sheet-shelf systems. It is derived from the Parallel Ice Sheet Model (Bueler and Brown, 2009). Velocities are calculated by superposition of two shallow stress balance approximations within the entire ice covered region: the shallow ice approximation (SIA) is dominant in grounded regions and accounts for shear deformation parallel to the geoid. The plug-flow type shallow shelf approximation (SSA) dominates the velocity field in ice shelf regions and serves as a basal sliding velocity in grounded regions. Ice streams can be identified diagnostically as regions with a significant contribution of membrane stresses to the local momentum balance. All lateral boundaries in PISM-PIK are free to evolve, including the grounding line and ice fronts. Ice shelf margins in particular are modeled using Neumann boundary conditions for the SSA equations, reflecting a hydrostatic stress imbalance along the vertical calving face. The ice front position is modeled using a subgrid-scale representation of calving front motion (Albrecht et al., 2011) and a physically-motivated calving law based on horizontal spreading rates. The model is tested in experiments from the Marine Ice Sheet Model Intercomparison Project (MISMIP). A dynamic equilibrium simulation of Antarctica under present-day conditions is presented in Martin et al. (2011).
- ItemSimilitude of ice dynamics against scaling of geometry and physical parameters(München : European Geopyhsical Union, 2016) Feldmann, Johannes; Levermann, AndersThe concept of similitude is commonly employed in the fields of fluid dynamics and engineering but rarely used in cryospheric research. Here we apply this method to the problem of ice flow to examine the dynamic similitude of isothermal ice sheets in shallow-shelf approximation against the scaling of their geometry and physical parameters. Carrying out a dimensional analysis of the stress balance we obtain dimensionless numbers that characterize the flow. Requiring that these numbers remain the same under scaling we obtain conditions that relate the geometric scaling factors, the parameters for the ice softness, surface mass balance and basal friction as well as the ice-sheet intrinsic response time to each other. We demonstrate that these scaling laws are the same for both the (two-dimensional) flow-line case and the three-dimensional case. The theoretically predicted ice-sheet scaling behavior agrees with results from numerical simulations that we conduct in flow-line and three-dimensional conceptual setups. We further investigate analytically the implications of geometric scaling of ice sheets for their response time. With this study we provide a framework which, under several assumptions, allows for a fundamental comparison of the ice-dynamic behavior across different scales. It proves to be useful in the design of conceptual numerical model setups and could also be helpful for designing laboratory glacier experiments. The concept might also be applied to real-world systems, e.g., to examine the response times of glaciers, ice streams or ice sheets to climatic perturbations.