Biogenic volatile organic compound (BVOC) emissions from European vegetation are a major precursor of tropospheric ozone and remain a key uncertainty in regional air-quality modelling. We present two high-resolution (0.1° × 0.1°) European BVOC emission datasets developed within the EU SEEDS project aimed at supporting scientific development within Copernicus Atmospheric Monitoring Service (CAMS). The datasets include BVOC species consistent with the RACM chemical mechanism and are generated by coupling the SURFEX land surface model with the MEGAN3.0 emission model.Emissions based on two land surface model simulations were analysed: (i) an open-loop SURFEX simulation available for 2018–2022, and (ii) a data-assimilation simulation in which satellite leaf area index (LAI) observations are assimilated, available for 2018–2020. In both cases, SURFEX is configured to allow vegetation phenological responses to meteorological variability, enabling a realistic representation of phenology. Evaluation against independent datasets shows that both simulations capture temporal variability in LAI and root-zone soil moisture, with improved skill in the analysis configuration.Given its importance for atmospheric chemistry, we focus on isoprene emissions. Interannual and seasonal variability in isoprene emissions is shown to be primarily driven by LAI variability, with specific events (e.g. summer 2019) linked to drought-induced vegetation stress simulated by SURFEX. Daily variability in isoprene emissions is evaluated against in-situ online isoprene concentration measurements at eight western European sites, revealing moderate to strong correlations across most site-year combinations. Comparisons with other bottom-up European isoprene inventories show that SURFEX-MEGAN3.0 emissions lie between the lower CAMS-GLOB-BIOv3.1 and higher MEGAN-MACC estimates, with differences in seasonality attributable largely to the underlying LAI datasets.These results highlight the important role of vegetation phenology, particularly LAI variability, in controlling BVOC emissions on monthly to interannual timescales, and demonstrate the added value of an Earth-system approach for BVOC emission modelling in support of air-quality assessments.ReferencesHamer, . D., Markelj, M., Rojas-Munoz, O., Bonan, B., Calvet, J.-C., Marécal, V., Guenther, A., Trimmel, H., Vallejo, I., Eckhardt, S., Sousa Santos, G., Sindelarova, K., Simpson, D., Schmidbauer, N., and Tarrasón, L.: Two Biogenic Volatile Organic Compound Emission Datasets over Europe Based on Land Surface Modelling and Satellite Data Assimilation, Earth Syst. Sci. Data Discuss. [preprint], https://doi.org/10.5194/essd-2025-442, in review, 2025.

Abstract Background Cloud water harbors diverse microbial communities despite its extreme oligotrophic conditions. However, the ecological and evolutionary dynamics of viruses in these transient atmospheric habitats remain poorly understood. Clouds have traditionally been regarded primarily as passive carriers of microorganisms rather than as active ecological environments supporting microbial interactions. In this study, cloud water was sampled at Mount Verde, Cape Verde Islands (744 m a.s.l.). We performed metagenomic analyses of iron-flocculated cloud water alongside genome analyses of a bacterial isolate and metagenome-assembled genomes using established bioinformatic approaches. Viral diversity, virus-host interactions, metabolic functions, genetic adaptations, and viral population dynamics across cloud events were investigated. In addition, UV-B resistance experiments were conducted for a novel cloud-water isolate. Results We isolated 24 cloud water bacteria, including four novel species lineages, and recovered 62 high-quality metagenome-assembled genomes, including 10 novel species lineages. We identified 458 viral operational taxonomic units and 237 virus-host linkages across diverse prokaryotic hosts, revealing active viral predation across diverse bacterial taxa. In addition, CRISPR spacer matches from isolates of novel bacterial lineages such as Deinococcus nubigenus MPC36 were found. Viruses carried genes involved in host adaptation to environmental stressors, including cold-shock response, UV radiation resistance, and osmotic stress. In addition, viral populations exhibited SNP-level microdiversity and shifts in single-nucleotide variant composition across temporally proximate cloud events, indicating rapid population turnover. Experimental characterization of the cloud isolate Curtobacterium nubigenum MPC39 further revealed pronounced resistance to UV-B radiation and the presence of an inducible prophage, Curtobacterium phage vB_CnuS_Cirrus1 assigned to the new viral family Nebulaviridae , which could be validated in transmission electron microscopy. Reconstructed genomes from cloud-associated bacteria encoded carbon monoxide dehydrogenase genes and UV resistance genes, suggesting trace gas metabolism and enhanced UV protection as survival strategies in oligotrophic cloud droplets. In silico replication rates estimated using iRep were consistent with active bacterial replication at the time of sampling. Conclusions Together, these findings demonstrate that clouds are not merely passive carriers of microorganisms, but dynamic atmospheric ecosystems in which virus-host interactions shape microbial diversity and contribute to microbial turnover, atmospheric dispersal, and cloud-water biogeochemistry.

Abstract. Closed-to-open cell mixed-phase cloud transitions within marine cold air outbreaks subjected to strong turbulent surface fluxes remain poorly understood despite their importance to high-latitude climate. The Cold-Air outbreak Experiment in the Sub-Arctic Region (CAESAR) research aircraft sampled closed-cells with cloud condensation nuclei concentrations surpassing 680 cm-3, decreasing to 90 cm-3 across a transition to open-cells. The aerosol likely originated from Siberian industrial emissions. With fetch, liquid water paths (LWPs) increase from 120 g m-2 to 270 g m-2 and cloud-top effective diameters increase from 10 μm to 16 μm, coincident with more riming. Ice particle number concentrations (Ni) are generally 2 L-1 or less, but exceed ice nucleating particle number concentrations by 100x. As the cloud-top inversion weakens and the boundary layer deepens further, ice precipitation co-exists with lidar-observed surface cold pools, modulated by entrainment events, juxtaposed with surface-based plumes of warm moist air. Open-cells contain isolated LWP peaks surpassing 500 g m-2 collocated with strong updrafts, adjacent to glaciated cloud. Ni surpasses 10 L-1 at cloud temperatures < -15 °C. Precipitation shafts contain abundant large graupel (> 5 mm diameter) with liquid-equivalent precipitation intensities reaching 3 mm hr-1 developing cold pools with virtual potential temperature depressions reaching 1.3 K. Nonetheless, buoyancy fluxes of 200-250 W m-2 prevent sub-cloud decoupling. The updrafts supporting liquid water production occur at the upwind edge of the cold pools. This case expands the observations needed to better understand mixed-phase Arctic cloud processes.

Biogenic volatile organic compound (BVOC) emissions from vegetation represent a major source of volatile compounds globally and play an important role as precursors for tropospheric ozone. Understanding their emissions is therefore crucial for quantifying the impact of ozone on air quality. We present two datasets of biogenic volatile organic compound emissions that cover the European modelling domain of the Copernicus Atmospheric Monitoring Service at a resolution of 0.1 degrees & times; 0.1 degrees to support the study of European scale air quality. The compounds included in the dataset follow the VOCs included in the regional atmospheric chemistry model mechanism (RACM). The datasets were produced within the framework of the EU's SEEDS project. We produced each dataset by coupling modelling output variables from the SURFEX land surface model with the MEGAN3.0 BVOC emission model. In one instance, the SURFEX model was run in free-running mode, which we term the open-loop (OL) and in the other case we assimilated satellite observations of leaf area index (LAI), which we term the analysis. The OL and analysis land surface model outputs form the basis for each emission dataset that are called SURFEX-MEGAN3.0 OL (10.7910/DVN/LAUVTU, Hamer et al., 2025a) and SURFEX-MEGAN3.0 analysis (10.7910/DVN/69G1FX, Hamer et al., 2025b), respectively. The OL dataset is available over a five-year period from 2018-2022 and the analysis dataset is available over the three-year period 2018-2020. SURFEX was run for both the OL and analysis simulations in a configuration that allowed simulated vegetation to respond to variations in meteorology over time to more realistically track vegetation phenology. Evaluation of the land surface model output LAI and root-zone soil moisture (RZSM) showed that the OL and analysis simulations had good skill at tracking temporal changes in both variables, with the analysis performing better in each instance. We perform a variety of evaluations on the isoprene emissions specifically given the importance of this compound for atmospheric chemistry. We evaluated the temporal variability of isoprene emissions in both datasets and found that the majority of the interannual and monthly variability was linked to variability in LAI that in specific cases, like the summer of 2019, could be linked to drought impacts on vegetation growth simulated by SURFEX. We evaluated the daily temporal variability of the OL and analysis isoprene emission datasets against in-situ online observations of isoprene concentrations at 8 sites in western Europe and found moderate to strong correlation between the emissions and observations in almost all location-year pairings. We also evaluated the OL and analysis emission datasets against other published bottom-up isoprene emission datasets over the same European domain used in this study. We found that the SURFEX-MEGAN3.0 OL and analysis isoprene emission datasets lie between the minimum (CAMS-GLOB-BIOv3.1) and maximum (MEGAN-MACC) published emission datasets based on bottom-up approaches. Furthermore, we were able to attribute differences in seasonality between SURFEX-MEGAN3.0 and other emission inventories to differences in the temporal variability of the underlying LAI dataset used to compile them. Overall, our findings show the importance of variability in LAI in controlling isoprene emissions on monthly to annual timescales. Combining this with the demonstrated skill of the emissions in evaluation with independent data, this points towards the value of an Earth-system approach to BVOC emission modelling.

Methane is a powerful greenhouse gas with a shorter lifetime than carbon dioxide (CO2), making it an important target for near-term climate action. The Global Methane Pledge (GMP) aims to cut anthropogenic methane emissions by 30% from 2020 levels by 2030. Using an Earth system model with interactive CH4 sources and sinks, we assess the Pledge's impact through 2050. Results show that current GMP commitments deliver only a 10% cut by 2030-well below the target. Only the maximum technically feasible reduction (MTFR) pathway can achieve the 30% goal. By 2050, current GMP commitments lowers methane concentrations by 3% relative to 2025, while MTFR achieves 8%. Both pathways slow warming slightly, avoiding about 0.1 degrees C of global temperature rise, with the Arctic seeing the greatest benefits (up to 2 degrees C less warming). Without wider participation, the GMP with current signatories will fall short of its targets and Paris Agreement goals.

Abstract Plastic pollution is an ever-growing concern around the globe, with current research portraying the presence of contaminants even in previously pristine regions. Especially small particles, micro- and nanoplastics, can be dispersed over long stretches by atmospheric transport, reaching even the most remote areas of the planet. This study aims to portray current levels of nanoplastic contamination (< 1 μm) by analyzing snow samples from high-altitude glaciers in the European Alps (2300–3800 m) via TD-PTR-MS. We were able to detect six common plastic pollutants, with PE and PP contributing over 60% of the total mass concentration measured across all sites. On average, nanoplastic concentrations of 85 ng mL − 1 (range: 3.6 ± 12 ng mL − 1 – 470 ± 76 ng mL − 1 ) were detected at sampled glacier surfaces. Access to these remote regions was enabled through a citizen-science initiative, involving trained mountaineers for sample collection. The results were applied to atmospheric modelling, highlighting possible point sources of contamination. While agriculture and local plastic factories were revealed as potential sources of nanoplastics, spatial and temporal limitations complicated comparison inside the sample cluster.

Abstract. Long-term observations of black carbon (BC) are essential for assessing the effectiveness of air-quality and climate policies and for understanding the potential feedback of climate change. We present an 18-year dataset (2007–2024) of equivalent black carbon (eBC) mass concentration derived from the aerosol absorption coefficient measured at the Monte Cimone WMO/GAW Global Station (CMN, 2165 m a.s.l.), a key high-altitude observatory in the Mediterranean climate-change hotspot. The dataset is accompanied by detailed description of: i) the infrastructural evolution, including instrumentation and sampling system; ii) the definition of quality assurance protocol and data treatment; iii) the definition of the uncertainties. This documentation provides a comprehensive description and quantification of the data reliability over 18 years of measurements, necessary for the scientific valorisation of the dataset. Applications of the dataset demonstrate: i) a statistically significant decline in eBC concentrations over the last two decades ii) a strong linkage between eBC variability and boundary-layer dynamics when combined with ERA5 reanalysis; iii) a variable agreement with FLEXPART simulations across temporal scales, supporting its use for model evaluation. The dataset, openly available through the ITINERIS HUB, provides one of the longest continuous eBC records in the Mediterranean troposphere and represents a valuable resource for climate studies, trend assessments, and multi-station integration within ACTRIS and GAW.

The implementation of air pollution mitigation strategies requires not only high-quality continuous measurements of pollutants but also proper definitions of ways to differentiate between transported and locally produced contributions, as only the latter can be effectively reduced by authorities. To address this issue, we propose a new approach for partitioning monitored black carbon (BC) concentrations into city-produced (urban) and transported fractions using a combination of measured and modeled data. Two simultaneous measurement campaigns (warm season 2022 and cold season 2022/23) were conducted in two urban environments: Vilnius (Lithuania) and Warsaw (Poland). In the cold season in Warsaw, BC mass concentration was 90% higher than in the warm season, while in Vilnius, an increase of 44% was observed, as compared to the warm season. Aerosol optical properties showed more complex aerosol mixtures of dust, BC and brown carbon (BrC) during the cold season, forming larger particles. Single scattering albedo (SSA) anti-correlated with BCFF, proving that fossil fuel (FF) combustion contributes to the warming effect in both cities. A positive correlation between the population density of the emission areas of transported BC and the BC mass concentrations in Vilnius and Warsaw was found. The impact of transported BC on the local BC levels in the cities was of 53 +/- 8% and 36 +/- 28% in the cold season and of 23 +/- 7% and 13 +/- 21% in the warm season for Warsaw and Vilnius, respectively. Thus, the approach of BC partitioning showed that in the cold season, the two cities suffered from worse air quality, in part due to more transported BC.

Plastic contamination has become a global concern, with evidence even in remote regions like Antarctica. While macro- and microplastics have been documented in Antarctic marine ecosystems, their presence in soils - particularly submicro- and nanoplastics - remains poorly studied. This study analyses soil samples from the McMurdo Dry Valleys collected on January 8th to 28th, 2023, and reports the first detection of nanoplastics - including polypropylene, polyethylene, polyethylene terephthalate, polystyrene, polyvinyl chloride, and tyre wear particles - using thermal desorption proton transfer reaction mass spectrometry. These plastics were detected at multiple topsoil sampling sites (n = 13), with concentrations reaching up to 295 ng g⁻¹ with nanoplastics detected above polymer-specific method detection limits at 54% of sites (median: 26.6 ng g⁻¹). They were also detected at lower concentrations in deeper soil layers (> 20 cm; n = 4), where nanoplastics were present at 50% of the sampled sites (median: 1.95 ng g⁻¹). Lagrangian particle dispersion model FLEXPART suggested seasonal deposition patterns, with inputs from both local sources and long-range atmospheric transport. This evidence shows that soils in one of Earth's most pristine environments are not exempt from plastic contamination, with the reported concentrations providing a crucial baseline for global pollution assessments. These findings also highlight the urgent need to study plastic fate, transport, and ecological impacts in polar regions.

Black carbon (BC) is a product of incomplete combustion, is climate relevant and has negative impacts on human health. It absorbs radiation as an aerosol in the atmosphere, but also changes the albedo of snow covered surfaces and leads to earlier melting. The origins are either anthropogenic or natural, with different annual cycles. While anthropogenic sources like domestic burning peak in the winter, natural sources like wild fires and agricultural burning peak during spring and summer. Due to its short lifetime are the global atmospheric concentrations highly variable.We use ground based observational data from different global networks, and the atmospheric transport model FLEXPART driven by ERA5 meteorological analysis, emission inventories and the inversion framework FLEXINVERT. By minimizing the mismatch modelled and observed BC concentration we improve the emission inventories for natural emissions (GFAS) and anthropogenic emissions (LRTAP). For every year up to 50 stations are used and each observation is matched with a 50 day FLEXPART backward calculation.We discuss the distribution and sources of global BC aerosols over the period 2015 to 2022 and compare existing emission inventories with the improved constraints of the global BC emissions derived with the FLEXINVERT inversion framework.

Summary Volcanic ash represents an extreme and dynamic habitat, yet it hosts diverse microbial communities with largely unexplored viral diversity. This study investigated bacterial and viral populations in volcanic ash from Mount Etna (Italy) collected during the eruption, focusing on microbial novelty, activity, and virus-host interactions. Taxonomic profiling revealed that Pseudomonas and Telluria were the dominant bacterial genera, both frequently detected in airborne environments. In contrast, enrichment cultures with volcanic ash were dominated by spore-forming members of the phylum Bacillota, highlighting their resilience under harsh conditions. Metagenomic analysis recovered 19 high-quality metagenome-assembled genomes, including four previously undescribed bacterial species. Replication rate estimates showed that certain taxa were metabolically active, particularly at one sampling site. The presence of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) systems with spacers matching viral sequences suggested viral predation pressure on volcanic ash. A total of 1139 viral operational taxonomic units (vOTUs) were identified, of only around half (660 vOTUs) showed similarities to known phages, underscoring the presence of novel viruses. Shared vOTUs across sites revealed the presence of both a core virome and site-specific viral populations. Virus-host predictions indicated frequent interactions with hosts from multiple Gammaproteobacterial genera. Additionally, a 336 kb jumbo phage genome exhibited extensive metabolic capabilities and genetic autonomy. Experimental work identified a unique lytic Bacillus -infecting phage (“Phoenix”) with limited propagation capacity. Furthermore, prophage induction experiments revealed active, morphologically diverse temperate phages across multiple bacterial host strains. Overall, these findings highlight volcanic ash as a reservoir of microbial and viral diversity, shaped by environmental extremes and dynamic ecological interactions. Highlights Metagenomic and cultivation experiments were used to study viruses on volcanic ash Novelty of viral and bacterial species was detected Viral-bacterial interactions in metagenomes from volcanic ash were detected Cultivatable bacteria were mainly spore-forming Bacilli species and harbored inducible prophages

In the summer of 2019, the Arctic region registered exceptionally high aerosol optical depth (AOD) values over Svalbard, linked to intense biomass burning (BB) and volcanic activity across the Northern Hemisphere. This study presents a comprehensive, multi-instrumental analysis of the aerosol conditions in and around Ny-& Aring;lesund (Spitsbergen, Norway), combining data from ground-based sun-photometry, in-situ observations, active remote sensing (ground-based and on satellite), and atmospheric dispersion modelling (FLEXPART). Despite high AOD was observed during all the period, three different aerosol events are identified in the atmospheric column (6-10 July, 25-28 July, and 6-17 August). In contrast, in-situ surface stations only recorded significant aerosol load during 5-9 July, 30 August, and 12 September, suggesting that most of the aerosol particles remained above the boundary layer. Lidar and photometric observations revealed the presence of spherical, weakly absorbing Accumulation-mode particles (with effective radii between 0.1 and 0.2 mu m) in both the troposphere and stratosphere, with persistent layers extending above 10 km. Simulations carried out with FLEXPART correlate well with the measurements, attributing the observed aerosol events to multiple sources, including Siberian and North American wildfires, the Raikoke (Russia) volcanic eruption, and anthropogenic pollution. While the simulations show a contribution from volcanic aerosols, the contribution from biomass-burning aerosols in the upper troposphere and lower stratosphere were likely more significant under the atmospheric conditions of summer 2019. Overall, the aerosol radiative impact during this long-lasting event was substantial, with a mean reduction in direct solar radiation of approximately -74 Wm-2 during July and August. This work shows how the use of dispersion modelling together with multiple observation sources allows to achieve a more complete description of the atmospheric aerosol events and contributes to a better understanding of the overall picture.

Abstract. Biomass burning (BB) aerosols emitted over Southern Africa (SAF) and South America (SA) represent a major seasonal perturbation to the usually pristine atmosphere of the South-West Indian Ocean (SWIO) but remain understudied in this remote region. Following a multi-instrumental approach, we characterize BB plumes reaching Reunion Island (21° S, 55° E) during September 2017, combining ground-based measurements (sun-photometer, lidars, Fourier Transform Infrared spectrometer), spaceborne observations, CAMS EAC4 reanalysis, and the Lagrangian transport model FLEXPART. Aerosol optical depth at 550 nm over Reunion reached unusually high values (0.16–0.42) during the second half of September, with organic matter contributing up to 60 %. MODIS imagery revealed two large-scale smoke plumes originating from SAF and SA transported toward the SWIO, and CALIOP profiles showed smoke layers extending from 4 to 9 km of altitude above Madagascar. On September 19th, a single layer was identified over Reunion between 2.8 and 4.7 km with an Ångström exponent (Å₃₅₅/₅₃₂) of 1.32±0.23, consistent with moderately aged BB particles essentially originating from SAF. On September 25th, two vertically decoupled layers were identified: a lower layer (3.3–5 km, Å=1.45±0.12) associated with mixed aged aerosols of SAF and SA origin, and an upper, drier layer (5–9 km, Å=1.60±0.06) of potentially fresher SAF smoke, consistent with rapid convective uplift into the mid-troposphere. This study offers new insights concerning the dynamical processes that govern aerosol variability over Reunion Island and highlights its value as a strategic long-term observational site in the SWIO.

Abstract. Carbonaceous aerosol measurements are scarce in the Western Balkans. We present year-round observations (2023 to 2025) from the urban-background site Ada Marina, Belgrade (44.790166° N, 20.4169° E; 71 m a.s.l.), combining three source apportionment (SA) approaches (organic tracers, the aethalometer model, and a novel application of non-negative matrix factorization, NNMF, to aethalometer data) with FLEXPART modelling and gridded emissions to identify source regions and transport pathways. The annual mean equivalent black carbon (eBC) concentration was 1.42 ± 1.40 µg m⁻³, indicating a moderate eBC burden, while PM₁₀ (27.7 µg m⁻³) remained below current EU and Serbian limits but above WHO guidelines and upcoming EU limits. All SA methods consistently show strong winter to summer seasonality driven by residential wood combustion (RWC). The aethalometer model attributes 64 % of eBC to solid fuel (eBCSF) in winter versus 18 % in summer, NNMF yields 51 % versus 6 %, and offline analysis of tracers show elemental carbon from biomass burning (ECBB) contributions of 48 % versus 3 %. Organic carbon (OC) shows similar seasonality, with biomass burning contributing 59 % in winter and 4 % in summer. A substantial wintertime secondary OC component linked to RWC is evident, with secondary OC (SOC) strongly correlated with ECBB and eBCSF (R² = 0.78 to 0.86), alongside enhanced brown carbon absorption (29 % annually; 44 % in winter). FLEXPART indicates pollution episodes are dominated by regional transport within Serbia and the Pannonian Basin, while correlations between modelled residential emissions and observations demonstrate strong model skill.

Coastal aerosols are formed through the complex mixing between marine air masses and continental emissions, which originate from both natural and anthropogenic sources. The properties of coastal aerosols are decisive for their interaction with sunlight and their influence on clouds, as well as the potential health implications for the population in these areas. In this study, the aerosol properties and sources at Aarhus Bay, Denmark, were investigated by combining in situ aerosol light scattering and absorption with size distribution measurements and footprint analysis by FLEXPART. Our analysis demonstrates a considerable contribution of anthropogenic aerosols from both fossil fuel combustion and biomass burning, as well as periods with highly scattering aerosols. Furthermore, good agreement was found between in situ and modelled black-carbon data. Combining in situ measurements and FLEXPART analysis further evidenced a major impact of local emissions, as well as a few long-range transport intrusions.

This study presents a comprehensive assessment of unique parallel measurements of surface airborne and deposited microplastics (AMPs) across urban and remote sites in Norway, employing pyrolysis-GC/MS for polymer-specific analysis. MPs were detected in nearly all samples, with significantly higher concentrations and fluxes observed in urban areas like Oslo, where tire wear particles (TWP) dominated (>90% of AMP mass). Seasonal peaks in TWP coincided with the transition to winter tires, while remote sites showed consistent but lower AMP levels, indicating long-range transport (LRT) from European source regions. Parallel measurements of suspended and deposited AMPs revealed consistent polymer signatures, highlighting common sources and transport pathways. Although urban TWP contributions to PM2.5 were generally low, episodic events reached up to 30%, raising concerns about human exposure. The dual dataset enabled a robust cross-validation of atmospheric loading estimates and facilitated integration into advanced transport models for remote sites. Our findings confirm AMPs as significant components of urban air pollution and subsequent carriers of chemical and biological contaminants to remote regions, emphasizing the need for targeted monitoring and mitigation strategies.

Abstract. Biogenic volatile organic compound (BVOC) emissions from vegetation represent a major source of volatile compounds globally and play an important role as precursors for tropospheric ozone. Understanding their emissions is therefore crucial for quantifying the impact of ozone on air quality. We present two datasets of biogenic volatile organic compound emissions that cover the European modelling domain of the Copernicus Atmospheric Monitoring Service at a resolution of 0.1° × 0.1° to support the study of European scale air quality. The compounds included in the dataset follow the VOCs included in the regional atmospheric chemistry model mechanism (RACM). The datasets were produced within the framework of the EU's SEEDS project. We produced each dataset by coupling modelling output variables from the SURFEX land surface model with the MEGAN3.0 BVOC emission model. In one instance, the SURFEX model was run in free-running mode, which we term the open-loop (OL) and in the other case we assimilated satellite observations of leaf area index (LAI), which we term the analysis. The OL and analysis land surface model outputs form the basis for each emission dataset that are called SURFEX-MEGAN3.0 OL and SURFEX-MEGAN3.0 analysis, respectively. The OL dataset is available over a five-year period from 2018–2022 and the analysis dataset is available over the three-year period 2018–2020. SURFEX was run for both the OL and analysis simulations in a configuration that allowed simulated vegetation to respond to variations in meteorology over time to more realistically track vegetation phenology. Evaluation of the land surface model output LAI and root-zone soil moisture (RZSM) showed that the OL and analysis simulations had good skill at tracking temporal changes in both variables, with the analysis performing better in each instance. We perform a variety of evaluations on the isoprene emissions specifically given the importance of this compound for atmospheric chemistry. We evaluated the temporal variability of isoprene emissions in both datasets and found that the majority of the interannual and monthly variability was linked to variability in LAI that in specific cases, like the summer of 2019, could be linked to drought impacts on vegetation growth simulated by SURFEX. We evaluated the daily temporal variability of the OL and analysis isoprene emission datasets against in-situ online observations of isoprene concentrations at 8 sites in western Europe and found moderate to strong correlation between the emissions and observations in almost all location-year pairings. We also evaluated the OL and analysis emission datasets against other published bottom-up isoprene emission datasets over the same European domain used in this study. We found that the SURFEX-MEGAN3.0 OL and analysis isoprene emission datasets lie between the minimum (CAMS-GLOB-BIOv3.1) and maximum (MEGAN-MACC) published emission datasets based on bottom-up approaches. Furthermore, we were able to attribute differences in seasonality between SURFEX-MEGAN3.0 and other emission inventories to differences in the temporal variability of the underlying LAI dataset used to compile them. Overall, our findings show the importance of variability in LAI in controlling isoprene emissions on monthly to annual timescales. Combining this with the demonstrated skill of the emissions in evaluation with independent data, this points towards the value of an Earth-system approach to BVOC emission modelling.