An extreme wave event in eastern Yucatán,Mexico: Evidence of a palaeotsunami event during the Mayan times

The Yucatán Peninsula, Mexico, has typically been considered a tectonically stable region with little significant seismic activity. The region though, is one that is regularly affected by hurricanes. A detailed survey of ca 100 km of the eastern Yucatán and Cozumel coast identified the presence of ridges containing individual boulders measuring >1 m in length. The boulder ridges reach 5 m in height and their origin is associated with extreme wave event activity. Previously modelled tsunami waves from known seismically active zones in the region (Muertos Thrust Belt and South Caribbean Deformed Belt) are not of sufficient scale in the area of the Yucatán Peninsula to have produced the boulder ridges recorded in this study. The occurrence of hurricanes in this region is more common, but two of the most destructive (Hurricane Gilbert 1988 and Hurricane Wilma 2005) produced coastal waves too small to have created the ridges recorded here. In this paper, a new tsunami model with a source area located on the Motagua/Swan Island Fault System has been generated that indicates a tsunami event may have caused the extreme wave events that resulted in the deposition of the boulder ridges.     

Rapid shoreline flooding enhances water turbidity by sediment resuspension: An example in a large Tibetan lake

Rapid water level rise due to climate change has the potential to remobilize loose sediments along shorelines and increase the turbidity of nearshore waters, thereby impacting water quality and aquatic ecosystem health. Siling Lake is one of the largest and most rapidly expanding lakes on the Tibetan Plateau. Between 2000 and 2017, this lake experienced an increase in water level of about 8 m and a doubling in water turbidity. Here, using this lake as a study site, we used a wave model and high-resolution remote sensing of turbidity (Landsat-8) to assess the potential connection between water-level rise, enhanced wind-driven sediment resuspension and water turbidity. Our analysis revealed that strong bottom shear stresses triggered by wind-generated waves over newly flooded areas were related to an increase in water turbidity. The spatial variability of Siling Lake turbidity showed a strong dependence on local wind characteristics and fetch. Two factors combined to drive the increase in turbidity: (1) high wave energy leading to high bottom shear stresses, and (2) flooding of unvegetated shallow areas. Using a new relationship between wave energy and turbidity developed here, we expect the increase in turbidity of Siling Lake to taper off in the near future due to the steep landscape surrounding the lake that will prevent further flooding. Our results imply that rising water levels along the coast are not only expected to influence terrestrial ecosystems but could also change water quality. The methodology presented herein could be applied to other shorelines affected by a rapid increase in water level. © 2020 John Wiley & Sons, Ltd.     

On the morphology of radial sand ridges

We use a hydrodynamic model applied to an idealized fan-shaped basin to explore the morphology and dynamics of radial sand ridges in a convergent coastal system. A positive morphological feedback between channel incision and flow redistribution is responsible for the formation of the channel-ridge pattern. The selection mechanism of bottom wavelength is associated with flow concentration in the deeper part of the channels. Our results are compared to sediment and hydraulic dynamics in the radial sand ridges (RSRs) in China. In a convergent, sloping basin the tangentially averaged tidal velocity peaks at 47 km from the apex. This distance is similar to the arc distance, 62 km, where the RSRs are most incised. An offshore shift in tidal phase results in stronger flows near the north coastline, explaining the presence of asymmetric channel patterns. A numerical stability analysis indicates that small radial oscillations with a wavelength of 10° to 15° maximize the velocity in the troughs. This oscillation wavelength also emerges in the RSRs, which display a peak in spectral energy at a radial wavelength between 25° to 37.5°. High-resolution numerical simulations in the RSRs confirm that flow concentration occurs in the deeper part of the channels, keeping them flushed. We therefore conclude that the RSRs display morphometric characteristics similar to other tidal incisions, like tidal inlets and intertidal channels. This result further supports the dominant role of tidal prism and related peak velocities in incising coastal landscapes. © 2020 John Wiley & Sons, Ltd.     

Back‐barrier flooding by storm surges and overland flow

Low‐elevation areas within a sandy barrier island are subject to flooding via saturation overland flow following moderate storm surges and rainfall events. Using a high resolution topographic survey and simple hydrology models, we estimate the discharge and velocities from storm surge return flow and saturation overland flow. Results show that return flow velocities are of the same magnitude as the critical velocity necessary to mobilize sand when a hydraulic connection between the watershed and back‐barrier bay is present. Storms of moderate strength and rainfall intensity may be sufficient to keep the return channels open within the back‐barrier, thus providing natural conduits for water exchange from overwash events during extreme storm surges triggered by hurricanes. Copyright © 2011 John Wiley & Sons, Ltd.     

Earthquake dynamics on circular faults: a review 1970–2015

We review the development of earthquake dynamics taken from the point of view of the origin of seismic radiation instead of the detailed study of rupture propagation on complex surfaces. Many features of seismic radiation can be explained by simple models that serve as elementary canonical problems. Some of these properties are very well known like the fact that at low frequencies, the seismic spectrum is proportional to the seismic moment. At high frequencies, on the other hand, radiation is generated by the motion of the rupture front, in particular stopping phases and geometrical obstacles (barriers). A rupture front moving at constant speed does not produce far-field radiation. For many practical applications, for determining source size and stress drop, for example, it is not necessary to determine geometrical details of the rupture. For such cases, a simple circular crack model is quite sufficient. An improvement on this method is to do dynamic inversion on simple, elliptical-shaped sources and letting the rupture start arbitrarily from a point on the fault. This problem can be solved nowadays with finite differences and a variety of inversion techniques.     

Dynamics of a fringe mangrove forest detected by Landsat images in the Mekong River Delta,Vietnam

Mangrove forests dominate many tropical coastlines and are one of the most bio‐diverse and productive environments on Earth. However, little is known of the large‐scale dynamics of mangrove canopies and how they colonize intertidal areas. Here we focus on a fringe mangrove forest located in the Mekong River Delta, Vietnam; a fast prograding shoreline where mangroves are encroaching tidal flats. The spatial and temporal evolution of the mangrove canopy is studied using a time series of Landsat images spanning two decades as well as Shuttle Radar Topography Mission (SRTM) elevation data. Our results show that fast mangrove expansion is followed by an increase in Normalized Difference Vegetation Index (NDVI) in the newly established canopy. We observe three different dynamics of the mangrove fringe: in the southwest part of the fringe, near a deltaic distributary where the fringe boundary is linear, the canopy expands uniformly on the tidal flats with a high colonization rate and high NDVI values. In the northeast part of the fringe, near another distributary, the canopy expands at a much lower rate with low NDVI values. In the fringe center, far from the river mouths, the fringe boundary is highly irregular and mangroves expansion in characterized by sparse vegetated patches displaying low NDVI values. We ascribe these different dynamics to wave action and southwest longshore transport triggered by energetic northeasterly monsoons during winter. We further link the large‐scale dynamics of the fringe to small‐scale physical disturbances (waves, erosion and deposition) that might prevent the establishment of mangrove seedlings. Based on these results, we include mangrove encroachment in an already published conceptual model of progradation of the Mekong River Delta. We conclude that high NDVI values and a constantly linear vegetation–water interface are indicative of stable mangrove canopies undergoing fast expansion, probably triggered by sediment availability at the shore. Our results can be applied more generally to mangrove forests growing in minerogenic and high tidal range environments with high sediment inputs. Copyright © 2016 John Wiley & Sons, Ltd.     

Montecarlo DLA-type simulations of wetting effects in fluid displacement in porous media

In this work, we report diffusion-limited aggregation (DLA)-type Montecarlo computations of a stochastic model of displacement of a viscous fluid by another that preferentially wets a porous medium, for the case when both fluids are immiscible in the absence of buoyancy forces. The model has the aim to simulate cooperative invasion processes found in experiments of immiscible wetting displacement. The model considers the nonlocal effects of the Laplacian pressure field and the capillary forces via hydrodynamic equations in the Darcy regime with a boundary condition for the pressure at the interface. The boundary condition contains two different types of disorder: the capillary term, which constitutes an additive random disorder, and a term containing an effective random surface tension, which couples to a curvature (it constitutes a multiplicative random term that carries nonlocal information of the whole pressure). We generate different displacement patterns for different setting of the parameters of the model. We analyze these patterns by studying the scaling properties of the interface that separate the two fluids and calculating the fractal dimension of the interface. The results show the existence of three distinct regimes of scaling. One regime at the smallest-length scales is due to the multiplicative random disorder together with the nonlocal coupling; it reveals itself in a roughness exponent α ≈ 0.80. Additionally, we find a DLA-type scaling regime with a roughness exponent α ≈ 0.60 at the largest scales and intermediate scaling regime with α ≈ 0.70 corresponding to invasion percolation with trapping. Each regime has definite scaling ranges that depend on the capillary number and the relative wetting tendency of the fluids. The behavior of the fractal dimensions of the interfaces of the aggregates constitutes a further confirmation of the existence of three scaling regimes and the multi-self-affinity of the perimeter of the interface boundaries.