Controls on basalt terrace formation in the eastern Lesotho highlands

The objective of this study is to ascertain the relative importance of lithological controls and geomorphological processes in the development of Drakensberg basalt terraces. Various hypotheses for terrace formation are considered, including geological controls, macroscale geomorphology, and climatic–geomorphological controls. The variations in strength and relative age differences for scarp surface exposures on two slopes of varying aspect are determined, so that a comparison can be made between various scarp outcrops and relative rates of weathering. Scarp outcrops were measured for their rock mass strength, rock surface roughness, and percentage lichen cover. The ethylene glycol test was performed on prepared rock samples to determine susceptibility to tensional breakup. A satellite image depicting the distribution of late-lying contemporary snow was used to assist in the construction of topo-climatic linkages with scarp terrace localities. Findings show that terraces are most common on southeast-facing aspects, and coincide with the major joint strike direction. Lithological factors are thus considered the primary control to such terrace development. It is suggested that a different set of geomorphological processes operates on various slope-altitudinal and slope-orientational positions at any given time.     

A multigrid solver for 3D electromagnetic diffusion

The performance of a multigrid solver for the time‐harmonic electromagnetic problem in geophysical settings is investigated. The frequencies are sufficiently small for waves travelling at the speed of light to be negligible, so that a diffusive problem remains. The discretization of the governing equations is obtained by the finite‐integration technique, which can be viewed as a finite‐volume generalization of Yee's staggered grid scheme. The resulting set of discrete equations is solved by a multigrid method. The convergence rate of the multigrid method decreased when the grid was stretched. The slower convergence rate of the multigrid method can be compensated by using bicgstab2 , a conjugate‐gradient‐type method for non‐symmetric problems. In that case, the multigrid solver acts as a preconditioner. However, whereas the multigrid method provides excellent convergence with constant grid spacings, it performs less than satisfactorily when substantial grid stretching is used.     

Exploring some issues in acoustic full waveform inversion

The least‐squares error measures the difference between observed and modelled seismic data. Because it suffers from local minima, a good initial velocity model is required to avoid convergence to the wrong model when using a gradient‐based minimization method. If a data set mainly contains reflection events, it is difficult to update the velocity model with the least‐squares error because the minimization method easily ends up in the nearest local minimum without ever reaching the global minimum. Several authors observed that the model could be updated by diving waves, requiring a wide‐angle or large‐offset data set. This full waveform tomography is limited to a maximum depth. Here, we use a linear velocity model to obtain estimates for the maximum depth. In addition, we investigate how frequencies should be selected if the seismic data are modelled in the frequency domain. In the presence of noise, the condition to avoid local minima requires more frequencies than needed for sufficient spectral coverage. We also considered acoustic inversion of a synthetic marine data set created by an elastic time‐domain finite‐difference code. This allowed us to validate the estimates made for the linear velocity model. The acoustic approximation leads to a number of problems when using long‐offset data. Nevertheless, we obtained reasonable results. The use of a variable density in the acoustic inversion helped to match the data at the expense of accuracy in the inversion result for the density.     

Initial formation of channel–shoal patterns in double–inlet systems

Ocean Dynamics - Channel–shoal patterns are often observed in the back–barrier basins of inlet systems and are important from both an economical and ecological point of view. Focussing...     

Research Note: A workaround for the corner problem in numerically exact non-reflecting boundary conditions

Simulations of wave propagation in the Earth usually require truncation of a larger domain to the region of interest to keep computational cost acceptable. This introduces artificial boundaries that should not generate reflected waves. Most existing boundary conditions are not able to completely suppress all the reflected energy, but suffice in practice except when modelling subtle events such as interbed multiples. Exact boundary conditions promise better performance but are usually formulated in terms of the governing wave equation and, after discretization, still may produce unwanted artefacts. Numerically exact non-reflecting boundary conditions are instead formulated in terms of the discretized wave equation. They have the property that the numerical solution computed on a given domain is the same as one on a domain enlarged to the extent that waves reflected from the boundary do not have the time to reach the original truncated domain. With a second- or higher-order finite-difference scheme for the one-dimensional wave equation, these boundary conditions follow from a recurrence relation. In its generalization to two or three dimensions, a recurrence relation was only found for a single non-reflecting boundary on one side of the domain or two of them at opposing ends. The other boundaries should then be zero Dirichlet or Neumann. If two non-reflecting boundaries meet at a corner, translation invariance is lost and a simple recurrence relation could not be found. Here, a workaround is presented that restores translation invariance by imposing classic, approximately non-reflecting boundary conditions on the other sides and numerically exact ones on the two opposing sides that otherwise would create the strongest reflected waves with the classic condition. The exact ones can also be applied independently. As a proof of principle, the method is applied to the two-dimensional acoustic wave equation, discretized on a rectangular domain with a second-order finite-difference scheme and first-order Enquist–Majda boundary conditions as approximate ones. The method is computationally costly but has the advantage that it can be reused on a sequence of problems as long as the time step and the sound speed values next to the boundary are kept fixed.     

Assessment of the geometry and volumes of small surface water reservoirs by remote sensing in a semi-arid region with high reservoir density

Water fluxes in highly impounded regions are heavily dependent on reservoir properties. However, for large and remote areas, this information is often unavailable. In this study, the geometry and volume of small surface reservoirs in the semi-arid region of Brazil were estimated using terrain and shape attributes extracted by remote sensing. Regression models and data classification were used to predict the volumes, at different water stages, of 312 reservoirs for which topographic information is available. The power function used to describe the reservoir shapes tends to overestimate the volumes; therefore, a modified shape equation was proposed. Among the methods tested, four were recommended based on performance and simplicity, for which the mean absolute percentage errors varied from 24 to 39%, in contrast to the 94% error achieved with the traditional method. Despite the challenge of precisely deriving the flooded areas of reservoirs, water management in highly reservoir-dense environments should benefit from volume prediction based on remote sensing.     

Transatlantic latest Cretaceous mosasaurs (Reptilia,Lacertilia) from the Maastrichtian type area and New Jersey

Cranial and postcranial elements of the mosasaurid reptiles Mosasaurus hoffmanni Mantell 1829 and Plioplatecarpus marshi Dollo 1882 from the Maastrichtian type area, and of Mosasaurus maximus Cope 1869 and Plioplatecarpus depressus (Cope 1869) from New Jersey are compared. Views held by previous authors are discussed. It is concluded that these European and North American taxa are conspecific.     

The influence of geological, geochemical, and biogenic habitat heterogeneity on seep biodiversity

Cold seeps are among the most heterogeneous of all continental margin habitats. Abiotic sources of heterogeneity in these systems include local variability in fluid flow, geochemistry, and substrate type, which give rise to different sets of microbial communities, microbial symbiont-bearing foundation species, and associated heterotrophic species. Biogenic habitats created by microbial mats and the symbiotic species including vesicomyid clams, bathymodiolin mussels, and siboglinid tubeworms add an additional layer of complexity to seep habitats. These forms of habitat heterogeneity result in a variety of macrofaunal and meiofaunal communities that respond to changes in structural complexity, habitat geochemistry, nutrient sources, and interspecific interactions in different ways and at different scales. These responses are predicted by a set of theoretical metacommunity models, the most appropriate of which for seep systems appears to be the 'species sorting' concept, an extension of niche theory. This concept is demonstrated through predictable patterns of community assembly, succession, and beta-level diversity. These processes are described using a newly developed analytical technique examining the change in the slope of the species accumulation curve with the number of habitats examined. The diversity response to heterogeneity has a consistent form, but quantitatively changes at different seep sites around the world as the types of habitats present and the size-classes of fauna analyzed change. The increase in beta diversity across seep habitat types demonstrates that cold seeps and associated biogenic habitats are significant sources of heterogeneity on continental margins globally.