A filtering approach for estimating lake water quality from remote sensing data

In this paper we consider the estimation of lake water quality constituent distributions from hyperspectral remote sensing data. We present a computational approach that can be used to assimilate information from mathematical evolution models into data processing. The method is based on a reduced order iterated extended Kalman filter, and a convection–diffusion model is used to describe the movement of the water quality constituents. The performance of the technique is evaluated in a simulation study. The results show that the filter approach with an appropriate evolution model yields estimates that have better spatial and temporal resolutions than those obtained with conventional methods. Furthermore, the use of a feasible evolution model may make it possible to obtain information also on the concentrations in the lower parts of the lake.     

Spatio-temporal characteristics of small-scale wave–current ripples on the Ameland ebb-tidal delta

Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Early Cretaceous sedimentary provenance and structural evolution of the central Sanandaj–Sirjan Zone,Iran: implications for palaeogeographic reconstructions of the northern Neo-Tethyan margin

ABSTRACT

The Early Cretaceous was an important epoch in the evolution of the Earth system in which major tectonic episodes occurred, especially along the Alpine–Himalayan belt. The paucity of reliable palaeogeographic data from the central segment of this geological puzzle, however, hampers the reconstruction of a panoramic view of its Early Cretaceous palaeogeography and geodynamic setting. Here we present multidisciplinary provenance data from Lower Cretaceous strata of the overriding plate of the Neo-Tethyan subduction zone (the Sanandaj–Sirjan Zone; SSZ, of central Iran), including structural, basin-fill evolution, petrographic and geochemical analyses. Sandstone provenance analysis of Lower Cretaceous red beds suggests the occurrence of sub-mature litho-quartzose sandstones attributed to an active continental arc margin in convergent setting predominantly derived from plutonic, quartzose sedimentary and metamorphic rocks exposed in the central SSZ. Weathering indices indicate moderate chemical weathering in the source area which may be related to close source-to-sink relationships or arid climate. Our palaeogeographic reconstructions and original geological mapping indicate that the erosion of uplifted basement rocks exposed in horst blocks provided the sediment sources for the syn-extensional deposition of uppermost Jurassic–lowermost Cretaceous conglomerates and Lower Cretaceous siliciclastic red beds within a continental retro-arc basin during initiation of the ‘Neo-Tethys 2?. The polyphase tectonic reactivation along the principal fault of the study area controlled the syn- and post-extensional tectonostratigraphic evolution that reflect the corresponding mechanical decoupling/coupling along the northern Neo-Tethyan plate margin.     

Geology of the Acropolis prospect,South Australia,constrained by high-precision CA-TIMS ages

Abstract

Acropolis is an Fe-oxide–copper–gold prospect ～20?km from Olympic Dam, South Australia, and marked by near-coincident gravity and magnetic anomalies. Prospective Fe-oxide–apatite?±?sulfide veins occur in Mesoproterozoic and Paleoproterozoic volcanic and granitoid host units beneath unmineralised sedimentary formations. We have produced a geological map and history of the prospect using data from 16 diamond drill holes, including LA-ICPMS and high-precision CA-TIMS ages. The oldest unit is megacrystic granite of the Donington Suite (ca 1850?Ma). A non-conformity spanning ca 250 My separates the Donington Suite and felsic lavas and ignimbrites of the Gawler Range Volcanics (GRV; 1594.03?±?0.68?Ma). The GRV were intruded by granite of the Hiltaba Suite (1594.88?±?0.50?Ma) and felsic dykes (1593.88?±?0.56?Ma; same age as the Roxby Downs Granite at Olympic Dam). The felsic dykes are weakly altered and lack Fe-oxide–apatite–sulfide veins, suggesting that they post-date the main hydrothermal event. If correct, this relationship implies that the main hydrothermal event at Acropolis was ca 1594?Ma and pre-dated the main hydrothermal event at Olympic Dam. The GRV at Acropolis are the same age as the GRV at Olympic Dam and ca 3–7 My older than the GRV exposed in the Gawler Ranges. The gravity and magnetic anomalies coincide with sections through the GRV, Hiltaba Suite and Donington Suite that contain abundant, wide, Fe-oxide veins. The GRV, Hiltaba Suite and Donington Suite are unconformably overlain by the Mesoproterozoic Pandurra Formation or Neoproterozoic Stuart Shelf sedimentary formations. The Pandurra Formation shows marked lateral variations in thickness related to paleotopography on the underlying units and post-Pandurra Formation pre-Neoproterozoic faults. The Stuart Shelf sedimentary formations have uniform thicknesses.

1. KEY POINTS

2. Fe-oxide–apatite?±?sulfide veins are hosted by the Gawler Range Volcanics (1594.03?±?0.68?Ma), the Hiltaba Suite granite (1594.88?±?0.50?Ma) and Donington Suite granite (ca 1850?Ma).

3. The age of felsic dykes (1593.88?±?0.56?Ma) interpreted to be post-mineralisation implies that the main hydrothermal event at Acropolis was ca 1594?Ma.

4. The Gawler Range Volcanics at Acropolis are the same age as the Gawler Range Volcanics at Olympic Dam and ca 3 to 7 My older than the Gawler Range Volcanics exposed in the Gawler Ranges.

     

An Approach to a Visual Language for Convective Fluid Flows: Visualization of Geophysical Simulations

It was known that deep within numbers and binary data from simulations of geophysical convective flows resided various patterns. Two models of convective fluid flows were being considered. One was a model of two-dimensional (768 × 256) air convection with finite Prandtl number of one and Rayleigh number of 10⁸−10¹⁰, and another was a model of three-dimensional (up to 120 × 120 × 90) mantle convection with infinite Prandtl number and Rayleigh number of 10⁶−10⁸. Clearly, phenomena existed which superceded each individual dimensionless computer model to provide a piece of information regarding actual fluid flows. The problem was how to find, prove, and communicate these patterns and phenomena for convection simulations with gigabytes of data. In a search for such an analytical and communicative tool, the alternative of visualization was considered. The need for visualization was recognized and discussed. Then, utilizing both two- and three-dimensional models of high Rayleigh number convection, basic techniques of style and content were developed. Applications of the visualization techniques were designed utilizing IBM’s Data Explorer in order to create communicative images and movies, and after the applications, the problems of data storage and transfer became apparent. Throughout the process though, it became clear how important the language of vision actually could be in the geophysics community. In a field in which words such as plumes and internal waves have in ways replaced mathematics as the basic language for science, there is a need for another resource, another language-the visualization of convective fluid flows.     

A discussion on models of thermals

For a thermal starting from rest, the scales of motion consistent with the initial conditions are given. An alternative time scale based on the motion of the thermal is derived. The anticipated similarity solutions for thermals are summarised and possible qualitative differences between solutions are given. Within this consistent framework previously published laboratory and numerical models of thermals are discussed. Reasons why numerical models have not rigorously demonstrated the existence of a self-similarity solution are considered. Comparisons of all available results show that a single similarity solution valid for all thermals does not exist.     

Non-Darcian groundwater flow in leaky aquifers

Abstract

A simplified method has been developed for solving leaky aquifer non-Darcian flow hydraulics. The principle of volumetric approach is combined with the confined-aquifer, time-dependent drawdown equation in an observation well. The groundwater flow in the leaky aquifer is assumed to obey a non-Darcian flow law of exponential type. The results are obtained in the form of type-curve expressions from which the necessary bundles of curves are drawn for a set of selective non-Darcian flow aquifer parameters. Although application of the methodology appears as rather limited but it provides a scientific contribution and extension of leaky aquifer theory towards nonlinear flow conditions. The methodology developed herein is applied to some actual field data from the eastern sedimentary basin in the Kingdom of Saudi Arabia.     

Application of linear and nonlinear techniques in river flow forecasting in the Kilombero River basin,Tanzania / Application de techniques linéaires et non-linéaires à la prévision des débits dans le bassin de la Rivière Kilombero,en Tanzanie

Abstract

Modelling of the rainfall–runoff transformation process and routing of river flows in the Kilombero River basin and its five sub-catchments within the Rufiji River basin in Tanzania was undertaken using three system (black-box) models—a simple linear model, a linear perturbation model and a linear varying gain factor model—in their linear transfer function forms. A lumped conceptual model—the soil moisture accounting and routing model—was also applied to the sub-catchments and the basin. The HEC-HMS model, which is a distributed model, was applied only to the entire Kilombero River basin. River discharge, rainfall and potential evaporation data were used as inputs to the appropriate models and it was observed that sometimes the system models performed better than complex hydrological models, especially in large catchments, illustrating the usefulness of using simple black-box models in datascarce situations.     

Ages of sediment-hosted Himalayan Pb–Zn–Cu–Ag polymetallic deposits in the Lanping basin,China: Re–Os geochronology of molybdenite and Sm–Nd dating of calcite

The Lanping basin is a significant Pb–Zn–Cu–Ag mineralization belt of the Sanjiang Tethyan metallogenic province in China. Over 100 thrust-controlled, sediment-hosted, Himalayan base metal deposits have been discovered in this basin, including the largest sandstone-hosted Pb–Zn deposit in the world (Jinding), and several Cu ± Ag ± Co deposits (Baiyangping, Baiyangchang and Jinman). These deposits, with total reserves of over 16.0 Mt Pb + Zn, 0.6 Mt Cu, and 7000 t Ag, are mainly hosted in Meso-Cenozoic mottled clastic rocks, and strictly controlled by two Cenozoic thrust systems developed in the western and eastern segments of the Lanping basin.To define the metallogenic history of the study area, we dated nine calcite samples associated with copper sulfides from the Jinman Cu deposit by the Sm–Nd method and five molybdenite samples from the Liancheng Cu–Mo deposit by the Re–Os method. The calcite Sm–Nd age for the Jinman deposit (58 ± 5 Ma) and the molybdenite Re–Os age for the Liancheng deposit (48 ± 2 Ma), together with previously published chronological data, demonstrate (1) the Cu–Ag mineralization in the western Lanping basin mainly occurred in three episodes (i.e., ∼56–54, 51–48, and 31–29 Ma), corresponding to the main- and late-collisional stages of the Indo–Asian orogeny; and (2) the Pb–Zn–Ag (±Cu) mineralization in the eastern Lanping basin lacked precise and direct dating, however, the apatite fission track ages of several representative deposits (21 ± 4 Ma to 32 ± 5 Ma) may offer some constraints on the mineralization age.     

Simulation of the Magnetothermal Instability

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (∂ P/∂z) (∂ ln T/∂ z) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows. J. M. Stone: Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544