CN111487216B - 一种二氧化碳通量反演方法、系统 - Google Patents
一种二氧化碳通量反演方法、系统 Download PDFInfo
- Publication number
- CN111487216B CN111487216B CN202010373517.0A CN202010373517A CN111487216B CN 111487216 B CN111487216 B CN 111487216B CN 202010373517 A CN202010373517 A CN 202010373517A CN 111487216 B CN111487216 B CN 111487216B
- Authority
- CN
- China
- Prior art keywords
- carbon dioxide
- flux
- value
- rams
- assimilation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 title claims abstract description 374
- 239000001569 carbon dioxide Substances 0.000 title claims abstract description 187
- 229910002092 carbon dioxide Inorganic materials 0.000 title claims abstract description 187
- 230000004907 flux Effects 0.000 title claims abstract description 130
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000004088 simulation Methods 0.000 claims abstract description 30
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 13
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 24
- 238000002485 combustion reaction Methods 0.000 claims description 19
- 238000010672 photosynthesis Methods 0.000 claims description 18
- 230000029553 photosynthesis Effects 0.000 claims description 18
- 239000002803 fossil fuel Substances 0.000 claims description 14
- 239000002028 Biomass Substances 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 11
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 9
- 239000010908 plant waste Substances 0.000 claims description 8
- 239000000446 fuel Substances 0.000 claims description 6
- 238000012950 reanalysis Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 26
- 238000009826 distribution Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 10
- 238000011160 research Methods 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 9
- 238000005070 sampling Methods 0.000 description 8
- 241000894007 species Species 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000005855 radiation Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004177 carbon cycle Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000003344 environmental pollutant Substances 0.000 description 3
- 239000005431 greenhouse gas Substances 0.000 description 3
- 231100000719 pollutant Toxicity 0.000 description 3
- 239000000443 aerosol Substances 0.000 description 2
- 239000000926 atmospheric chemistry Substances 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 101100495256 Caenorhabditis elegans mat-3 gene Proteins 0.000 description 1
- 244000025254 Cannabis sativa Species 0.000 description 1
- 241000272165 Charadriidae Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 239000000809 air pollutant Substances 0.000 description 1
- 231100001243 air pollutant Toxicity 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000013626 chemical specie Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000006552 photochemical reaction Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种二氧化碳通量反演方法、系统，解决现有方法和系统推广性差且不适用GMI载荷的问题。所述方法包含：获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值；基于RAMS‑CMAQ模拟模型，模拟区域大气陆地碳通量交换模式；将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS‑CMAQ模型的排放清单中二氧化碳通量之和作为二氧化碳通量的先验值，采用集合四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。所述系统使用所述方法。本发明实现了应用GMI载荷计算二氧化碳通量。
Description
技术领域
本发明涉及卫星遥感领域，尤其涉及一种二氧化碳通量反演方法、系统。
背景技术
在全球碳循环研究中全球大气化学模式在研究全球CO2(二氧化碳)时空分布及传输机制模拟分析和诊断中发挥了重要作用。现有的通过资料同化技术结合大气化学传输模式和CO2地面站点观测资料，是反演地表CO2通量的一种有效方法，但由于大部分观测结果本身还只是代表观测站周边的生态系统CO2收支特征，将站点的观测结果直接外推到更大的区域尺度会导致较大的不确定性。另外，由于地基观测站点数量较稀少且在全球的分布极不均匀，在很大程度上降低了同化系统在很多区域优化CO2通量估算的可行性。同时，我国高分五号卫星搭载的GMI(大气温室气体监测仪)载荷还未用于区域二氧化碳通量的估算。
发明内容
本发明提供一种二氧化碳通量反演方法、系统，解决现有方法和系统推广性差且不适用GMI载荷的问题。
为解决上述问题，本发明是这样实现的：
第一方面，本发明实施例指出一种二氧化碳通量反演方法，包含以下步骤：获取GMI(Greenhouse-gas Measurement Instrument)数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值；基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在RAMS-CMAQ模型中引入气态物种二氧化碳，驱动场使用RAMS提供的三维气象输入资料，排放清单中的地气交换通量取自陆地生态系统光合作用和呼吸作用模型的输出结果；根据所述二氧化碳平均柱浓度观测值、二氧化碳廓线观测值、RAMS-CMAQ模型，计算得到二氧化碳平均柱浓度模拟值：
其中，
为所述二氧化碳平均柱浓度模拟值，
为所述二氧化碳平均柱浓度观测值，h0为气压权重函数，A为GMI的平均核函数，xa为所述二氧化碳廓线观测值，xh为所述RAMS-CMAQ模型模拟的垂直方向二氧化碳廓线插值到GMI分层高度处的浓度值；将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为模式模拟的二氧化碳通量的先验值，采用集合四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
优选地，所述陆地生态系统光合作用模型的输出结果为利用植被生态系统的光能利用率模型估算的陆地生态系统总初级生产力。
优选地，所述RAMS-CMAQ模型的排放清单包含化石燃料排放，森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源，二氧化碳海气交换通量和二氧化碳地气交换通量。
优选地，所述化石燃料排放采用亚洲区域排放清单(Regional Emissioninventory in Asia)0.25°×0.25°月平均排放清单，所述森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源使用全球火情排放数据库(Global Fire EmissionsDatabase)0.5°×0.5°月平均排放源清单，所述二氧化碳海气交换通量采用3°经度×3°纬度二氧化碳通量。
优选地，所述RAMS-CMAQ模型的区域气象模式驱动场取自美国NECP再分析资料。
优选地，所述GMI数据的获取周期为24小时。
优选地，所述RAMS-CMAQ模拟模型中的区域化学模式与区域大气模式的投影中心点重合、网格距相同，模拟区域为6654×5440km2，水平分辨率64km×64km，水平方向网格点为105×86，垂直方向采用地形追随坐标，从地表至21km不等距分为15层，最低7层高度与RAMS重合。
第二方面，本发明实施例还指出一种二氧化碳通量反演系统，使用所述方法，所述系统包含：卫星数据处理模块、大气处理模块、同化处理模块；所述卫星数据处理模块，用于获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值；所述大气处理模块，用于基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在RAMS-CMAQ模型中引入气态物种二氧化碳，驱动场使用RAMS提供的三维气象输入资料，排放清单中的地气交换通量取自陆地生态系统光合作用和呼吸作用模型的输出结果；所述同化处理模块，用于根据所述二氧化碳平均柱浓度观测值、二氧化碳廓线观测值、RAMS-CMAQ模型，计算得到二氧化碳柱平均浓度模拟值；还用于将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为二氧化碳通量的先验值，采用集合四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
本发明有益效果包括：本发明适用于我国高分五号卫星搭载的GMI载荷的二氧化碳通量反演方法，在CMAQ模式中引入并耦合了陆地生态系统光合作用和呼吸作用的二氧化碳模块，使该模式成为污染气体和碳循环为一个大气的处理系统，从而得到的二氧化碳通量更准确、提高了反演精度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为一种二氧化碳通量反演方法流程实施例；
图2为一种二氧化碳通量反演系统实施例。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在全球碳循环研究中全球大气化学模式在研究全球CO2时空分布及传输机制模拟分析和诊断中发挥了重要作用。近年来由于城市化石燃料排放研究，需要提高区域碳收支的认识和验证，对排放源的估算对更精细尺度上CO2的传输和变化的模拟需求也越来越迫切。由于排放源分布与地表性质的高度不均匀性，CO2浓度差异主要表现在大气边界层内，所以如何获得地面CO2通量分布就成了大气CO2“净排放”估算的核心问题。
目前研究成果已证明使用区域模式模拟CO2的可行性并发现了高分辨率区域CTM(chemical transport models)在模拟方面的优势。区域模式研究表明有各种因素影响CO2时空分布的，包括地形、地气交换通量的日变化及水平分布不均一性、输送过程与CO2通量之间的相关性等。如将Weather Research and Forecast Model(WRF)和生态系统光合与呼吸模型(Vegetation Photosynthesis and Respiration Model，VPRM)耦合。用WRF模式和STILT(拉格朗日型)耦合了VPRM模型模拟高山和海岸复杂地形下垫面CO2浓度的分布。在WRF-Chem中引入CO2气态物种(作为示踪物传输)，并加入CO2人为源、生物源及海洋源，考察不同水平分辨率的CO2通量在不同时间尺度上(季节、天气、日变化)的模拟情况。对比全球模式GEOS-Chem模拟和高塔观测数据(北京站与合肥站)表明合肥站的CO2浓度日变化大于北京站。
近年来国内很多研究试图通过实测陆地碳汇的技术途径来确定CO2地气交换，在典型区域的陆地生态系统开展了CO2通量长期观测与试验工作，获得了高时间分辨率的CO2通量和环境变化数据，为开展不同时间尺度的CO2通量变化及其环境响应机理研究提供了方便。
通过资料同化技术结合大气化学传输模式和CO2地面站点观测资料，是反演地表CO2通量的一种有效方法。由于卫星资料相比地基台站观测具有更高的空间覆盖率，星载CO2观测有热红外TIR(thermal infrared wavelengths)和短波近红外SWIR(short-waveinfrared wavelengths)两种探测方式，热红外遥感对于大气低层的CO2变化特征的表征不确定性比较大。日本的GOSAT(Greenhourse gases Observing SATellite)、NASA的OCO-2(Orbiting Carbon Observatory-2)和我国的高分五号卫星GMI载荷，分别已于2008年1月23日、2014年7月2日和2018年5月9日成功发射，其搭载的短波近红外探测器接收由地表反射的短波太阳辐射信息。但卫星资料需要严格控制反演误差、大气模式传输误差以及CO2先验通量的误差。
以往的许多研究已经通过资料同化方法给出了更精确的CO2通量估算。但由于大部分观测结果本身还只是代表观测站周边的生态系统CO2收支特征，将站点的观测结果直接外推到更大的区域尺度会导致较大的不确定性。另外，由于地基观测站点数量较稀少且在全球的分布极不均匀，在很大程度上降低了同化系统在很多区域优化CO2通量估算的可行性。由于卫星资料相比地基台站观测具有更高的空间覆盖率，研究表明这些观测资料对于区域CO2地表源汇的优化提供新的支持。
日本的GOSAT、NASA的OCO-2和我国的高分五号卫星GMI载荷，其搭载的短波近红外探测器接收由地表反射的短波太阳辐射信息。本提案采用我国自己研制的高分五号卫星的GMI载荷的二氧化碳数据。但卫星资料需要严格控制反演误差、大气模式传输误差以及CO2先验通量的误差。到目前为止，还没有利用我国高分五号卫星GMI载荷探测资料来同化，并同化二氧化碳通量的研究。
本发明创新点如下：第一，本发明针对我国高分五号卫星温室气体探测载荷GMI，首次实现了基于GMI载荷的二氧化碳柱浓度反演方法和系统，可提高GMI载荷的工程应用范围；第二，本发明在CMAQ模式中引入并耦合了陆地生态系统光合作用和呼吸作用的二氧化碳模块，使该模式成为污染气体和碳循环为一个大气的处理系统，完善了CMAQ模型；第三，本发明基于卫星观测误差与反演误差、RAMS-CMAQ模式模拟误差的同化技术，开展二氧化碳浓度与通量的同时同化反演，反演精度更高。
以下结合附图，详细说明本发明各实施例提供的技术方案。
图1为一种二氧化碳通量反演方法流程实施例，通过高分五号卫星的GMI载荷数据实现二氧化碳通量反演，作为本发明实施例，一种二氧化碳通量反演方法，具体包含以下步骤：
步骤101，获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值。
在步骤101中，获取高分五号卫星GMI载荷数据，GMI载荷的光谱通道分别是：0.759-0.769μm、1.568-1.583μm、1.642-1.658μm和2.043-2.058μm，其中用于CO2反演的两个通道1.642-1.658μm和2.043-2.058μm的光谱分辨率为0.27波数，即利用0.759-0.769μm和1.568-1.583μm波段数据反演二氧化碳柱浓度，所用的反演方法是基于最优化理论的非线性迭代反演方法，为修正的阻尼牛顿方法。
在步骤101中，所述二氧化碳廓线观测值为模式分层内插得到的数据。
需要说明的是，二氧化碳柱浓度反演方法为现有技术，这里不具体说明。
在步骤101中，所述GMI载荷数据的获取周期为24小时。
步骤102，基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在RAMS-CMAQ模型中引入气态物种二氧化碳，驱动场使用RAMS提供的三维气象输入资料，排放清单中地气交换通量取自陆地生态系统光合作用和呼吸作用模型的输出结果。
在步骤102中，基于现有的RAMS-CMAQ模型耦合陆地生态系统光合作用模型和呼吸作用模型，模拟了区域大气陆地碳通量交换模式。
需要说明的是，RAMS-CMAQ由区域大气模式RAMS(Regional AtmosphericModeling System)和区域大气化学模式CMAQ(Models-3Community Multiscale AirQuality)构成，可模拟多种污染物的浓度分布情况。RAMS-CMAQ是基于美国EPA开发的第三代多尺度空气质量模式CMAQ并由区域大气模式RAMS提供三维气象场数据。除了气象模块，CMAQ包含的主要模块有：排放源处理模块，用于对自然和人为源排放的读入进行处理；化学传输模块(Chemical-Transport Modeling System,CCTM)，对污染物在大气中的传输、演变进行模拟，包括的机制有水平扩散、垂直对流、气相化学反应机制(包括光化学反应速率算法)、液相化学反应机制(包括云参数化算法)、气溶胶动力学、气体和颗粒物沉降算法等。
在步骤102中，在RAMS-CMAQ模型中引入气态物种二氧化碳：在CMAQ模式系统中的化学传输模块使用碳键机制CB05(Carbon Bond mechanism)，化学传输模块中含有化学物种52种，核心反应156个。CMAQ模块中引入气态物种CO2并保留模式中的其他物种，可以实现CO2和所有传统监测的空气污染物的同时模拟和监测。需要指出的是，RAMS-CMAQ的结合使用现有技术，这里不再赘述。
在步骤102中，RAMS-CMAQ模拟的驱动场使用RAMS提供的三维气象输入资料，考虑到大气CO2区别于其他污染物的一些独特性质，如在大气中具有较长的寿命，较高的背景浓度，较强的地气交换影响，在运用CMAQ区域模式模拟CO2时必须注意一些必要条件，包括来自全球模式的真实初始场和侧边界条件，其大尺度天气背景场取自美国NCEP 6小时间隔的1°×1°再分析资料。
CO2不参与化学反应，但其浓度受地气交换通量、海气交换通量、化石燃料排放和生物质燃烧等四类源的影响。在本发明实施例中，大气CO2的地表源汇主要包括化石燃料排放、地气交换、海气交换和生物质燃烧，其中，化石燃料和生物质燃烧部分排放的CO2通量来自现有的排放源清单。本发明实施例采用真实观测资料同化，所述CMAQ模型的排放清单包含：化石燃料排放，森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源，二氧化碳海气交换和二氧化碳地气交换。
进一步地，所述化石燃料排放，采用亚洲区域排放清单(Regional Emissioninventory in Asia)2015年亚洲0.25°×0.25°月平均排放清单；所述森林大火、草原燃烧及农作物废弃物燃烧释放的生物质燃烧排放源使用全球火点排放数据库(Global FireEmissions Database)0.5°×0.5°月平均排放源清单；所述二氧化碳地气交换清单的网格为经度3°×纬度2°，3小时步长的二氧化碳通量。
需要说明的是，所述RAMS-CMAQ模型的排放清单可以是本发明实施例中的排放清单，也可以是其他排放清单，这里不做特别限定。
在步骤102中，所述RAMS-CMAQ模拟模型中的CMAQ投影中心点与RAMS重合(例如，35.0°N,116.0°E)，网格距与RAMS相同，模拟区域为6654×5440km2，水平分辨率64km×64km，水平方向网格点为105×86，覆盖整个东亚地区；垂直方向采用地形追随坐标，从地表至大约21km不等距分为15层，最低7层高度与RAMS重合，且有近一半分布于2km以下，以便更好地模拟大气边界层结构。
在步骤102中，所述地气交换通量取自陆地生态系统光合作用模型的输出结果，所述陆地生态系统光合作用模型的输出结果为利用植被生态系统的光能利用率模型估算的陆地生态系统总初级生产力。
具体地，近地面CO2浓度是CO2地气交换的结果，是指陆地生态系统和大气之间的CO2通量交换形成的浓度分布。利用植被生态系统的光能利用率模型来估算陆地生态系统的总初级生产力GPP(Gross Primary Production)：
GPP＝εn×FPAR×QPAR (1)
其中，这里εn是光能利用率，FPAR是植被冠层吸收的光合有效辐射比率(theFraction of Photosynthetically Active Radiution Absorbed by vegetation,FPAR)，QPAR是光合有效辐射(Photosynthetically Active Radiation,PAR)，是波长400～700nm范围内的地表太阳有效辐射，在本申请方案中主要从区域气候模式RAMS的模拟中获得。
公式(1)中的FPAR可以从归一化植被指数(Normalized Difference VegetationIndex,NDVI),或增强植被指数(Enhanced Vegetation Index,EVI)获得，即FPAR＝EVI。
其中，ρblue、ρred和ρnir分别是MODIS蓝光(459～479nm)、红光(620～670nm)和近红外(841～875nm)波段的地表反射率。
公式(1)中的εn可以用下式表示：
εn＝ε0×f(T)×f(W) (3)
其中，ε0是最佳环境下的最大光能利用率，f(T)和f(W)是光合作用中的大气温度与植被水分限制函数，ε0是重要参数，与不同植被类型有关，植被类型的分布从现有的卫星遥感的植被分类图中获取，式中的T是日平均温度,光合作用的温度限制函数f(T)：
其中，Topt是植被生长的最佳温度25℃，Tmax是植被光合作用的最高温度,一般为40℃。
公式(3)中f(W)是水分对光合作用的限制函数，f(W)可以利用卫星探测的植被冠层水分指数来计算为：
其中，LSWImax是植被生长期间的最大水分指数，LSWI是植被生长期间的水分指数，可以从下式计算获得：
其中，ρnir和ρswir是卫星数据MODIS传感器近红外841～875nm和1628～1652nm波段的反射率。
步骤103，根据所述RAMS-CMAQ模型、二氧化碳平均柱浓度观测值、二氧化碳廓线观测值，计算得到二氧化碳平均柱浓度模拟值：
其中，
为所述二氧化碳平均柱浓度模拟值，
为所述二氧化碳平均柱浓度观测值，为GMI卫星反演过程中使用的CO2平均柱浓度先验值，h0为气压权重函数，A为GMI的平均核函数，xa为所述二氧化碳廓线观测值，xh为所述RAMS-CMAQ模型模拟的垂直方向二氧化碳廓线插值到GMI分层高度处的浓度值。
需要说明的是，A可以从卫星反演结果中读取，为卫星监测结果的平均核函数，xa为二氧化碳反演结果按模式分层要求内插的廓线观测值。
在步骤103中，卫星资料同化的关键技术是如何构造观测算子，即如何将CMAQ模式输出的CO2浓度与卫星资料观测的CO2浓度转换到同一个尺度，本发明实施例利用Rodgers提供的方法构造观测算子，具体公式为公式(7)。
步骤104，将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为二氧化碳通量的先验值，采用集合四维变分同化(POD-4DVar)方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
在步骤104中，所述二氧化碳通量的先验值为所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和，例如，化石燃料排放，森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源，二氧化碳海气交换通量和二氧化碳地气交换通量这四种通量之和。
在步骤104中，集合四维变分同化算法为现有方法，具体可参考文献“xiangjunTian,Zhenghui Xie and Qin Sun.A POD-based ensemble four-dimensionalvariational assimilation method，Tellus(2011),63A,805–816”。
在步骤104中，如何结合集合四维变分同化方法计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值为本发明创新内容，本发明将所述二氧化碳平均柱浓度观测值作为观测值，所述二氧化碳平均柱浓度模拟值作为观测算子，利用两次CMAQ模式产生二氧化碳浓度背景场与采样场，得到二氧化碳浓度背景场的微小增量、先验通量线性比例因子的微小增量、二氧化碳柱浓度模拟观测值增量以及二氧化碳柱浓度观测值增量，基于POD-4DVar核心同化算法，计算得到二氧化碳浓度与先验通量的优化值。
具体地，集合四维变分同化算法中，在一个同化的时间窗口内，通过构建代价函数描述状态量分析值和真值之间的差异，把数据同化问题转换为一个极值求解问题。在满足动态约束的条件下，最小化状态预测值和观测值之间的“距离”，使得这种“距离”最小的状态量即为最优状态估计量。
所构建的代价函数为：
其中，x是模式初始状态值，即需要极小化的控制变量，xb为背景场变量，B为背景误差协方差矩阵，H为观测算子，R为观测误差协方差矩阵，上标T表示转置，上标b为背景场，S为同化窗口中的总观测时次，xk为在强约束条件下从初始时刻t0到积分时刻tk对模式积分得到，
这里M为模式预报算子，yk为第k个观测值。
2011年基于蒙特卡罗方法和本征正交分解技术(Monte Carlo method and theproper orthogonal decomposition，POD)提出了集合四维变分同化算法POD-4DVar算法，该算法通过极小化4DVar的增量代价函数，得到起始条件的最优增量，增量代价函数具体可参考文献“xiangjun Tian,Zhenghui Xie and Qin Sun.A POD-based ensemble four-dimensional variational assimilation method，Tellus(2011),63A,805–816”，即为：
其中，x'＝x-xb是初始时刻背景场的微小增量，y'(x')模拟的观测增量，y'obs为观测增量，xa'为控制变量x相对于背景场xb的微小增量x'优化后的结果。可以根据POD-4DVar算法可求解出分析值的增量xa'，根据x'的定义，可求出控制变量x优化后的最终分析解为x＝xa'+xb。
对于浓度同化，微小增量x'是由两次运行CMAQ模式产生，分别是CO2浓度背景模拟和采样模拟。在背景模拟中，背景二氧化碳浓度xb由CMAQ在同化窗口内积分模拟产生。采样模拟中，采样场二氧化碳浓度xm由CMAQ在采样窗口内积分产生。由此可求得控制变量x的微小增量为xm-xb。由步骤103中所述观测算子，求解的二氧化碳背景场对应的柱浓度模拟值为H(xb)，二氧化碳采样场对应柱浓度模拟值为H(xm)，据此可求出模拟的二氧化碳柱浓度增量y'(x')＝H(xm)-H(xb)。对应的观测增量y′obs可由步骤101中所述GMI反演方法所获取的二氧化碳柱浓度值yobs(所述二氧化碳平均柱浓度观测值)与二氧化碳背景场柱浓度值模拟求得，即y'obs＝yobs-H(xb)。
利用公式(9)求解出分析值的增量，即可求出优化后的二氧化碳浓度分析值为：xa'+xb，即为所述二氧化碳柱浓度同化值。
对于通量同化，在步骤104中，同化系统假设地表CO2先验通量为F*(x,y,t)，本研究通过引入CO2观测资料来优化先验通量F*(x,y,t)的估算，通过同化表征通量的线性比例因子来实现通量的同化。
F(x,y,t)＝λ(x,y,t)×F*(x,y,t) (10)
其中，F(x,y,t)为所述二氧化碳通量同化值，F*(x,y,t)为所述地表CO2先验通量，λ(x,y,t)为先验通量线性比例因子。
与浓度同化不同，通量同化中，待优化的控制变量x是地表二氧化碳先验通量的线性比例因子λ(x,y,t)，xb为线性比例因子的背景λb，λm为通量线性比例因子的采样场。因此，控制变量x的微小增量为λm-λb。将此微小增量，与浓度同化中所求得的模拟的二氧化碳柱浓度增量y'(x')、及对应的观测增量y′obs，带入公式(9)求解出分析值的增量λa′，即可求出优化后的线性比例因子分析值λa为λa′+λb，据此可计算出优化后的通量，即所述二氧化碳通量同化值。
将优化后的二氧化碳浓度分析值作为二氧化碳浓度的初始场，与优化后的二氧化碳通量，同时驱动预报模式CMAQ，预报得到下一时刻的二氧化碳浓度值，作为下一个同化窗口的初始状态值，进入下个同化窗口进行同化。
本发明实施例提供了一种二氧化碳通量反演方法，实现了通过GMI载荷数据对近地面二氧化碳通量的反演，结合RAMS-CMAQ模型，引入了大气光合作用效应的影响，并采用同化方法，使得二氧化碳浓度和通量的反演结果更准确。
图2为一种二氧化碳通量反演系统实施例，使用了二氧化碳通量反演方法，作为本发明实施例，一种二氧化碳通量反演系统，包含：卫星数据处理模块11、大气处理模块12、同化处理模块13。
所述卫星数据处理模块，用于获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值。
所述大气处理模块，用于基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在CMAQ模型中引入气态物种二氧化碳，CMAQ模型的驱动场使用RAMS提供的三维气象输入资料，CMAQ模型的排放清单中陆地生态系统和大气的地气交换通量取自陆地生态系统光合作用模型的输出结果。
所述同化处理模块，用于，根据所述二氧化碳平均柱浓度观测值、二氧化碳廓线观测值、RAMS-CMAQ模型，计算得到二氧化碳柱平均浓度模拟值；
还用于，将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为二氧化碳通量的先验值，采用四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
在本发明实施例中，所述同化处理模块，利用两次CMAQ模式产生二氧化碳浓度背景场与采样场，得到二氧化碳浓度背景场的微小增量、先验通量线性比例因子的微小增量、二氧化碳柱浓度模拟观测值增量以及二氧化碳柱浓度观测值增量，基于POD-4DVar核心同化算法，计算得到二氧化碳浓度与先验通量的优化值。
本发明实施例中，各模块进行处理的具体方法在第1实施例中已经详细论述，这里不重复论述。
本发明实施例的观测数据为2019年1月，现有源清单中包含化石燃料排放，森林大火、草原燃烧及农作物废弃物燃烧释放的生物质燃烧排放源，CO2海气交换和地气交换，将该四类通量之和作为CO2通量的先验值F1 *，以GMI提供的CO2平均柱浓度值作为观测资料，设置同化系统参数，集合数N＝145，局地化半径R＝320km。考虑到卫星资料的时间密度相对地基资料较低，将优化窗口和观测窗口设为24h，以保证每次同化中由适量的观测资料进入。
为了对同化结果进行检验，进行一组由先验通量驱动的CMAQ模拟的对照组实验，给出同化结果和对照模拟结果。以初步检验该同化系统对卫星资料的同化效果，探讨该区域同化系统在实际应用方面的潜力。
在卫星资料反演过程中，为控制数据质量大多数探测(如太阳天顶角大于85°，受云和气溶胶等影响不满足晴空条件，探测点下垫面为海洋但未开启glint模态等)被排除，每条卫星轨道上反演数据产品中质量为“Good”大约仅包含10～100个探测，对应到本研究中模式的东亚区域内每条卫星轨道上大约有0～60个探测符合反演条件。
在本发明实施例中，首先通过分析CO2浓度的同化值来探讨该区域同化系统在加入卫星资料后的同化效果，此外也进行了一组由先验通量驱动的CMAQ模拟，从而与同化结果进行对比。
本发明实施例提供了通过二氧化碳通量反演方法计算得到的二氧化碳浓度同化值，说明了可以同时同化二氧化碳浓度值和通量值。
可以看到，加入CO2柱浓度信息对与地表CO2通量的改进效果是非常明显。特别是在日本、中国珠三角以及京津冀地区等人为活动聚集的区域，二者之间的差异在0.05～0.2μmole m-2s-1之间，表明了先验排放源清单中可能存在对CO2通量的低估，由于这些地区在冬季主要受化石燃料排放的影响，那么这里估算的不确定性可能主要由化石燃料排放清单引起。因此，相比以往在基于生态系统尺度上的仅关注CO2地气和海气通量同化，该区域同化系统考虑在格点尺度上将CO2通量作为一个整体在模式网格内估算的同化方案表现出了一些优势。同化前后CO2通量差异与近地面CO2浓度具有类似的水平分布特征。
需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。
Claims (8)
1.一种二氧化碳通量反演方法，其特征在于，包含以下步骤：
获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值；
基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在RAMS-CMAQ模型中引入气态物种二氧化碳，驱动场使用RAMS提供的三维气象输入资料，排放清单中地气交换通量取自陆地生态系统光合作用和呼吸作用模型的输出结果；
根据所述二氧化碳平均柱浓度观测值、二氧化碳廓线观测值、RAMS-CMAQ模型，计算得到二氧化碳平均柱浓度模拟值：
其中，
为所述二氧化碳平均柱浓度模拟值，
为所述二氧化碳平均柱浓度观测值，h0为气压权重函数，A为GMI的平均核函数，xa为所述二氧化碳廓线观测值，xh为所述RAMS-CMAQ模型模拟的垂直方向二氧化碳廓线插值到GMI分层高度处的浓度值；
将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为模式模拟的二氧化碳通量的先验值，采用集合四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
2.如权利要求1所述的二氧化碳通量反演方法，其特征在于，陆地生态系统光合作用模型的输出结果为利用植被生态系统的光能利用率模型估算的陆地生态系统总初级生产力。
3.如权利要求1所述的二氧化碳通量反演方法，其特征在于，所述RAMS-CMAQ模型的排放清单包含化石燃料排放，森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源，二氧化碳海气交换通量和二氧化碳地气交换通量。
4.如权利要求1所述的二氧化碳通量反演方法，其特征在于，所述RAMS-CMAQ模型的区域气象模式驱动场取自美国NECP再分析资料。
5.如权利要求1所述的二氧化碳通量反演方法，其特征在于，所述GMI数据的获取周期为24小时。
6.如权利要求1所述的二氧化碳通量反演方法，其特征在于，所述RAMS-CMAQ模型中的区域化学模式与区域大气模式的投影中心点重合、网格距相同，设立的模拟区域为6654×5440km2，水平分辨率64km×64km，水平方向网格点为105×86，垂直方向采用地形追随坐标，从地表至21km不等距分为15层，最低7层高度与RAMS重合。
7.如权利要求3所述的二氧化碳通量反演方法，其特征在于，所述化石燃料排放采用亚洲区域排放清单0.25°×0.25°月平均排放清单，所述森林大火、草原燃料及农作物废弃物燃烧释放的生物质燃烧排放源使用全球火情排放数据库0.5°×0.5°月平均排放源清单，所述二氧化碳海气交换通量采用3°经度×3°纬度二氧化碳通量。
8.一种二氧化碳通量反演系统，使用权利要求1～7任一项所述方法，其特征在于，所述系统包含：卫星数据处理模块、大气处理模块、同化处理模块；
所述卫星数据处理模块，用于获取GMI数据，进行二氧化碳柱浓度反演、得到二氧化碳平均柱浓度观测值，进行高度内插、得到二氧化碳廓线观测值；
所述大气处理模块，用于基于RAMS-CMAQ模型，模拟区域大气陆地碳通量交换模式：在CMAQ模型中引入气态物种二氧化碳，CMAQ模型的驱动场使用RAMS提供的三维气象输入资料，CMAQ模型的排放清单中陆地生态系统和大气的地气交换通量取自陆地生态系统光合作用和呼吸作用模型的输出结果；
所述同化处理模块，用于，
根据所述二氧化碳平均柱浓度观测值、二氧化碳廓线观测值、RAMS-CMAQ模型，计算得到二氧化碳柱平均浓度模拟值；
将所述二氧化碳平均柱浓度模拟值作为观测算子，将所述RAMS-CMAQ模型的排放清单中二氧化碳通量之和作为二氧化碳通量的先验值，采用集合四维变分同化方法，计算得到二氧化碳通量同化值和二氧化碳柱浓度同化值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010373517.0A CN111487216B (zh) | 2020-05-06 | 2020-05-06 | 一种二氧化碳通量反演方法、系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010373517.0A CN111487216B (zh) | 2020-05-06 | 2020-05-06 | 一种二氧化碳通量反演方法、系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111487216A CN111487216A (zh) | 2020-08-04 |
| CN111487216B true CN111487216B (zh) | 2023-01-24 |
Family
ID=71795433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010373517.0A Active CN111487216B (zh) | 2020-05-06 | 2020-05-06 | 一种二氧化碳通量反演方法、系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111487216B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112161943B (zh) * | 2020-09-29 | 2021-04-16 | 中国科学院地理科学与资源研究所 | TanSat卫星XCO2反演数据的偏差校正方法和系统 |
| CN112949085B (zh) * | 2021-03-25 | 2024-06-07 | 中国空间技术研究院 | 基于天基co2通量数据的陆地碳循环同化方法及系统 |
| CN113834902B (zh) * | 2021-08-16 | 2024-06-14 | 中国人民解放军国防科技大学 | 一种基于四维变分同化的二氧化硫排放源反演方法 |
| CN114357894B (zh) * | 2022-03-11 | 2022-07-12 | 阿里巴巴达摩院(杭州)科技有限公司 | 大气污染物的处理方法、存储介质以及计算机终端 |
| CN115565063B (zh) * | 2022-03-24 | 2023-04-11 | 中国矿业大学(北京) | 基于气候潜在补偿的矿区植被碳汇贡献量计算与分析方法 |
| CN116205022A (zh) * | 2022-06-07 | 2023-06-02 | 中国科学院地理科学与资源研究所 | 基于co2/co比率反演核算高分辨率人为co2排放的方法及系统 |
| CN114970184B (zh) * | 2022-06-07 | 2024-04-02 | 中国科学院地理科学与资源研究所 | 同步反演高分辨率人为co2排放与自然co2通量的同化方法及系统 |
| CN114724647B (zh) * | 2022-06-09 | 2023-01-31 | 南京大学 | 一种高时空分辨率co2通量反演系统及方法 |
| CN117217027B (zh) * | 2023-11-02 | 2024-01-30 | 中国人民解放军国防科技大学 | 基于四维变分同化的污染物点源廓线排放估算方法和装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107273995A (zh) * | 2016-04-08 | 2017-10-20 | 株式会社日立制作所 | 空气质量预报方法 |
| CN109858686A (zh) * | 2019-01-14 | 2019-06-07 | 南京大学 | 一种基于EnKF的地面排放清单反演优化方法 |
-
2020
- 2020-05-06 CN CN202010373517.0A patent/CN111487216B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107273995A (zh) * | 2016-04-08 | 2017-10-20 | 株式会社日立制作所 | 空气质量预报方法 |
| CN109858686A (zh) * | 2019-01-14 | 2019-06-07 | 南京大学 | 一种基于EnKF的地面排放清单反演优化方法 |
Non-Patent Citations (2)
| Title |
|---|
| Simulation of Non-Homogeneous CO2 and its impact on regional temperature in east asia;xiaodong xie et.al;《 Special collection on aerosol-cloud-radiation interactions》;20180709;全文 * |
| 基于多源卫星遥感的高分辨率全球碳同化系统研究;居为民等;《地球科学进展》;20161130;第31卷(第11期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111487216A (zh) | 2020-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111487216B (zh) | 一种二氧化碳通量反演方法、系统 | |
| Ryu et al. | MODIS-derived global land products of shortwave radiation and diffuse and total photosynthetically active radiation at 5 km resolution from 2000 | |
| Ye et al. | Constraining fossil fuel CO2 emissions from urban area using OCO‐2 observations of total column CO2 | |
| Deng et al. | Inferring regional sources and sinks of atmospheric CO 2 from GOSAT XCO 2 data | |
| Monteil et al. | The regional EUROpean atmospheric transport inversion COMparison, EUROCOM: first results on European wide terrestrial carbon fluxes for the period 2006–2015 | |
| Schuh et al. | Evaluating atmospheric CO2 inversions at multiple scales over a highly inventoried agricultural landscape | |
| McCallum et al. | Satellite-based terrestrial production efficiency modeling | |
| CN106372730B (zh) | 利用机器学习的植被净初级生产力遥感估算方法 | |
| CN112884079A (zh) | 一种基于Stacking集成模型的近地面二氧化氮浓度估算方法 | |
| Beringer et al. | SPECIAL—Savanna patterns of energy and carbon integrated across the landscape | |
| Palmer et al. | A measurement-based verification framework for UK greenhouse gas emissions: an overview of the Greenhouse gAs Uk and Global Emissions (GAUGE) project | |
| Matross et al. | Estimating regional carbon exchange in New England and Quebec by combining atmospheric, ground-based and satellite data | |
| CN111579504A (zh) | 基于光学遥感的大气污染成分垂直分布反演方法 | |
| CN116050612A (zh) | 基于多技术融合的大气温室气体监测站点选址方法及系统及存储介质 | |
| Lian et al. | Sensitivity to the sources of uncertainties in the modeling of atmospheric CO 2 concentration within and in the vicinity of Paris | |
| CN114970184B (zh) | 同步反演高分辨率人为co2排放与自然co2通量的同化方法及系统 | |
| CN115186437B (zh) | 碳同位素联合同化模型及甄别人为碳排放与自然碳通量区域同化系统的构建方法 | |
| Messerschmidt et al. | Evaluation of seasonal atmosphere–biosphere exchange estimations with TCCON measurements | |
| CN114724647B (zh) | 一种高时空分辨率co2通量反演系统及方法 | |
| Kaminski et al. | Assessing the Impact of Atmospheric CO2 and NO2 Measurements From Space on Estimating City-Scale Fossil Fuel CO2 Emissions in a Data Assimilation System | |
| Zhao et al. | Understanding greenhouse gas (GHG) column concentrations in Munich using WRF | |
| Wang et al. | Error correlation between CO 2 and CO as constraint for CO 2 flux inversions using satellite data | |
| SAITO et al. | Overview of model systems for global carbon dioxide and methane flux estimates using GOSAT and GOSAT-2 observations | |
| Maselli et al. | Validating an integrated strategy to model net land carbon exchange against aircraft flux measurements | |
| CN115876948A (zh) | 基于卫星柱浓度和4d-letkf混合同化算法的碳卫星同化系统及其构建方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |