中国沿海小黄鱼种群动力学及其空间异质性研究

小黄鱼是我国近海四种最重要的经济鱼类之一，在过去的几十年中小黄鱼种群及其两个地理亚种群经历了巨大的变化。小黄鱼的种群动力学研究，对于维持这一重要渔业的可持续管理是至关重要的。目前仅有的两个小黄鱼种群动力学模型只涵盖了较短的时间，且没有关注其空间差异。本文研究了1968年至2015年间黄渤海和东海小黄鱼的种群动力学模型，包含和不包含空间分层结构的两种贝叶斯模型被用于分析其种群动力学的大尺度空间异质性。本文分析了不同的假设，来研究小黄鱼种群动力学潜在的变化趋势。研究结果表明小黄鱼种群动力学特征具有明显的时间和空间变化。种群的增长速度从20世纪八十年代开始增加，而可捕系数从1981年到2015年增加了两倍多。与黄渤海亚种群相比，东海的小黄鱼亚种群生长速度更快，遭受的捕捞压力也更大。基于最大可持续产量MSY的参考点表明，无论是整个小黄鱼种群还是两个亚种群，近年来都有非常高的过度捕捞风险。因此我国小黄鱼的渔业管理急需更加保守的管理策略，同时考虑其地域差异。本文所用的方法可以应用于其他种类的资源评估和渔业管理，尤其是具有空间异质性和数据有限的种类。     

Pricing Weather Derivatives Index based on Temperature: The Case of Bahir Dar,Ethiopia

In this paper we present a stochastic model for daily average temperature to calculate the temperature indices upon which temperature-based derivatives are written. We propose a seasonal mean and volatility model that describes the daily average temperature behavior using the mean-reverting Ornstein-Uhlenbeck process. We also use higher order continuous-time autoregressive process with lag 3 for modeling the time evolution of the temperatures after removing trend and seasonality. Our model is fitted to 11 years of data recorded, in the period 1 January 2005 to 31 December 2015, Bahir Dar, Ethiopia, obtained from Ethiopia National Meteorological Services Agency. The analytical approximation formulas are used to price heating degree days (HDD) and cooling degree days (CDD) futures. The suggested model is analytically tractable for derivation of explicit prices for CDD and HDD futures and option. The price of the CDD future is calculated, using analytical approximation formulas. Numerical examples are presented to indicate the accuracy of the method. The results show that our model performs better to predict CDD indices.     

西太平洋暖池1990s中期的年代际扩张对夏季登陆中国台风的影响

应用美国联合预警中心(Joint Typhoon Warning Center,JTWC)的台风最佳路径资料、美国国家海洋大气局(National Oceanic and Atmospheric Administration,NOAA)的扩展海表面温度资料以及美国国家环境预报中心(National Centers for Environmental Prediction,NCEP)和美国国家大气科学研究中心(National Center for Atmospheric Research,NCAR)的大气环流场资料,研究了20世纪90年代西太平洋暖池(简称暖池)年代际扩张对西北太平洋台风和登陆中国沿岸台风的影响。研究发现,相比于暖池扩张前期(1965—1992),后期(1993—2013)台风生成在西北太平洋中部区域(10°—20°N,135°—145°E)显著减弱,在10°—20°N,145°—160°E区域和南海北部区域则表现出增多的特点。台风移动路径变异特征呈现为移动进入南海和登陆中国东部沿岸的西行和西北行路径减少,登陆日本的转向型路径增多,同时登陆我国海南岛和东南部沿岸的台风增多。进一步探查这种影响的可能原因发现,与暖池扩张密切相关的太平洋年代际变化引起的纬向环流的变异是西北太平洋中部台风生成减少的主要原因;而南海北部台风生成增多则归因于南海区域局地环流特征的变异。同时,南海北部台风生成增多是登陆我国海南岛和东南沿岸台风增多的主要决定因素。     

SWAN模型对实验室小尺度浅水波浪模拟研究

基于浅水斜坡地形的物理模型试验数据,考察SWAN模型对实验室小尺度浅水波浪的模拟效果,进而检验其浅水项的模拟精度。模拟中采用直接输入初始测点的实测海浪谱进行造波,重点考察浅水中三波相互作用和变浅破碎两个源项,对不同工况下,SWAN模式在水深条件变化下的有效波高、谱平均周期、海浪谱演化的模拟能力进行研究。研究表明:模拟的有效波高较符合实测波浪的增长和衰减,但谱平均周期计算值明显偏小;海浪谱的能量转移机制同实测有较大区别,频谱模拟结果出现高频高估、低频低估现象。对两个源项进行对比分析得出三波相互作用对海浪谱的能量转换影响远大于变浅破碎耗散。想要提高近岸区谱平均周期和海浪谱的模拟精度则SWAN模型中三波非线性项的计算精确度仍需更多研究和改进。     

基于散射计风场数据的台风强度诊断方法——以海洋二号卫星数据为例

本文提出了一种基于散射计风场数据的台风定强方法。该方法定义了台风所在海区的一个圆形区域上的风速平均值来对台风进行定强，并使用海洋二号卫星散射计历史数据对该方法进行了验证。验证结果表明，该台风定强方法可以克服散射计反演台风风速过低的问题。与台风强度监测的历史记录比较的结果显示，在台风强度低于"强台风"这一强度等级时，该方法能够有效估算台风强度。     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of rectangular piles subjected to lateral loads in nonhomogeneous soil using a Winkler model approach

An analytical approach using a Winkler model based on two lateral soil displacement components in a three‐dimensional soil is investigated to provide analytical solutions of horizontal response of a rectangular pile subjected to lateral loads in nonhomogeneous soil. The two lateral displacement components of a soil surrounding the rectangular pile are represented by the Fourier series of displacement potential functions in the elastic three‐dimensional analysis. The lateral stiffness coefficient of the rectangular pile shaft in nonhomogeneous soil is derived from the rocking stiffness coefficient taking into account rocking rotation of a rigid pile shaft. The relationship between horizontal displacement, rotation, moment, and shear force for the rectangular pile subjected to horizontal loads in nonhomogeneous soil is obtainable in the form of the recurrence equation. The formulation of lateral displacement and rotation for a rectangular pile subjected to lateral loads on the pile base in nonhomogeneous soil is proposed by taking into account Mindlin's equation and the equivalent thickness for soil layers in the equivalent elastic method. The difference of lateral behavior between square and circular piles subjected to lateral loads is insignificant. The effect of aspect ratio of the rectangular pile on the lateral behavior is great for the lower stiffness ratio between pile and soil and the larger length–equivalent diameter ratio. The effect of the value of Poisson's ratio of soil on lateral stiffness coefficient is relatively small except Poisson's ratio close to 0.5. The comparison of the results calculated by the current method for a rectangular pile subjected to lateral loads in nonhomogeneous soil has shown good agreement with those obtained from the analytical methods and the finite element method. Copyright © 2015 John Wiley & Sons, Ltd.     

Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment

Assessing catchment runoff response remains a key research frontier because of limitations in current observational techniques to fully characterize water source areas and transit times in diverse geographical environments. Here, we report a study that combines empirical data with modelling to identify dominant runoff processes in a sparsely monitored humid tropical catchment. The analysis integrated isotope tracers into conceptual rainfall–runoff models of varying complexity (from 5 to 11 calibrated parameters) that are able to simulate discharge and tracer concentrations and track the evolving age of stream water exiting the catchment. The model structures can be seen as competing hypotheses of catchment functioning and were simultaneously calibrated against uncertain streamflow gaugings and a 2‐year daily isotope rainfall–runoff record. Comparison of the models was facilitated using global parameter sensitivity analysis and the resulting effect on calibration. We show that a variety of tested model structures reproduced water and tracer dynamics in stream, but the simpler models failed to adequately reproduce both. The resulting water age distributions of the tested models varied significantly with little similarity between the stream water age and stored water age distributions. The sensitivity analysis revealed that only some of the more complex models (from eight parameters) could be better constrained to infer more plausible water age distributions and catchment storage estimates. These models indicated that the age of water stored in the catchment is generally older compared with the age of water fluxes, with evapotranspiration age being younger compared with streamflow. However, the water age distributions followed a similar temporal behaviour dominated by climatic seasonality. Stream water ages increased during the dry season (greater than 1 year) and decreased with increased streamflow (a few weeks old) during the wet season. We further show that the ratios of the streamwater age to stored water age distribution and the water age distribution of actual evapotranspiration to the stored water age distribution from constrained models could potentially serve as useful hydrological indicators of catchment functioning. Copyright © 2016 John Wiley & Sons, Ltd.     

A two‐dimensional coupled hydromechanical discontinuum model for simulating rock hydraulic fracturing

Within the framework of our discontinuous deformation analysis for rock failure algorithm, this paper presents a two‐dimensional coupled hydromechanical discontinuum model for simulating the rock hydraulic fracturing process. In the proposed approach, based on the generated joint network, the calculation of fluid mechanics is performed first to obtain the seepage pressure near the tips of existing cracks, and then the fluid pressure is treated as linearly distributed loads on corresponding block boundaries. The contribution of the hydraulic pressure to the initiation/propagation of the cracks is considered by adding the components of these blocks into the force matrix of the global equilibrium equation. Finally, failure criteria are applied at the crack tips to determine the occurrence of cracking events. Several verification examples are simulated, and the results show that this newly proposed numerical model can simulate the hydraulic fracturing process correctly and effectively. Although the numerical and experimental verifications focus on one unique preexisting crack, because of the capability of discontinuous deformation analysis in simulating block‐like structures, the proposed approach is capable of modeling rock hydraulic fracturing processes. Copyright © 2014 John Wiley & Sons, Ltd.