the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Theoretical assessment of the ability of the MicroCarb satellite city-scan observing mode to estimate urban CO2 emissions
Kai Wu
Denis Jouglet
Liang Feng
Download
- Final revised paper (published on 30 Jan 2023)
- Preprint (discussion started on 08 Sep 2022)
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2022-749', Anonymous Referee #1, 21 Oct 2022
This publication investigates the potential of the MicroCarb satellite mission planned to be launched in the 2023/2024 timeframe to estimate urban CO2 emissions. The potential is analyzed for Paris and London using synthetic model experiments for a total of 15 days distributed over three representative months of the year with varying illumination conditions, biospheric activity, and cloud cover.
The MicroCarb mission will have three different observation configurations including a "city-mode", which allows observing a target with a two times higher spatial resolution (2.25 km x 2.25) and a two or three times wider swath than in the nominal nadir mode by performing two or three pitch maneuvers when flying over the target. The study is important as it provides a scientific basis for the decision on whether a two-sweep or three-sweep mode should be implemented.
The manuscript is very well written and concise, the model simulations are state-of-the-art, and the figures are informative. I have only a few minor points and one main concern that actually not only applies to this study but to other similar studies as well. I consider the manuscript acceptable after addressing these (primarily minor) points.
Main concern:
Like other studies employing an OSSE approach with synthetically generated observations, the conclusions on the ability of the satellite to constrain (urban) emissions are quite optimistic, and I would argue too optimistic. Some of the limitations are addressed in the manuscript, e.g. in the "Conclusions and discussion" section where the results are presented as "best-case" scenario. I am worried, however, that papers like this generate unrealistic high expectations of city governments or other stakeholders, which the scientific community, once the satellites are in orbit, will not be able to meet. The study is very carefully designed and executed, but in my view the discussion on the limitations of the study induced by the different assumptions is falling short. Probably the easiest solution would be to add a paragraph or two to the "conclusions and discussion" section, which so far actually contains very little discussion.
Here is a (incomplete) list of my concerns:
- My general concern is that systematic errors in any of the components (e.g. atmospheric transport, observations, biospheric fluxes) are much more critical than random errors. The latter can always be reduced by increasing the number of observations, and since most OSSEs only consider random errors, the results tend to be overly optimistic. In the present study, the three-sweep mode simply leads to better results because of the larger number of observations, but this is not the only point that should be considered when comparing the two modes. There is only a short discussion on the critical issue of determining the background against which the urban enhancements are measured. Reliably determining the background concentration levels will likely be essential, especially in the presence of systematic observation biases, and will need a sufficient number of cloud-free pixels outside of the urban plume. This may be another important advantage of the three-sweep mode that is not addressed in this study.
- The closed-loop experiment assumes unbiased transport errors. Although this limitiation is acknowledged by the authors, there is essentially no discussion on how such errors could affect the emission estimation. Is there no information available from the meteorological community about biases and errors in their analyzed wind products? Relative errors in wind speeds are almost certainly largest under low wind conditions, but these are the conditions that produce the strongest signals in an OSSE experiment as presented here. Thus, the situations working best in this OSSE may actually be the most problematic in reality. Probably there is an optimal range between low wind speeds, where relative errors are large, and large wind speeds, where the urban enhancements are too small to be detected. Unfortunately, the present study presents no insights into this question.
- Cloud cover is also treated in a random fashion. Depending on cloud cover, a probability is estimated for a pixel to be cloud-free or not, which results in a random pattern of cloud-free observations. However, clouds tend to be organized and to obscure one part of the image more than another part. As a result, the satellite will likely be blind to emissions in some parts of the city, an effect that is not captured by a purely random pattern. Furthermore, with increasing cloudiness, also cloud-free pixels tend to show biases in XCO2 due to 3D cloud radiation effects (e.g. Massie et al., 2017; https://doi.org/10.1002/2016JD026111). As a result, the ability to constrain urban emissions will likely drop faster with increasing cloud cover than suggested by Fig. 6. Figure 6 would tell the operators of Microcarb that a 50% cloud cover would reduce the error reduction for Paris only from 13% to 10% as compared to a cloud-free case (for a 1 ppm observation uncertainty). I doubt that this will be true in reality.
- Point sources are excluded in the OSSE, but point sources are an important part of reality, not only outside but also inside cities. In many cities, emissions from waste incinerators, combined heat and power plants or industrial stacks make up a large fraction of the total CO2 emissions and the emissions are usually not constant in time. X-STILT, to my knowledge, only computes the sensitivity to emissions at the surface but not to elevated emissions. The mixture of emissions released near the surface and from stacks will pose a major challenge to quantify urban emissions.
- Microcarb as other satellites will only measure at a single time of the day and will observe a given city only very rarely in a given year (actually, it would be good to know how many samples can be obtained over a city under optimal conditions). Emissions have a diurnal cycle and change from day-to-day. However, cities would like to better know their annual mean emissions and how they evolve over the years. Because of the uncertainties in temporal variability of emissions, observations on a few days per year (if possible at all) will provide only small constraints on annual mean emissions. Some of these challenges were addressed in Broquet et al. (2018; https://doi.org/10.5194/amt-11-681-2018), which also performed an OSSE for Paris.
Minor points and corrections:
- I had to look up other documentation of MicroCarb to fully understand the city-mode. If the satellite would only perform pitch-maneuvers (as I thought after reading section 2.1), it would miss many cities or observe them only partially. However, as stated on https://microcarb.cnes.fr/en/MICROCARB/GP_mission.htm, the pointing-mirror allows also movements along the roll-axis by pointing the line of sight on either side of the satellite ground track up to ±200 km. This should be explained more clearly, because otherwise it is difficult to understand why the cities are exactly in the center of the scans.
Furthermore, I assume that the sampling pattern shown in Figure 2 only holds when the city is exactly in nadir. How does the pattern look like when the satellite points some 100-200 km to the side (would be good to show e.g. in the Annex)? This will likely be necessary in many cases to have good coverage of a given target. Is the distortion of the pattern and the degradation of the resolution important or not? Would that affect the conclusions of the study? - It is probably correct that urban areas are responsible for ~70% of global emissions, but at the same time it is clear that only a fraction of this is actually emitted within the city boundaries, since a large part of the power consumed by cities is generated by power plants outside. I always find this number of 70% misleading when it is presented without further explanation.
- The introduction section misses a number of references that I consider important for this paper:
Other CO2 satellite OSSE studies for urban areas: Broquet et al. 2018 (see above) and Kuhlmann et al. (2020; https://doi.org/10.5194/amt-13-6733-2020).
Study on CO2 emissions of the city of London based on aircraft observations: Pitt et al. (2019; https://doi.org/10.5194/acp-19-8931-2019).
Study on large CO2 emitters (including cities) as observed from OCO-2: Chevallier et al. (2022; https://doi.org/10.1029/2021GL097540). - Page 3, Line 60: Change "pass over" to "passes over"
- P3, L68: Change "analyses" to "analysis" (singular form needed here)
- P3, L87: The ACT swath seems much closer to 30 km than to 40 km, and similarly for the 3-sweep mode seems closer to 45 km than to 60 km.
- P4, L98: Change "where g defined" to "where g is defined"
- P4, L110: I was confused by the usage of the term "scene", which apparently refers to a single observation/pixel rather than a full "scene" of a given area.
- P5, L143: The high correlation coefficient of 0.98 gives a false impression of accuracy. It only tells us that the diurnal cycle is on average well represented, which seems comparatively trivial as it follows the mean diurnal cycle of radiation. The more important message is that the fluxes are in a realistic range.
- P7, L187: Please provide more details on the "realistic noise model". Which parameters determine the noise and how was it applied in the present study? Do you explicitly consider, for example, surface reflectance, and if so, based on which data set?
- P7, L198: What is the assumption of a 25% random error for the biological fluxes based on? I can't track this value from the material presented in the paper or in Figure A2. Is this uncertainty assumed to be purely random and spatially uncorrelated?
- P7, L200: The eigenvalue decomposition sounds interesting. Please provide more information on how it works.
Citation: https://doi.org/10.5194/egusphere-2022-749-RC1 - AC1: 'Reply on RC1', Kai Wu, 12 Dec 2022
-
RC2: 'Comment on egusphere-2022-749', Anonymous Referee #2, 21 Oct 2022
GENERAL COMMENTS
This study evaluates the ability of the upcoming MicroCarb satellite to retrieve anthropogenic CO2 emissions based on synthetic measurements. The study follows a well-established approach for such observation system simulation experiments: generating perfect observations based on flux inventories and a model of atmospheric transport, perturbing them based on a model of instrument error, using the perturbed observations to retrieve fluxes in an inverse model of atmospheric transport, and analyzing these fluxes.
While the modeling and analysis methods are sound, I have some concerns on the setup of the experiments and, consequently, applicability of the results to real measurements. In addition, I believe that the interpretation of a key result is inaccurate.
The manuscript is mostly easy to read. A few paragraphs, especially on methods and experimental setup, need clarifications and additional information to better understand the results.
In summary, the manuscript is valuable for the CO2 flux estimation community to understand the capabilities of the upcoming MicroCarb satellite. However, especially the major issues summarized above need to be addressed before publication, which might require additional simulations or better communication of the limitations in abstract and conclusions.
SPECIFIC COMMENTSIn my opinion, the major issues to address are:
- The authors omit large point sources in their analyses. Therefore, I'm not sure how the results relate to real measurements, which of course see the integrated signal of all sources (see comments on Sect. 2.3).
- I believe that the authors' interpretation of their results on biogenic fluxes is not accurate (see comments on Sect. 3.3).
- Information on the optimization method is missing, namely how the state vector is set up and, related, how the posterior anthropogenic flux component is obtained (see comments on Sect. 2.7)
- The manuscript needs more information on the observation error model (see comments on Sect. 2.7)Below I elaborate on these major and some other points in detailed comments arranged by section.
Abstract or Introduction:
State somewhere early on that the study analyses only one out of many sources of uncertainty in emission estimation with MicroCarb, i.e. random measurement errors (as discussed in Sect. 4).2.2. Cloud screening
This is probably a minor issue, but please comment:
The authors decide cloudiness with a random number generator while in reality there should be some spatial correlation when clouds are bigger than satellite pixels. So e.g. for cloud cover 50%, the analysis is valid for clouds that are smaller or similar in size as the satellite pixel. I don't know whether this is a realistic assumption or what impact larger clouds would have.
2.3 Anthropogenic emissions
The authors exclude point sources from the analysis, citing smaller uncertainty than the areal anthropogenic sources they do include. I don't understand that rationale because anthropogenic fluxes are the main focus of the study and the real measurements will of course include all sources. I think that this omission is problematic, because it alters some results of the study:
- As mentioned by the authors, the prior uncertainty of the areally diffuse sources is higher than that of the point sources. Therefore, the random prior uncertainty used in this study is higher than that of the total anthropogenic emissions. Therefore, the reduction of random uncertainty (the "skill" of the observation system) of city-wide integrated anthropogenic emissions is overestimated.
- The atmospheric signal is reduced due to the omission of point sources, which could imply that (relative) biases in posterior emissions are overestimated in this study.
- Omitting point sources means that this study is about an artificial problem that includes only areally diffuse anthropogenic emissions, whereas the abstract includes statements on "total emissions".Every OSSE underestimates the skill of an observation system because they cannot account for/quantify all sources of uncertainty, and one may argue that it is fine to only focus on the harder flux component. However, as explained above, some uncertainties are overestimated while some are underestimated w.r.t. total anthropogenic emissions, so how the results relate to the "real" problem of integrated signals is complex. One way out could be to argue that point source plumes may be removed from the real MicroCarb signal (i.e. reducing the real problem to the one studied here), but this would introduce new errors, and their analysis would fill its own dedicated paper.
I think that the best way forward would be to repeat the analyses but with the strong point sources included. Alternatively, the limitations arising from the omission need to be clearly stated and discussed in the manuscript (but in my opinion this would diminish the value of the manuscript because only a partial problem is studied).
2.7 ExperimentsThe section on measurement errors needs to be expanded: please provide a reference or add a dedicated, if short, section on synthetic measurement errors, including how they were estimated and possible limitations (e.g. if there are error sources that were not considered).
Please clarify:
Are the mean prior fluxes the true fluxes + bias or the true fluxes + bias + random flux error realization? The text sounds like the latter, but in figures 4 and 5, the variations in the prior look very close to those of the true fluxes. Given that the random error in the fossil fluxes is the same as the bias, and the correlation length is small compared to the satellite swath/sweep, I would expect to see differences in the structures more clearly, i.e. on the order of the bias between prior and truth.More information on the state vector is needed:
- How are anthropogenic emissions inferred from the posterior fluxes? Are individual flux components (e.g. anthropogenic, GPP, respiration) optimized as separate state vector elements? Or are total emissions optimized and then partitioned (how?) into the individual flux components?
- On what resolution are fluxes estimated, i.e. what is the size of the state vector?
- Is there a temporal resolution?
A new section dedicated to the state vector could be helpful for structuring this additional information.
3.3 Biogenic fluxesUpon first reading the manuscript, only here did I realize that results so far only considered the anthropogenic flux component. In hindsight it's obvious, because figures 4 and 5 don't show negative fluxes/enhancements. Nonetheless, this approach is unexpected and needs to be made clearer throughout the manuscript, e.g. in the introductory paragraph of Sect. 3 (for example, there are 18 scenarios, not 9), in subsection headings, and in figure captions (and perhaps in Sect. 2.7).
I believe that the interpretations of the results in this paragraph (lines 253ff) are not accurate. The Bayesian posterior random flux error is independent of enhancements (see Eq. (4)), so smaller enhancements do not, as the authors state, limit the ability to reduce the random flux errors (note that the posterior flux is of course not independent of the random error or enhancements). What actually limits the ability to retrieve the anthropogenic component should be uncertainty in the biogenic fluxes, not, as the authors state, the fact that they are negative in summer. Since the authors set the uncertainty of the biogenic fluxes at 25% of the flux, it's higher in summer (according to Fig. A2), so the inversion dumps more corrections into this component then (assuming the different flux components are optimized separately, which, as mentioned above, is not clear). This should explain why the posterior anthropogenic fluxes get closer to the truth (Fig. 8) and why the random uncertainty reduction is better without bio fluxes (Fig. 9), especially in summer.
TECHNICAL CORRECTIONS
Recurring
- Please include the DOI in all references that have oneLine 10: I think it would be clearer if the structure of the sentence were switched around, i.e.:
The three-sweep observing strategy, which generally outperforms the two-sweep mode by virtue of its wider scan area that typically yields more cloud-free scenes, can retrieve the total emissions of the truth within 7% over Paris and 21% over London.Line 14: See comments on Sect. 3.3. I think it's the random error, not the fact that the biospheric signal is negative.
Line 19: Is there a better peer-reviewed reference for this statement? Bulkeley 2013 is not peer-reviewed.
Line 29: I'd remove "integrated".
Line 35: "subject to" doesn't really fit here - perhaps "but they are subject to"
Line 45: "region thereby" -> "region, thereby"
Line 52: "continuity of collecting ..." -> "continuity with other satellites collecting ..."
Line 61: I'd explicitly write out "TCCON"
Line 68: "analyses" -> "analysis"
Line 92: "still exist" -> maybe you mean "may not be predicted"?
Line 98: "is defined as"
Lines 90-106: Add somewhere in the beginning something along the lines of "We follow the approach by Palmer et al. (2011) and briefly outline the method here" - the section is copied (partly verbatim) from there.
Line 143: In the comparison of Eddy data and SMUrF results, please include the cumulative NEE, over one year or perhaps during the months analyzed in the study, alongside r and slope.
Lines 150ff: State how long back in time the backtrajectories are.
Line 189: "Realistic 20% random error": I don't find this number in a quick scan of either of the cited references. Please elaborate a bit on how you arrive at 20% to represent transport error.
Line 197: Please also provide the bias in terms of percentages of the fluxes.
Line 200: "closer to one" -> "close to one"
Line 202: Similar to line 189: Please explain briefly the reasons behind the choice of 10km as correlation length. A quick look at the cited references might not make it obvious to the reader.
Line 215: "Figure 4" -> "Figures 4a, b, e and f"
Lines 281f: Would the results be better if the position of the city scan were adjusted to the wind direction? I.e. scan further South when the wind is from the North? You could add this to the conclusions.
Line 282: See comments on Sect. 3.3. I think it's the random error, not the fact that the biospheric signal is negative.
Lines 285ff: This should be a major conclusion, alongside the call for measuring additional species.
Lines 296f: "Quantifying sub-city scale emissions requires improving the accuracy and precision of satellite-based CO2 measurements." - Does this refer to the MicroCarb skill in this study, or is it a general statement?
Fig. 2:
- Indicate the overpass direction
- Caption: "halve the cloud cover from one to 0.5" is unclear to me. Sounds like cloud cover up to 50% gets a pass, which shouldn't be the case
- Update to this comment after reading a bit further: It seems that this phrase refers to the artificial halving described in line 109. If so, the phrase is confusing in the caption of Fig. 2 and I would either leave it out or refer the reader to Sect. 2.2.Fig. 5:
- Just a comment - the figure is really helpful!Fig. A1:
- Caption: "red star marks" -> "red stars mark"Fig. A6:
- Do these figures only include anthropogenic or also biogenic fluxes?Citation: https://doi.org/10.5194/egusphere-2022-749-RC2 - AC2: 'Reply on RC2', Kai Wu, 12 Dec 2022