The MONARCH high-resolution reanalysis of desert dust aerosol over Northern Africa, the Middle East and Europe (2007–2016)
- 1Barcelona Supercomputing Center (BSC), Barcelona, Spain
- 2Geophysical Fluid Dynamics Laboratory (GFDL), Princeton, NJ, USA
- 3Consiglio Nazionale delle Ricerche-Istituto di Scienze dell’Atmosfera e del Clima (CNR-ISAC), Italy
- 4Izaña Atmospheric Research Center, AEMET, Santa Cruz de Tenerife, Spain
- 5LISA, UMR CNRS 7583, Université Paris-Est-Créteil, Université de Paris, Institut Pierre-Simon Laplace (IPSL), Créteil, France
- 6Universitat Politècnica de Catalunya - BarcelonaTech (UPC), Department of Project and Construction Engineering, Terrassa, Spain
- 7Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK-TRO), Department Troposphere Research, Karlsruhe, Germany
- 8Consiglio Nazionale delle Ricerche-Istituto di Metodologie per l’Analisi Ambientale (CNR-IMAA), Italy
- 9NASA Goddard Institute for Space Studies (GISS), New York, NY, USA
- 10Department of Earth Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, the Netherlands
- 11University of Twente, Department of Governance and Technology for Sustainability (BMS-CSTM), the Netherlands
- 12Finnish Meteorological Institute (FMI), Weather and Climate Change Impact Research, Finland
- 13State Meteorological Agency (AEMET), Spain
- 14ICREA, Catalan Institution for Research and Advanced Studies, Barcelona, Spain
- 1Barcelona Supercomputing Center (BSC), Barcelona, Spain
- 2Geophysical Fluid Dynamics Laboratory (GFDL), Princeton, NJ, USA
- 3Consiglio Nazionale delle Ricerche-Istituto di Scienze dell’Atmosfera e del Clima (CNR-ISAC), Italy
- 4Izaña Atmospheric Research Center, AEMET, Santa Cruz de Tenerife, Spain
- 5LISA, UMR CNRS 7583, Université Paris-Est-Créteil, Université de Paris, Institut Pierre-Simon Laplace (IPSL), Créteil, France
- 6Universitat Politècnica de Catalunya - BarcelonaTech (UPC), Department of Project and Construction Engineering, Terrassa, Spain
- 7Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK-TRO), Department Troposphere Research, Karlsruhe, Germany
- 8Consiglio Nazionale delle Ricerche-Istituto di Metodologie per l’Analisi Ambientale (CNR-IMAA), Italy
- 9NASA Goddard Institute for Space Studies (GISS), New York, NY, USA
- 10Department of Earth Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, the Netherlands
- 11University of Twente, Department of Governance and Technology for Sustainability (BMS-CSTM), the Netherlands
- 12Finnish Meteorological Institute (FMI), Weather and Climate Change Impact Research, Finland
- 13State Meteorological Agency (AEMET), Spain
- 14ICREA, Catalan Institution for Research and Advanced Studies, Barcelona, Spain
Abstract. One of the challenges in studying desert dust aerosol along with its numerous interactions and impacts is the paucity of direct in-situ measurements, particularly in the areas most affected by dust storms. Satellites typically provide columnintegrated aerosol measurements, but observationally-constrained continuous 3D dust fields are needed to assess dust variability, climate effects and impacts upon a variety of socio-economic sectors. Here, we present a high resolution regional reanalysis data set of desert dust aerosols that covers Northern Africa, the Middle East and Europe along with the Mediterranean sea and parts of Central Asia, and the Atlantic and Indian Oceans between 2007 and 2016. The horizontal resolution is 0.1° latitude × 0.1° longitude, and the temporal resolution is 3 hours. The reanalysis was produced using Local Ensemble Transform Kalman Filter (LETKF) data assimilation in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) developed at the Barcelona Supercomputing Center (BSC). The assimilated data are coarse-mode dust optical depth retrieved from the Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue Level 2 products. The reanalysis data set consists of upper air (dust mass concentrations and extinction coefficient), surface (dust deposition and solar irradiance fields, among them) and total column (e.g., dust optical depth and load) variables. Some dust variables, such as concentrations and wet and dry deposition, are expressed for a binned size distribution that ranges from 0.2 to 20 μm in particle diameter. Both analysis and first-guess (analysis-initialized simulation) fields are available for the variables that are diagnosed from the state vector. A set of ensemble statistics is archived for each output variable, namely the ensemble mean, standard deviation, maximum and median. The spatial and temporal distribution of the dust fields follows well-known dust cycle features controlled by seasonal changes in meteorology and vegetation cover. The analysis is statistically closer to the assimilated retrievals than the first-guess, which proves the consistency of the data assimilation method. Independent evaluation using AERONET dust-filtered optical depth retrievals indicates that the reanalysis data set is highly accurate (mean bias = −0.05, RMSE = 0.12, r = 0.81 when compared to retrievals from the spectral de-convolution algorithm on a 3-hourly basis). Verification statistics are broadly homogeneous in space and time with regional differences that can be partly attributed to model limitations (e.g., poor representation of small-scale emission processes), presence of aerosols other than dust in the observations used in the evaluation, and differences in the number of observations among seasons. Such a reliable high-resolution historical record of atmospheric desert dust will allow a better quantification of dust impacts upon key sectors of society and economy, including health, solar energy production and transportation. The reanalysis data set (Di Tomaso et al., 2021) is distributed via a Thematic Real-time Environmental Distributed Data Service (THREDDS) at BSC and freely available at http://hdl.handle.net/21.12146/c6d4a608-5de3-47f6-a004-67cb1d498d98.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(8330 KB)
-
Supplement
(5360 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(8330 KB) -
Supplement
(5360 KB) - BibTeX
- EndNote
Journal article(s) based on this preprint
Enza Di Tomaso et al.
Interactive discussion
Status: closed
-
RC1: 'Comment on essd-2021-358', Anonymous Referee #1, 26 Nov 2021
The MONARCH high-resolution reanalysis of desert dust aerosol over Northern Africa, the Middle East and Europe (2007-2016)
Enza Di Tomaso et al.
SUMMARY
The manuscript describes a data set of assimilated dust aerosol concentrations and optical properties. The data covers an area including Europe, North Africa, and Middle East, at a resolution of 0.1x0.1 degrees. The time period is 2007-2016, and data is available every 3 hour. Simulations were done using the MONARCH atmospheric model. AOD observations of MODIS were assimilated using an ensemble approach (LETKF), the product consists of ensemble mean, std.dev., maximum, and median.
The production of the data set is extensively described. Some extra clarification could be made (see SPECIFIC COMMENTS below), but in general the assimilation procedure is described sufficiently detailed to understand what has been done, including the generation of the ensemble, calibration, data selection, error statistics used, and assimilation sequence. Where necessary references are present where details could be found. With this the paper could serve as a reference for studies that actually use the data set, as indicated in section 7.
With some minor clarifications the manuscript could be published in this journal.
GENERAL COMMENTS
The manuscript describes a validation of the data set in terms of AOD (or specific, DOD, Dust Optical Depth). This is also the quantity that is assimilated, and it therefore makes sense to use this as first validation. For a data set related to dust, it would however be useful to have also an idea on the dust concentrations themselves, and how accurate these are. The only dust-related results are shown in Figure 4, but no comparison with observations has been made. Will there be a validation of the dust concentrations included in the follow-up papers mentioned in Section 7? It would be useful to have that clearly mentioned. Also, some remarks could be made already on the dust concentrations themselves and how they are changed by the analysis.
For example, what is the impact of the calibration described in Section 6.1 on the dust load in the ensemble members? If my interpretation is correct, the calibration factors for the dust emissions range from 0.004 for the K14 emission scheme, to 2.65 for the MB95 scheme. This is a huge difference; does it mean that the K14 scheme by default has a huge over-estimation? After calibration, do the ensemble members have dust concentrations that are more or less in the same range?
The adjustment of the dust concentrations depends strongly on how DOD is calculated, thus on the optical properties and the radiance computations. Is there any idea on how accurate these computations are? With incorrect optical properties computed, the dust concentrations might require unrealistic perturbations to obtain the correct DOD’s. The meteorological data is also relevant for this computation I guess; since this comes in the ensemble from two different models (MERRA2 and ERA-Interim), is there a clear difference seen between the DOD’s computed for different meteo?
Some clarification on the ensemble generation would be useful. Section 3 describes that a 12-member ensemble is used, with each member choosing one-of-two meteo sets, one-of-three emission schemes, and a random value for (among others) the friction threshold; is that indeed what is done? I guess that each member then keeps it’s choice for meteo and emission scheme, but are the emission parameters changing in time or per grid cell?
SPECIFIC COMMENTS
Lines 145-146: What does a Frequency-of-Occurrence of 0.20 mean? That in 20% of the days dust is observed over a location?
Line 230: What is meant with a time slot centered around 12 UTC? Aren’t more MODIS orbits assimilated then, with different time slots?
The analysis weights in Eq. (2) might provide some information on a preference of the analysis for certain ensemble members, for example the members with a specific emission scheme. Is that indeed possible, and has some information been deduced from them?
Section 6 describes that an observation screening is applied. Is it kept which fraction of the observations has been rejected, and whether that is especially in certain regions?
Table 6: What is exactly done for “DOD-mixed2”? How could AE be <0.75 and >1.2 ?
Line 492: The “DOD-mixed2” leads to more zero values; should that be visible in Figure 10 then as a a higher density?
Line 502. Is the change in statistics associated with changed conditions? Or could it be related to a degradation of the data?
Line 510: Which complexity of the topography is relevant here?
SPELL AND GRAMMAR
Line 78: chemical formula should not in math mode
Line 79: “additionally”
Lines 197, 199: shouldn’t it be: ‘friction velocity “for” wind velocity’ ?
Subsection 6.1: shouldn’t this be a section?
-
RC2: 'Comment on essd-2021-358', Julie Letertre-Danczak, 24 Feb 2022
This article on the data set of reanalysis from MONARCH over Nothern Africa, Middle East and Europe for 2007-2016 is well documented.
The explanation of the data set is rich of information and the metadata are clear.
However I have few comments on the data set:
Some confusion are made on between Modis Collection 6 and collection 6.1 (some clarifications are needed).
In the table 1 which is the overview of the experiment that has generated the data is missing the data assimilation window (it is mentionned at the earliest in section 5, adding it in the table will be beneficial).
The dust bins description is referred but again a table making an easy finding of the information will be helpful for someone who would like to use the data.
The organisation of the sections should be revisited (the section 6 should come directly after section 3 as it is more details on the model.2 more general comments:
1) The output of this experiment should be compared to a denial experiment where no data assimilation will be performed to understand why the emission scheme is over producing (issue with the scheme and the climatology behind? scheme loaded with dust transport? ...) That will be interesting to conduct maybe for an other article.
2) The comparison with an independent set of data: Modis is evaluated against AERONET, so if your system is converging toward your data by your data assimilation, it seems logical that it will also converge toward AERONET even if it is not a guarantee. I would suggest that you add another independant det of data in the evaluation process for the eventual article that I have suggested before (lidar, dry deposition measurement, ...).The bibliography has not been checked.
Final comment:
This article on the description of a dataset correspond to the criteria expected for a data journal. -
AC1: 'Comment on essd-2021-358', Enza Di Tomaso, 09 Apr 2022
The comment was uploaded in the form of a supplement: https://essd.copernicus.org/preprints/essd-2021-358/essd-2021-358-AC1-supplement.pdf
Peer review completion




Interactive discussion
Status: closed
-
RC1: 'Comment on essd-2021-358', Anonymous Referee #1, 26 Nov 2021
The MONARCH high-resolution reanalysis of desert dust aerosol over Northern Africa, the Middle East and Europe (2007-2016)
Enza Di Tomaso et al.
SUMMARY
The manuscript describes a data set of assimilated dust aerosol concentrations and optical properties. The data covers an area including Europe, North Africa, and Middle East, at a resolution of 0.1x0.1 degrees. The time period is 2007-2016, and data is available every 3 hour. Simulations were done using the MONARCH atmospheric model. AOD observations of MODIS were assimilated using an ensemble approach (LETKF), the product consists of ensemble mean, std.dev., maximum, and median.
The production of the data set is extensively described. Some extra clarification could be made (see SPECIFIC COMMENTS below), but in general the assimilation procedure is described sufficiently detailed to understand what has been done, including the generation of the ensemble, calibration, data selection, error statistics used, and assimilation sequence. Where necessary references are present where details could be found. With this the paper could serve as a reference for studies that actually use the data set, as indicated in section 7.
With some minor clarifications the manuscript could be published in this journal.
GENERAL COMMENTS
The manuscript describes a validation of the data set in terms of AOD (or specific, DOD, Dust Optical Depth). This is also the quantity that is assimilated, and it therefore makes sense to use this as first validation. For a data set related to dust, it would however be useful to have also an idea on the dust concentrations themselves, and how accurate these are. The only dust-related results are shown in Figure 4, but no comparison with observations has been made. Will there be a validation of the dust concentrations included in the follow-up papers mentioned in Section 7? It would be useful to have that clearly mentioned. Also, some remarks could be made already on the dust concentrations themselves and how they are changed by the analysis.
For example, what is the impact of the calibration described in Section 6.1 on the dust load in the ensemble members? If my interpretation is correct, the calibration factors for the dust emissions range from 0.004 for the K14 emission scheme, to 2.65 for the MB95 scheme. This is a huge difference; does it mean that the K14 scheme by default has a huge over-estimation? After calibration, do the ensemble members have dust concentrations that are more or less in the same range?
The adjustment of the dust concentrations depends strongly on how DOD is calculated, thus on the optical properties and the radiance computations. Is there any idea on how accurate these computations are? With incorrect optical properties computed, the dust concentrations might require unrealistic perturbations to obtain the correct DOD’s. The meteorological data is also relevant for this computation I guess; since this comes in the ensemble from two different models (MERRA2 and ERA-Interim), is there a clear difference seen between the DOD’s computed for different meteo?
Some clarification on the ensemble generation would be useful. Section 3 describes that a 12-member ensemble is used, with each member choosing one-of-two meteo sets, one-of-three emission schemes, and a random value for (among others) the friction threshold; is that indeed what is done? I guess that each member then keeps it’s choice for meteo and emission scheme, but are the emission parameters changing in time or per grid cell?
SPECIFIC COMMENTS
Lines 145-146: What does a Frequency-of-Occurrence of 0.20 mean? That in 20% of the days dust is observed over a location?
Line 230: What is meant with a time slot centered around 12 UTC? Aren’t more MODIS orbits assimilated then, with different time slots?
The analysis weights in Eq. (2) might provide some information on a preference of the analysis for certain ensemble members, for example the members with a specific emission scheme. Is that indeed possible, and has some information been deduced from them?
Section 6 describes that an observation screening is applied. Is it kept which fraction of the observations has been rejected, and whether that is especially in certain regions?
Table 6: What is exactly done for “DOD-mixed2”? How could AE be <0.75 and >1.2 ?
Line 492: The “DOD-mixed2” leads to more zero values; should that be visible in Figure 10 then as a a higher density?
Line 502. Is the change in statistics associated with changed conditions? Or could it be related to a degradation of the data?
Line 510: Which complexity of the topography is relevant here?
SPELL AND GRAMMAR
Line 78: chemical formula should not in math mode
Line 79: “additionally”
Lines 197, 199: shouldn’t it be: ‘friction velocity “for” wind velocity’ ?
Subsection 6.1: shouldn’t this be a section?
-
RC2: 'Comment on essd-2021-358', Julie Letertre-Danczak, 24 Feb 2022
This article on the data set of reanalysis from MONARCH over Nothern Africa, Middle East and Europe for 2007-2016 is well documented.
The explanation of the data set is rich of information and the metadata are clear.
However I have few comments on the data set:
Some confusion are made on between Modis Collection 6 and collection 6.1 (some clarifications are needed).
In the table 1 which is the overview of the experiment that has generated the data is missing the data assimilation window (it is mentionned at the earliest in section 5, adding it in the table will be beneficial).
The dust bins description is referred but again a table making an easy finding of the information will be helpful for someone who would like to use the data.
The organisation of the sections should be revisited (the section 6 should come directly after section 3 as it is more details on the model.2 more general comments:
1) The output of this experiment should be compared to a denial experiment where no data assimilation will be performed to understand why the emission scheme is over producing (issue with the scheme and the climatology behind? scheme loaded with dust transport? ...) That will be interesting to conduct maybe for an other article.
2) The comparison with an independent set of data: Modis is evaluated against AERONET, so if your system is converging toward your data by your data assimilation, it seems logical that it will also converge toward AERONET even if it is not a guarantee. I would suggest that you add another independant det of data in the evaluation process for the eventual article that I have suggested before (lidar, dry deposition measurement, ...).The bibliography has not been checked.
Final comment:
This article on the description of a dataset correspond to the criteria expected for a data journal. -
AC1: 'Comment on essd-2021-358', Enza Di Tomaso, 09 Apr 2022
The comment was uploaded in the form of a supplement: https://essd.copernicus.org/preprints/essd-2021-358/essd-2021-358-AC1-supplement.pdf
Peer review completion




Journal article(s) based on this preprint
Enza Di Tomaso et al.
Data sets
MONARCH high-resolution reanalysis data set of desert dust aerosol over Northern Africa, the Middle East and Europe Di Tomaso, E., Escribano, J., Basart, S., Macchia, F., Benincasa, F., Bretonnière, P.-A., Buñuel, A., Castrillo, M., Gonçalves, M., Jorba, O., Klose, M., Montané, G., Olid, M., Pérez García-Pando, C. http://hdl.handle.net/21.12146/c6d4a608-5de3-47f6-a004-67cb1d498d98
Enza Di Tomaso et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
789 | 268 | 38 | 1,095 | 75 | 22 | 23 |
- HTML: 789
- PDF: 268
- XML: 38
- Total: 1,095
- Supplement: 75
- BibTeX: 22
- EndNote: 23
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Cited
1 citations as recorded by crossref.
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(8330 KB) - Metadata XML
-
Supplement
(5360 KB) - BibTeX
- EndNote