550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (902)
- Doctoral Thesis (194)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (21)
- Conference Proceeding (18)
- Part of a Book (9)
- Diploma Thesis (8)
- diplomthesis (7)
- Other (5)
- Report (5)
- Master's Thesis (4)
- Review (4)
- Periodical (3)
- Bachelor Thesis (2)
- Habilitation (2)
- Preprint (2)
Language
Keywords
- climate change (11)
- Klima (8)
- Climate change (7)
- Klimaänderung (7)
- Modellierung (7)
- COSMO-CLM (6)
- Klimawandel (6)
- Palaeoclimate (6)
- precipitation (6)
- Atmospheric chemistry (5)
Institute
- Geowissenschaften (778)
- Geowissenschaften / Geographie (144)
- Biodiversität und Klima Forschungszentrum (BiK-F) (62)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (56)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (10)
- Physik (8)
Semi-arid African ecosystems influence trends and variability in global terrestrial carbon dynamics. However, there are uncertainties in potential effects of future climates for semi-arid ecosystems, especially for niche ecosystems. At the same time, African ecosystems provide the livelihoods and ecosystem services for around 1.4 billion people. Future population growth and associated changes in land use pose a challenge for the protection of African biodiversity. Therefore, this work focussed on future impacts of climate change on African ecosystems and carbon dynamics and also for African protected areas (PAs), where they may cooccur with other global change factors. Another focus was on uncertainties associated with future projections and with modelling the Nama Karoo, as an example of a semi-arid niche ecosystem. Dynamic vegetation models (DVMs) were the main research tool.
In Chapter 2, we analysed climate change impacts on African ecosystems and carbon pools until the end of the 21st century and associated uncertainties based on an ensemble of vegetation simulations with the DVM adaptive dynamic vegetation model (aDGVM). We investigated the impact of increased atmospheric CO2 concentrations and two climate change scenarios (medium (RCP4.5) and high emissions (RCP8.5); RCP - representative concentration pathway) on vegetation changes. Differences in the simulated vegetation were primarily driven by assumptions about the influence of CO2 on plants. Elevated CO2 concentrations led to increased total aboveground vegetation biomass and shrub encroachment into grasslands and savannas for both climate scenarios. In simulations without the direct influence of CO2 on plants, there was hardly any shrub encroachment and vegetation biomass decreased or varied between a slight decrease in some cases and a slight increase in others. Based on these results, biome changes due to climate change are likely in Africa in the future. Due to the large uncertainties in future projections, strategies to adapt to climate change must be flexible.
The simulated vegetation in Chapter 2 represented potential, natural vegetation and is particularly suitable to investigate PAs. However, PAs do not exist isolated from their environment and social developments. In Chapter 3, the vegetation projections with CO2 effect from Chapter 2 were combined with projections for population density and land use. Except for many PAs in North Africa, most PAs were adversely affected by at least one of the three drivers by the end of the 21st century in both investigated scenarios ("middle-of-the-road" and "fossil-fuelled development"). Cooccurrence of the drivers varied by region and scenario for PAs. Both scenarios implied increasing challenges for the conservation of African biodiversity in PAs. The impact of climate change on vegetation is likely to be exacerbated by socio-economic change for most African PAs. Strong mitigation of future climate change together with equitable societal development may facilitate successful ecosystem conservation.
The simulations in Chapters 2 and 3 showed large-scale patterns of vegetation change, but their low resolution makes them unsuitable for local analyses. In Chapter 4, the challenges of simulating smaller scale, semi-arid ecosystems and their carbon cycle were analysed for the Nama Karoo with the aDGVM2 and its shrub module. The aDGVM2 is based on the aDGVM, but represents plants more flexibly. In all tested aDGVM2 configurations, the carbon fluxes improved compared to initial simulations but still overestimated them. The measured morphology of the dwarf shrubs and soil water dynamics were not reproduced in aDGVM2. Semi-arid soil water dynamics and coping strategies of semi-arid dwarf shrubs under drought stress are not adequately implemented in the aDGVM2. Further field research on semi-arid water and carbon dynamics of vegetation is necessary to parameterise the aDGVM2 for dwarf shrubs. If these challenges are overcome, DVMs can be a powerful tool for much-needed research on the impacts of climate change on the Nama Karoo.
The analyses have shown that climate change under medium to high emission scenarios is likely to lead to large-scale changes in ecosystems and the carbon balance in Africa. Because lower emissions scenarios come with less uncertainty, climate change adaptation strategies likely need to be less complex or extensive if climate change is minimised. For African PAs, the challenges of climate change may be exacerbated by socio-economic factors to a regionally varying extent. This research suggests that successful ecosystem conservation depends on climate change mitigation measures and ensuring equitable, sustainable development. The shown uncertainties, e.g., in the implementation of the CO2 effect on plants or vegetation dynamics in more niche ecosystems, help to focus future research efforts and increase our understanding of the range of plausible futures we may need to adapt to.
First-principles modeling techniques offer the ability to simulate a wide range of systems under different physical conditions, such as temperature, pressure, and composition, without relying on empirical knowledge. Density functional theory (DFT), a quantum mechanical method, has become an exceptionally successful framework for materials science modeling. Employing DFT makes it possible to gain valuable insights into the fundamental state of a system, enabling the reliable determination of equilibrium crystal structures. Over time, DFT has become an essential tool that can be incorporated into various schemes for predicting the properties of a material related to its structure, insulating/metallic behavior, magnetism, and optics. DFT is regularly applied in numerous fields, spanning from fundamental subjects in condensed matter physics to the study of large-scale phenomena in geosciences. In the latter, the effectiveness of DFT stems from its ability to simulate the properties found on the Earth, other planets, and meteorites, which may pose challenges for their direct study or laboratory investigation.
In this thesis, a comprehensive examination of a family of monosulfides and a perovskite heterostructure was conducted. These materials are relevant for their potential applications in technology, energy harvesting, and in the case of monosulfides, their speculated abundance on the planet Mercury.
Firstly, a DFT approach was used to analyze two non-magnetic monosulfides, CaS and MgS. We determined their structural properties and then focused on the modeling of their reflectivity in the infrared region. The calculation of the reflectivity considered both harmonic and anharmonic contributions. In the harmonic limit, the non-analytic correction was employed to accurately determine the LO/TO splitting, which is necessary to delimit the retstrahlend band, that is, the maximum of the reflectivity. The anharmonic effects given by up to three-phonon and isotopic scatterings, which were included using perturbation theory, primarily smeared the reflectivity spectra edges in the high-wave region.
Secondly, four polymorphs of MnS were studied using a combination of first-principles methods to simulate their antiferromagnetic (AFM) and paramagnetic (PM) states. The integration of DFT+$U$ with special quasirandom structures (SQS) supercells, and occupation matrix control techniques was crucial for achieving convergence, structural optimization accuracy, and obtaining finite energy band gaps and local magnetic moments in the PM phases. The addition of the Hubbard $U$ correction was necessary to treat the highly-correlated Mn $d$-electrons. The success of our approach was clear based on our electronic structure predictions for the PM rock-salt B1-MnS polymorph. Experimentally this phase has been observed to be an insulator, but multiple \emph{ab initio} works resulted previously in metallic behavior. Our computations, on the other hand, predicted insulating and magnetic properties that compare well with available measurements. Additionally, the pressure-field stability of the four MnS polymorphs was studied. In the case of the PM phases, B1-MnS was identified to be the most stable up to about 21 GPa, then transforming into the B31-MnS polymorph. This finding was in close agreement with high-pressure experiments reporting a similar phase transformation. The optical properties of B1-, B4-, and B31-MnS were also simulated. The SQS technique was used to obtain soft-mode-free phonon band structures within the harmonic approximation. Then, the anharmonic effects were included, and the reflectivity was calculated for B1-MnS and B4-MnS. In both cases, a good agreement for the LO/TO splitting was achieved in comparison to experimental results.
Lastly, the oxygen-deficient heterostructure of LaAlO$_{3-\delta}$ /SrTiO$_{3-\delta}$ was investigated also employing DFT+$U$, with a particular emphasis on the potential impact of vacancy clustering at the interface. Six distinct configurations of pairs of vacancies were studied and their energies were compared to find the most stable one. The orbital reconstruction of Ti orbitals was also examined based on their location with respect to the vacancies and the local magnetic moments were calculated. The final results showed that linearly arranged vacancies located opposite to Ti ions give the most energetically stable configuration.
Climatology of morphology and cloud-radiative properties of marine low-level mixed-phase clouds
(2023)
Marine stratocumuli cover about 40 - 60% of the ocean surface. They self-organize into different morphological regimes. The two organized cellular regimes are called open and closed mesoscale-cellular convective (MCC) clouds. In mid-to-high latitudes, open and closed cells are the two most frequent types of MCC clouds. In particular, many MCC clouds consist of a mixture of vapor, liquid droplets, and ice particles, referred to as mixed-phase clouds (MPCs). Even for the same cloud fraction, the albedo of open cells is, on average, lower than that of closed MCC clouds. Cloud phase and morphology individually influence the cloud radiative effect. Thus, this thesis investigates the relationships between the cloud phase, MCC organization, cell size, and differences regarding the cloud-radiative effect.
This thesis focuses on space-borne retrievals to achieve extensive temporal and spatial coverage. The liDAR-raDAR (DARDAR) version 2 product collocates two active and one passive satellite: CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Moderate Resolution Imaging Spectroradiometer (MODIS). The cloud phase of DARDAR is vertically integrated to establish a single cloud phase at each data point. The MCC classification data set based on the liquid water path (LWP) of MODIS scenes is collocated with the DARDAR product to determine the MCC organization. Cell-size statistics of both MCC clouds are obtained using a marker-based image segmentation method on MODIS reflectance scenes. In addition, based on MODIS reflectance scenes, a convolutional neural network (CNN) is developed to classify open and closed MCC scenes to avoid missing mature MPCs with a low LWP.
The first part of this thesis explores the relationships between cloud phase, morphology, and cloud albedo in the Southern Ocean (SO). At a given cloud-top temperature (CTT), seasonal changes in the mixed-phase fraction, defined as the number of MPCs divided by the sum of MPC and supercooled liquid cloud (SLC) pixels, are stronger than the morphological changes. Therefore, external factors seem to influence these changes instead of morphology. The dependence of cloud phase on cloud-top height (CTH) is more substantial than on CTT in clouds with CTHs below 2.5 km. The previously observed acceleration of closed-to-open transition in MPCs, known as preconditioning, is not the primary driver of climatological cloud morphology statistics in the SO. The morphological differences in cloud albedo are more pronounced in SLCs than in MPCs. This change in albedo alters the cloud radiative effect in the SO by 21Wm−2 to 39Wm−2 depending onseason and cloud phase.
Open and closed MCC clouds exhibit larger equivalent cell diameters in the MPCs than in SLCs in austral summer, whereas, in austral winter, the SLCs are larger. The cell’s aspect ratio accounts for varying CTHs. Closed cells have smaller aspect ratios than open cells, so their cell diameter is smaller, independent of CTH. While the seasonal differences in closed cells are due to changes in CTH, the seasonal aspect ratio differences in open cells are mainly caused by MPCs. With increasing aspect ratios, the cloud albedo decreases in both open and closed MCC clouds, with the most substantial decrease in open MPCs clouds. This leads to cloud-radiative changes of 60 - 75Wm−2 in the SO, depending on cloud phase and aspect ratio.
The established CNN exhibits a good accuracy of 80.6%, with even higher accuracies in the Open (85.5%) and Closed (87.3%) categories. The global MCC climatology based on the CNN generally agrees well with previous MCC distributions. The most notable difference occurs in the Northern Hemisphere (NH) in boreal winter, with a higher occurrence frequency of closed and open MCC clouds. This might indicate missing MPCs in previous studies based on the LWP and some restricted to warm cloud scenes. Thus, the developed CNN seems to better represent the different morphologies in MPCs than in previous classifications.
In conclusion, this thesis shows that understanding the dependencies of cloud phase, cloud morphology, and cell size is important to enhance predictions of the cloud-radiative effect and thus, it is important to evaluate how cloud phase, cloud morphology, and cellsize change in a warming climate.
Due to massive energetic investments in woody support structures, trees are subject to unique physiological, mechanical, and ecological pressures not experienced by herbaceous plants. Despite a wealth of studies exploring trait relationships across the entire plant kingdom, the dominant traits underpinning these unique aspects of tree form and function remain unclear. Here, by considering 18 functional traits, encompassing leaf, seed, bark, wood, crown, and root characteristics, we quantify the multidimensional relationships in tree trait expression. We find that nearly half of trait variation is captured by two axes: one reflecting leaf economics, the other reflecting tree size and competition for light. Yet these orthogonal axes reveal strong environmental convergence, exhibiting correlated responses to temperature, moisture, and elevation. By subsequently exploring multidimensional trait relationships, we show that the full dimensionality of trait space is captured by eight distinct clusters, each reflecting a unique aspect of tree form and function. Collectively, this work identifies a core set of traits needed to quantify global patterns in functional biodiversity, and it contributes to our fundamental understanding of the functioning of forests worldwide.
Die nachfolgende Dissertation wurde an der Goethe-Universität Frankfurt am Institut für Geowissenschaften (FB 11) in der Arbeitsgruppe Kristallographie und Mineralogie (AG Winkler) verfasst. Die Betreuung der hier durchgeführten Arbeiten erfolgte hauptsächlich durch Prof. B. Winkler in Zusammenarbeit mit Dr. L. Bayarjargal, PD Dr. E. Haussühl und PD Dr. V. Vinograd. Bei dem vorliegenden Manuskript handelt es sich um eine kumulative bzw. publikationsbasierte Dissertation, welche die Forschungsergebnisse verschiedener bereits veröffentlichter wissenschaftlicher Fachartikel zusammenfasst.
Die Arbeit beschreibt verschiedene Synthesen und Untersuchungen an Carbonaten und teilt sich im Wesentlichen in zwei Abschnitte. Zum einen wurden Experimente mit Carbonaten bei Extrembedingungen bzw. unter hohen Drücken und hohen Temperaturen durchgeführt, wie sie auch im Inneren der Erde zu finden sind. Im zweiten Teil wurden Carbonate bei Raumbedingungen synthetisiert und der Einbau von Seltenerdelementen untersucht. Grundsätzlich werden jedoch in beiden Teilen dieser Arbeit die Strukturen und Eigenschaften verschiedener Carbonate und eine mögliche Kationensubstitution bzw. die Synthese isostruktureller Verbindungen erforscht.
The Izu–Bonin–Mariana volcanic arc is situated at a convergent plate margin where subduction initiation triggered the formation of MORB-like forearc basalts as a result of decompression melting and near-trench spreading. International Ocean Discovery Program (IODP) Expedition 352 recovered samples within the forearc basalt stratigraphy that contained unusual macroscopic globular textures hosted in andesitic glass (Unit 6, Hole 1440B). It is unclear how these andesites, which are unique in a stratigraphic sequence dominated by forearc basalts, and the globular textures therein may have formed. Here, we present detailed textural evidence, major and trace element analysis, as well as B and Sr isotope compositions, to investigate the genesis of these globular andesites. Samples consist of K2O-rich basaltic globules set in a glassy groundmass of andesitic composition. Between these two textural domains a likely hydrated interface of devitrified glass occurs, which, based on textural evidence, seems to be genetically linked to the formation of the globules. The andesitic groundmass is Cl rich (ca. 3000 µg/g ), whereas globules and the interface are Cl poor (ca. 300 µg/g ). Concentrations of fluid-mobile trace elements also appear to be fractionated in that globules and show enrichments in B, K, Rb, Cs, and Tl, but not in Ba and W relative to the andesitic groundmass, whereas the interface shows depletions in the latter, but is enriched in the former. Interestingly, globules and andesitic groundmass have identical Sr isotopic composition within analytical uncertainty ( 87Sr∕86Sr of 0.70580 ± 10 ), indicating that they likely formed from the same source. However, globules show high δ11 B (ca. + 7 ‰ ), whereas their host andesites are isotopically lighter (ca. – 1 ‰ ), potentially indicating that whatever process led to their formation either introduced heavier B isotopes to the globules, or induced stable isotope fractionation of B between globules and their groundmass. Based on the bulk of the textural information and geochemical data obtained from these samples, we conclude that these andesites likely formed as a result of the assimilation of shallowly altered oceanic crust (AOC) during forearc basaltic magmatism. Assimilation likely introduced radiogenic Sr, as well as heavier B isotopes to comparatively unradiogenic and low δ11B forearc basalt parental magmas (average 87Sr∕86Sr of 0.703284). Moreover, the globular textures are consistent with their formation being the result of fluid-melt immiscibility that was potentially induced by the rapid release of water from assimilated AOC whose escape likely formed the interface. If the globular textures present in these samples are indeed the result of fluid-melt immiscibility, then this process led to significant trace element and stable isotope fractionation. The textures and chemical compositions of the globules highlight the need for future experimental studies aimed at investigating the exsolution process with respect to potential trace element and isotopic fractionation in arc magmas that have perhaps not been previously considered.
Knowledge about the initial tectonic and depositional dynamics, as well as the influence of early rifting on climate and environmental evolution remains speculative to a large extent, because sediments are usually deeply buried. Within the East African Rift System, inversion tectonics uplifted a few of these successions to the surface hence presenting rare windows into the pre-rift depositional history. One such example, an exceptional 700 m long and up to 60 m high fresh road cut provided the opportunity to study in detail initial rift successions of the southern Albertine Rift (Western Uganda). This focusses on the basal and poorly known Middle to Late Miocene in order to unravel the climatic, environmental, hydrological and tectonic evolution of the initial Albertine Rift. A large and robust multi-proxy dataset was gathered comprising 169 m of stratigraphic thickness, which spans from 14.5 to 4.9 Ma according to a revised lithostratigraphic model. Fieldwork comprised logging of the sedimentary record, spectral gamma ray, magnetic susceptibility and 2D wall mapping with photomosaics. Additionally, the sections were sampled for bulk mineral and clay mineral analysis. The succession exposes a suite of lithofacies and architectural elements detailing the evolution of a fluvio-lacustrine system. Five depositional environments were identified which show an overall back-stepping trend from an alluvial plain to a delta plain and finally palustrine/shallow lacustrine conditions. Mesoscale base-level cycles, preservation potential of architectural elements, and stacking pattern exhibit limited accommodation space. However, it increases over time. This overall trend indicates increasing tectonic subsidence, which can be explained by flexural downwarp within the pre-rift phase and in the upper part grading into fault-controlled crustal extension of the syn-rift phase, which more and more disrupted a large-scale river system. From the Middle Miocene up to the early Pliocene, this study revealed that palaeoclimate trends become marked by increasing and more fluctuating Th concentrations, loss of feldspar, intercalated lenses of hydroxosulphate minerals, and a shift from smectite-dominated to kaolinite-dominated clays. These signals are all interpreted as detrital except for the hydroxosulphates, and they mirror the increasing intensity of chemical weathering and stripping of soils in the catchment. A trend towards increasing humidity is supported by an increase in lacustrine sediment facies and a lake-level rise. Nevertheless, intercalation of hydroxosulphate, ferricretes and pedogenised horizons prove ongoing seasonality and dry intervals. Finally, based on a revised stratigraphic model a sequence stratigraphic correlation of the outcrop's depositional cycles with basin-scale cycles is presented. According to these cycles, transition from the pre-rift to the syn-rift stage is marked by an unconformity and a tectonic pulse in the latest Miocene. However, the response of fluvial supply, the depositional system as well as climate conditions are less punctuated and characterised by gradual trends and temporal delays. The long pre-rift phase (ca 10 Myr) and the gradual transition to the syn-rift phase is in accordance with the active rifting model, which is based on thermal thinning of the lithosphere by asthenospheric upwelling.
Herein, the high-pressure/high-temperature synthesis (11 GPa, 650 °C) of Tb3B10O17(OH)5 in a modified Walker-type multianvil device is presented. The structure of this rare-earth borate was determined by single-crystal X-ray diffraction methods and was found to crystallize orthorhombically in the space group Pmn21 (no. 31) with the unit cell parameters a = 16.2527(4), b = 4.4373(1), and c = 8.8174(2) Å. The new compound was further characterized using infrared spectroscopy, energy-dispersive X-ray spectroscopy, second harmonic generation (SHG) measurements, and temperature-dependent X-ray powder diffraction. Tb3B10O17(OH)5 decomposes to β-Tb(BO2)3 at temperatures higher than 460 °C. With increasing temperatures, the formation of μ-TbBO3 was observed, which transforms to π-TbBO3 upon cooling.
High-pressure/high-temperature synthesis of the new boron-rich terbium hydroxyborate Tb3B12O19(OH)7
(2023)
Monoclinic Tb3B12O19(OH)7 was obtained by multianvil high-pressure/high-temperature syntheses at 6 GPa and 650 °C. The crystal structure was investigated by single-crystal X-ray diffraction methods and space group C2 (no. 5) with the unit cell parameters a = 24.2299(5) Å, b = 4.4667(1) Å, c = 7.0964(2) Å, β = 94.58(1)°, and two formula units per cell were revealed. Powder X-ray diffraction, infrared spectroscopy and the investigation of its second harmonic generation properties support the proposed structural model.
In this dissertation, different aspects of turbulent transport and thermally driven flows over complex terrain are investigated. Two publications concentrate on the vertical heat and moisture exchange in the convective boundary layer over mountainous terrain. To study this, Large-Eddy Simulation (LES) is used. Both turbulent and advective transport mechanisms are evaluated over the simple orography of a quasi-two-dimensional, periodic valley with prescribed surface fluxes. Here, terrain elevation varies along only one of the horizontal coordinate axes. Even a relatively shallow orography, possibly unresolved in existing numerical weather prediction models, modifies the domain-averaged moisture and temperature profiles. For the analysis, the flow is decomposed into a local turbulent part, a local mean circulation, and a large-scale part. An analysis of the turbulent kinetic energy and turbulent heat and moisture flux budgets shows that the thermal circulation significantly contributes to the vertical transport. It is found that thermal upslope winds are important for the moisture transport from the valley to the mountain tops. In total, moisture export out of the valley is mostly accomplished by the mean circulation. On the temperature distribution, which is horizontally relatively homogeneous, the thermal circulation has a weaker impact. If an upper-level wind is present, it interacts with the thermal circulation. This weakens the vertical transport of moisture and thus reduces its export out of the valley. The heat transport is less affected by the upper-level wind because of its weaker dependence on the thermal circulation. These findings were corroborated in a more realistic experiment simulating the full diurnal cycle using radiation forcing and an interactive land surface model.
Based on these results, coherent turbulent structures in the convective boundary layer over non-flat terrain are studied in further detail. A conditional sampling method based on the concentration of a decaying passive tracer is implemented in order to identify the boundary-layer plumes objectively. Conditional sampling allows to quantify the contribution of plume structures to the vertical transport of heat and moisture. In case of the idealized valley, vertical transport by coherent structures is the dominant contribution to the turbulent components of both heat and moisture flux. It is comparable in magnitude to the advective transport by the mean slope-wind circulation, although it is more important for heat than for moisture transport. A set of less idealized simulations considers the flow over three-dimensional terrain. In this case, conditional sampling is carried out by using a simple domain-decomposition approach. We demonstrate that thermal updrafts are generally more frequent on hill tops than over the surroundings, but they are less persistent on the windward sides when large-scale winds are present in the free atmosphere.
The tools for flow decomposition and budget analysis are also applied in another idealized case with a quasi-two-dimensional valley featuring the stable boundary layer. Here, the formation of a low stratus cloud is investigated. The main driver for the cloud formation is radiative cooling due to outgoing longwave radiation. Despite a purely horizontal flow, the advection terms in the prognostic equations for heat and moisture produce vertical mixing across the upper cloud edge leading to a loss of cloud water content. However, this behavior is not due to any kind of thermally-driven circulation. Instead, this spurious mixing is caused by the diffusive error of the advection scheme in regions where the sloping surfaces of the terrain-following vertical coordinate intersect the cloud top. It is shown that the intensity of the (spurious) numerical diffusion strongly depends on the horizontal resolution, the order of advection, and the choice of the scalar advection scheme. A LES with 4 m horizontal resolution serves as a reference. For horizontal resolutions of a few hundred meters, carried out with a model setup as it is used in Numerical Weather Prediction, a strong reduction of the simulated liquid-water path is observed. In order to keep the (spurious) numerical diffusion at coarser resolutions small, at least a fifth-order advection scheme should be used. In the present case, a WENO scalar advection scheme turns out to increase the numerical diffusion along a sharp cloud edge compared to an upwind scheme. Furthermore, the choice of the vertical coordinate has a strong impact on the simulated liquid-water path over orography. With a modified definition of the terrain-following sigma coordinate, it is possible to produce cloud water where the classical sigma coordinate does not allow any cloud formation.