Intercomparison of coastal area morphodynamic models

Five different coastal area morphodynamic models have been set up to run on the same offshore breakwater layout and an intercomparison carried out on the hydrodynamic and morphodynamic output produced by each scheme. In addition, the predicted morphodynamics was checked against available laboratory and field data.It is concluded that the models are capable of producing realistic estimates for the dominant morphodynamic features associated with offshore breakwaters. Coupling of the wave, current and sediment transport components of each scheme is shown to yield bathymetry which attains a state of equilibrium, unlike models which are based on the initial transport field only.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

遥感蚀变信息检测中背景与干扰问题研究

在不同自然景观区的遥感蚀变信息检测中,背景与干扰是2个复杂而又容易混淆的研究对象。本文结合光谱混合模型理论与随机变量概率分布相关理论,深入分析了遥感数据的光谱空间几何结构特征(重点剖析二维散点图),系统研究了背景与干扰的关系、复杂度以及干扰地物类型的划分等问题。这些问题的探讨为遥感蚀变信息检测的最佳方案设计和蚀变信息潜在性评价与准确提取提供了技术思路。     

曲流点坝内部剩余油形成与分布规律物理模拟

建立了曲流点坝物理模型,模拟地下曲流点坝储层进行驱替实验。观察注入剂的运动规律,同时,在模型上布置了测量电极,以双电极法测量驱替过程中垂向电阻变化,反映注入剂波及高度的变化规律。分析了侧积泥岩夹层的建筑结构对注入剂在点坝砂体中的流动规律,以及对剩余油的形成分布的影响。通过实验研究证实,在曲流点坝内部,侧积泥岩夹层对注入剂具有强烈的遮挡作用,受其影响形成侧积体差异型剩余油、重力型剩余油以及压力异常型剩余油三种类型剩余油,主要分布在点坝的中上部。在实际生产中,利用水平井钻遇点坝上部进行开采是主要的手段。     

Glacio-lacustrine stratigraphy,aquifer characterization and contaminant transport: a case study in South Lake Tahoe,California, USA

A hydrogeologic model that has been used by many researchers and consultants to describe an area of South Lake Tahoe, California, USA impacted by MTBE contamination describes a relatively homogeneous unconfined aquifer comprised of poorly sorted glacial outwash deposits, within which water-supply wells are able to exert significant alteration in natural groundwater flow. A re-examination of the area’s hydrogeology is presented, which supports a layered heterogeneous aquifer system constructed of alternating fine and coarser-grained glacio-lacustrine depositional units. This re-evaluation was accomplished through a review of lithologic logs across an area of approximately 1 km², combined with observations of significant hydraulic head differences and knowledge of the depositional environments controlled by Pleistocene Lake Tahoe high stands. Many of the fine-grained units observed at depths from 6 to 15 m, although relatively thin, are generally continuous and serve as significant barriers to groundwater flow. The vertical migration of contamination across these fine-grained units to deeper groundwaters was facilitated by cross-screened monitoring wells installed as part of site investigation activities. This conclusion highlights the importance of geologic characterization and proper monitoring well construction at contaminated site investigations.

---

Electronic supplementary material   The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Comparison of numerical models of two debris flows in the Cortina d’ Ampezzo area,Dolomites, Italy

The accurate prediction of runout distances, velocities and the knowledge of flow rheology can reduce the casualties and property damage produced by debris flows, providing a means to delineate hazard areas, to estimate hazard intensities for input into risk studies and to provide parameters for the design of protective measures. The application of most of models that describe the propagation and deposition of debris flow requires detailed topography, rheological and hydrological data that are not always available for the debris-flow hazard delineation and estimation. In the Cortina d’Ampezzo area, Eastern Dolomites, Italy, most of the slope instabilities are represented by debris flows; 325 debris-flow prone watersheds have been mapped in the geomorphological hazard map of this area. We compared the results of simulations of two well-documented debris flows in the Cortina d’Ampezzo area, carried on with two different single-phase, non-Newtonian models, the one-dimensional DAN-W and the two-dimensional FLO-2D, to test the possibility to simulate the dynamic behaviour of a debris flow with a model using a limited range of input parameters. FLO-2D model creates a more accurate representation of the hazard area in terms of flooded area, but the results in terms of runout distances and deposits thickness are similar to DAN-W results. Using DAN-W, the most appropriate rheology to describe the debris-flow behaviour is the Voellmy model. When detailed topographical, rheological and hydrological data are not available, DAN-W, which requires less detailed data, is a valuable tool to predict debris-flow hazard. Parameters obtained through back-analysis with both models can be applied to predict hazard in other areas characterized by similar geology, morphology and climate.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

The “Zero-order Model” revisited: Paul Schindler's influence on the development of trace metal scavenging models

In 1975 Paul Schindler produced the first oceanic trace metal scavenging model to explicitly include the role of surface chemistry as a control on trace metal water column residence times. The eighteen years that have elapsed since the publication of Schindler's seminal paper have seen the development of a variety of oceanic scavenging models; yet, the fundamental insight of his Zero-order Model remains the benchmark. This paper describes the role of Paul Schindler's work on surface chemistry in providing a framework for the current generation of trace element scavenging models.