Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

Tsunami generation,propagation, and run-up with a high-order Boussinesq model

In this work we extend a high-order Boussinesq-type (finite difference) model, capable of simulating waves out to wavenumber times depth kh < 25, to include a moving sea-bed, for the simulation of earthquake- and landslide-induced tsunamis. The extension is straight forward, requiring only an additional term within the kinematic bottom condition. As first test cases we simulate linear and nonlinear surface waves generated from both positive and negative impulsive bottom movements. The computed results compare well against earlier theoretical, numerical, and experimental values. Additionally, we show that the long-time (fully nonlinear) evolution of waves resulting from an upthrusted bottom can eventually result in true solitary waves, consistent with theoretical predictions. It is stressed, however, that the nonlinearity used far exceeds that typical of geophysical tsunamis in the open ocean. The Boussinesq-type model is then used to simulate numerous tsunami-type events generated from submerged landslides, in both one and two horizontal dimensions. The results again compare well against previous experiments and/or numerical simulations. The new extension compliments recently developed run-up capabilities within this approach, and as demonstrated, the model can therefore treat tsunami events from their initial generation, through their later propagation, and final run-up phases. The developed model is shown to maintain reasonable computational efficiency, and is therefore attractive for the simulation of such events, especially in cases where dispersion is important.     

Swimming behavior of Dungeness crab,Cancer magister Dana,megalopae in still and moving water

In the Grays Harbor estuary, juvenile Dungeness crab (Cancer magister Dana) are found at higher densities in epibenthic shell deposits compared to open mud flat. Differences in predation rate between habitats have been suggested to be due to habitat preference and differential survival. Megalopae preferred shell over open space in still-water conditions. However, it is not known whether megalopae are able to select shell in flowing water since larval preference is known to differ between still and flowing water. Here we report the first experimental study of swimming behavior of Dungeness crab megalopae in a range of current velocities (0–40 cm s^?1) equivalent to natural flow in Grays Harbor estuary. Experiments were conducted in daylight using a recirculating flume. Megalopae swimming speeds ranged from 8.5 cm s^?1 (8 body lengths s^?1) in still water to 44.8 cm s^?1 (44 body lengths s^?1) at flow speeds of 40 cm s^?1, Neither swimming behavior nor sheltering behavior in shells showed any flow-related pattern. Megalopae spent a large proportion of time swimming against the current and made headway upstream against all current velocities tested. The results suggest that Dungeness crab megalopae are able to maneuver and actively search for settlement sites under current velocities found in natural habitats, including intertidal shell deposits, and support the hypothesis of active selection of shell by megalopae.     

Studies of compaction and subsidence in the Orinoco Oil Belt

To control vertical movements in the Orinoco Oil Belt, a high-precision levelling network was installed in an area of about 100 km². Links to the National Geodetic Network were established.

A total of 190 BMs with deep pile foundations is already installed and the first measurement campaign was concluded. The accuracies obtained are those of first-order levelling works.     

A ∼6100 14C yr record of El Niño activity from the Galápagos Islands

Lithostratigrahic and mineralogic analyses of sediments from hypersaline Bainbridge Crater Lake, Galápagos Islands, provide evidence of past El Niño frequency and intensity. Laminated sediments indicate that at least 435 moderate to very strong El Niño events have occurred since 6100 ¹⁴C yr BP (7130 cal yr BP), and that frequency and intensity of events increased at about 3000 ¹⁴C yr BP (3100 cal yr BP). El Niño activity was present between 6100 and 4000 ¹⁴C yr BP (4600 cal yr BP) but infrequent. The Bainbridge record indicates that there has been considerable millennial-scale variability in El Niño since the mid-Holocene.     

Diatoms on earth,comets, europa and in interstellar space

There exists a close correspondence between the measured infrared properties of diatoms and the infrared spectrum of interstellar dust as observed in the Trapezium nebula and toward the galactic center source GC-IRS 7. Diatoms and bacteria also exhibit an absorbance peak near 2200 », which is found to agree with the observed ultraviolet absorbance properties of interstellar grains. We review the observational data and consider the known properties of diatoms and bacteria. It is suggested that these characteristics are consistent with the concept of a cosmic microbiological system in which these or similar microorganisms might exist on comets, Europa and in interstellar space.     

WSO and the winds of massive stars: the gate to the metal-poor Local Universe

The spectrographs on-board the World Space Observatory (WSO) will provide access to the 1020–1800 Å wavelength range with unprecedented sensitivity. Previous observatories operating in the 1150–2000 Å range (such as IUE and HST-STIS) have proved extremely useful to study the winds of OB type stars, which leave their most prominent imprints in the far ultraviolet range. The addition of the λ < 1200 Å wavelengths is critical as it contains important diagnostic lines for mass loss and shocks in the wind, as found by FUSE-based analyses.WSO will enable quantitative spectroscopic analyses of blue massive stars in the Local Group beyond the Magellanic Clouds. The results will lead to the characterization of their winds as a function of metallicity, and shed new light on current urging questions regarding radiation driven winds.     

Mercury surface composition: Integrating petrologic modeling and remote sensing data to place constraints on FeO abundance

A significant opaque component in Mercury’s crust is inferred based on albedo and spectral observations. Previous workers have favored iron-titanium bearing oxide minerals as the spectrally neutral opaque. A consequence of this hypothesis is that Mercury’s surface would have a high FeO content. An array of remote sensing techniques have not provided definitive constraints on the FeO content of Mercury’s surface. However, spectral observations have not detected a diagnostic 1 μm absorption band and have thus limited the FeO in coexisting silicates to <2 wt.% FeO. In this paper, we assess equilibrium among oxide and silicate minerals to constrain the distribution of iron between opaque oxides and silicates under a variety of environmental conditions. Equilibrium modeling is favored here because the geologic process that produced Mercury’s low-albedo intermediate terrain must have occurred globally, which favors a common widespread igneous process. Based on our modeling, we find that iron-rich ilmenite cannot occur with silicates that do not display a 1 μm absorption feature unless plagioclase abundances are high. However, such high plagioclase abundances are precluded by Mercury’s low albedo. Incorporating equilibrium crystallization modeling with spectral and albedo constraints we find the iron abundance of Mercury’s intermediate terrain is ?10 wt.% FeO. This intermediate iron composition matches constraints provided by visible albedo and total neutron absorption observed by MESSENGER. In fact, the total neutron absorption of mixtures of oxide, plagioclase, olivine and pyroxene for the oxide abundances estimated for Mercury, favor Mg-rich members of the ilmenite-geikielite solid-solution series. This work offers compositional constraints for Fe, Ti, and Mg that will be testable by various MESSENGER instrument data sets after it begins its orbital mission.