An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

A robust algorithm for ellipse-based discrete element modelling of granular materials

With recent research indicating the importance of the rolling mechanism of deformation in granular systems consisting of perfectly round particles, it has become popular to use ellipse-shaped particles in the Discrete Element Method (DEM) numerical model. Inherent in this technique is the need for accurately computing ellipse to ellipse intersection, in order to properly detect contact formation and compute relative contact velocities. However, the commonly used algorithms for computing ellipse-ellipse intersection are generally poorly conditioned and can be inaccurate. An alternate method for computing ellipse-ellipse intersection is developed and presented which results in a well-conditioned, stable and accurate contact detection method. These modification are incorporated into the general DEM algorithm.     

From Celsius to Chaos to Cyclones: Using Temperature-Conversion Equations to Introduce Advanced Mathematical Concepts in Earth Science Courses

The Fahrenheit-to-Celsius temperature-conversion equation is a basic component of many introductory earth science courses. Despite its simplicity, it presents a challenge to students and instructors alike because residents of the United States are unfamiliar with the Celsius scale. By solving for the point at which these two temperature scales are equal, it is possible to use the equations for temperature conversion as a springboard to more advanced topics. It is demonstrated that temperature-conversion equations and chaotic equations can be solved using identical numerical and graphical techniques. As a result, the fundamental concepts of chaos theory and numerical methods can be introduced to students in the context of the simplest equations in the earth sciences. These solution methods are applied to the quantitative theory of the extratropical cyclone as an example of the utility and broad scope of this educational approach.     

A crustal suture and lineament in North Africa

A possible suture between an old craton in West Africa and the Arabian—Nubian shield is marked by a persistent zone of Phanerozoic weakness. This zone extends between the Nile and Niger deltas and passes through the Qattara embayment, along the Bahariya fault, through the Kufra basin, between the Tibesti and Jabal al Uwaynat massifs, through the Bahr el Ghazal and Lake Chad areas, and along the Benoue trough.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

A comparison of the photospheric and interplanetary magnetic fields during the flight of Mariner IV

The polarity of the interplanetary magnetic field observed by the Mariner IV spacecraft is compared to the polarity of the photospheric magnetic field observed with the solar magnetograph at Mt. Wilson Observatory. Unlike the results obtained from observations during the flight of IMP-I, these polarities are not well correlated when the photospheric polarity is determined from data along a narrow latitudinal strip. It is suggested that the structure of the interplanetary field is often related to major features in the photospheric field that are observed over a broad range of solar latitudes.     

Testing optically stimulated luminescence dating on sand-sized quartz of deltaic deposits from the Sperchios delta plain,central Greece

This study reports on the first investigation into the potential of luminescence dating to establish a chronological framework for the depositional sequences of the Sperchios delta plain, central Greece. A series of three borehole cores(20 m deep) and two shallow cores(4 m deep), from across the delta plain, were extracted, and samples were collected for luminescence dating. The luminescence ages of sand-sized quartz grains were obtained from small aliquots of quartz, using the Single-Aliquot Regenerative-dose(SAR) protocol.The equivalent dose determination included a series of tests and the selection of the Minimum Age Model(MAM) as the most appropriate statistical model. This made it possible to confirm the applicability of quartz Optically Stimulated Luminescence(OSL) dating to establish absolute chronology for deltaic sediments from the Sperchios delta plain.Testing age results of the five cores showed that the deltaic sediments were deposited during the Holocene.A relatively rapid deposition is implied for the top ~14 m possibly as a result of the deceleration in the rate of the sea-level rise and the transition to terrestrial conditions, while on the deeper parts, the reduced sedimentation rate may indicate a lagoonal or coastal environment.     

MoisturEC: A New R Program for Moisture Content Estimation from Electrical Conductivity Data

Noninvasive geophysical estimation of soil moisture has potential to improve understanding of flow in the unsaturated zone for problems involving agricultural management, aquifer recharge, and optimization of landfill design and operations. In principle, several geophysical techniques (e.g., electrical resistivity, electromagnetic induction, and nuclear magnetic resonance) offer insight into soil moisture, but data‐analysis tools are needed to “translate” geophysical results into estimates of soil moisture, consistent with (1) the uncertainty of this translation and (2) direct measurements of moisture. Although geostatistical frameworks exist for this purpose, straightforward and user‐friendly tools are required to fully capitalize on the potential of geophysical information for soil‐moisture estimation. Here, we present MoisturEC, a simple R program with a graphical user interface to convert measurements or images of electrical conductivity (EC) to soil moisture. Input includes EC values, point moisture estimates, and definition of either Archie parameters (based on experimental or literature values) or empirical data of moisture vs. EC. The program produces two‐ and three‐dimensional images of moisture based on available EC and direct measurements of moisture, interpolating between measurement locations using a Tikhonov regularization approach.