南极纳尔逊冰帽的雪层及其成冰作用研究

纳尔逊冰帽雪—粒雪的演化依赖于融水渗浸冻结作用下的暖变质过程。积雪的密实化过程的快慢,取决于温度条件和融水的参与程度,以及自身的物理状况。粒雪的密实化过程表现为均匀且变幅小。纳尔逊冰帽成冰深度在23～25m,成冰历时17～19年。成冰带分为暖渗浸—重结晶带,渗浸—冻结带,消融带。     

南极洲纳尔逊冰帽浅层粒雪／冰的晶体组构特征

对纳尔逊冰帽３个地点浅层粒雪／冰芯的晶体组构测量表明，冰晶ｃ轴在冰帽表面具有优势方位，２０ｍ深度内普遍出现双扇形组构，在２０－３０ｍ深度具有向单极大型变化的趋势。这与所报道的中低纬度山地温冰川上的晶体组构特征有明显不同。估计大量融水渗浸和再冻结作用对冰晶组构有重要影响，但详细机制尚不清楚     

珠穆朗玛峰北坡东绒布冰川成冰作用的新认识

冰川成冰作用的研究对于选择冰芯钻取点具有重要的科学意义。前人对珠穆朗玛峰北坡冰川成冰作用的研究,由于缺少高海拔区域的实测资料而具有一定的局限性。文章通过1998年东绒布冰川垭口处(6 500 m a. s. l.)11 m冰芯和海拔6 450 m处20 m冰芯剖面的成冰作用过程研究,认识到由于水、热条件的逐年波动,冰川成冰作用也处于变化之中。珠穆朗玛峰北坡东绒布冰川高海拔区域,在一定的水、热条件下(如气温较低和降水量较大等),再冻结-重结晶作用依然占主导地位,该成冰作用至少在垭口部位是有分布的。而一般在气温较高或降水量较少等条件下,冰川的成冰作用则以冷渗浸-重结晶作用为主。     

天山东部冰芯pH值和电导率的大气环境空间差异

冰芯电导率是大气环境变化的替代性指标.对我国天山东部三个研究点奎屯哈希勒根48号冰川、乌鲁木齐河源1号冰川以及哈密庙儿沟平顶冰川粒雪芯中pH值和电导率的近期环境记录进行了分析研究.结果显示,三个研究点粒雪芯反映的近期pH值和电导率变化趋势不同:奎屯电导率随冰芯深度增加表现出升高趋势,哈密庙尔沟平顶冰川和乌鲁木齐河源1号冰川随冰芯深度增加表现出一定下降趋势;冰芯电导率平均值在哈密庙儿沟冰川最大,而其他两个点相对较小,这与矿物粉尘浓度和离子浓度的分布一致;雪冰电导率与粉尘及离子相关性分析表明,电导率主要受中亚粉尘活动影响,同时与Ca2+、Mg2+、Na+等相关性很高.通过比较显示,天山东部与我国西部其他研究点以及极地雪冰电导率的差异可以很好地反映大气环境的空间差异.     

南极洲乔治王岛柯林斯冰帽冰芯层位和密度变化的一般特征

柯林斯冰帽两支冰芯层位记录证实了该冰帽主冰穹顶部（海拔约７００ｍ）属暖渗浸带，小冰穹顶部（海拔约２５０ｍ）属渗浸带。雪、冰的层位分布和密度变化包含了一定的测年信息。主冰穹顶冰芯密度－深度曲线在表层呈现韵律性变化，与层位记录中的融化冻结现象相对应，据此粗略划分年层，断定当地年积累雪层厚度为３－３．５ｍ，折合水当量１６５０－１９２５ｋｇ／ｍ２ａ，年平均积累速率约为２．０ｍ／ａ（冰当量）。主冰穹顶成冰深度为３８－３９ｍ，此深度以上密度自上而下缓慢增加，但以下由于含水层的出现，密度迅速升高，在５－６ｍ区间达到９００ｋｇ／ｍ３。小冰穹冰芯除表层外，密度基本在８００－９００ｋｇ／ｍ３之间，冰芯中雪冰互层，存在污化面，４０ｍ以下发现很厚的火山灰沉积物。小冰穹平均年积累率约为０．７ｍ／ａ（冰当量），成冰深度７－８ｍ，成冰年限为１０年左右     

南极威德尔海周边雪冰和大气气溶胶中的MSA,nssSO_4~(2-)浓度及其比率研究

分析了威德尔海周边雪冰和大气中ＭＳＡ和ｎｓｓＳＯ２－４浓度资料。结果表明,威德尔海是造成这一地区大气中高含量生物硫化物的“源”,其周边岛屿及冰架表层雪冰记载和反映了大气中生物硫化物的空间分布和季节变化特点。源区产出硫化物的强度存在显著差异：靠近南极半岛一侧高于靠近Ｆｉｌｃｈｎｅｒ－Ｒｏｎｎｅ冰架一侧。尽管诸多因素对输运和沉降过程发生影响,表层积雪ＭＳＡ和ｎｓＳＯ２－４浓度仍呈现出很有规律性的空间分布：随沉降地点离海距离和海拔高度的增加,浓度值递减;但在一定海拔之下,“高程效应”不显著。冰芯中ＭＳＡ浓度垂向分布显示出季节分配“滞后”现象,在表层表现为受天气条件制约的输运过程造成的相位差,在深层则归因于有待验证的“迁移”机制作用的结果。大气和雪冰样品对比研究结果表明,在类似威德尔海这样的高纬地区,大气中的ＳＯ２－４和ＭＳＡ“信号”在雪冰样品中会遭到不同程度的减弱。但就空间变化趋势和季节分布而言,表层雪冰仍是大气中组分（在本文是ＳＯ２－４和ＭＳＡ）的良好载体。十分接近的大气和雪冰ＭＳＡ对ｎｓｓＳＯ２－４（或ＳＯ２－４）比率证实,这两种大气组分在由大气洗净和沉降到冰雪过程中只发生微弱的分离作用。这也就是大气和雪冰中     

南极威德尔海海冰晶体类型分析

为了解南极威德尔海不同区域海冰形成过程,分析了2006年9—10月威德尔海海冰调查中的27支冰芯晶体并描述了冰芯反映的冰层形成过程;给出2支完整冰芯晶体照片和多边形晶体雪冰照片,以及对应照片的冰层描述和冰盐度、密度、粒径。根据冰芯晶体揭示的海冰形成过程,发现尽管现场考察均选择较大表面平整的浮冰作业,但冰芯晶体揭示出实际冰层可以是受动力和热力联合作用形成的重叠冰、固结冰脊、二年冰;而从表面到底面属于纯热力学生长形成的冰层所占比例较小。积雪经过密实和水分参与,再冻结形成多边形晶体雪冰也是南极海冰主要组成之一。粒状晶体、过渡区混合晶体和柱状晶体海冰分别占27支冰芯总长度的28.7%、14.4%和55.2%。海冰的区域分布表现出海冰边缘区以纯热力学生长的一年冰为主;在冰架前缘,受动力作用形成的重叠冰和固结冰脊的比例增加,并且存在二年冰和多边形晶体雪冰;冰间湖内生长有大面积纯热力学生长的海冰并向外输送。     

南半球环状模与南极冰芯中气候信息关联研究综述

冰芯作为气候环境信息的良好载体,在研究过去大气环流异常中具有其他气候代用指标难以比拟的优势。本文重点综述了南极冰芯气候记录与南半球环状模(SAM-Southern Annular Mode)的关联特征:SAM与冰芯记录为同期变化,冰芯中NO-3、海盐组分(以Na+为代表)和水同位素比值δ18O同SAM指数具有较好的关联性;冰芯中SO2-4和MSA同SAM存在一定的关联,但相关性不显著;西南极冰芯记录的积累率同SAM具有较强的关联。需要说明的是,冰芯中各参数与SAM的相关性存在较大的空间差异。     

极地气溶胶及雪冰中非海盐硫酸根与甲基磺酸时空分布及环境指示

非海盐硫酸根(nssSO_4~(2-))和甲基磺酸(MSA)是极地气溶胶和雪冰介质中主要的含硫化合物。空间尺度上,南半球nssSO_4~(2-)和MSA具有明显的纬度效应特征;时间尺度上, nssSO_4~(2-)和MSA具有季节变化特征(夏高冬低),南极冰芯中MSA和nssSO_4~(2-)浓度均在冰期增大,北极冰芯记录出现MSA浓度在冰期减小、但是nssSO_4~(2-)浓度在冰期仍然增大的现象。利用极地nssSO_4~(2-)和MSA可以指示海洋初级生产力、大气温度、海冰范围变化以及人类活动影响等气候环境信息,并对利用nssSO_4~(2-)和MSA开展全球变化研究进行了展望。     

中国天山冰川区降水、积雪pH和电导率季节变化特征分析

根据乌鲁木齐河源区2003-2004年度的大气降水样品、1号冰川一个完整年周期的连续雪坑样品及哈密庙儿沟平顶冰川和奎屯51号冰川雪坑样品的pH和电导率资料,系统探讨大气降水、表层雪、雪坑样品中pH和电导率的季节变化特征和沉积过程中pH和电导率的变化特点,以及不同冰川区雪坑中pH值和电导率的空间差异。结果表明:大气降水及表层雪中pH和电导率有着相似的季节变化趋势,主要受冰川所在区域大气中可溶性离子浓度和沉积后过程的影响;对1号冰川大气降水、表层雪和雪坑样品中pH和电导率的对比分析发现,从降水到表层雪,最后到雪坑样品的演化过程中,pH和电导率均呈减小趋势,而且pH和电导率的最高值均出现在春季;自西向东随着地理环境的变化,雪坑中样品的pH值呈升高趋势,哈密庙尔沟平顶冰川最大,乌鲁木齐河源1号冰川次之,奎屯哈希勒根51号冰川最小,但电导率却呈现出了不同的变化规律。     

利用地球重力场模型确定南极大地水准面

利用我国最新地球重力场模型 WDM94,给出了南极 (纬度范围为 - 60°～ - 90°)大地水准面高和平均空间重力异常。为了全面总结分析南极大地水准面特征 ,收集了国外最新地球重力场模型 OSU91 ( 360阶次 )和 JGMOSU( 360阶次 ) ,计算了相应的大地水准面高和平均重力异常。其结果分别与 WDM94的结果作了比较 ,WDM94与 OSU91和 JGMOSU的大地水准面高标准差分别为± 1 .90 m和± 2 .0 9m,平均空间重力异常标准差分别为± 8.97mgal和± 9.32 mgal     

珠穆朗玛峰重力值的推估──兼论推估方法的基本原理

在山区 ,尤其是在有全球第三极之称的喜马拉雅山区 ,当相邻点间距不大时 ,如何利用这些点上的重力与地形 (高程 )数据推估待求点的重力值 ,这对难以攀登和不能用仪器观测的山峰很有意义。研究指出 ,在地形负荷的波长很短时 ,具有一定强度的地壳足以能够支撑这种负荷 ,因此 ,不能用Airy Heiskanen和Pratt Hayford局部补偿模型作重力推估 ;由于空间异常主要受地形起伏制约 ,因此借助于邻近重力点的地形 (高程 )作推估会得到满意的结果。基于这一思路 ,我们采用了 4种有关公式 ,有效地推估了第三极之巅珠穆朗玛峰顶上的重力值 ,该值为(976 970± 7)× 1 0 - 5m·s- 2 。这一结果为精确推求珠峰大地水准面和正高提供了必要的数据 ,若用均衡的方法来推估 ,则可能相差近 1 0 0× 1 0 - 5m·s- 2 。     

南极重力场特征及界面的计算与解释

本文根据ＯＳＵ９１Ａ地球重力场模型和地形资料，计算了南极洲的空间重力异常和布格重力异常，分析了空间重力异常变化剧烈的原因及其与高程的相关性，同时根据布格重力异常用两层界面的反演方法计算了冰盖厚度和地壳厚度，冰盖厚度较大的地区位于极区的东南部，而极区周围和西南部地区厚度较小。地壳较厚的地区位于极区的东南部，最大达５６ｋｍ，西南部地区地壳较薄，最小值为８ｋｍ。     

中国南极长城站测绘基准系统的建立

经1985年12月～1987年3月,三年度夏测绘科学考察,在中国南极长城站地区建立了一套完整的、精确的测绘基准系统。它包括大地坐标系统、高程系统和重力参考系统。本文叙述了建立这些基准系统的方法和数学模型,并分析了测量精度。     

Validation and fusion of different databases in preparation of high-resolution geoid determination

An up to date determination of a high-resolution geoid requires the use of best available databases concerning digital terrain model (DTM), bathymetry, global geopotential model and gravity field. The occasion to revisit methods to validate and merge different data sets has been created by a new project for the determination of a new European Geoid.
Since the computation of the latest European geoid and quasi-geoid model (EGG97), significant new or improved data sets have become available, such as new global geopotential models from CHAMP and GRACE missions, new national and global DTMs and new or upgraded gravity data sets.
In the context of the new European Gravity and Geoid Project (EGGP), within the IAG Commission 2, some data validation tests have been performed in the Italian zone.
In the area 19°× 17° wide, covering Italy, three kinds of tests have been performed: comparison among different DTMs in order to choose the best one to be used; comparisons in terms of geoid computation in some coastal areas, to evaluate bathymetry effects, and the validation of the EIGEN-CG01C and EIGEN-CG03C new global models up to degree and order 360.
These preliminary tests lead to the choice of SRTM DTM (integrated in no-data holes), with an added bathymetry derived by the Italian 1:25 000 official cartography near the coasts and the NOAA bathymetry in high seas. The validation of the new global models and the comparison with EGM96 model show that, in terms of geoid computation, the EGM96 yields better results. Moreover, the validation of new available land gravity data and the cross-validation of two sets of gravity data on sea have been completed.     

The reduction of aliasing in gravity anomalies and geoid heights using digital terrain data

Observations of gravity can be aliased by virtue of the logistics involved in collecting these data in the field. For instance, gravity measurements are often made in more accessible lowland areas where there are roads and tracks, thus omitting areas of higher relief in between. The gravimetric determination of the geoid requires mean terrain-corrected free-air anomalies; however, anomalies based only on the observations in lowland regions are not necessarily representative of the true mean value over the topography. A five-stage approach is taken that uses a digital elevation model, which provides a more accurate representation of the topography than the gravity observation elevations, to reduce the unrepresentative sampling in the gravity observations. When using this approach with the Australian digital elevation model, the terrain-corrected free-air anomalies generated from the Australian gravity data base change by between 77.075 and −84.335 mgal (−0.193 mgal mean and 2.687 mgal standard deviation). Subsequent gravimetric geoid computations are used to illustrate the effect of aliasing in the Australian gravity data upon the geoid. The difference between 'aliased' and 'non-aliased' gravimetric geoid solutions varies by between 0.732 and −1.816 m (−0.058 m mean and 0.122 m standard deviation). Based on these conceptual arguments and numerical results, it is recommended that supplementary digital elevation information be included during the estimation of mean gravity anomalies prior to the computation of a gravimetric geoid model.     

Detailed 1×1° gravimetric Indian Ocean geoid and comparison with GEOS-3 radar altimeter geoid profiles

Summary. A new set of 1×1° mean free-air anomalies in the Indian Ocean is determined on the basis of previously published free-air anomaly maps (Talwani & Kahle) and the most recent Lamont surface ship gravity measurements. The data are then used to compute a (total) 1×1° gravimetric Indian Ocean geoid. The computation is carried out by combining the Goddard Space Flight Center (GSFC) GEM-6 geoid and a difference geoid that corresponds to the differences between the set of 1×1° surface gravity values and the GEM-6 gravity anomalies. The difference geoid is highest over the Madagascar Ridge (+ 20 m) and lowest over the Timor Trough (-30 m). The total geoid is compared with GEOS-3 radar altimeter derived geoid profiles and geophysical implications are discussed.     

The relationship between gravity and bathymetry in the Pacific Ocean

Summary. Surface-ship and satellite derived data have been compiled in new free-air gravity anomaly, bathymetry and geoid anomaly maps of the Pacific Ocean basin and its margin. The maps are based on smoothed values of the gravity anomaly, bathymetry and geoid interpolated on to a 90 × 90 km grid. Each smoothed value was obtained by Gaussian filtering measurements along individual ship and subsatellite tracks. The resulting maps resolve features in the gravity, bathymetry and geoid with wavelengths that range from a few hundred to a few thousand kilometres. The smoothed values of bathymetry and geoid anomaly have been corrected for age. The resulting maps show the Pacific ocean basin is associated with a number of ENE–WSW-trending geoid anomaly highs with amplitudes of about ± 5 m and wavelengths of about 3000 km. The most prominent of these highs correlate with the Magellan seamounts–Marshall Gilbert Islands–Magellan rise and the Hess rise–Hawaiian ridge regions. The correlation between geoid anomaly and bathymetry cannot be explained by models of static compensation, but is consistent with a model in which the geoid anomaly and bathymetry are supported by some form of dynamic compensation. We suggest that the dynamic compensation, which characterizes oceanic lithosphere older than 80 Myr, is the result of mantle convection on scales that are smaller than the lithospheric plates themselves.     

Three-dimensional gravity inversion for Moho depth at rifted continental margins incorporating a lithosphere thermal gravity anomaly correction

This paper describes a method for determining Moho depth, lithosphere thinning factor (γ= 1 − 1/β) and the location of the ocean–continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere. The lateral density changes caused by the elevated geotherm in thinned continental margin and adjacent ocean basin lithosphere produce a significant lithosphere thermal gravity anomaly which may be in excess of −100 mGal, and for which a correction must be made in order to determine Moho depth accurately from gravity inversion. We describe a method of iteratively calculating the lithosphere thermal gravity anomaly using a lithosphere thermal model to give the present-day temperature field from which we calculate the lithosphere thermal density and gravity anomalies. For continental margin lithosphere, the lithosphere thermal perturbation is calculated from the lithosphere thinning factor (γ= 1 − 1/β) obtained from crustal thinning determined by gravity inversion and breakup age for thermal re-equilibration time. For oceanic lithosphere, the lithosphere thermal model used to predict the lithosphere thermal gravity anomaly may be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere. A correction is made for crustal melt addition due to decompression melting during continental breakup and seafloor spreading. We investigate the sensitivity of the lithosphere thermal gravity anomaly and the predicted Moho depth from gravity inversion at continental rifted margins to the methods used to calculate and condition the lithosphere thermal model using both synthetic models and examples from the North Atlantic.     

The absolute gravity measurement by FG5 gravimeter at Great Wall Station,Antarctica

Gravity measurement is of great importance to the height datum in Antarctica.The absolute gravity measurement was carried out at Great Wall Station,Antarctica,using FG5 absolute gravity instrument.The gravity data was processed with corrections of earth tide,ocean tide,polar motion and the atmospher,and the RMS is within ±3×10-8 ms-2.The vertical and horizontal gravity gradients were measured using 2 LaCoaste & Romberg(LCR) gravimeters.The absolute gravity measurement provides the fundamental data for the validation and calibration of the satellite gravity projects such as CHAMP,GRACE and GOCE,and for the high accuracy geoid model.