阿拉斯加的 39 个站点和加拿大西北部的 12 个站点)的季节性沉降、活动层厚度 (ALT)等数据

举报
此星光明 发表于 2024/03/28 09:43:04 2024/03/28
【摘要】 ​ ABoVE: Active Layer Thickness from Airborne L- and P- band SAR, Alaska, 2017, Ver. 3简介文件修订日期:2022-09-15数据集版本: 3摘要本数据集以 30 米的分辨率提供了整个 ABoVE 域 51 个站点(包括阿拉斯加的 39 个站点和加拿大西北部的 12 个站点)的季节性沉降、活动层厚度 (ALT...

 ABoVE: Active Layer Thickness from Airborne L- and P- band SAR, Alaska, 2017, Ver. 3

简介

文件修订日期:2022-09-15

数据集版本: 3

摘要

本数据集以 30 米的分辨率提供了整个 ABoVE 域 51 个站点(包括阿拉斯加的 39 个站点和加拿大西北部的 12 个站点)的季节性沉降、活动层厚度 (ALT)、垂直土壤水分剖面和不确定性的估计值。ALT和土壤水分剖面检索同时使用L波段和P波段合成孔径雷达(SAR)数据,这些数据由NASA/JPL无人驾驶飞行器合成孔径雷达(UAVSAR)仪器在2017年北极极地脆弱性实验(ABoVE)机载活动中获取。数据以 NetCDF 第 4 版格式提供,并附有一个用于根据数据估算土壤体积含水量的 python 脚本。
该产品由永冻土动力学观测站(PDO)项目创建,用于估算 2017 年 6 月和 9 月获取的 L 波段干涉合成孔径雷达(InSAR)对中由于活动层融化而导致的季节性下沉。它还通过机载下冠层和下表面微波观测站(AirMOSS)在 8 月份获取的 P 波段极坐标合成孔径雷达(PolSAR)反向散射来估算土壤体积含水量(VWC)的垂直剖面。联合检索同时使用 L 波段和 P 波段反向散射得出的季节性沉降来估算 ALT 和垂直土壤水分剖面以及不确定性。

北极-北方脆弱性实验(ABoVE)是美国国家航空航天局(NASA)陆地生态计划的一项实地活动,从 2015 年开始在阿拉斯加和加拿大西部进行,为期 8 到 10 年。ABoVE 的研究将基于实地的过程级研究与机载和卫星传感器获得的地理空间数据产品联系起来,为提高分析和建模能力奠定了基础,这些能力是了解和预测北极和北方地区生态系统对气候变化的反应及其对社会的影响所必需的。

空间覆盖范围:整个 ABoVE 域的 51 个站点,包括阿拉斯加的 39 个站点和加拿大的 12 个站点

ABoVE 参考地点

域:ABoVE 核心区域和扩展区域

州/地区:阿拉斯加州、育空地区、西北地区

网格单元:见表 1

空间分辨率30 m

时间覆盖范围:2017-06-19 至 2017-09-16

时间分辨率:估算值为一次性估算值。地表变形估算值代表 2017 年解冻季节发生的沉降。解冻深度和土壤水分估算值被认为代表最大解冻时的活动层土壤剖面,在这种情况下,解冻深度相当于 2017 年的 ALT。

表 1.研究区域。纬度和经度以十进制度表示。

Site Site Code UAVSAR Flight Lines ABoVE Grid Cells North Latitude South Latitude East Longitude West Longitude
Alaska
Ambler ambler ambler Bh005v002, Bh006v002 67.25 66.33 -157.82 -159.41
Anaktuvuk anaktu 70.27 68.38 -148.04 -153.49
Atqasuk atqasu atqasu Bh007v001, Bh008v001 71.05 69.67 -156.79 -158.30
Barrow barrow barrow Bh008v001 71.49 70.25 -154.44 -156.86
Bonanza Creek bonanz 65.63 64.40 -146.12 -149.14
Chevak chevak 62.00 61.07 -164.63 -166.27
Coldfoot coldfo coldfo Bh007v003, Bh006v004, Bh007v004 67.43 66.20 -149.67 -151.12
Council counci counci Bh003v002, Bh004v002, Bh004v001 65.55 64.51 -163.25 -165.15
Delta Junction deltaj deltjA, deltjB, deltjC, djNEON Bh006v005, Bh006v006, Bh006v007, Bh007v006, Bh007v007 64.50 62.85 -142.08 -147.61
Denali denali 64.70 62.97 -148.36 -152.40
Deadhorse dhorse 70.52 68.69 -145.97 -151.23
Fort Smith N fsmitN 61.66 59.77 -111.92 -113.91
Fort Smith S fsmitS 60.27 59.50 -111.20 -112.28
ftreso 61.75 60.91 -112.37 -114.33
Good Hope gdhope 67.98 65.98 -128.20 -131.92
Huslia huslia huslia Bh005v003 65.54 65.18 -154.82 -156.80
Inigok inigok inigok Bh007v002, Bh008v002 70.68 69.71 -151.23 -154.20
Ivotuk ivotuk 69.02 67.91 -154.52 -157.26
Katmai National Park katmai katmaA, katmaB Bh001v006, Bh002v006 58.49 58.12 -154.58 -157.06
Kluane A kluanA 61.58 60.24 -136.78 -139.26
Kluane B kluanB 61.77 60.70 -137.21 -139.58
Kougarok kougar kougar Bh004v002, Bh005v002 65.71 65.45 -159.96 -163.03
Koyuk koyukk koyukk Bh004v002, Bh004v003 65.06 64.71 -159.19 -161.12
Lake Clark lclark 61.44 59.70 -152.19 -155.72
Mcpherson mcpher 67.74 66.85 -133.81 -136.87
Noatak noatak 68.54 67.34 -159.39 -163.10
Norman Wells nwells 66.43 64.26 -124.92 -129.59
Poorman poorma poorma Bh004v004, Bh005v004 64.42 64.26 -153.12 -156.23
Scoaoi scoaoi 60.91 59.91 -119.78 -121.83
Scotty scotty 61.78 60.44 -120.21 -122.00
SnagyK snagyk 63.06 62.00 -139.88 -142.01
Snake River sriver sriver Bh003v004, Bh003v005 61.68 60.77 -156.05 -157.48
Teller teller 65.98 64.25 -163.37 -166.73
Toolik toolik toolik Bh008v003 68.84 68.51 -148.97 -150.03
Watson watson 60.84 59.22 -127.26 -130.91
Wolf Creek wolfcr 60.75 60.25 -134.47 -135.74
WrigLN wrigLN 63.71 62.85 -122.73 -123.96
Yukon Flats yflats yflatE, yflatW Bh007v004, Bh008v004, Bh007v005, Bh008v005 67.21 65.73 -144.74 -147.42
Yukon-Kuskokwin Delta ykdelt ykdelA, ykdelB Bh002v003 61.32 61.05 -161.86 -163.89
Canada
Aklavik Highway aklavi aklavi Bh010v005 68.37 68.01 -133.11 -135.55
Behchoko behcho behcho Bh013v011, Bh013v010 62.67 62.08 -116.41 -116.88
Daring Lake daring daring Bh014v009, Bh015v009, Bh014v010 65.00 63.65 -111.09 -113.33
Faber Lake faberl faberl Bh013v009, Bh013v010 64.39 63.86 -116.95 -118.21
Kakisa Lake A kakisA kakisA Bh012v011, Bh012v012 61.30 60.80 -116.42 -118.02
Kakisa Lake B kakisB kakisB Bh012v011 61.40 60.67 -117.38 -117.99
Old Crow Airport oldcrB oldcrB Bh008v005, Bh009v005 67.63 66.85 -139.58 -143.17
Old Crow Flats oldcrA oldcrA Bh009v005 68.22 67.66 -139.05 -140.05
Fort Providence provid provid Bh013v011, Bh012v011 62.20 61.29 -116.21 -117.90
Snare River snarer snarer Bh012v010, Bh013v010 63.96 62.72 -116.56 -118.30
Inuvik to Tuk Highway tukhwy tukhwy Bh010v005, Bh011v005 69.51 68.18 -132.92 -134.41
Yellowknife yellow yellow Bh010v005, Bh011v005 62.58 61.94 -112.34 -113.82

数据文件信息

 

这里有 51 个 netCDF 第 4 版 (*.nc4) 格式的文件(每个站点一个),以及一个根据数据估算土壤容积含水量的 python 脚本。

netCDF 文件命名为 PDO_ReSALT_site_2017_03.nc4(例如,PDO_ReSALT_aklavi_2017_03.nc4),其中站点作为站点代码提供在表 1 中。

为用户提供了 Python 脚本文件 generate_pdo_soil_vwc.py,用于生成平均到感兴趣深度的土壤容积含水量(VWC)。该脚本需要使用 numpy 和 gdalr 库。执行脚本的步骤如下

python generate_pdoo_soil_vwc.py -s site_code -a depth_avg_in_cm

其中,site_code 是表 1 中列出的站点,depth_avg_in_cm 是以厘米为单位的平均深度。

数据文件详情

缺失值用 -9999 表示。

使用的投影为 "Canada_Albers_Equal_Area_Conic"(EPSG:102001);Proj.4 字符串:"+proj=aea +lat_1=50 +lat_2=70 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m no_defs"。

应用与推导

土壤湿度和活动层厚度(ALT)是了解高纬度地区永久冻土和活动层动态如何应对气候变暖的关键变量。InSAR 等遥感技术提供了测量 ALT 的方法(Liu 等,2012 年)。

质量评估

InSAR 相关数据用于掩盖不良数据和不确定性量化。相干性小于 0.35 的差相关像素被屏蔽,因为它们不能提供可靠的相位估计值。利用克拉默-拉奥约束(Tough 等人,1995 年)对 InSAR 相关性的相位不确定性进行估计,并将其转换为厘米级的下沉不确定性。由于在一个解冻季节内缺乏重复的时间机载观测,Cramer-Rao 边界估计值被认为是下限估计值,实际的不确定性可能比我们报告的要大。尽管如此,对相位校准不确定性的经验估算得出的点值约为 0.5-1.0 厘米。不确定值的合理范围在克拉默-拉奥边界和 1.0-2.0 厘米之间。

根据绝对辐射校准,P 波段 PolSAR 后向散射的不确定性假定为 0.5 dB(Chapin 等人,2015 年)。根据前向模型,构建了一个包含所有可能解决方案的数据立方体。所有解决方案的成本函数值均小于 P 波段后向散射不确定性。这些解决方案的平均值被报告为所需地球物理变量的最佳估计值,标准偏差被报告为不确定性。

数据采集、材料和方法
InSAR 处理

InSAR 利用不同时间获取的合成孔径雷达图像之间的相位差,生成雷达视线方向上的地表变形干涉图。对 NASA/JPL UAVSAR (https://uavsar.jpl.nasa.gov) 提供的干涉相位采用 SNAPHU 算法(Chen 和 Zebker,2002 年)进行解包。从视线(LOS)观测几何图形解包的干涉相位被转换为垂直运动(Chen 等人,2016 年)。大气和对流层含水量的变化会给干涉相位值带来偏差(Jolivet 等人,2011 年)。对于有明显大气噪声的扫面,采用高通滤波器去除与大气相关的长度尺度上的信号(Lohman 和 Simons,2005 年)。

日测温数据(Thornton 等人,2016 年)用于生成每个站点的累计解冻度日(ADDT)时间序列。利用 ADDT 时间序列对测量的下沉量进行调整,以估算解冻季节期间活动层的季节性总下沉量。对 Schaefer 等人(2015 年)的技术进行了修改,以根据整个解冻季节与 L 波段合成孔径雷达图像采集日期之间的 ADDT 差值来调整测得的下沉量。有些站点未包含 ADDT 修正,在第三版处理过程中添加了 ADDT 修正。

参考站点的已知沉降值用于将相对阶段转换为绝对季节沉降。如果有可靠的 ALT 原位实地测量数据,则使用土壤膨胀模型将其转换为季节下沉值(Schaefer 等人,2015 年;Michaelides 等人,2019 年)。对于没有可靠原位数据的地点,根据全场景相对变形值直方图,将变形的第 5 百分位数选为零。所有季节性形变、相关性和形变不确定性的估计值都投影到 30 米 ABoVE 参考网格上。

数据版本

Version Release Date Revision Notes
3 2022-08-31 Current dataset. For this version, new flight lines were added, reference points were updated, a new seasonal calibration was used, atmospheric and numerical noise was removed, additional uncertainty was quantified, and discontinuities were corrected. See Section 5 for details. 
2 2021-04-01 In Version 2, retrieval of P- and L-band data was performed in tandem. For this version, the seasonal subsidence model and the radar backscattering model were integrated using a joint inversion process. Additional sites across the ABoVE domain were added. This version is now superseded and available only upon request. ABoVE: Active Layer Thickness Derived from Airborne L- and P-band SAR, Alaska, 2017, https://doi.org/10.3334/ORNLDAAC/1796
1 2019-05-08 Version 1 of this dataset. This version is now superseded and available only upon request. ABoVE: Active Layer Thickness Derived from Airborne L- and P-band SAR, Alaska, 2017, https://doi.org/10.3334/ORNLDAAC/1676

代码

!pip install leafmap
!pip install pandas
!pip install folium
!pip install matplotlib
!pip install mapclassify
 
import pandas as pd
import leafmap
 
url = "https://github.com/opengeos/NASA-Earth-Data/raw/main/nasa_earth_data.tsv"
df = pd.read_csv(url, sep="\t")
df
 
leafmap.nasa_data_login()
 
 
results, gdf = leafmap.nasa_data_search(
    short_name="ABoVE_ReSALT_InSAR_PolSAR_V3_2004",
    cloud_hosted=True,
    bounding_box=(-166.73, 57.83, -110.42, 71.52),
    temporal=("2017-07-20", "2017-08-08"),
    count=-1,  # use -1 to return all datasets
    return_gdf=True,
)
 
 
gdf.explore()
 
#leafmap.nasa_data_download(results[:5], out_dir="data")

引用

Chen, R.H., R.J. Michaelides, J. Chen, A.C. Chen, L.K. Clayton, K. Bakian-Dogaheh, L. Huang, E. Jafarov, L. Liu, M. Moghaddam, A.D. Parsekian, T.D. Sullivan, A. Tabatabaeenejad, E. Wig, H.A. Zebker, and Y. Zhao. 2022. ABoVE: Active Layer Thickness from Airborne L- and P- band SAR, Alaska, 2017, Ver. 3. ORNL DAAC, Oak Ridge, Tennessee, USA. ABoVE: Active Layer Thickness from Airborne L- and P- band SAR, Alaska, 2017, Ver. 3, https://doi.org/10.3334/ORNLDAAC/2004

网址推荐

0代码在线构建地图应用

https://sso.mapmost.com/#/login?source_inviter=nClSZANO

机器学习

https://www.cbedai.net/xg 


【版权声明】本文为华为云社区用户原创内容,转载时必须标注文章的来源(华为云社区)、文章链接、文章作者等基本信息, 否则作者和本社区有权追究责任。如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。