Numerical models of convection in the earth's mantle

Calculations on convection in Newtonian fluids at large Rayleigh numbers show behaviour in general agreement with geophysical observations. The calculated values of surface velocities, regional gravity and topographic anomalies all lie within the range of values found at the earth's surface. The sign of the gravity anomaly over rising fluid is positive because the gravitational effect of the surface deformation is greater than that of the density deficit. This result appears to agree with geophysical observations and should permit the flow within the mantle to be followed in some detail.     

Is Bangkok at risk of marine flooding? Evidence relating to the historical floods of AD 1785 and 1983

Natural Hazards - Thailand’s low-lying capital city Bangkok is a dense metropolis, one of Asia’s rapidly growing ‘megacities’, and home to over 10 million...     

A Synthetic Haematite Reference Material for LA‐ICP‐MS U‐Pb Geochronology and Application to Iron Oxide‐Cu‐Au Systems

In both nature and synthetic experiments, the common iron oxide haematite (α‐Fe₂O₃) can incorporate significant amounts of U into its crystal structure and retain radiogenic Pb over geological time. Haematite is a ubiquitous component of many ore deposit types and, therefore, represents a valuable hydrothermal mineral geochronometer, allowing direct constraints to be placed on the timing of ore formation and upgrading. However, to date, no suitable natural haematite reference material has been identified. Here, a synthetic haematite U‐Pb reference material (MR‐HFO) is characterised using LA‐ICP‐MS and ID‐TIMS. Centimetre‐scale ‘chips’ of synthesised α‐Fe₂O₃ were randomly microsampled via laser ablation‐extraction and analysed using ID‐TIMS. Reproducible U/Pb and Pb/Pb measurements were obtained across four separate chips (n = 13). Subsequently, an evaluation of the suitability MR‐HFO in constraining U‐Pb data via LA‐ICP‐MS is presented using a selection of natural samples ranging from Cenozoic to Proterozoic in age. The MR‐HFO normalised U‐Pb ratios are more concordant and ages more accurate versus the same LA‐ICP‐MS spot analyses normalised to zircon reference material, when compared with independently acquired ID‐TIMS data from the same natural haematite grains. Results establish MR‐HFO as a suitable reference material for LA‐ICP‐MS haematite U‐Pb geochronology.     

Shallow velocity structure along the Hirapur–Mandla profile using traveltime inversion of wide-angle seismic data, and its tectonic implications

In order to investigate the velocity structure, and hence shed light on the related tectonics, across the Narmada–Son lineament, traveltimes of wide-angle seismic data along the 240 km long Hirapur–Mandla profile in central India have been inverted. A blocky, laterally heterogeneous, three-layer velocity model down to a depth of 10 km has been derived. The first layer shows a maximum thickness of the upper Vindhyans (4.5 km s⁻¹ ) of about 1.35 km and rests on top of normal crystalline basement, represented by the 5.9 km s⁻¹ velocity layer. The anomalous feature of the study is the absence of normal granitic basement in the great Vindhyan Graben, where lower Vindhyan sediments (5.3 km s⁻¹ ) were deposited during the Precambrian on high-velocity (6.3 km s⁻¹ ) metamorphic rock. The block beneath the Narmada–Son lineament represents a horst feature in which high-velocity (6.5 km s⁻¹ ) lower crustal material has risen to a depth of less than 2 km. South of the lineament, the Deccan Traps were deposited on normal basement during the upper Cretaceous period and attained a maximum thickness of about 800 m.     

Unstructured mesh generation and landcover-based resistance for hydrodynamic modeling of urban flooding

Urban flood inundation modeling with a hydrodynamic flow solver is addressed in this paper, focusing on strategies to effectively integrate geospatial data for unstructured mesh generation, building representation and flow resistance parameterization. Data considered include Light Detection and Ranging (LiDAR) terrain height surveys, aerial imagery and vector datasets such as building footprint polygons. First, a unstructured mesh-generation technique we term the building-hole method (BH) is developed whereby building footprint data define interior domain boundaries or mesh holes. A wall boundary condition depicts the impact of buildings on flood hydrodynamics. BH provides an alternative to the more commonly used method of raising terrain heights where buildings coincide with the mesh. We term this the building-block method (BB). Application of BH and BB to a flooding site in Glasgow, Scotland identifies a number of tradeoffs to consider at resolutions ranging from 1 to 5 m. At fine resolution, BH is shown to be similarly accurate but execute faster than BB. And at coarse resolution, BH is shown to preserve the geometry of buildings and maintain better accuracy than BB, but requires a longer run time. Meshes that ignore buildings completely (no-building method or NB) also support surprisingly good flood inundation predictions at coarse resolution compared to BH and BB. NB also supports faster execution times than BH at coarse resolution because the latter uses localized refinements that mandate a greater number of computational cells. However, with mesh refinement, NB converges to a different (and presumably less-accurate) solution compared to BH and BB. Using the same test conditions, Hunter et al. [Hunter NM, Bates PD, Neelz S, Pender G, Villanueva I, Wright NG, Liang D, et al. Benchmarking 2D hydraulic models for urban flood simulations. ICE J Water Manage 2008;161(1):13–30] compared the performance of dynamic-wave and diffusive-wave models and reported that diffusive-wave models under-predicted the longitudinal penetration of the flood zone due to important inertial effects. Here, we find that a relatively coarse-mesh implementation of a dynamic-wave model suffers from the same drawback because of numerical diffusion. This shows that whether diffusion is achieved through the mathematics or numerics, the effect on flood extent is similar. Finally, several methods of distributing resistance parameters (e.g., Manning n) across the Glasgow site were evaluated including methods that utilize aerial imagery-based landcover classification data, MasterMap^® landcover classification data and LiDAR-based feature height data (e.g., height of shrubs or hedges). Results show that landcover data is more important than feature height data in this urban site, that shadows in aerial imagery can cause errors in landcover classification which degrade flood predictions, and that aerial imagery offers a more detailed mapping of trees and bushes than MasterMap^® which can locally impact depth predictions but has little impact on flood extent.     

The electromagnetic response of a horizontal electric dipole buried in a multi‐layered earth

The electromagnetic response of a horizontal electric dipole transmitter in the presence of a conductive, layered earth is important in a number of geophysical applications, ranging from controlled‐source audio‐frequency magnetotellurics to borehole geophysics to marine electromagnetics. The problem has been thoroughly studied for more than a century, starting from a dipole resting on the surface of a half‐space and subsequently advancing all the way to a transmitter buried within a stack of anisotropic layers. The solution is still relevant today. For example, it is useful for one‐dimensional modelling and interpretation, as well as to provide background fields for two‐ and three‐dimensional modelling methods such as integral equation or primary–secondary field formulations. This tutorial borrows elements from the many texts and papers on the topic and combines them into what we believe is a helpful guide to performing layered earth electromagnetic field calculations. It is not intended to replace any of the existing work on the subject. However, we have found that this combination of elements is particularly effective in teaching electromagnetic theory and providing a basis for algorithmic development. Readers will be able to calculate electric and magnetic fields at any point in or above the earth, produced by a transmitter at any location. As an illustrative example, we calculate the fields of a dipole buried in a multi‐layered anisotropic earth to demonstrate how the theory that developed in this tutorial can be implemented in practice; we then use the example to examine the diffusion of volume charge density within anisotropic media—a rarely visualised process. The algorithm is internally validated by comparing the response of many thin layers with alternating high and low conductivity values to the theoretically equivalent (yet algorithmically simpler) anisotropic solution, as well as externally validated against an independent algorithm.     

Insights into the Microwave Instruments Onboard the Fengyun 3DSatellite: Data Quality and Assimilation in theMet Office NWP System

正1Met Office, Exeter EX1 3PB, UK2National Satellite Meteorological Center, China Meteorological Administration, Beijing 100081, China