364 results
Search Results
Now showing 1 - 10 of 364
- 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.
- ItemASAMgpu V1.0 - A moist fully compressible atmospheric model using graphics processing units (GPUs)(München : European Geopyhsical Union, 2012) Horn, S.In this work the three dimensional compressible moist atmospheric model ASAMgpu is presented. The calculations are done using graphics processing units (GPUs). To ensure platform independence OpenGL and GLSL are used, with that the model runs on any hardware supporting fragment shaders. The MPICH2 library enables interprocess communication allowing the usage of more than one GPU through domain decomposition. Time integration is done with an explicit three step Runge-Kutta scheme with a time-splitting algorithm for the acoustic waves. The results for four test cases are shown in this paper. A rising dry heat bubble, a cold bubble induced density flow, a rising moist heat bubble in a saturated environment, and a DYCOMS-II case.
- 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.
- ItemOn the influence of wind and waves on underwater irradiance fluctuations(München : European Geopyhsical Union, 2012) Hieronymi, M.; Macke, A.The influence of various wind and wave conditions on the variability of downwelling irradiance Ed (490 nm) in water is subject of this study. The work is based on a two-dimensional Monte Carlo radiative transfer model with high spatial resolution. The model assumes conditions that are ideal for wave focusing, thus simulation results reveal the upper limit for light fluctuations. Local wind primarily determines the steepness of capillary-gravity waves which in turn dominate the irradiance variability near the surface. Down to 3 m depth, maximum irradiance peaks that exceed the mean irradiance Ed by a factor of more than 7 can be observed at low wind speeds up to 5 m s−1. The strength of irradiance fluctuations can be even amplified under the influence of higher ultra-gravity waves; thereby peaks can exceed 11 Ed. Sea states influence the light field much deeper; gravity waves can cause considerable irradiance variability even at 100 m depth. The simulation results show that under realistic conditions 50% radiative enhancements compared to the mean can still occur at 30 m depth. At greater depths, the underwater light variability depends on the wave steepness of the characteristic wave of a sea state; steeper waves cause stronger light fluctuations.
- 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.
- ItemNew developments in the representation of Saharan dust sources in the aerosol-climate model ECHAM6-HAM2(München : European Geopyhsical Union, 2016) Heinold, Bernd; Tegen, Ina; Schepanski, Kerstin; Banks, Jamie R.In the aerosol-climate model ECHAM6-HAM2, dust source activation (DSA) observations from Meteosat Second Generation (MSG) satellite are proposed to replace the original source area parameterization over the Sahara Desert. The new setup is tested in nudged simulations for the period 2007 to 2008. The evaluation is based on comparisons to dust emission events inferred from MSG dust index imagery, Aerosol Robotic Network (AERONET) sun photometer observations, and satellite retrievals of aerosol optical thickness (AOT). The model results agree well with AERONET measurements especially in terms of seasonal variability, and a good spatial correlation was found between model results and MSG-SEVIRI (Spinning-Enhanced Visible and InfraRed Imager) dust AOT as well as Multi-angle Imaging SpectroRadiometer (MISR) AOT. ECHAM6-HAM2 computes a more realistic geographical distribution and up to 20 % higher annual Saharan dust emissions, using the MSG-based source map. The representation of dust AOT is partly improved in the southern Sahara and Sahel. In addition, the spatial variability is increased towards a better agreement with observations depending on the season. Thus, using the MSG DSA map can help to circumvent the issue of uncertain soil input parameters. An important issue remains the need to improve the model representation of moist convection and stable nighttime conditions. Compared to sub-daily DSA information from MSG-SEVIRI and results from a regional model, ECHAM6-HAM2 notably underestimates the important fraction of morning dust events by the breakdown of the nocturnal low-level jet, while a major contribution is from afternoon-to-evening emissions.
- ItemHigh-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.
- ItemThe regional aerosol-climate model REMO-HAM(München : European Geopyhsical Union, 2012) Pietikäinen, J.-P.; O'Donnell, D.; Teichmann, C.; Karstens, U.; Pfeifer, S.; Kazil, J.; Podzun, R.; Fiedler, S.; Kokkola, H.; Birmili, W.; O'Dowd, C.; Baltensperger, U.; Weingartner, E.; Gehrig, R.; Spindler, G.; Kulmala, M.; Feichter, J.; Jacob, D.; Laaksonen, A.REMO-HAM is a new regional aerosol-climate model. It is based on the REMO regional climate model and includes most of the major aerosol processes. The structure for aerosol is similar to the global aerosol-climate model ECHAM5-HAM, for example the aerosol module HAM is coupled with a two-moment stratiform cloud scheme. On the other hand, REMO-HAM does not include an online coupled aerosol-radiation nor a secondary organic aerosol module. In this work, we evaluate the model and compare the results against ECHAM5-HAM and measurements. Four different measurement sites were chosen for the comparison of total number concentrations, size distributions and gas phase sulfur dioxide concentrations: Hyytiälä in Finland, Melpitz in Germany, Mace Head in Ireland and Jungfraujoch in Switzerland. REMO-HAM is run with two different resolutions: 50 × 50 km2 and 10 × 10 km2. Based on our simulations, REMO-HAM is in reasonable agreement with the measured values. The differences in the total number concentrations between REMO-HAM and ECHAM5-HAM can be mainly explained by the difference in the nucleation mode. Since we did not use activation nor kinetic nucleation for the boundary layer, the total number concentrations are somewhat underestimated. From the meteorological point of view, REMO-HAM represents the precipitation fields and 2 m temperature profile very well compared to measurement. Overall, we show that REMO-HAM is a functional aerosol-climate model, which will be used in further studies.
- ItemModeling of wave-induced irradiance variability in the upper ocean mixed layer(München : European Geopyhsical Union, 2012) Hieronymi, M.; Macke, A.; Zielinski, O.A Monte Carlo based radiative transfer model has been developed for calculating the availability of solar radiation within the top 100 m of the ocean. The model is optimized for simulations of spatial high resolution downwelling irradiance Ed fluctuations that arise from the lensing effect of waves at the water surface. In a first step the accuracy of simulation results has been verified by measurements of the oceanic underwater light field and through intercomparison with an established radiative transfer model. Secondly the potential depth-impact of nonlinear shaped single waves, from capillary to swell waves, is assessed by considering the most favorable conditions for light focusing, i.e. monochromatic light at 490 nm, very clear oceanic water with a low chlorophyll a content of 0.1 mg m−3 and high sun elevation. Finally light fields below irregular wave profiles accounting for realistic sea states were simulated. Our simulation results suggest that under open ocean conditions light flashes with 50% irradiance enhancements can appear down to 35 m depth, and light variability in the range of ±10% compared to the mean Ed is still possible in 100 m depth.
- ItemASAM v2.7: A compressible atmospheric model with a Cartesian cut cell approach(München : European Geopyhsical Union, 2015) Jähn, M.; Knoth, O.; König, M.; Vogelsberg, U.In this work, the fully compressible, three-dimensional, nonhydrostatic atmospheric model called All Scale Atmospheric Model (ASAM) is presented. A cut cell approach is used to include obstacles and orography into the Cartesian grid. Discretization is realized by a mixture of finite differences and finite volumes and a state limiting is applied. Necessary shifting and interpolation techniques are outlined. The method can be generalized to any other orthogonal grids, e.g., a lat–long grid. A linear implicit Rosenbrock time integration scheme ensures numerical stability in the presence of fast sound waves and around small cells. Analyses of five two-dimensional benchmark test cases from the literature are carried out to show that the described method produces meaningful results with respect to conservation properties and model accuracy. The test cases are partly modified in a way that the flow field or scalars interact with cut cells. To make the model applicable for atmospheric problems, physical parameterizations like a Smagorinsky subgrid-scale model, a two-moment bulk microphysics scheme, and precipitation and surface fluxes using a sophisticated multi-layer soil model are implemented and described. Results of an idealized three-dimensional simulation are shown, where the flow field around an idealized mountain with subsequent gravity wave generation, latent heat release, orographic clouds and precipitation are modeled.