Search Results

Now showing 1 - 10 of 27
  • Item
    3+2 + X : what is the most useful depolarization input for retrieving microphysical properties of non-spherical particles from lidar measurements using the spheroid model of Dubovik et al. (2006)?
    (Katlenburg-Lindau : Copernicus, 2019) Tesche, Matthias; Kolgotin, Alexei; Haarig, Moritz; Burton, Sharon P.; Ferrare, Richard A.; Hostetler, Chris A.; Müller, Detlef
    The typical multiwavelength aerosol lidar data set for inversion of optical to microphysical parameters is composed of three backscatter coefficients (β) at 355, 532, and 1064 nm and two extinction coefficients (α) at 355 and 532 nm. This data combination is referred to as a 3β C 2α or 3 + 2 data set. This set of data is sufficient for retrieving some important microphysical particle parameters if the particles have spherical shape. Here, we investigate the effect of including the particle linear depolarization ratio (δ) as a third input parameter for the inversion of lidar data. The inversion algorithm is generally not used if measurements show values of d that exceed 0.10 at 532 nm, i.e. in the presence of nonspherical particles such as desert dust, volcanic ash, and, under special circumstances, biomass-burning smoke. We use experimental data collected with instruments that are capable of measuring d at all three lidar wavelengths with an inversion routine that applies the spheroidal light-scattering model of Dubovik et al. (2006) with a fixed axis-ratio distribution to replicate scattering properties of non-spherical particles. The inversion gives the fraction of spheroids required to replicate the optical data as an additional output parameter. This is the first systematic test of the effect of using all theoretically possible combinations of d taken at 355, 532, and 1064 nm as input in the lidar data inversion. We find that depolarization information of at least one wavelength already provides useful information for the inversion of optical data that have been collected in the presence of non-spherical mineral dust particles. However, any choice of d will give lower values of the single-scattering albedo than the traditional 3 + 2 data set. We find that input data sets that include d355 give a spheroid fraction that closely resembles the dust ratio we obtain from using β532 and d532 in a methodology applied in aerosol-type separation. The use of d355 in data sets of two or three d? reduces the spheroid fraction that is retrieved when using d532 and d1064. Use of the latter two parameters without accounting for d355 generally leads to high spheroid fractions that we consider not trustworthy. The use of three d instead of two δ, including the constraint that one of these is measured at 355 nm does not provide any advantage over using 3 + 2 + d355 for the observations with varying contributions of mineral dust considered here. However, additional measurements at wavelengths different from 355 nm would be desirable for application to a wider range of aerosol scenarios that may include non-spherical smoke particles, which can have values of d355 that are indistinguishable from those found for mineral dust. We therefore conclude that - depending on measurement capability - the future standard input for inversion of lidar data taken in the presence of mineral dust particles and using the spheroid model of Dubovik et al. (2006) might be 3+2Cδ355 or 3 + 2 + δ355 + δ532. © 2019 The Author(s).
  • Item
    Modeling the global emission, transport and deposition of trace elements associated with mineral dust
    (München : European Geopyhsical Union, 2015) Zhang, Y.; Mahowald, N.; Scanza, R.A.; Journet, E.; Desboeufs, K.; Albani, S.; Kok, J.F.; Zhuang, G.; Chen, Y.; Cohen, D.D.; Paytan, A.; Patey, M.D.; Achterberg, E.P.; Engelbrecht, J.P.; Fomba, K.W.
    Trace element deposition from desert dust has important impacts on ocean primary productivity, the quantification of which could be useful in determining the magnitude and sign of the biogeochemical feedback on radiative forcing. However, the impact of elemental deposition to remote ocean regions is not well understood and is not currently included in global climate models. In this study, emission inventories for eight elements primarily of soil origin, Mg, P, Ca, Mn, Fe, K, Al, and Si are determined based on a global mineral data set and a soil data set. The resulting elemental fractions are used to drive the desert dust model in the Community Earth System Model (CESM) in order to simulate the elemental concentrations of atmospheric dust. Spatial variability of mineral dust elemental fractions is evident on a global scale, particularly for Ca. Simulations of global variations in the Ca / Al ratio, which typically range from around 0.1 to 5.0 in soils, are consistent with observations, suggesting that this ratio is a good signature for dust source regions. The simulated variable fractions of chemical elements are sufficiently different; estimates of deposition should include elemental variations, especially for Ca, Al and Fe. The model results have been evaluated with observations of elemental aerosol concentrations from desert regions and dust events in non-dust regions, providing insights into uncertainties in the modeling approach. The ratios between modeled and observed elemental fractions range from 0.7 to 1.6, except for Mg and Mn (3.4 and 3.5, respectively). Using the soil database improves the correspondence of the spatial heterogeneity in the modeling of several elements (Ca, Al and Fe) compared to observations. Total and soluble dust element fluxes to different ocean basins and ice sheet regions have been estimated, based on the model results. The annual inputs of soluble Mg, P, Ca, Mn, Fe and K associated with dust using the mineral data set are 0.30 Tg, 16.89 Gg, 1.32 Tg, 22.84 Gg, 0.068 Tg, and 0.15 Tg to global oceans and ice sheets.
  • Item
    Year-round stratospheric aerosol backscatter ratios calculated from lidar measurements above northern Norway
    (Göttingen : Copernicus GmbH, 2019) Langenbach, A.; Baumgarten, G.; Fiedler, J.; Lübken, F.-J.; Von Savigny, C.; Zalach, J.
    We present a new method for calculating backscatter ratios of the stratospheric sulfate aerosol (SSA) layer from daytime and nighttime lidar measurements. Using this new method we show a first year-round dataset of stratospheric aerosol backscatter ratios at high latitudes. The SSA layer is located at altitudes between the tropopause and about 30 km. It is of fundamental importance for the radiative balance of the atmosphere. We use a state-of-the-art Rayleigh-Mie-Raman lidar at the Arctic Lidar Observatory for Middle Atmosphere Research (ALOMAR) station located in northern Norway (69N, 16E; 380ma.s.l.). For nighttime measurements the aerosol backscatter ratios are derived using elastic and inelastic backscatter of the emitted laser wavelengths 355, 532 and 1064nm. The setup of the lidar allows measurements with a resolution of about 5 min in time and 150 m in altitude to be performed in high quality, which enables the identification of multiple sub-layers in the stratospheric aerosol layer of less than 1 km vertical thickness. We introduce a method to extend the dataset throughout the summer when measurements need to be performed under permanent daytime conditions. For that purpose we approximate the backscatter ratios from color ratios of elastic scattering and apply a correction function. We calculate the correction function using the average backscatter ratio profile at 355nm from about 1700 h of nighttime measurements from the years 2000 to 2018. Using the new method we finally present a year-round dataset based on about 4100 h of measurements during the years 2014 to 2017. © Author(s) 2019.
  • Item
    Calibration of Raman lidar water vapor profiles by means of AERONET photometer observations and GDAS meteorological data
    (München : European Geopyhsical Union, 2018) Dai, Guangyao; Althausen, Dietrich; Hofer, Julian; Engelmann, Ronny; Seifert, Patric; Bühl, Johannes; Mamouri, Rodanthi-Elisavet; Wu, Songhua; Ansmann, Albert
    We present a practical method to continuously calibrate Raman lidar observations of water vapor mixing ratio profiles. The water vapor profile measured with the multiwavelength polarization Raman lidar PollyXT is calibrated by means of co-located AErosol RObotic NETwork (AERONET) sun photometer observations and Global Data Assimilation System (GDAS) temperature and pressure profiles. This method is applied to lidar observations conducted during the Cyprus Cloud Aerosol and Rain Experiment (CyCARE) in Limassol, Cyprus. We use the GDAS temperature and pressure profiles to retrieve the water vapor density. In the next step, the precipitable water vapor from the lidar observations is used for the calibration of the lidar measurements with the sun photometer measurements. The retrieved calibrated water vapor mixing ratio from the lidar measurements has a relative uncertainty of 11 % in which the error is mainly caused by the error of the sun photometer measurements. During CyCARE, nine measurement cases with cloud-free and stable meteorological conditions are selected to calculate the precipitable water vapor from the lidar and the sun photometer observations. The ratio of these two precipitable water vapor values yields the water vapor calibration constant. The calibration constant for the PollyXT Raman lidar is 6.56 g kg−1 ± 0.72 g kg−1 (with a statistical uncertainty of 0.08 g kg−1 and an instrumental uncertainty of 0.72 g kg−1). To check the quality of the water vapor calibration, the water vapor mixing ratio profiles from the simultaneous nighttime observations with Raman lidar and Vaisala radiosonde sounding are compared. The correlation of the water vapor mixing ratios from these two instruments is determined by using all of the 19 simultaneous nighttime measurements during CyCARE. Excellent agreement with the slope of 1.01 and the R2 of 0.99 is found. One example is presented to demonstrate the full potential of a well-calibrated Raman lidar. The relative humidity profiles from lidar, GDAS (simulation) and radiosonde are compared, too. It is found that the combination of water vapor mixing ratio and GDAS temperature profiles allow us to derive relative humidity profiles with the relative uncertainty of 10–20 %.
  • Item
    Target categorization of aerosol and clouds by continuous multiwavelength-polarization lidar measurements
    (Katlenburg-Lindau : Copernicus, 2017) Baars, Holger; Seifert, Patric; Engelmann, Ronny; Wandinger, Ulla
    Absolute calibrated signals at 532 and 1064 nm and the depolarization ratio from a multiwavelength lidar are used to categorize primary aerosol but also clouds in high temporal and spatial resolution. Automatically derived particle backscatter coefficient profiles in low temporal resolution (30 min) are applied to calibrate the lidar signals. From these calibrated lidar signals, new atmospheric parameters in temporally high resolution (quasi-particle-backscatter coefficients) are derived. By using thresholds obtained from multiyear, multisite EARLINET (European Aerosol Research Lidar Network) measurements, four aerosol classes (small; large, spherical; large, non-spherical; mixed, partly nonspherical) and several cloud classes (liquid, ice) are defined. Thus, particles are classified by their physical features (shape and size) instead of by source. The methodology is applied to 2 months of continuous observations (24 h a day, 7 days a week) with the multiwavelength-Raman-polarization lidar PollyXT during the High-Definition Clouds and Precipitation for advancing Climate Prediction (HD(CP)2) Observational Prototype Experiment (HOPE) in spring 2013. Cloudnet equipment was operated continuously directly next to the lidar and is used for comparison. By discussing three 24 h case studies, it is shown that the aerosol discrimination is very feasible and informative and gives a good complement to the Cloudnet target categorization. Performing the categorization for the 2-month data set of the entire HOPE campaign, almost 1 million pixel (5 min×30 m) could be analysed with the newly developed tool. We find that the majority of the aerosol trapped in the planetary boundary layer (PBL) was composed of small particles as expected for a heavily populated and industrialized area. Large, spherical aerosol was observed mostly at the top of the PBL and close to the identified cloud bases, indicating the importance of hygroscopic growth of the particles at high relative humidity. Interestingly, it is found that on several days non-spherical particles were dispersed from the ground into the atmosphere.
  • Item
    Schneefernerhaus as a mountain research station for clouds and turbulence
    (München : European Geopyhsical Union, 2015) Risius, S.; Xu, H.; Di Lorenzo, F.; Xi, H.; Siebert, H.; Shaw, R.A.; Bodenschatz, E.
    Cloud measurements are usually carried out with airborne campaigns, which are expensive and are limited by temporal duration and weather conditions. Ground-based measurements at high-altitude research stations therefore play a complementary role in cloud study. Using the meteorological data (wind speed, direction, temperature, humidity, visibility, etc.) collected by the German Weather Service (DWD) from 2000 to 2012 and turbulence measurements recorded by multiple ultrasonic sensors (sampled at 10 Hz) in 2010, we show that the Umweltforschungsstation Schneefernerhaus (UFS) located just below the peak of Zugspitze in the German Alps, at a height of 2650 m, is a well-suited station for cloud–turbulence research. The wind at UFS is dominantly in the east–west direction and nearly horizontal. During the summertime (July and August) the UFS is immersed in warm clouds about 25 % of the time. The clouds are either from convection originating in the valley in the east, or associated with synoptic-scale weather systems typically advected from the west. Air turbulence, as measured from the second- and third-order velocity structure functions that exhibit well-developed inertial ranges, possesses Taylor microscale Reynolds numbers up to 104, with the most probable value at ~ 3000. In spite of the complex topography, the turbulence appears to be nearly as isotropic as many laboratory flows when evaluated on the "Lumley triangle".
  • Item
    The Fifth International Workshop on Ice Nucleation phase 2 (FIN-02): Laboratory intercomparison of ice nucleation measurements
    (Katlenburg-Lindau : Copernicus, 2018) DeMott, Paul J.; Möhler, Ottmar; Cziczo, Daniel J.; Hiranuma, Naruki; Petters, Markus D.; Petters, Sarah S.; Belosi, Franco; Bingemer, Heinz G.; Brooks, Sarah D.; Budke, Carsten; Burkert-Kohn, Monika; Collier, Kristen N.; Danielczok, Anja; Eppers, Oliver; Felgitsch, Laura; Garimella, Sarvesh; Grothe, Hinrich; Herenz, Paul; Hill, Thomas C. J.; Höhler, Kristina; Kanji, Zamin A.; Kiselev, Alexei; Koop, Thomas; Kristensen, Thomas B.; Krüger, Konstantin; Kulkarni, Gourihar; Levin, Ezra J. T.; Murray, Benjamin J.; Nicosia, Alessia; O'Sullivan, Daniel; Peckhaus, Andreas; Polen, Michael J.; Price, Hannah C.; Reicher, Naama; Rothenberg, Daniel A.; Rudich, Yinon; Santachiara, Gianni; Schiebel, Thea; Schrod, Jann; Seifried, Teresa M.; Stratmann, Frank; Sullivan, Ryan C.; Suski, Kaitlyn J.; Szakáll, Miklós; Taylor, Hans P.; Ullrich, Romy; Vergara-Temprado, Jesus; Wagner, Robert; Whale, Thomas F.; Weber, Daniel; Welti, André; Wilson, Theodore W.; Wolf, Martin J.; Zenker, Jake
    The second phase of the Fifth International Ice Nucleation Workshop (FIN-02) involved the gathering of a large number of researchers at the Karlsruhe Institute of Technology's Aerosol Interactions and Dynamics of the Atmosphere (AIDA) facility to promote characterization and understanding of ice nucleation measurements made by a variety of methods used worldwide. Compared to the previous workshop in 2007, participation was doubled, reflecting a vibrant research area. Experimental methods involved sampling of aerosol particles by direct processing ice nucleation measuring systems from the same volume of air in separate experiments using different ice nucleating particle (INP) types, and collections of aerosol particle samples onto filters or into liquid for sharing amongst measurement techniques that post-process these samples. In this manner, any errors introduced by differences in generation methods when samples are shared across laboratories were mitigated. Furthermore, as much as possible, aerosol particle size distribution was controlled so that the size limitations of different methods were minimized. The results presented here use data from the workshop to assess the comparability of immersion freezing measurement methods activating INPs in bulk suspensions, methods that activate INPs in condensation and/or immersion freezing modes as single particles on a substrate, continuous flow diffusion chambers (CFDCs) directly sampling and processing particles well above water saturation to maximize immersion and subsequent freezing of aerosol particles, and expansion cloud chamber simulations in which liquid cloud droplets were first activated on aerosol particles prior to freezing. The AIDA expansion chamber measurements are expected to be the closest representation to INP activation in atmospheric cloud parcels in these comparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX and potassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam (Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively, and a bacterial INP (Snomax®). Considered together, the agreement among post-processed immersion freezing measurements of the numbers and fractions of particles active at different temperatures following bulk collection of particles into liquid was excellent, with possible temperature uncertainties inferred to be a key factor in determining INP uncertainties. Collection onto filters for rinsing versus directly into liquid in impingers made little difference. For methods that activated collected single particles on a substrate at a controlled humidity at or above water saturation, agreement with immersion freezing methods was good in most cases, but was biased low in a few others for reasons that have not been resolved, but could relate to water vapor competition effects. Amongst CFDC-style instruments, various factors requiring (variable) higher supersaturations to achieve equivalent immersion freezing activation dominate the uncertainty between these measurements, and for comparison with bulk immersion freezing methods. When operated above water saturation to include assessment of immersion freezing, CFDC measurements often measured at or above the upper bound of immersion freezing device measurements, but often underestimated INP concentration in comparison to an immersion freezing method that first activates all particles into liquid droplets prior to cooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber-Portable Ice Nucleation Chamber), and typically slightly underestimated INP number concentrations in comparison to cloud parcel expansions in the AIDA chamber; this can be largely mitigated when it is possible to raise the relative humidity to sufficiently high values in the CFDCs, although this is not always possible operationally. Correspondence of measurements of INPs among direct sampling and post-processing systems varied depending on the INP type. Agreement was best for Snomax® particles in the temperature regime colder than -10°C, where their ice nucleation activity is nearly maximized and changes very little with temperature. At temperatures warmer than -10°C, Snomax® INP measurements (all via freezing of suspensions) demonstrated discrepancies consistent with previous reports of the instability of its protein aggregates that appear to make it less suitable as a calibration INP at these temperatures. For Argentinian soil dust particles, there was excellent agreement across all measurement methods; measures ranged within 1 order of magnitude for INP number concentrations, active fractions and calculated active site densities over a 25 to 30°C range and 5 to 8 orders of corresponding magnitude change in number concentrations. This was also the case for all temperatures warmer than -25°C in Tunisian dust experiments. In contrast, discrepancies in measurements of INP concentrations or active site densities that exceeded 2 orders of magnitude across a broad range of temperature measurements found at temperatures warmer than -25°C in a previous study were replicated for illite NX. Discrepancies also exceeded 2 orders of magnitude at temperatures of -20 to -25°C for potassium feldspar (K-feldspar), but these coincided with the range of temperatures at which INP concentrations increase rapidly at approximately an order of magnitude per 2°C cooling for K-feldspar. These few discrepancies did not outweigh the overall positive outcomes of the workshop activity, nor the future utility of this data set or future similar efforts for resolving remaining measurement issues. Measurements of the same materials were repeatable over the time of the workshop and demonstrated strong consistency with prior studies, as reflected by agreement of data broadly with parameterizations of different specific or general (e.g., soil dust) aerosol types. The divergent measurements of the INP activity of illite NX by direct versus post-processing methods were not repeated for other particle types, and the Snomax° data demonstrated that, at least for a biological INP type, there is no expected measurement bias between bulk collection and direct immediately processed freezing methods to as warm as -10°C. Since particle size ranges were limited for this workshop, it can be expected that for atmospheric populations of INPs, measurement discrepancies will appear due to the different capabilities of methods for sampling the full aerosol size distribution, or due to limitations on achieving sufficient water supersaturations to fully capture immersion freezing in direct processing instruments. Overall, this workshop presents an improved picture of present capabilities for measuring INPs than in past workshops, and provides direction toward addressing remaining measurement issues.
  • Item
    Concentration and variability of ice nuclei in the subtropical maritime boundary layer
    (Katlenburg-Lindau : EGU, 2018) Welti, André; Müller, Konrad; Fleming, Zoë L.; Stratmann, Frank
    Measurements of the concentration and variability of ice nucleating particles in the subtropical maritime boundary layer are reported. Filter samples collected in Cabo Verde over the period 2009-2013 are analyzed with a drop freezing experiment with sensitivity to detect the few rare ice nuclei active at low supercooling. The data set is augmented with continuous flow diffusion chamber measurements at temperatures below -24 °C from a 2-month field campaign in Cabo Verde in 2016. The data set is used to address the following questions: what are typical concentrations of ice nucleating particles active at a certain temperature? What affects their concentration and where are their sources? Concentration of ice nucleating particles is found to increase exponentially by 7 orders of magnitude from -5 to -38 °C. Sample-to-sample variation in the steepness of the increase indicates that particles of different origin, with different ice nucleation properties (size, composition), contribute to the ice nuclei concentration at different temperatures. The concentration of ice nuclei active at a specific temperature varies over a range of up to 4 orders of magnitude. The frequency with which a certain ice nuclei concentration is measured within this range is found to follow a lognormal distribution, which can be explained by random dilution during transport. To investigate the geographic origin of ice nuclei, source attribution of air masses from dispersion modeling is used to classify the data into seven typical conditions. While no source could be attributed to the ice nuclei active at temperatures higher than -12 °C, concentrations at lower temperatures tend to be elevated in air masses originating from the Sahara.
  • Item
    Three-dimensional evolution of Saharan dust transport towards Europe based on a 9-year EARLINET-optimized CALIPSO dataset
    (Katlenburg-Lindau : EGU, 2017) Marinou, Eleni; Amiridis, Vassilis; Binietoglou, Ioannis; Tsikerdekis, Athanasios; Solomos, Stavros; Proestakis, Emannouil; Konsta, Dimitra; Papagiannopoulos, Nikolaos; Tsekeri, Alexandra; Vlastou, Georgia; Zanis, Prodromos; Balis, Dimitrios; Wandinger, Ulla; Ansmann, Albert
    In this study we use a new dust product developed using CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) observations and EARLINET (European Aerosol Research Lidar Network) measurements and methods to provide a 3-D multiyear analysis on the evolution of Saharan dust over North Africa and Europe. The product uses a CALIPSO L2 backscatter product corrected with a depolarization-based method to separate pure dust in external aerosol mixtures and a Saharan dust lidar ratio (LR) based on long-term EARLINET measurements to calculate the dust extinction profiles. The methodology is applied on a 9-year CALIPSO dataset (2007-2015) and the results are analyzed here to reveal for the first time the 3-D dust evolution and the seasonal patterns of dust over its transportation paths from the Sahara towards the Mediterranean and Continental Europe. During spring, the spatial distribution of dust shows a uniform pattern over the Sahara desert. The dust transport over the Mediterranean Sea results in mean dust optical depth (DOD) values up to 0.1. During summer, the dust activity is mostly shifted to the western part of the desert where mean DOD near the source is up to 0.6. Elevated dust plumes with mean extinction values between 10 and 75 Mm-1 are observed throughout the year at various heights between 2 and 6 km, extending up to latitudes of 40° N. Dust advection is identified even at latitudes of about 60° N, but this is due to rare events of episodic nature. Dust plumes of high DOD are also observed above the Balkans during the winter period and above northwest Europe during autumn at heights between 2 and 4 km, reaching mean extinction values up to 50 Mm-1. The dataset is considered unique with respect to its potential applications, including the evaluation of dust transport models and the estimation of cloud condensation nuclei (CCN) and ice nuclei (IN) concentration profiles. Finally, the product can be used to study dust dynamics during transportation, since it is capable of revealing even fine dynamical features such as the particle uplifting and deposition on European mountainous ridges such as the Alps and Carpathian Mountains.
  • Item
    Improving the LPJmL4-SPITFIRE vegetation–fire model for South America using satellite data
    (Katlenburg-Lindau : Copernicus, 2019) Drüke, Markus; Forkel, Matthias; von Bloh, Werner; Sakschewski, Boris; Cardoso, Manoel; Bustamante, Mercedes; Kurths, Jürgen; Thonicke, Kirsten
    Vegetation fires influence global vegetation distribution, ecosystem functioning, and global carbon cycling. Specifically in South America, changes in fire occurrence together with land-use change accelerate ecosystem fragmentation and increase the vulnerability of tropical forests and savannas to climate change. Dynamic global vegetation models (DGVMs) are valuable tools to estimate the effects of fire on ecosystem functioning and carbon cycling under future climate changes. However, most fire-enabled DGVMs have problems in capturing the magnitude, spatial patterns, and temporal dynamics of burned area as observed by satellites. As fire is controlled by the interplay of weather conditions, vegetation properties, and human activities, fire modules in DGVMs can be improved in various aspects. In this study we focus on improving the controls of climate and hence fuel moisture content on fire danger in the LPJmL4-SPITFIRE DGVM in South America, especially for the Brazilian fire-prone biomes of Caatinga and Cerrado. We therefore test two alternative model formulations (standard Nesterov Index and a newly implemented water vapor pressure deficit) for climate effects on fire danger within a formal model–data integration setup where we estimate model parameters against satellite datasets of burned area (GFED4) and aboveground biomass of trees. Our results show that the optimized model improves the representation of spatial patterns and the seasonal to interannual dynamics of burned area especially in the Cerrado and Caatinga regions. In addition, the model improves the simulation of aboveground biomass and the spatial distribution of plant functional types (PFTs). We obtained the best results by using the water vapor pressure deficit (VPD) for the calculation of fire danger. The VPD includes, in comparison to the Nesterov Index, a representation of the air humidity and the vegetation density. This work shows the successful application of a systematic model–data integration setup, as well as the integration of a new fire danger formulation, in order to optimize a process-based fire-enabled DGVM. It further highlights the potential of this approach to achieve a new level of accuracy in comprehensive global fire modeling and prediction.