freeboard2thickness文档
freeboard2thickness是南极的一部分映射工具Matlab(格林et al ., 2017)。点击这里在AMT功能的完整列表。
freeboard2thickness估计冰层厚度从表面海拔高度,假设hyrostatic平衡。
内容
语法
T = freeboard2thickness (F) T = freeboard2thickness (…, rhoi, iceDensity) T = freeboard2thickness (…, rhow, waterDensity) T = freeboard2thickness(…,“罗斯”,snowDensity) T = freeboard2thickness (…, Ts, snowThickness)
描述
T = freeboard2thickness (F)估计冰层厚度T在米高度大地水准面F的米。
T = freeboard2thickness (…, rhoi, iceDensity)指定冰密度公斤/米^ 3。默认冰密度是917公斤/米^ 3。
T = freeboard2thickness (…, rhow, waterDensity)指定水密度公斤/米^ 3。默认的水密度是1027公斤/米^ 3。
T = freeboard2thickness(…,“罗斯”,snowDensity)指定雪密度公斤/米^ 3。默认雪密度350 kg / m ^ 3,但是默认积雪厚度为0 m,所以雪密度值只会影响计算如果指定积雪厚度。
T = freeboard2thickness (…, Ts, snowThickness)指定在米雪厚度。默认的积雪厚度为0 m。
示例1:
冰山表面的海拔是80厘米,它是纯冰。冰山有多厚?
T = freeboard2thickness (0.8)
T = 7.4691
你可以回到姐姐的其他方式功能:
F = thickness2freeboard (7.47)
F = 0.8001
示例2:
冰山表面的海拔高度是80厘米,包括40厘米的积雪的密度为300公斤/米^ 3。冰山有多厚?
T = freeboard2thickness (0.8,“t”,0.4,“罗斯”,300)
T = 4.8255
示例3:当猪浮动
这个例子看起来在一个流线松岛冰川。得到一个实际的基于表面流线速度的坐标可以使用流线功能,但是现在让我们输入一些手动点和他们的阴谋南极洲的MODIS马赛克的背景图片。我沿着流线致密化点100米的间隔pspath功能:
(纬度、经度)= pspath ((-75.3585 - -75.2041 -75.0089 - -74.8032),…(-98.0220 - -99.4048 -100.7342 - -101.5879),…100);%(100米间距以及跟踪)plotps(纬度,经度,“b -”,“线宽”,2)modismoaps轴从
如果你有Bedmap2插件AMT沿着线高程资料非常简单:
图bedmap2_profile(纬度、经度);
这似乎有点reduntant,但情节Bedmap2表面为红色。我们还需要绘制它的函数距离沿着纬度,经度路径相同的方式bedmap2_profile画上面的距离。所以使用pathdistps得到累积距离跟踪:
%得到Bedmap2表面高度:sfz = bedmap2_interp(纬度,经度,“表面”);%计算沿流线距离公里:d = pathdistps(纬度,经度,“公里”);%画出插值Bedmap2表面作为粗红线:情节(d、sfz“r”,“线宽”,2)
现在我们可以使用水压估计冰的基础假设。首先假设冰架是由纯冰密度为917公斤/米^ 3。
从干舷%计算冰层厚度:厚度= freeboard2thickness (sfz);%底部高程表面-厚度:基础= sfz-thickness;%画出静水冰基厚红线:情节(d,基地,“r”,“线宽”,2)%的指数最后有限的地表高程测量:印第安纳州=找到(isfinite(基地),1,“最后一次”);%一个文本标签:文本(d(印第安纳州),基地(印第安纳州),“0 m FAC”,“fontweight”,“大胆”,“颜色”,“r”)
这并不是一个了不起的比赛!我记得,冰架基地海拔,Bedmap2使用水压的假设后占一定的粒雪冰列中的空气含量。我们可以做同样的通过设置“雪”厚度20 m和积雪的密度设置为0公斤/立方米。我认为Ligtenberg数据显示,积雪空气含量是猪20米。让我们再试试这个20 m FAC:
thickness2 = freeboard2thickness (sfz,“t”,20岁,“罗斯”,0);base2 = sfz-thickness2;情节(d、base2“b”,“线宽”2)文本(d(印第安纳州),base2(印第安纳州),“20 m FAC”,“fontweight”,“大胆”,“颜色”,“b”)
冰川变薄和接地冰线撤退
床上= bedmap2_interp(纬度,经度,“床上”);基础= sfz - freeboard2thickness (sfz,“t”15岁的“罗斯”,0);%找到冰基地电梯从床上的索引:印第安纳州=找到(基床> = 1“第一”);%画一个红色圆圈静压接地冰线:情节(d(印第安纳州),基地(印第安纳州),“罗”,“markersize”,10)
这是接地冰线撤退的粗略近似值的函数冰川变薄。我们将循环,发现上面的接地线位置一样,但对于一个ever-thinning猪:
% 10厘米表面高程变化的步骤从0到40 m:dh = 0:0.1:40;% Preallocate接地冰线位置向量:gl =南(大小(dh));%新建一个概要图图bedmap2_profile(纬度、经度);% colormap:颜色=喷气(长度(dh));%遍历每个稀疏值:为k = 1:长度(dh);%表面高程与k的稀释:sfzi = sfz - dh (k);%液压基础高度稀释:basei = sfzi - freeboard2thickness (sfzi“t”15岁的“罗斯”,0);%绘制表面和冰基地海拔高度:情节(d、sfzi“- - -”,“颜色”、颜色(k,:))情节(d basei“- - -”,“颜色”颜色(k,:))%找到冰基地电梯从床上的索引:印第安纳州=找到(床basei > = 1“第一”);%日志接地线位置:gl (k) = d(印第安纳州);结束
现在我们可以把接地线位置作为表面高程变化的函数。
图绘制(gl, dh)轴紧盒子从包含({接地线位置的;沿着流线的距离(公里)})ylabel“冰架表面高程变化(m)”
我应该强调的是,上面的分析是强大的原油。我不会相信这些确切的数字,因为我们完全忽视了压力在冰,接地区域动力学,床上的不确定性,冰层厚度的不确定性,等等,但冰川学的基本原则是:从右上角开始的情节,随着冰层变薄,接地线撤退。有时候,一点点的变薄会导致很多接地冰线撤退。
援引AMT
如果这个函数或其他任何部分的南极地图工具对你很有用,请引用描述了AMT的纸。
格林,c。,Gwyther, D. E., & Blankenship, D. D. Antarctic Mapping Tools for Matlab.电脑与地球科学。104 (2017)pp.151 - 157。doi: 10.1016 / j.cageo.2016.08.003。
作者信息
这个函数是由乍得a。格林德克萨斯大学的地球物理研究所(UTIG), 2017年4月,包含在南极映射工具软件包。