Landslide dam failure and flood hydraulics. Part II: coupled mathematical modelling

A coupled 2D mathematical modelling study of landslide dam failure and flood is presented, complementing our experimental investigation presented in the companion paper. The model is built upon the shallow water hydrodynamic equations. The governing equations are numerically solved using the total-variation-diminishing version of the second-order weighted-average-flux method along with the HLLC (Harten, Lax and van Leer with Contact wave restored) approximate Riemann solver. Two parameters related to bed-load sediment transport and critical slope stability are calibrated using the measured stage hydrographs from two runs of the flume experiments. The calibrated model is then applied to other independent runs of the experiments featuring different inflow discharges, dam geometry, dam composition and initial breach dimensions. It is found to be able to satisfactorily reproduce the measured stage hydrographs and the widening of initial breach. The experimental observation of the prime role of the inflow discharge and initial breach in dictating the dam failure process and flood is unequivocally resolved, along with the impacts of dam geometry as well the content of cohesive clay and gravel in the dam. Interestingly, the downstream peak discharge and stage of the flood are substantially reduced by initial breach, which clearly exemplifies its role in modulating the flooding.     

A preliminary seismic model for the region of the west-Bohemian earthquake swarms

Summary A layered model of the crustal and upper mantle structure for the region of the West-Bohemian earthquake swarms is proposed. This model represents a compilation of the results of deep seismic soundings, surface-wave investigations, and some local seismic measurements in the western regions of the Bohemian Massif.     

Hydrocode Simulations of Chicxulub Crater Collapse and Peak-Ring Formation

We use hydrocode modeling to investigate dynamic models for the collapse of the Chicxulub impact crater. Our aim is to integrate the results from numerical simulations with kinematic models derived from seismic reflection and wide-angle velocity data to further our understanding of the formation of large impact craters. In our simulations, we model the collapse of a 100-km diameter, bowl-shaped cavity formed in comprehensively fractured crustal material. To facilitate wholesale collapse, we require that the strength of the target be significantly weakened. In the present model, we achieve this using acoustic fluidization, where strong vibrations produced by the expanding shock wave cause extreme pressure fluctuations in the target. At times and positions where the overburden pressure is sufficiently counteracted, the frictional resistance is reduced, enabling the rock debris to flow. Our simulations produce a collapsed crater that contains most of the features that we observe in the seismic data at Chicxulub. In particular, we observe a topographic peak ring, formed as material that is originally part of the central uplift collapses outward and is thrust over the inwardly collapsing transient crater rim. This model for peak-ring generation has not been previously demonstrated by numerical simulations and predicts that the peak ring is composed of deeply derived material and that the stratigraphy within the peak ring is overturned.     

Methods of computing the partial derivatives of dispersion curves

Summary Three methods of computing the partial derivatives of dispersion curves are briefly described and compared: numerical differentiation, variational methods and an implicit function method.     

Forces at the Kerr singularity

In a classical-like language introduced in a previous paper, the fact that the ring singularity of the Kerr spacetime cannot be left/reached by timelike geodesic which is not bound to the equatorial plane is interpreted as a result of the interplay of the gravitational and – mainly – the centrifugal forces.     

Axisymmetric annulus convection at unit Prandtl number

Abstract

We examine the role played in annulus flows by mechanisms dependent upon the Prandtl number, σ. Solutions are obtained at σ = 1 for both the real annulus system and for the hypothetical “free annulus” system (free slip lateral boundaries). These solutions are compared with previously obtained solutions at σ = 7.

In the free annulus, the solution at σ = 1 differs radically from that at σ = 7. The σ = 1 solution appears to be essentially a finite amplitude mode due to Solberg instability whereas the solution at σ = 7 manifests a flow caused by the diffusive overturning mechanism.

The variation with σ of the real annulus flow is not so fundamental but some differences in the dynamical structures are noted.     

Evidence for an impact‐induced magnetic fabric in Allende,and exogenous alternatives to the core dynamo theory for Allende magnetization

We conducted a paleomagnetic study of the matrix of Allende CV3 chondritic meteorite, isolating the matrix's primary remanent magnetization, measuring its magnetic fabric and estimating the ancient magnetic field intensity. A strong planar magnetic fabric was identified; the remanent magnetization of the matrix was aligned within this plane, suggesting a mechanism relating the magnetic fabric and remanence. The intensity of the matrix's remanent magnetization was found to be consistent and low (~6 μT). The primary magnetic mineral was found to be pyrrhotite. Given the thermal history of Allende, we conclude that the remanent magnetization was formed during or after an impact event. Recent mesoscale impact modeling, where chondrules and matrix are resolved, has shown that low‐velocity collisions can generate significant matrix temperatures, as pore‐space compaction attenuates shock energy and dramatically increases the amount of heating. Nonporous chondrules are unaffected, and act as heat‐sinks, so matrix temperature excursions are brief. We extend this work to model Allende, and show that a 1 km/s planar impact generates bulk porosity, matrix porosity, and fabric in our target that match the observed values. Bimodal mixtures of a highly porous matrix and nominally zero‐porosity chondrules make chondrites uniquely capable of recording transient or unstable fields. Targets that have uniform porosity, e.g., terrestrial impact craters, will not record transient or unstable fields. Rather than a core dynamo, it is therefore possible that the origin of the magnetic field in Allende was the impact itself, or a nebula field recorded during transient impact heating.     

Tidal amplification and along‐strike process variability in a mixed‐energy paralic system prograding onto a low accommodation shelf,Edgeøya,Svalbard

The study describes the depositional development and sediment partitioning in a prograding paralic Triassic succession. The deposits are associated with the advance of large prism‐scale clinoforms across a shallower platform area. Approaching the platform, the limited accommodation and associated relative higher rates of deposition generated straighter clinoforms with lower foreset angles. The vertical restriction across the platform is interpreted to have amplified the tidal signature. Sediment was redistributed from the coast into increasingly sandy delta‐front deposits, compared to offshore equivalents. The deposits comprise extensive compound dune fields of amalgamated and increasingly clean sandbodies up‐section. Rapid deposition of significant amounts of sand led to differential subsidence and growth‐faulting in the delta front, with downthrown fault blocks further amplifying the tidal energy through funnelling. A mixed‐energy environment created along‐strike variability along the delta front with sedimentation governing process‐regime. Areas of lower sedimentation were reworked by wave and storm‐action, whereas high sedimentation rates preserved fluvially dominated mouth bars. A major transgression, however, favoured tidally dominated deposits also in these areas, attributed to increasing rugosity of the coastline. Formation of an extensive subaqueous platform between the coast and delta front dampened incoming wave energy, and tidally dominated deposits dominate the near‐shore successions. Meanwhile formation of wave‐built sand‐bars atop the platform attest to continued wave influence. The strong tidal regime led to the development of a heterolithic near‐shore tidally dominated channel system, and sandier fluvial channels up‐river. The highly meandering tidal channels incising the subaqueous platform form kilometre wide successions of inclined heterolithic stratification. The fluvially dominated channels which govern deposition on the delta plain are narrower and slightly less deep, straighter, generally symmetric and filled with cleaner sands. This study provides important insight into tidal amplification and sand redistribution during shallowing on a wide shelf, along with along‐strike process‐regime variability resulting from variations in sediment influx.     

Stability of co-orbital ring material with applications to the Janus-Epimetheus system

An analytical model that describes the evolution of ring particles that are co-orbital with two larger bodies on near-circular and near-planar orbits has been formulated. This can be used to estimate the lifetime of the particles within the ring. All the examples investigated, such as the Janus-Epimetheus (JE) system, indicate that the particles should be removed from the co-orbital region within half a synodic period (∼4 years for JE). Numerical modelling confirms the predictions of the model. When the masses are on eccentric orbits the particles remain within the co-orbital system for longer. Our results suggest that the ring associated with Janus and Epimetheus must be continually fed with material, probably by meteoroid impacts on the two satellites.