Change Analysis of Three Major Antarctic Ice Shelves Based on Multi-source Remote Sensing Data

In recent years, melting and calving happen in the Antarctic ice shelves. In this paper, seven periods of coastlines were applied to provide an analysis of the ice front changes of Ross, Filchner-Ronne and Amery ice shelves with the inclusion of 1997 and 2000 Radarsat products, 2003/2004 and 2008/2009 MODIS products and 2006, 2012, 2015 coastline which were extracted from MODIS images. Change area, SCE (Shoreline Change Envelope) and NSM (Net Shoreline Movement) were applied to analyze the variation of the ice shelf front. The results shows that, the ice front of Amery ice shelf has advanced since 1997 and the total outward extension distance of the ice front was about 20 to 25 km while the advance area reached 3.03×10³ km². Ross ice shelf and Filchner-Ronne ice shelf continued to advance after ice calving events under the driver action of glacier. However, the advance area was less than the retreat area and the net change area is respectively -9.39×10³ km² and -5.86×10³ km². The retreat distance of the collapse area were up to 53 km and 39 km in the two biggest ice shelves.     

小冰期时中国南方地区降水模式的差异研究

小冰期是过去一千年中全球气候变化的重要事件之一。关于小冰期时中国季风区和西风影响区气候变化的对比研究众多,但是缺乏中国南方地区降水模式时空差异的研究,难以了解中国南方地区降水变化规律。为了系统地了解小冰期时中国南方地区降水的复杂性,本文将中国南方地区划分为东南—华南沿海地区、中部地区以及西南地区三个区域,总共选取了19条高分辨率的古气候记录进行对比研究,主要有以下几点认识:(1)相对于中世纪暖期而言,小冰期期间中国南方东南—华南沿海地区的气候偏湿,这可能与雨带在中国南方的滞留时间延长和沿海地区受台风的影响增强有关。(2)中国中部地区秦岭南麓和神农架高山林区在小冰期时期主要呈"冷湿"的模式,差异在于秦岭南麓区域主要在小冰期中后期偏湿,这与中部其他区域偏"冷干"的模式不同。这种区域差异可能是由于地形地势和大气环流的复杂性导致。(3)中国西南地区受印度夏季风和东亚夏季风的共同影响,且该区域地形复杂,其气候变化在小冰期时期存在更加明显的空间差异,没有呈现出比较一致的降水模式。与小冰期期间的降水变化不同的是,近30年东南—华南沿海地区除了台湾和雷州半岛,其他区域降水明显减少,可能受气温和人类活动等因素的影响。通过结合高分辨率的古气候记录,我们系统分析了中国南方小冰期的干湿模式在时空上的差异及其可能的影响因子,这对于认识小冰期时中国南方不同区域降水的复杂性及未来旱涝灾害的防控具有一定意义。     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground-penetrating radar

Across the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large-scale storms, while providing important habitats and recreational opportunities. Understanding their long-term evolution is vital in managing their condition in a changing climate. Recently, ground-penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre-scale resolution images of near-surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo-environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground-truth control. Nonetheless, GPR data interpretation can be non-intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi-automated GPR feature-extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo-landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty-organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR-interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.     

Tectonic Movement and Global Climate Change

Glaciation between northern hemisphere and southern hemisphere were synchronous, the ice age occurred not in high but in low value of the eccentricity of the earth's orbit. Such facts went against the precession principle of the astronomical theory of ice age. The inhomogeneous distribution of climate consisted with the inhomogeneous distribution of ocean and continent. The north/south antisymmetry may be attributed to southward deviation of the thermal center and northward deviation of the mass center within the mantle demonstrated by seismic tomography. The core - mantle angular momentum makes rotational energy into thermal energy and mantle plumes erupt in the ocean bottom. The earth's deformation by tidal force makes the eruption of mantle plumes strong. They are the reason that glaciation between the Northern Hemisphere and Southern Hemisphere are synchronous and the ice age occurred in low value of the eccentricity of the earth' s orbit. The tectonic movement is playing a most important part in global climate change.     

沟道偏转地形对滑坡碎屑流运动的影响研究

在沟谷地形中,滑坡碎屑流的运动常受到地形的影响,导致其运动方向发生改变,进而影响到滑坡运动速度和堆积特征。本文利用三维离散元素法,对四川都江堰三溪村高速远程滑坡进行模拟,研究滑坡体不同部位的块体失稳后,在沟道偏转地形主导下的滑坡碎屑流前缘的运动速度、各部位滑体的速度变化过程和堆积特征,并提出沟道偏转地形耗能模型分析了地形偏转造成的动能消耗。研究结果显示:滑坡前缘在地形偏转位置运动方向发生变化,导致运动速度突降;由于滑坡不同部位的滑块相对于地形偏转点具有不同的撞击角度,导致其撞击后产生不同的偏转角度,滑块的偏转角度越大,速度变化越大;由沟道偏转地形导致的滑坡运动速度减小反映了偏转地形对滑坡的动能产生的耗散,动能耗散率与cos2θ(θ为偏转地形在水平面上的偏转角度)成反比;不同部位滑块的堆积长度随偏转角度的增大而减小。本研究分析了沟谷地形偏转对滑坡碎屑流运动速度作用机制及不同部位岩土体堆积范围的影响,可为该类地形条件下滑坡的运动机制研究和防灾减灾工作提供参考。     

利用大地水准面模型计算垂线偏差的方法及精度分析

根据大地水准面与垂线偏差的关系,设计合理的计算方案,给出利用大地水准面模型计算垂线偏差的简化公式,并通过模拟计算探讨大地水准面相对精度、取点间距和已知点选取及个数对计算结果的影响。利用GEOID12B模型分别计算GSVS2011、GSVS2014项目中各测站点和美国西部区域（40°~45°N,100°~105°W,分辨率为1′）的垂线偏差,并与GSVS项目垂线偏差实测值和DEFLEC12B模型值进行比较。结果表明,垂线偏差南北分量和东西分量的计算精度均优于±0.5″,说明利用相对精度为cm甚至亚cm级的大地水准面模型可获取较高精度的垂线偏差。     

2017-08-09新疆精河Ms6.6地震附近区域相对应变的变化分析

2017年精河M_S6.6地震震前,尼勒克、巴伦台和小泉沟台钻孔应变仪记录到显著的应变变化,经现场落实确认,异常是可靠的。根据这3个台的记录资料计算相对应变的变化,结果表明,地震前最大（小）主应变大小和方向分别出现明显加速和急剧偏转变化。此外还发现,相对应变场异常期间主方向与震源机制解P轴方位较为一致。     

Body motion in liquid under ice plate with snow cover

The motion of a submarine in liquid under an ice plate covered with flooded snow is considered. The ice is modelled as an elastic plate and the snow cover is modelled as a viscous layer on the top of the plate. The submarine is modelled as a slender solid of revolution with scale 1:300. The experimental and theoretical study of the influence of the viscous snow layer on deflections of the floating ice plate is conducted. The viscous layer reduces the amplitudes of flexural-gravity waves. The greatest influence of the viscous layer on the plate deflections is achieved for velocities of the submarine, where the waves of maximum amplitude are generated. Theoretical results are in good qualitative and quantitative agreement with the model experiments.     

Modelling of the passive turning of a turret moored vessel in conditions of compact ice

The concept of the offshore oil and gas field development using floating ship-shaped platforms, frequently named floating production units (FPU further), with turret mooring is widely used in the world, including regions with harsh environment. The direct transfer of this concept to the Arctic seas is not possible, mostly due to the difficulties of maintaining a weathervaning, or passive turning regime in heavy ice conditions. The main danger relates to the fact that the expected FPU rotation under the action of high ice loads can be accompanied by a translational displacement of the FPU away from the mooring point. As a result, the mooring system may reach an overloaded state until the FPU turns to a favourable position relative to the ice drift direction. In the paper, we are focused on the investigation of a mathematical model of the passive FPU turning on a spot under the assumption that the ice cover is described by a rigid-plastic continuum. The study is performed both analytically and by numerical simulations. A number of specific FPU motion patterns are analytically derived from the model in quasi-static approximation in the form of successive limit states of the system FPU – ice continuum. Some results of the corresponding numerical simulation are presented that confirm the existence of similar solutions in the full dynamic setting of the problem. A partial parametric analysis of the problem is also performed.