You are here: Home » Institute » Image of the Month

Image of the Month

April 2023:

Figure from Karmouche et al. (2022):
Figure from Karmouche et al. (2022): (a) Whisker plot showing the distribution of F1-scores across the CMIP6 LEs for the causal analysis for: the complete period (light blue boxes), the In-Phase regime (dark blue boxes) and the Out-of-Phase regime (green boxes). White scatter points denote the mean LE F1-scores. (b) Reference causal network estimated from reanalysis during the Out-of-Phase regime (left, with low-pass AMV and PDV time series below) compared to networks and time series from three CMIP6 simulations (right, with simulated low-pass AMV and PDV time series below each network) with the best network similarity i.e. highest F1-score.

Regime-oriented causal analysis of Atlantic-Pacific teleconnections in CMIP6

Modes of natural climate variability from interannual to multidecadal timescales have large effects on regional and global climate with important socio-economic impacts (Eyring et al., 2021). Despite their importance, systematic evaluation of climate models and their simulation of internal variability remains a challenging task (Eyring et al., 2019). Here, causal discovery provides a way to estimate dynamical climate dependencies from timeseries data. In Karmouche et al. (2022), we perform a regime-oriented causal model evaluation to assess the ability of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6, Eyring et al., 2016) to represent the observed changing interactions between the Pacific Decadal Variability (PDV) and the Atlantic Multidecadal Variability (AMV) and their extra-tropical teleconnections to the Pacific North- and Pacific South- American Patterns (PNA and PSA1). The indices for these modes of climate variability are calculated by the Climate Variability Diagnostic Package for Large Ensembles (CVDP-LE, Phillips et al., 2020).

Our figure of the month shows a comparison of causal networks for reanalysis datasets (ERSSTv5 for PDV and AMV, ERA20C_ERA5 for PNA and PSA1) and CMIP6 Large Ensemble (LE) historical simulations during the period 1900-2014. The interactions between the aforementioned modes of climate variability are hypothesized to follow different regimes depending on the physical state of the Atlantic or Pacific (i.e. AMV or PDV phase). To identify phase-dependent sets of connections (causal networks), we perform the analysis on shorter periods (regimes) that represent different phases based on the low-pass filtered AMV and PDV indices, thus identifying the regime-oriented analysis. Here, we refer to the time periods where AMV and PDV are on the same phase (both positive or both negative) by “In-Phase regime”, whereas for the periods where the two indices are on opposite phases (one negative while the other is positive) we use the term “Out-of-Phase regime”.

We run the PCMCI+ (Runge et al., 2019; Runge, 2020) causal discovery method, available as part of the Tigramite package, to estimate causal networks from CMIP6 Large Ensemble models and compare results to those from reanalysis data during different specified periods. Blue links in the resulting causal graphs translate an opposite-sign response while red links a same-sign response. The labels on the links, when shown, denote the time-delay of the response in years. The figure shows (a) a whisker plot for the distribution of F1-scores, a measure of similarity between causal networks following methods from Nowack et al. (2020), across the CMIP6 LEs during the complete period, In-Phase regime, and Out-of-Phase regime with light blue, dark blue, and green colored boxes, respectively. The white scatter points indicate that on average, CESM2, CanESM5, MIROC6, and MPI-ESM1-2-LR LEs exhibit better network similarity with observations during the Out-of-Phase regime. CanESM5 and MIROC6 LEs have the highest scores (0.92) during this regime, as indicated by the location markers on the whisker plot. In (b) we compare the causal graphs (three panels on the right) of these top-scoring realizations and their low-pass filtered AMV and PDV time series to those generated from reference reanalysis datasets (left). The networks in (b) agree on the 1-year lagged AMV→PDV link, illustrating an opposite-sign response. The positive contemporaneous PDV—PNA link is directed differently in CanESM5 r11i1p2f1 compared to reanalysis, and remains unoriented in CanESM5 r17i1p2f1 and MIROC6 r20i1p1f1. The Out-of-Phase graphs of these realizations also demonstrate an agreement on a same-sign contemporaneous AMV—PNA connection, which is weaker (lower cross-MCI value) than the PDV—PNA link. Overall, the regime-oriented causal model evaluation followed in Karmouche et al. (2022) has the potential of a powerful methodology that can be applied in a number of environment-related topics, offering tremendous insight to improve the understanding of the complex earth system and the state-of-the-art of climate modeling.


Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geoscientific Model Development, 9, 1937–1958, 2016.

Eyring, V., Cox, P. M., Flato, G. M., Gleckler, P. J., Abramowitz, G., Caldwell, P., Collins, W. D., Gier, B. K., Hall, A. D., Hoffman, F. M., et al.: Taking climate model evaluation to the next level, Nature Climate Change, 9, 102–110, 2019

Eyring, V., Gillett, N., Achutarao, K., Barimalala, R., Barreiro Parrillo, M., Bellouin, N., Cassou, C., Durack, P., Kosaka, Y., McGregor, S., et al.: Human Influence on the Climate System: Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, IPCC Sixth Assessment Report, 2021

Karmouche, S., Galytska, E., Runge, J., Meehl, G. A., Phillips, A. S., Weigel, K., and Eyring, V.: Regime-oriented causal model evaluation of Atlantic-Pacific teleconnections in CMIP6, EGUsphere [preprint],, 2022.

Nowack, P., Runge, J., Eyring, V., and Haigh, J. D.: Causal networks for climate model evaluation and constrained projections, Nature communications, 11, 1–11, 2020.

Phillips, A. S., Deser, C., Fasullo, J., Schneider, D., and Simpson, I.: Assessing Climate Variability and Change in Model Large Ensembles: A User’s Guide to the “Climate Variability Diagnostics Package for Large Ensembles”, version 1, version, 1, 0, 2020

Runge, J.: Discovering contemporaneous and lagged causal relations in autocorrelated nonlinear time series datasets, in: Conference on Uncertainty in Artificial Intelligence, pp. 1388–1397, PMLR, 2020.

Runge, J., Nowack, P., Kretschmer, M., Flaxman, S., and Sejdinovic, D.: Detecting and quantifying causal associations in large nonlinear time series datasets, Science advances, 5, eaau4996, 2019

March 2023:

Figure: ROxCOMP22 campaign in the SAPHIR chamber, a) View of SAPHIR and the instruments participating in the campaign, b) IUP participants, c) inside the SAPHIR chamber, d) detail of the PeRCEAS inlet inside SAPHIR, e) PeRCEAS signal during the measurement of peroxy radicals. The red arrows indicate the position of the inlet in the chamber and the container for the PeRCEAS instrument.

Workshop on the ROxCOMP22 campaign

Radical chemistry is of great importance for understanding the atmospheric chemical processing and formation of secondary pollutants. The TROLAS (Tropospheric Radical Observations and Laser Absorption Spectroscopy) group at the Institute of Environmental Physics has wide experience in the investigation of the distribution of peroxy radicals and their precursors, observed within airborne and ground based field measurements in different environments.

Scientists from IUP-TROLAS successfully participated with the airborne PeRCEAS
(Peroxy Radical Chemical Enhancement and Absorption Spectrometer) instrument in the 1-month long international comparison study ROxCOMP22. ROxCOMP22 aimed at improving the quality of the peroxy radical detection by bringing together instruments using different techniques to measure organic peroxy radicals worldwide. The 16 chamber experiments made in August 2022 in the SAPHIR chamber at the Forschungszentrum Jülich provided a unique opportunity to identify potential interferences or systematic errors while investigating relevant radical reactions under controlled chemical conditions.

After a calibration phase in the home laboratories, the first workshop on the discussion of the obtained ROxCOMP22 data has taken place In February 2023. The data analysis is still on-going.

Contact: Lola Andrés Hernández

February 2023:

(click to enlarge)
(click to enlarge) Trends of equivalent water height (EWH) within the Indian subcontinent based on gravimetric data from the GRACE and GRACE-FO satellites (Boergens et al., 2020). The global map (top-left) shows linear EWH trends for the 20-year period of 2002-2021 over all continental regions except Greenland and Antarctica. The blue rectangle indicates the Ganges-Brahmaputra-Meghna region. The EWH time series within this region (bottom-left) is based on area weighed averages of the monthly data. The panel also includes EWH data after correction for river sediment transport as well as linear least-squares fits for both original (grey) and corrected (colored) EWH data. The regional map of EWH trends in the Indian subcontinent (top-right) shows the Ganges and the Brahmaputra catchments in blue and red, respectively. The time series of EWH within these catchments (bottom-right) iterate between the Ganges, the Brahmaputra as well as the mountain regions within these catchments. The shaded area in the local map indicates what region the individual time series are based on.

Gravity fields measured by geodetic satellites yield information on global mass variations and have proven crucial to monitor changes in water storage and fluxes. Thus, satellite gravimetry is a key component in the investigation of the massive groundwater depletion observed on the Indian subcontinent. In a current study, we look into the impact of river sediment transport on terrestrial EWH data measured by satellites of the Gravity recovery and Climate Experiment (GRACE) and its follow-on mission (GRACE-FO). While the mass loss by sediment transport in rivers is assumed to be below the detection limit of those satellites, it potentially impacts long term EWH trends derived from their data. The Indian subcontinent is drained by the Ganges and Brahmaputra rivers which constitute one of the the world's most sediment rich river systems. In this month's figure, we show the impact of sediment mass loss within the Ganges and Brahmaputra catchments on gravimetric estimates of EWH trends. Forthe whole Ganges-Brahmaputra-Meghna region, sediment transport accounts for roughly 4% of the gravity decrease currently attributed to groundwater depletion. Since most sediment is eroded from mountain regions, the impact of sediment transport on EWH trends in these regions is high (average of 14%).

Contact: Alexandra Klemme (

Boergens, Eva; Dobslaw, Henryk; Dill, Robert (2020): COST-G GravIS RL01 Continental Water Storage Anomalies. V. 0004. GFZ Data Services.

January 2023:

Seasonal climatology of tropospheric ozone column (TrOC) derived from OMPS aboard the Suomi-NPP satellite using limb-nadir matching (LNM). From left to right and top to bottom: 2012-2018 seasonal means in boreal winter, spring, summer, and autumn. SSA stands for South Atlantic Anomaly, which is a region where satellites are impacted by electronic disturbances due to an earth’s magnetic field anomaly. This anomaly lowers the quality of the derived ozone measurements (Orfanoz-Cheuquelaf et al., 2022).

The SUOMI-NPP platform carries three instruments that measure ozone in limb and nadir viewing geometries. The OMPS-Nadir Napper (NM) provides total column ozone from the ground to the upper atmosphere. From spectral data of the OMPS-Limb profiler (LP), vertical ozone profiles are derived. Using the ozone data from both instruments, the tropospheric ozone column (from the ground to the tropopause, located between 7 and 17 km altitude, dependent on latitude) is calculated by subtracting the stratospheric column from the total. The figure shows the global seasonal mean of tropospheric ozone from 2012 to 2018 as observed by OMPS. The seasonal maps show the variable outflow of tropospheric ozone e.g. from Africa (biomass-burning regions), North-East America, and Eastern Asia (both industrial regions), where ozone is produced from chemical precursors, which are mainly air pollutants.

In addition to the OMPS tropospheric ozone algorithm, we used our advanced algorithms to retrieve the total column and profiles of ozone from both OMPS instruments. Similar algorithms have been applied to SCIAMACHY satellite data (2002-2012). In the next step, we will combine LMN tropospheric ozone data from SCIAMACHY and OMPS to obtain a long-term data record spanning at least 16 years for investigating changes in decadal timescales.


Arosio, C., Rozanov, A., Malinina, E., Eichmann, K.-U., von Clarmann, T., and Burrows, J. P.: Retrieval of ozone profiles from OMPS limb scattering observations, Atmos. Meas. Tech., 11, 2135–2149,, 2018.

Ebojie, F., von Savigny, C., Ladstätter-Weißenmayer, A., Rozanov, A., Weber, M., Eichmann, K.-U., Bötel, S., Rahpoe, N., Bovensmann, H., and Burrows, J. P., Tropospheric column amount of ozone retrieved from SCIAMACHY limb–nadir-matching observations, Atmos. Meas. Tech., 7, 2073-2096,, 2014.

Ebojie, F., Burrows, J. P., Gebhardt, C., von Savigny, C., Rozanov, A., Weber, M., and Bovensmann, H.: Global and zonal tropospheric ozone variations from 2003–2011 as seen by SCIAMACHY, Atmos. Chem. Phys., 16, 417-436,, 2016.

Orfanoz-Cheuquelaf, A., Rozanov, A., Weber, M., Arosio, C., Ladstätter-Weißenmayer, A., and Burrows, J.P.: Total ozone column from Ozone Mapping and Profiler Suite Nadir Mapper (OMPS-NM) measurements using the broadband weighting function fitting approach (WFFA), Atmos. Meas. Tech., 14, 5771–5789,, 2021.

Orfanoz-Cheuquelaf, A, Arosio, C., Rozanov, A., Weber, M., Ladstätter-Weißenmayer, A., and Burrows, J.P., Tropospheric ozone column dataset from OMPS-LP/OMPS-NM limb-nadir matching, manuscript in preparation, 2022.

December 2022:

A priori [panel (a)] and optimized a posteriori [panel (b)] estimate for the secondary CO production in the atmosphere for September 2018. Panel (c) shows the absolute difference between the a posteriori and a priori emissions. The relative difference between a posteriori and a priori CO emissions using both TROPOMI observations and flask measurements is presented in panel (d). The black dots denote the location of the NOAA surface flask measurements.

Johann Rasmus Nüß1, Nikos Daskalakis1, Fabian Günther Piwowarczyk1, Angelos Gkouvousis2,3, Oliver Schneising1, Michael Buchwitz1, Maria Kanakidou1,2,3, Maarten C. Krol4,5, Mihalis Vrekoussis1,6,7

1 Institute of Environmental Physics (IUP-UB), University of Bremen, Bremen, Germany
2 Environmental Chemical Processes Laboratory (ECPL), Department of Chemistry, University of Crete, Heraklion, Greece
3 Center for the Study of Air Quality and Climate Change (C-STACC), Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas, Patras, Greece
4 Meteorology and Air Quality, Wageningen University and Research, Wageningen, the Netherlands
5 Institute for Marine and Atmospheric Research, Utrecht University, Utrecht, the Netherlands
6 Center of Marine Environmental Science (MARUM), University of Bremen, Germany
7 Climate and Atmosphere Research Center (CARE-C), The Cyprus Institute, Nicosia, Cyprus

Carbon monoxide (CO) in the atmosphere adversely affects air quality and climate, making knowledge about its sources crucial. However, current global bottom-up emission estimates retain significant uncertainties. We attempt to reduce these uncertainties by optimizing emission estimates for the second half of the year 2018 on a global scale with a focus on the Northern Hemisphere through the top-down approach of inverse modeling. Specifically, we introduce observations from the TROPOspheric Monitoring Instrument (TROPOMI) into the TM5-4DVAR model. The emissions are further constrained using NOAA surface flask measurements, with their locations denoted as black dots in the figures below.

Figure (a) shows the vertically integrated bottom-up estimate for the secondary, i.e. chemical, production of CO in the atmosphere for September 2018 as taken from the full-chemistry model TM5-MP. This production field is used as part of the initial emissions, the a priori, for our inversions. The resulting optimized field, the a posteriori, is shown in Figure (b). These emission fields lead to the smallest difference between the simulated mixing ratios and those observed by both TROPOMI and the NOAA flasks. Due to its strong dependence on photochemistry, secondary CO production rates vary by multiple orders of magnitude globally, making the absolute increments between a posteriori and a priori emissions hard to interpret (panel (c)). Therefore, we show their relative difference in Figure (d), i.e. the factor by which the emissions changed at each location.

Compared to the bottom-up estimates, we observe strong broad-scale emission reductions (up to 75 %) in China and India. In part, these reductions can be attributed to policy changes in China. However, the OH climatology we use to simulate chemical loss appears to be underestimated in that region, which skews the inversion towards lower emissions as well. Conversely, we find strong local emission increments over Europe and the Sahara. These are likely artifacts caused by the model’s limited capabilities to capture the surface flask measurements at these specific stations and cannot be confirmed by satellite observations.

Notably, the mixing ratios calculated using optimized emissions from the inversion driven solely by satellite observations agree very well with the flask measurements (comparison not shown). While this holds only south of 55° N, due to model limitations, it could potentially allow for near real-time inversions purely based on satellite observations. These would then be validated against and adjusted to the surface flasks as soon as they are available.


Nüß, J. R., Daskalakis, N., Piwowarczyk, F. G., Gkouvousis, A., Schneising, O., Buchwitz, M., Kanakidou, M., Krol, M. C., Vrekoussis, M: Efficacy of high-resolution satellite observations in inverse modeling of carbon monoxide emissions, submitted to Journal of Advances in Modeling Earth Systems (JAMES), 2022

Contact: Rasmus Nüß (

November 2022:

Fig. 1
Fig. 1
Fig. 2
Fig. 2
Fig. 3
Fig. 3

Methane anomalies in Greenland

Methane (CH4) represents the second most important greenhouse gas after carbon-dioxide (CO2) and is roughly 28-times more potent than CO2. It is released through human activities (agriculture, waste, fossil fuels,...) and natural processes, for example bio-chemical processes in wetlands or in thawing permafrost.

That is why we are investigating the atmospheric methane concentration (XCH4) within the framework of the junior-research group „Greenhouse gases in the Arctic“ using satellite data. Over Greenland strong methane anomalies could be observed (Fig. 1a), i.e. regions with very high or low concentrations. However these anomalies are not related to emissions from the permafrost but to inaccuracies in the digital elevation model (DEM) used in the calculations (Fig. 2), as our research has shown.
In addition to Greenland further regions show inaccuracies in older DEMs (Fig. 3). Following this finding a more recent DEM was implemented for the WFMD methane product which is developed at the IUP. The anomalies vanish in the new product both for Greenland (Fig 1b) and the other regions which improved the general quality of the product.


Hachmeister, J., Schneising, O., Buchwitz, M., Lorente, A., Borsdorff, T., Burrows, J. P., Notholt, J., and Buschmann, M.: On the influence of underlying elevation data on Sentinel-5 Precursor TROPOMI satellite methane retrievals over Greenland, Atmos. Meas. Tech., 15, 4063–4074,, 2022.


October 2022:

Difference in number (a), duration (b), average SPI (c), and severity index (d) of drought events between the RCP8.5 (2050–2100) and historic (1950 to 2000) multi-model mean of 15 CMIP5 models. Here, a drought event is defined as any number of consecutive months with an SPI <-2. For the SPI calculation a gamma distribution and a representative timescale of 6 months are used. Figure 5 from Weigel et al., 2021 produced with the ESMValTool recipe recipe_martin18grl.yml.

Future changes in droughts in the CMIP5 ensemble analyzed with the Earth System Model Evaluation Tool (ESMValTool)

Motivated by the dry summer of this year, our image of the month shows projections of future changes in number, duration, average standard precipitation index (SPI), and severity of droughts as simulated by the climate models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) in the RCP8.5 future scenario. The SPI describes local precipitation anomalies and can be used to identify meteorological droughts. By the end of the 21stcentury the projections show an increase of drought number, duration and severity mainly in the already dry, subtropical areas. A shift towards more frequent droughts is found for most parts of the globe for drought indices which in addition to precipitation consider the potential evapotranspiration, like the Standardized Precipitation Evapotranspiration Index (SPEI) or are based on water balance calculation, like the scaled Palmer Drought Severity Index (SC-PDSI), see e.g. Ruhe (2022).

The diagnostics to produce this figure were processed with the Earth System Model Evaluation Tool (ESMValTool), a community diagnostics and performance metrics tool designed to improve comprehensive and routine evaluation of Earth System Models (ESMs). It provides a highly optimized preprocessor written in Python (Righi et al. 2020) centralizing standard operations on input data such as regridding or computation of multi-model statistics to allow for efficient and user-friendly data processing. A large number of diagnostics is provided with the ESMValTool, covering quasi-operational and comprehensive evaluation of ESMs, like the calculation and plotting of mean state, trends, and variability for atmospheric, oceanic, and terrestrial variables (Eyring et al. 2020), as well as diagnostics for emergent constraints and analysis of future projections from ESMs in CMIP (Lauer et al. 2020). Our image of the month is the result of an ESMValTool recipe (recipe_martin18grl.yml) that has been implemented in the ESMValTool version 2.0 and publishedin Weigel et al. (2020) that summarizes diagnostics for extreme events, regional and impact evaluation.


Eyring, V., Bock, L., Lauer, A., Righi, M., Schlund, M., Andela, B., Arnone, E., Bellprat, O., Brötz, B., Caron, L.-P., Carvalhais, N., Cionni, I., Cortesi, N., Crezee, B., Davin, E.L., Davini, P., Debeire, K., de Mora, L., Deser, C., Docquier, D., Earnshaw, P., Ehbrecht, C., Gier, B.K., Gonzalez-Reviriego, N., Goodman, P., Hagemann, S., Hardiman, S., Hassler, B., Hunter, A., Kadow, C., Kindermann, S., Koirala, S., Koldunov, N., Lejeune, Q., Lembo, V., Lovato, T., Lucarini, V., Massonnet, F., Müller, B., Pandde, A., Pérez-Zanón, N., Phillips, A., Predoi, V., Russell, J., Sellar, A., Serva, F., Stacke, T., Swaminathan, R., Torralba, V., Vegas-Regidor, J., von Hardenberg, J., Weigel, K., & Zimmermann, K.: Earth System Model Evaluation Tool (ESMValTool) v2.0 – an extended set of large-scale diagnostics for quasi-operational and comprehensive evaluation of Earth system models in CMIP. Geoscientific Model Development, 13, 3383-3438, 2020.

Lauer, A., Eyring, V., Bellprat, O., Bock, L., Gier, B.K., Hunter, A., Lorenz, R., Pérez-Zanón, N., Righi, M., Schlund, M., Senftleben, D., Weigel, K., & Zechlau, S.: Earth System Model Evaluation Tool (ESMValTool) v2.0 – diagnostics for emergent constraints and future projections from Earth system models in CMIP. Geoscientific Model Development, 13, 4205-4228, 2020.
Righi, M., Andela, B., Eyring, V., Lauer, A., Predoi, V., Schlund, M., Vegas-Regidor, J., Bock, L., Brötz, B., de Mora, L., Diblen, F., Dreyer, L., Drost, N., Earnshaw, P., Hassler, B., Koldunov, N., Little, B., Loosveldt Tomas, S., & Zimmermann, K.: Earth System Model Evaluation Tool (ESMValTool) v2.0 – technical overview. Geosci. Model Dev., 13, 1179-1199, 2020.

Ruhe, L.: A Comparison of Drought Indices in CMIP6 Climate Projections, M.Sc. Thesis, University of Bremen, Bremen, Germany,, 2022.

Weigel, K., Bock, L., Gier, B. K., Lauer, A., Righi, M., Schlund, M., Adeniyi, K., Andela, B., Arnone, E., Berg, P., Caron, L.-P., Cionni, I., Corti, S., Drost, N., Hunter, A., Lledó, L., Mohr, C. W., Paçal, A., Pérez-Zanón, N., Predoi, V., Sandstad, M., Sillmann, J., Sterl, A., Vegas-Regidor, J., von Hardenberg, J., and Eyring, V.: Earth System Model Evaluation Tool (ESMValTool) v2.0 – diagnostics for extreme events, regional and impact evaluation, and analysis of Earth system models in CMIP, Geosci. Model Dev., 14, 3159–3184,, 2021.


ESMValTool Website

September 2022:

Melt ponds on Arctic sea ice:
Melt ponds on Arctic sea ice: Acquisition of new in-situ measurements for improving satellite retrievals

Melt ponds are puddles of surface melt water that form on the Arctic sea ice during summer. They reduce the reflectivity and albedo of the surface and thus increase the energy absorption of incoming solar radiation. This self-enhancing process is called ice-albedo feedback and plays a significant role for the amplified warming of the Arctic compared to the mid-latitudes.

The polar remote sensing group at IUP works on a retrieval that extracts surface albedo and melt pond fraction from optical satellite measurements ( In the past two months we were on a RV Polarstern expedition in the Arctic ( to collect important in-situ data that will help to improve the forward model in the satellite retrieval. With a hyperspectral camera, called WALL-E, detailed spectra of melt pond and ice surfaces were measured to investigate their variability and evolution. Because of the spatial extension of the measurements also an estimate of the melt pond fraction is possible. To determine the driving forces of albedo change, these measurements were repeated on the same ice floes three times between 15th and 31st July 2022 (visit 1 to 3). The ice floe got “darker” during this time period (from 0.59 to 0.43 reflectance), which is mainly due to the increase of melt pond fraction but also a reduced reflectance of the ice surface (see bar plot in upper right).

Dr. Gunnar Spreen(
Hannah Niehaus (

August 2022:

Global satellite maps of N2O from GOSAT-2

Nitrous oxide (N2O) is after carbon dioxide (CO2) and methane (CH4) the most important anthropogenic greenhouse gas. N2O is mainly produced by natural processes (e.g. by bacteria), but is also released in agriculture (fertilisers) and biomass burning processes. Although the atmospheric concentration of N2O is about 1000 times smaller than for CO2, the global warming potential of N2O is almost 300 times larger, which makes it a relevant contributor to climate change.

However, global measurements of N2O especially from satellites are sparse. The map in the image (from Noël et. al, 2022) shows the atmospheric distribution of XN2O (i.e. the total dry column average N2O) in April 2019 derived at IUP from measurements of the Japanese GOSAT-2 satellite.

As can be seen, total column concentrations of N2O are larger in the tropics than at higher latitudes. This is because most of the N2O is produced and located in the lower atmosphere, the troposphere, which reaches higher in the tropics.


Noël et al., Retrieval of greenhouse gases from GOSAT and GOSAT-2 using the FOCAL algorithm, Atmos. Meas. Tech., 15(11), 3401–3437, 2022, doi: rm10.5194/amt-15-3401-2022.


Stefan Noël (

July 2022:

A recent study conducted by the IUP "Cloud Aerosol Surface PArameter Retrieval" group as part of the Transregio (AC)³ investigates clouds, their properties, and the magnitude of the radiative effect of clouds in the Arctic based on satellite data.

Arctic clouds are an important factor in influencing Arctic climate. They have both warming and cooling effects. The reasons whether they warm or cool depend on different conditions. These include the microphysical properties of the clouds, the illumination, the thermodynamic phase of the clouds (i.e., to what degree they contain liquid or frozen water), and the reflectivity of the ground.

Although when viewed over the entire Arctic the cloudiness has not changed appreciably over time, some cloud properties have. For example, the optical thickness of liquid water clouds and that of pure ice clouds. Partially significant changes in cloud brightness (cloud albedo) are also obvious.

Reference: The aerosol, cloud and surface property group "Pan-Arctic and regional trends of reflectance, clouds and fluxes: implications for Arctic Amplification" (2022)

Kontakt/contact: Marco Vountas
or Luca Lelli

June 2022:

Figure 1:
Figure 1: (a) Diamond HK36TTC-ECO Super Dimona D-KWHV aircraft operated by Jade Hochschule Wilhelmshaven and the team from the University Bremen und Jade Hochschule. (b) University of Bremen measurement suite comprising an in-situ greenhouse gas analyser from LRG and a 5-hole turbulence probe. Both are mounted underneath the right wing in an underwing pod. (c) Fligh-by at the ICOS measurement tower in Steinkimmen/Ganderkesee operated be the DWD during the second calibration flight. (d) Vertical wind profile measured during the very first calibration flight of the turbulence probe.

New measurement suite to improve Greenhouse Gases observations from aircraft,
University Bremen, Institute of Environmental Physics

Scientists at the Institute of Environmental Physics are specialised in the development and deployment of various types of airborne sensors to locate and quantify emissions of anthropogenic greenhouse gases by measuring atmospheric Methane (CH4) and Carbon dioxide (CO2) distributions. To accurately estimate emissions from those distributions, precise knowledge of the local wind field is necessary. Therefore, in spring 2022, a newly acquired 5-hole turbulence probe together with a in-situ greenhouse gas analyser, were successfully mounted in one of the wing pods of the motor glider of Jade Hochschule Wilhelmshaven. The Diamond HK36TTC-ECO is a flexible platform for different remote sensing tasks, as it needs no certification process to install sensors. First test flights in cooperation with Jade Hochschule Wilhelmshaven were successfully conducted in northern Germany to calibrate the new instrument. Additionally, fly-bys at the ICOS measurement tower in Steinkimmen/Ganderkesee operated by the German Weather Service DWD, which observes the vertical wind profile at 5 different altitudes till an altitude of 250m, were carried out. Those measurements are used for comparisons and will improve the calibration quality of the wind probe further. Additional, calibration flights comparing measurements from airborne, LIDAR, drone, and tower measurements are planned.


May 2022:

Example Figure:
Example Figure: TROPOMI (A) Kd-UVAB, (B) Kd-UVA, and (C) Kd-blue in the Atlantic Ocean gridded at 0.083° as mean for 11 May to 9 June 2018. Accordingly, Kd-UVAB, Kd-UVA, and Kd-blue measured in-situ at 19 stations during expedition PS113 are overlayed as diamonds (match-ups) and circles (unmatched stations). Adapted from Oelker et al. 2022, Fig. 5.

TROPOMI-retrieved underwater light attenuation in three spectral regions in the ultraviolet to blue. Frontiers in Marine Science 9: 787992. doi:10.3389/fmars.2022.787992

Sunlight plays an important role for biological, chemical, and physical processes in the ocean. High-energetic ultraviolet (UV) radiation can have damaging and beneficial effects for aquatic organisms and its interaction with the ocean is generally complex. Most processes feedback with climate warming. Satellite-based observations of light penetration into the ocean in combination with modeling are used to understand these processes and make predictions for the future ocean and climate scenarios in general. Traditional satellite ocean color sensors don’t measure the ultraviolet range. Information on UV light penetration is only inferred indirectly from measurements in the visible wavelength range, naturally connected to lager uncertainties.

This study exploited backscattered UV to blue light data at continuous spectral resolution of 0.5 nm of the TROPOMI sensor onboard the Sentinel-5-Precursor satellite. We present the first direct satellite-based observations of shortwave penetration, in terms of the diffuse attenuation (Kd) into the ocean ranging from the ultraviolet to the blue spectral domain. Our approach is based on Differential Optical Absorption Spectroscopy to retrieve the vibrational Raman scattering (VRS) signal and then combined with coupled ocean-atmosphere radiative transfer modeling (RTM) to derive Kd in the UV range (312.5-338.5 nm and 356.5-390 nm), additionally to the blue Kd (390-423 nm). The VRS signal is well detected in TROPOMI measurements (fit errors <15%) and TROPOMI Kd retrievals exhibit low sensitivity to parametrization of oceanic and atmospheric effects and show good agreement to in situ Kd obtained from in situ measured underwater light spectra.

These products have high potential satisfying user needs in the modeling community which require spectral information on shortwave light penetration for improving estimates of the ocean’s heat budget, primary productivity, photochemical reaction rates of climatically important compounds, and the UV dose rates as an indicator for damaging effects on aquatic organisms.

April 2022:

Figure: Colour ratio gradients (a,b) as altitudinal changes of the radiance ratio L(869nm)/L(460nm), as well as retrieved extinction coefficients at 869 nm (c,d) from OMPS-LP limb observations of two selected orbits. Areas of NaN values are scratched out. The map in the top row shows the location of the orbits and the volcano (triangle).

Hunga Tonga-Hunga Ha'apai eruption

On 15th of January 2022, the undervolcano Hunga Tonga-Hunga Ha'apai (20.55°S, 175.40°W), ejected material consisting of gas, steam, and ash up to an altitude of 58 km. This plume height is exceptional and the highest known since satellite observations. Even the second largest volcanic eruption of Mount Pinatubo in the Philippines in 1991 „only“ reached a height of 35 km. Probably, the combination of volcanic heat and superheated moisture from the ocean pushes the aerosols in such unprecedented altitudes.

The uppermost part of the plume sublimated quickly due to the extreme dryness in the mesosphere. In around 30 km altitude, the volcanic plume formed an extensive umbrella carried westward by the strong stratospheric winds. As can be seen from the Figure, this umbrella rose and expanded due to thermal buoyancy and dispersion. It circled the globe within two weeks. The volcanic aerosols will remain in the stratosphere for a long time but have no significant impact on the global climate. The injected aerosol content of 0.4 teragram of sulphur dioxide was too low for that. For comparison: Mount Pinatubo emitted about 18.5 teragram of sulphur dioxide into the stratosphere, temporarily lowering the global temperature by about 0.6°C.


Malinina, E., Rozanov, A., Niemeier, U., Wallis, S., Arosio, C., Wrana, F., Timmreck, C., von Savigny, C., and Burrows, J. P.: Changes in stratospheric aerosol extinction coefficient after the 2018 Ambae eruption as seen by OMPS-LP and MAECHAM5-HAM, Atmos. Chem. Phys., 21, 14871–14891,, 2021.

March 2022:

Figure: Deployment of an acoustic Doppler current profiler (ADCP) during expedition SO283 with the research vessel SONNE in the South Atlantic, April 2021 (top). Internal wave energy flux (red) and the pathway of Agulhas Eddies (blue) in the South Atlantic. Stars and circles denote deployment positions of instruments (bottom left). Current measurements (arrows) along the ship track and sea surface height (SSH) from satellite altimetry (bottom right).

Expedition number 283 of the German research vessel SONNE, which was carried out in spring 2021, was rather special. Due to the Corona pandemic, many research cruises had been cancelled, and now several oceanographic moorings in the South Atlantic had to be recovered and new instruments deployed. Scientists and students from four institutes embarked on the long voyage, which startedin Emden for infection control reasons and ended there after 64 days at sea and more than 17,000 nautical miles. Our group from IUP deployed two tall moorings with current meters and thermistors as well as 5 bottom-mounted inverted echo sounders, which is part of the observational program of the Collaborative Research Center TRR 181 “Energy Transfers in Atmosphere and Ocean” (

The objective of the experiment is to get many-months-long time series of current velocities and temperature to study the interaction of internal gravity waves emanating from the Walvis Ridge and Agulhas Eddies that travel from south of Africa into the South Atlantic. These interactions between waves and eddies and their exchange of kinetic energy are not yet fully understood. Tidal currents at seamounts and continental shelves, for example, exciteinternal gravity waves that can travel hundreds of kilometers across ocean basins before eventually contributing to turbulent mixing in the water column. Their loss of energy is determined by a number of processes and interactions. In our project, we study the scattering and refraction of tidal waves in the interior of the oceans with the goal of better integrating these processes intoclimate models.

February 2022:

Figure: GMAP-2021 campaign in South Korea. Top left: Group picture of some participants. Top right: The IUP Bremen MAX-DOAS instrument on the roof of the NIER building. Bottom figures: Two days of NO2 observations in the Seoul Metropolitan Area. Background values are GEMS tropospheric NO2 columns retrieved by IUP Bremen. The overlaid line are matching car-DOAS measurements performed within +/- 2 hours of the satellite overpass.

The GMAP-2021 campaign
In October and November 2021, IUP Bremen participated in the GMAP-2021 (GEMS Map of Air Pollution) campaign in South Korea. This campaign brought together instruments from South Korea, the US, Belgium, the Netherlands, and Germany to collect data on air pollution in South Korea for the validation of the GEMS (Geostationary Environment Monitoring Spectrometer) satellite instrument. GEMS is the first geostationary satellite instrument dedicated to air quality, launched by the South Korean space agency in February 2020 and providing measurements of key pollutants over Asia in hourly resolution. Similar instruments will be launched by the US (TEMPO) and Europe (Sentinel-4) in the coming years.

During the GMAP campaign, IUP Bremen installed a multi-azimuth MAX-DOAS instrument on the rooftop of the NIER (National Institute of Environmental Research) building at Incheon. This instrument provides data on the abundance of NO2, HCHO, SO2and other pollutants valuable for long-term validation of GEMS retrievals. It will continue to operate at the location for the coming months. Similar instruments have been deployed by other groups in the Seoul Metropolitan Area and other parts of South Korea, and are expected to help in characterizing and improving the satellite data products.

In addition to the stationary measurements, a large number of mobile DOAS measurements was performed from cars to investigate the spatial variability of NO2, to monitor temporal changes and to estimate NOx emissions from Seoul. The observations show a large variability of NO2within individual GEMS satellite pixels, rapid changes over time and the impact of both, local emissions and transport to the observed NO2column amounts. These data are complemented by both in-situ and remote sensing observations from several aircraft overpasses over the same area.

Further reading:

Jhoon Kim et al., New Era of Air Quality Monitoring from Space: Geostationary Environment Monitoring Spectrometer (GEMS), BAMS, 2020, DOI:

January 2022:

Figure 1:
Figure 1: The Carbon Dioxide (CO2) plume emitted from the power plant Jänschwalde observed with the new airborne instrument MAMAP2D-Light could be detected up to 10 km away from the power plant itself. Reddish colors indicate high CO2 concentrations, while the white areas have been filtered out due to too low signal over water (Image: IUP, University Bremen).

Scientists of the Institute of Environmental Physics are developing new optical sensors to image atmospheric Methane (CH4) and CO2 distributions from aircraft in a similar way as future satellites, but with higher sensitivity to point source emissions. The prototype of the push broom imaging spectrometer – called MAMAP2D-Light – was successfully mounted in the wing pods of the motor glider of Jade Hochschule Wilhelmshaven (see IUP picture of the month June 2021) and successfully performed its first flight over the power plant Jänschwalde near Cottbus, Germany.

Subsequent data analysis resulted in the image shown above, where the CO2 plume of the power plant Jänschwalde is clearly visible in reddish colors. The emissions estimated from this data set matches within its uncertainty range the average emissions during the week of the overflight, and the results and additional performance characteristics have been presented at the fall meeting of the American Geophysical Union in December 2021 (1). In summer 2022 this instrument will be flown onboard the German high altitude research aircraft HALO as part of the international COMET 2.0 campaign targeting high latitude emissions of CH4 and CO2 from wetlands as well as geological seeps and oil/gas production in Canada.


(1) Jakob Borchardt, Konstantin Gerilowski, Sven Krautwurst, Wilke Thomssen, Jan Franke, Martin Kumm, Pascal Janßen, Jens Wellhausen, Heinrich Bovensmann, John P. Burrows(2021), The New Imaging Spectrometer MAMAP2D-Light– Initial Calibration and First Measurement Results, [A25G-1759] presented at 2021 Fall Meeting, AGU, 13-17 Dec.,