Numerical modeling of space-time wave extremes using WAVEWATCH III

A novel implementation of parameters estimating the space-time wave extremes within the spectral wave model WAVEWATCH III (WW3) is presented. The new output parameters, available in WW3 version 5.16, rely on the theoretical model of Fedele (J Phys Oceanogr 42(9):1601-1615, 2012) extended by Benetazzo et al. (J Phys Oceanogr 45(9):2261–2275, 2015) to estimate the maximum second-order nonlinear crest height over a given space-time region. In order to assess the wave height associated to the maximum crest height and the maximum wave height (generally different in a broad-band stormy sea state), the linear quasi-determinism theory of Boccotti (2000) is considered. The new WW3 implementation is tested by simulating sea states and space-time extremes over the Mediterranean Sea (forced by the wind fields produced by the COSMO-ME atmospheric model). Model simulations are compared to space-time wave maxima observed on March 10th, 2014, in the northern Adriatic Sea (Italy), by a stereo camera system installed on-board the “Acqua Alta” oceanographic tower. Results show that modeled space-time extremes are in general agreement with observations. Differences are mostly ascribed to the accuracy of the wind forcing and, to a lesser extent, to the approximations introduced in the space-time extremes parameterizations. Model estimates are expected to be even more accurate over areas larger than the mean wavelength (for instance, the model grid size).     

Seismic response of bridge piers on pile groups for different soil damping models and lumped parameter representations of the foundation

This paper presents a wide parametric study aimed at elucidating the influence, on the computed seismic response of bridge piers, of two related aspects of the model: (1) the adoption of the classical hysteretic or the causal Biot's damping models for the soil and (2) the use of two different lumped parameter models of different complexity and accuracy to approximate the impedances of the pile foundation. A total of 2072 cases, including different superstructures, pile foundations, soil deposits, and seismic input signals, are studied. The results are presented so that the influence of the different parameters involved in the analysis can be assessed. From an engineering point of view, both lumped parameter models provide, in general, sufficiently low errors. The choice of the most adequate model for each case will depend not only on the configuration of the structure and the soil-foundation system but also on the assumed soil damping model, whose influence on the computed seismic responses is relevant in many cases. The nonphysical behaviour provided by the classical hysteretic damping model for the soil at zero frequency generates issues in the process of fitting the impedance functions. It is also found that larger deck displacements are predicted by Biot's model due to the higher damping at low frequencies provided by the classical hysteretic damping model.     

Investigating the impact of surface wave breaking on modeling the trajectories of drifters in the northern Adriatic Sea during a wind-storm event

An accurate numerical prediction of the oceanic upper layer velocity is a demanding requirement for many applications at sea and is a function of several near-surface processes that need to be incorporated in a numerical model. Among them, we assess the effects of vertical resolution, different vertical mixing parameterization (the so-called Generic Length Scale –GLS– set of k–ε, k–ω, gen, and the Mellor–Yamada), and surface roughness values on turbulent kinetic energy (k) injection from breaking waves.First, we modified the GLS turbulence closure formulation in the Regional Ocean Modeling System (ROMS) to incorporate the surface flux of turbulent kinetic energy due to wave breaking. Then, we applied the model to idealized test cases, exploring the sensitivity to the above mentioned factors. Last, the model was applied to a realistic situation in the Adriatic Sea driven by numerical meteorological forcings and river discharges. In this case, numerical drifters were released during an intense episode of Bora winds that occurred in mid-February 2003, and their trajectories compared to the displacement of satellite-tracked drifters deployed during the ADRIA02-03 sea-truth campaign.Results indicted that the inclusion of the wave breaking process helps improve the accuracy of the numerical simulations, subject to an increase in the typical value of the surface roughness z₀. Specifically, the best performance was obtained using α_CH = 56,000 in the Charnok formula, the wave breaking parameterization activated, k–ε as the turbulence closure model. With these options, the relative error with respect to the average distance of the drifter was about 25% (5.5 km/day). The most sensitive factors in the model were found to be the value of α_CH enhanced with respect to a standard value, followed by the adoption of wave breaking parameterization and the particular turbulence closure model selected.     

The Radio and X-ray Mode-Switching Pulsar PSR B0943+10

Observations obtained in the last years challenged the widespread notion that rotation-powered neutron stars are steady X-ray emitters. Besides a few allegedly rotation-powered neutron stars that showed ‘magnetar-like’ variability, a particularly interesting case is that of PSR B0943+10. Recent observations have shown that this pulsar, well studied in the radio band where it alternates between a bright and a quiescent mode, displays significant X-ray variations, anticorrelated in flux with the radio emission. The study of such synchronous radio/X-ray mode switching opens a new window to investigate the processes responsible for the pulsar radio and high-energy emission. Here we review the main X-ray properties of PSR B0943+10 derived from recent coordinated X-ray and radio observations.     

Three-Dimensional Mapping of Air Flow at an Urban Canyon Intersection

In this experimental work both qualitative (flow visualisation) and quantitative (laser Doppler anemometry) methods were applied in a wind tunnel in order to describe the complex three-dimensional flow field in a real environment (a street canyon intersection). The main aim was an examination of the mean flow, turbulence and flow pathlines characterising a complex three-dimensional urban location. The experiments highlighted the complexity of the observed flows, particularly in the upwind region of the intersection. In this complex and realistic situation some details of the upwind flow, such as the presence of two tall towers, play an important role in defining the flow field within the intersection, particularly at roof level. This effect is likely to have a strong influence on the mass exchange mechanism between the canopy flow and the air aloft, and therefore the distribution of pollutants. This strong interaction between the flows inside and outside the urban canopy is currently neglected in most state-of-the-art local scale dispersion models.     

Seismic soil–structure interaction in multi‐span bridges: Application to a railway bridge

The effects of soil‐structure interaction on the seismic response of multi‐span bridges are investigated by means of a modelling strategy based on the domain decomposition technique. First, the analysis methodology is presented: kinematic interaction analysis is performed in the frequency domain by means of a procedure accounting for radiation damping, soil–pile and pile‐to‐pile interaction; the seismic response of the superstructure is evaluated in the time domain by means of user‐friendly finite element programs introducing suitable lumped parameter models take into account the frequency‐dependent impedances of the soil–foundation system. Second, a real multi‐span railway bridge longitudinally restrained at one abutment is analyzed. The input motion is represented by two sets of real accelerograms: one consistent with the Italian seismic code and the other constituted by five records characterized by different frequency contents. The seismic response of the compliant‐base model is compared with that obtained from a fixed‐base model. Pile stress resultants due to kinematic and inertial interactions are also evaluated. The application demonstrates the importance of performing a comprehensive analysis of the soil–foundation–structure system in the design process, in order to capture the effects of soil‐structure interaction in each structural element that may be beneficial or detrimental. Copyright © 2011 John Wiley & Sons, Ltd.     

Prediction of equilibrium Li isotope fractionation between minerals and aqueous solutions at high P and T: An efficient ab initio approach

The mass-dependent equilibrium stable isotope fractionation between different materials is an important geochemical process. Here we present an efficient method to compute the isotope fractionation between complex minerals and fluids at high pressure, P, and temperature, T, representative for the Earth’s crust and mantle. The method is tested by computation of the equilibrium fractionation of lithium isotopes between aqueous fluids and various Li bearing minerals such as staurolite, spodumene and mica. We are able to correctly predict the direction of the isotope fractionation as observed in the experiments. On the quantitative level the computed fractionation factors agree within 1.0‰ with the experimental values indicating predictive power of ab initio methods. We show that with ab initio methods we are able to investigate the underlying mechanisms driving the equilibrium isotope fractionation process, such as coordination of the fractionating elements, their bond strengths to the neighboring atoms, compression of fluids and thermal expansion of solids. This gives valuable insight into the processes governing the isotope fractionation mechanisms on the atomic scale. The method is applicable to any state and does not require different treatment of crystals and fluids.