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.     

非饱和土中水平入渗方程显式求解

在求解非饱和态土中水分入渗问题时,水力函数是体积含水率或者吸力的函数,致使其控制方程呈现出强非线性的特征,进而使得其求解变得十分困难。基于水分在土体介质中流动耗时取极值路径的选择这一假定,引入时间泛函,基于变分法原理将水平入渗问题转化为泛函极值问题。通过求解Euler–Lagrange方程,结合边界条件,得到非线性瞬态水平入渗问题的显式解析解。结合Brooks-Corey型水力函数,显式地求解出该类型非饱和态土的体积含水率发展分布规律。通过计算4种不同类型土体的水平入渗规律,将求解结果与已有结果以及数值结果进行对比,验证了该方法的有效性。结果表明:体积含水率分布与位置距离和湿润峰距离比值呈幂函数关系,指数取决于土-水特征曲线的形状参数;初始条件与边界条件会对体积含水率分布造成不同程度的影响。     

静压桩承载力时间效应的研究进展

由于静压桩沉桩后桩周土重塑,静压桩承载力表现出随着休止期的延长而增长的特性。本文从静压桩沉桩后桩周土体内孔隙水消散固结的角度出发,对静压桩承载力时间效应的理论和试验分别进行归纳,结合孔隙水消散路径及固结模型,对桩周土体初始超静孔隙水压力大小及其分布特征进行总结,分析承载力各种测试方法的优缺点,对静压桩承载力的时效性进行深化研究,并探讨了不同地质条件、不同桩的类型对休止期内静压桩承载力的影响,进一步对基于实测数据得出的经验公式进行总结。讨论了基于不同本构关系模型的应力场及位移场解答和沉桩后孔隙水压力消散解答,在此基础上总结了桩基极限承载力理论公式;探讨了黏性土、砂土条件下,考虑超固结比、不排水抗剪强度和塑性指数比对桩基极限承载力系数A的影响,在此基础上归纳了桩基极限承载力经验公式。建议在经验公式基础上设置多重参数,以提升经验公式的精确度,并完善对不同桩、土类型的参数解答;利用BP神经网络,导入静压桩承载力相关参数,以得到针对不同地质条件、桩型、休止期的承载力最优解。     

Integration of zircon and apatite U–Pb geochronology and geochemical mapping of the Wude basalts(Emeishan large igneous province):A tool for a better understanding of the tectonothermal and geodynamic evolution of the Emeishan LIP

Radiogenic isotopic dating and Lu–Hf isotopic composition using laser ablation-inductively coupled plasma-mass spectrometry(LA-ICP-MS)of the Wude basalt in Yunnan province from the Emeishan large igneous province(ELIP)yielded timing of formation and post-eruption tectonothermal event.Holistic lithogeochemistry and elements mapping of basaltic rocks were further reevaluated to provide insights into crustal contamination and formation of the ELIP.A zircon U–Pb age of 251.3±2.0 Ma of the Wude basalt recorded the youngest volcanic eruption event and was consistent with the age span of 251-263 Ma for the emplacement of the ELIP.Such zircons hadε_Hf(t)values ranging from7.3 to+2.2,identical to those of magmatic zircons from the intrusive rocks of the ELIP,suggesting that crust-mantle interaction occurred during magmatic emplacement,or crust-mantle mixing existed in the deep source region prior to deep melting.The apatite U–Pb age at 53.6±3.4 Ma recorded an early Eocene magmatic superimposition of a regional tectonothermal event,corresponding to the Indian–Eurasian plate collision.Negative Nb,Ta,Ti and P anomalies of the Emeishan basalt may reflect crustal contamination.The uneven Nb/La and Th/Ta values distribution throughout the ELIP supported a mantle plume model origin.Therefore,the ELIP was formed as a result of a mantle plume which was later superimposed by a regional tectonothermal event attributed to the Indian–Eurasian plate collision during early Eocene.     

Geochronological and geochemical evidence for a Late Ordovician to Silurian arc-back-arc system in the northern Great Xing'an Range,NE China

The early Paleozoic tectonic evolution of the Xing'an-Mongolian Orogenic Belt is dominated by two oceanic basins on the northwestern and southeastern sides of the Xing'an Block,i.e.,the Xinlin-Xiguitu Ocean and the Nenjiang Ocean.However,the early development of the Nenjiang Ocean remains unclear.Here,we present zircon U-Pb geochronology and whole-rock elemental and Sr-Nd isotopic data on the gabbros in the Xinglong area together with andesitic tuffs and basalts in the Duobaoshan area.LA-ICP-MS zircon U-Pb dating of gabbros and andesitic tuffs yielded crystallization ages of 443-436 Ma and 452-451 Ma,respectively.The Early Silurian Xinglong gabbros show calc-alkaline and E-MORB affinities but they are enriched in LILEs,and depleted in HFSEs,with relatively low U/Th ratios of 0.18-0.36 andεNd(t)values of-1.6 to+0.5.These geochemical features suggest that the gabbros might originate from a mantle wedge modified by pelagic sediment-derived melts,consistent with a back-arc basin setting.By contrast,the andesitic tuffs are characterized by high MgO(>5 wt.%),Cr(138-200 ppm),and Ni(65-110 ppm)contents,and can be termed as high-Mg andesites.Their low Sr/Y ratios of 15.98-17.15 and U/Th values of 0.24-0.25 and moderate(La/Sm)_n values of 3.07-3.26 are similar to those from the Setouchi Volcanic Belt(SW Japan),and are thought to be derived from partial melting of subducted sediments,and subsequent melt-mantle interaction.The Duobaoshan basalts have high Nb(8.44-10.30 ppm)and TiO2 contents(1.17-1.60 wt.%),typical of Nb-enriched basalts.They are slightly younger than regional adakitic rocks and have positiveεNd(t)values of+5.2 to+5.7 and are interpreted to be generated by partial melting of a depleted mantle source metasomatized by earlier adakitic melts.Synthesized with coeval arc-related igneous rocks from the southeastern Xing'an Block,we propose that the Duobaoshan high-Mg andesitic tuffs and Nbenriched basalts are parts of the Late Ordovician and Silurian Sonid Zuoqi-Duobaoshan arc belt,and they were formed by the northwestern subduction of the Nenjiang Ocean.Such a subduction beneath the integrated Xing'an-Erguna Block also gave rise to the East Ujimqin-Xinglong igneous belt in a continental back-arc basin setting.Our new data support an early Paleozoic arc-back-arc model in the northern Great Xing'an Range.     

Solution of the propagation of the waves in open channels by the transfer matrix method

Many problems in mechanics can be solved by the use of the transfer matrix method. The use of this method in hydraulics engineering is not widespread and only limited studies are available. In this study, linearized St. Venant equations were used and the use of transfer matrix in ocean engineering was investigated for long waves in open channels, and numerical application was carried out. The results obtained through the transfer matrix method, which is quite easy to use, program and comprehend, showed similar results obtained from the characteristics method and finite differences method.     

港域波浪数学模型的改进与验证

通过物理模型对改进的港内波浪传播变形数学模型进行验证。该数学模型以推广的时变缓坡方程为控制方程,采用含松弛因子的ADI法求解,并对波浪反射和透射边界模拟方法进行改进。先通过物理模型试验确定斜向浪入射条件下抛石防波堤前的波浪反射系数,作为数学模型中部分反射边界模拟的依据。然后进行了一个典型港口内波浪折射、绕射和反射的模型试验,测量港内波浪分布。对比模型试验和数学模型计算的结果表明,数学模型可较好地模拟港内复杂地形和边界条件下规则波和不规则波的传播变形。     

Rossby waves with linear topography in barotropic fluids

Rossby waves are the most important waves in the atmosphere and ocean, and are parts of a large-scale system in fluid. The theory and observation show that, they satisfy quasi-geostrophic and quasi-static equilibrium approximations. In this paper, solitary Rossby waves induced by linear topography in barotropic fluids with a shear flow are studied. In order to simplify the problem, the topography is taken as a linear function of latitude variable y, then employing a weakly nonlinear method and a perturbation method, a KdV （Korteweg-de Vries） equation describing evolution of the amplitude of solitary Rossby waves induced by linear topography is derived. The results show that the variation of linear topography can induce the solitary Rossby waves in barotropic fluids with a shear flow, and extend the classical geophysical theory of fluid dynamics.     

Experimental investigation of tsunami bore impact pressure on a perforated seawall

Catastrophic failures of many tsunami barriers along the affected coasts during the 2011 Tohoku earthquake tsunami has prompted extensive investigation into improving and revising design codes for tsunami defence structures. To date, researchers and coastal engineers are investigating to understand the failure mechanisms and to find solutions so that the structures merely remain intact in the extreme event such as tsunami. Thus, the present work is motivated to experimentally study tsunami-induced bore pressures exerted on vertical seawalls; a solid vertical wall and a porous vertical seawall that consisted of a perforated front wall and a solid rear wall. Bores with various heights and velocities were generated by using the dam-break method. A porous seawall with 20% porosity of perforated front wall was used in this study. Bore pressures exerted on the solid rear wall and chamber oscillations that occurred in the experiments were also discussed. The experimental results showed that multiple peak pressures were observed during bore run-up phase in the time series of bore impacts. A predictive equation to estimate the maximum bore pressure on a perforated seawall was developed using multiple regression analysis. The proposed equation was also compared with previous empirical formulas.     

Hydrodynamic performance of a T-shaped floating breakwater

The comprehensive utilization of floating breakwaters, specially acting as a supporting structure for offshore marine renewable energy explorations, has received more and more attention recently. Based on linear water-wave theory, the hydrodynamic performance of a T-shaped floating breakwater is semi-analytically investigated through the matched eigenfunction expansion method (MEEM). Auxiliary functions, to speed up the convergence and improve the accuracy in the numerical computations, are introduced to represent the singular behavior of fluid field near the lower salient corners of the structure. The effects of the height and installation position of the vertical screen on the reflection and transmission coefficients, dynamic response and wave forces are examined. It is found that the presence of the screen shifts the resonance frequency of RAO for both surge and pitch modes to the low-frequency area, while has no effect on heave mode. The identical added masses, damping and transmission coefficients can be obtained in the cases where the screen holds the same distance away from the longitudinal central axis of the upper box-type structure. Moreover, a relatively small pitch response can be achieved in a wide wave–frequency range, when the breakwater is Γ-shaped.