Holocene alpine glaciation inferred from lacustrine sediments on northeastern Baffin Island,Arctic Canada

With accelerated melting of alpine glaciers, understanding the future state of the cryosphere is critical. Because the observational record of glacier response to climate change is short, palaeo‐records of glacier change are needed. Using proglacial lake sediments, which contain continuous and datable records of past glacier activity, we investigate Holocene glacier fluctuations on northeastern Baffin Island. Basal radiocarbon ages from three lakes constrain Laurentide Ice Sheet retreat by ca. 10.5 ka. High sedimentation rates (0.03 cm a^?1) and continuous minerogenic sedimentation throughout the Holocene in proglacial lakes, in contrast to organic‐rich sediments and low sedimentation rates (0.005 cm a^?1) in neighbouring non‐glacial lakes, suggest that glaciers may have persisted in proglacial lake catchments since regional deglaciation. The presence of varves and relatively high magnetic susceptibility from 10 to 6 ka and since 2 ka in one proglacial lake suggest minimum Holocene glacier extent ca. 6–2 ka. Moraine evidence and proglacial and threshold lake sediments indicate that the maximum Holocene glacier extent occurred during the Little Ice Age. The finding that glaciers likely persisted through the Holocene is surprising, given that regional proxy records reveal summer temperatures several degrees warmer than today, and may be due to shorter ablation seasons and greater accumulation‐season precipitation. Copyright © 2009 John Wiley & Sons, Ltd.     

Assessing uncertainties in a conceptual water balance model using Bayesian methodology / Estimation bayésienne des incertitudes au sein d’une modélisation conceptuelle de bilan hydrologique

Abstract

Abstract The aim of this study was to estimate the uncertainties in the streamflow simulated by a rainfall–runoff model. Two sources of uncertainties in hydrological modelling were considered: the uncertainties in model parameters and those in model structure. The uncertainties were calculated by Bayesian statistics, and the Metropolis-Hastings algorithm was used to simulate the posterior parameter distribution. The parameter uncertainty calculated by the Metropolis-Hastings algorithm was compared to maximum likelihood estimates which assume that both the parameters and model residuals are normally distributed. The study was performed using the model WASMOD on 25 basins in central Sweden. Confidence intervals in the simulated discharge due to the parameter uncertainty and the total uncertainty were calculated. The results indicate that (a) the Metropolis-Hastings algorithm and the maximum likelihood method give almost identical estimates concerning the parameter uncertainty, and (b) the uncertainties in the simulated streamflow due to the parameter uncertainty are less important than uncertainties originating from other sources for this simple model with fewer parameters.     

Statistical estimation of a shock center: Slate Islands astrobleme

Statistical models are proposed for the distribution of directions in three dimensions that are thought to point toward a single source. These models are based on the Fisher distribution. The method of maximum likelihood is used to obtain an estimate of position of the source and to provide corre-sponding confidence regions. When applied to shatter cone data from the Slate Islands, northern Lake Superior, the method yields estimates comparable with those obtained by Stesky and Halls (1982), as well as statistically valid confidence regions.     

重磁解释中健全的极大似然法及三维反演

本文介绍了JOAO.B.C.SILVA提出的重磁解释中健全的(Robust)极大似然反演方法的基本原理,讨论了重力二度体非线性几何参数的反演情况,并将该方法扩展到三维重、磁非线性几何参数反演中。重点讨论了场值中有较大随机干扰、有不同类型的地质噪音源影响下的反演结果。理论反演结果表明,该方法有比最小二乘反演方法更接近实际的解,从而提高了解的可靠性,对于重磁局部异常或孤立异常有较好的反演效果。     

Seismic recognition of local inhomogeneities

A simple method to contour local inhomogeneities using seismic data is proposed. It formalizes an approximate inversion method which is based on the interpretation of local inhomogeneities as making the differences between an actual seismic data set and a previous reference model. It uses the optimal statistical criteria of parameter estimation and recognition and the ray representation of the waves spreading. Any combination of direct, reflected and/or other types of waves may be used as the database. Inhomogeneities, having a size two times above the wavelength of the seismic waves, can be resolved. Laboratory experiments, using ultrasonic waves and analysis of data from field experiments, confirmed the theoretical results. The method can be used to search for ore bodies, kimberlite cubes, oiltraps, etc.     

方差—协方差分量极大似然估计的通用公式

本文由概括平差函数模型出发，按极大似然做估计原则导出了适用于所有平差函数模型的方差分量估计的通用公式，由Ｋ．Ｋｕｂｉｋ和Ｃ．Ｒ．Ｋｏｃｈ所导出的两个公式都是它的特例。     

线性规划在徐州地下水资源开发中的应用

本文从地下水系统的观点出发,对徐州地区中寒武统—中奥陶统碳酸盐岩地下水系统进行了分析。并应用线性规划的原理和方法对徐州市西郊四个水源地地下水进行了水文地质计算。在规定水位允许降深的约束条件下,求得至2000年,该区地下水最大开采量为14万米~3/天及各水源地相应取水量的最佳分配。文中采用特定条件图解法计算多边界、多主井抽水干扰影响的方法,具有通用意义。     

Application of bivariate extreme value distribution to flood frequency analysis: a case study of Northwestern Mexico

In Mexico, poverty has forced people to live almost on the water of rivers. This situation along with the occurrence of floods is a serious problem for the local governments. In order to protect their lives and goods, it is very important to account with a mathematical tool that may reduce the uncertainties in computing the design events for different return periods. In this paper, the Logistic model for bivariate extreme value distribution with Weibull-2 and Mixed Weibull marginals is proposed for the case of flood frequency analysis. A procedure to estimate their parameters based on the maximum likelihood method is developed. A region in Northwestern Mexico with 16 gauging stations has been selected to apply the model and regional at-site quantiles were estimated. A significant improvement occurs, measured through the use of a goodness-of-fit test, when parameters are estimated using the bivariate distribution instead of its univariate counterpart. Results suggest that it is very important to consider the Mixed Weibull distribution and its bivariate option when analyzing floods generated by a␣mixture of two populations.     

Hybrid Estimation of Semivariogram Parameters

Two widely used methods of semivariogram estimation are weighted least squares estimation and maximum likelihood estimation. The former have certain computational advantages, whereas the latter are more statistically efficient. We introduce and study a “hybrid” semivariogram estimation procedure that combines weighted least squares estimation of the range parameter with maximum likelihood estimation of the sill (and nugget) assuming known range, in such a way that the sill-to-range ratio in an exponential semivariogram is estimated consistently under an infill asymptotic regime. We show empirically that such a procedure is nearly as efficient computationally, and more efficient statistically for some parameters, than weighted least squares estimation of all of the semivariogram’s parameters. Furthermore, we demonstrate that standard plug-in (or empirical) spatial predictors and prediction error variances, obtained by replacing the unknown semivariogram parameters with estimates in expressions for the ordinary kriging predictor and kriging variance, respectively, perform better when hybrid estimates are plugged in than when weighted least squares estimates are plugged in. In view of these results and the simplicity of computing the hybrid estimates from weighted least squares estimates, we suggest that software that currently estimates the semivariogram by weighted least squares methods be amended to include hybrid estimation as an option.     

Maximum-likelihood ambiguity resolution based on Bayesian principle

 Based on the Bayesian principle and the fact that GPS carrier-phase ambiguities are integers, the posterior distribution of the ambiguities and the position parameters is derived. This is then used to derive the maximum posterior likelihood solution of the ambiguities. The accuracy of the integer ambiguity solution and the position parameters is also studied according to the posterior distribution. It is found that the accuracy of the integer solution depends not only on the variance of the corresponding float ambiguity solution but also on its values. Received: 27 July 1999 / Accepted: 22 November 2000