Numerical study of tides in Ontario Lacus,a hydrocarbon lake on the surface of the Saturnian moon Titan

In the context of the emergence of extra-terrestrial oceanography, we adapted an existing oceanographic model, SLIM (www.climate.be/slim), to the conditions of Titan, a moon of Saturn. The tidal response of the largest southern lake at Titan’s surface, namely Ontario Lacus, is simulated. SLIM solves the 2D, depth-averaged shallow water equations on an unstructured mesh using the discontinuous Galerkin finite element method, which allows for high spatial resolution wherever needed. The impact of the wind forcing, the bathymetry, and the bottom friction is also discussed. The predicted maximum tidal range is about 0.56 m in the southern part of the lake, which is more than twice as large as the previous estimates (see Tokano, Ocean Dyn 60:(4) 803–817 10.1007/s10236-010-0285-3 (Tokano 2010)). The patterns and magnitude of the current are also markedly different from those of previous studies: the tidal motion is not aligned with the major axis of the lake and the speed is larger nearshore. Indeed, the main tidal component rotates clockwise in an exact period of one Titan day and the tidal currents can reach 0.046 ms ^?1 close to the shores depending on the geometry and the bathymetry. Except for these specific nearshore regions, the current speed is less than 0.02 ms ^?1. Circular patterns can be observed offshore, their rotational direction and size varying along the day.     

Building damage scenarios based on exploitation of Housner intensity derived from finite faults ground motion simulations

In this paper earthquake damage scenarios for residential buildings (about 4200 units) in Potenza (Southern Italy) have been estimated adopting a novel probabilistic approach that involves complex source models, site effects, building vulnerability assessment and damage estimation through Damage Probability Matrices. Several causative faults of single seismic events, with magnitude up to 7, are known to be close to the town. A seismic hazard approach based on finite faults ground motion simulation techniques has been used to identify the sources producing the maximum expected ground motion at Potenza and to generate a set of ground motion time histories to be adopted for building damage scenarios. Additionally, site effects, evaluated in a previous work through amplification factors of Housner intensity, have been combined with the bedrock values provided by hazard assessment. Furthermore, a new relationship between Housner and EMS-98 macroseismic intensity has been developed. This relationship has been used to convert the probability mass functions of Housner intensity obtained from synthetic seismograms amplified by the site effects coefficients into probability mass function of EMS-98 intensity. Finally, the Damage Probability Matrices have been applied to estimate the damage levels of the residential buildings located in the urban area of Potenza. The proposed methodology returns the full probabilistic distribution of expected damage, thus avoiding average damage index or uncertainties expressed in term of dispersion indexes.     

Concept design and hydrodynamic optimization of an innovative SWATH USV by CFD methods

The paper presents the main characteristics of an innovative platform which has been conceived and designed to extend the operational capabilities of current unmanned surface vehicles in terms of platform stability in waves and of powering requirement at a relatively high speed. The main idea which rules the project is the realization of a small autonomous surface unit (about 6 m in length) capable of undertaking several tasks in the marine environment even with moderate rough sea conditions. The designed vessel has the ability to locate, recover, and launch other members of the autonomous fleet (like AUVs or other underwater devices) and at the same time could carry out a surveillance service of the surrounding areas. To manage these tasks, the vehicle is designed to provide a fairly good autonomy which is needed to face intermediate-range missions (100 nautical miles). The choice of a small waterplane area twin hull (SWATH) form has been motivated by its excellent properties of seakeeping qualities, combined with a non-conventional low resistance underwater hull shape, currently under patenting process, which is able to reduce to a minimum the resistance of the vessel especially at higher speeds. To obtain the most efficient profile of the underwater bodies, a systematic optimization with an automatic procedure based on a parametric definition of the geometry, a state-of-the-art computational fluid dynamics (CFD) flow solver, and a differential evolution global minimization algorithm have been created and used. As expected, all the final CFD computations on the best design have demonstrated the superior efficiency of the developed unconventional SWATH technology with respect to different alternatives of current hull typologies.     

Flow and transport in channels with submerged vegetation

This paper reviews recent work on flow and transport in channels with submerged vegetation, including discussions of turbulence structure, mean velocity profiles, and dispersion. For submerged canopies of sufficient density, the dominant characteristic of the flow is the generation of a shear-layer at the top of the canopy. The shear-layer generates coherent vortices by Kelvin-Helmholtz (KH) instability. These vortices control the vertical exchange of mass and momentum, influencing both the mean velocity profile, as well as the turbulent diffusivity. For flexible canopies, the passage of the KH vortices generates a progressive wave along the canopy interface, termed monami. The KH vortices formed at the top of the canopy penetrate a distance δ _e into the canopy. This penetration scale segregates the canopy into an upper layer of rapid transport and a lower layer of slow transport. Flushing of the upper canopy is enhanced by the energetic shear-scale vortices. In the lower layer turbulence is limited to length-scales set by the stem geometry, and the resulting transport is significantly slower than that of the upper layer.     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

Multi-method global sensitivity analysis of flood inundation models

Global sensitivity analysis is a valuable tool in understanding flood inundation models and deriving decisions on strategies to reduce model uncertainty. In this paper, a sensitivity analysis of a one-dimensional flood inundation model (HEC-RAS) on the River Alzette, Luxembourg, is presented. It is impossible to define sensitivity in a unique way and different methods can lead to a difference in ranking of importance of model factors. In this paper five different methods (Sobol, Kullback–Leibler entropy, Morris, regionalised sensitivity analysis and regression) are applied and the outcomes on selected examples compared. It is demonstrated that the different methods lead to completely different ranking of importance of the parameter factors and that it is impossible to draw firm conclusions about the relative sensitivity of different factors. Moreover, the uncertainty inherent in the sensitivity methods is highlighted.     

Dissimilarity among Ocean Reanalyses in Equatorial Pacific Upper-Ocean Heat Content and Its Relationship with ENSO

This study focuses on the temporal variation of dissimilarity in heat content(HC)anomalies in the upper 300 m of ocean(HC300A)in the equatorial Pacific(±10°N)and its response to the El Ni?o-Southern Oscillation(ENSO).The HC300A anomalies are derived from four ocean reanalyses that are commonly used in ENSO studies and are compared using a simple differencing method.The dissimilarity in HC300A is found to vary closely with the magnitude of ENSO(regardless of phase),meaning that it tends to be greater during strong ENSO events.However,the dissimilarity among ocean reanalyses persists after the event decays.This effect is more pronounced after strong events.The persistence of the dissimilarity after ENSO events is a result of a late maturation of the ENSO signal,its persistence,and the interruption of the signal decay due to follow-up ENSO events.The combined effect of these three factors slows down the decay of HC300A in the region and hence results in the slow decay of dissimilarity.It is also found that areas with a significant spread in vertical temperature profiles collocate with the ENSO signal during warm ENSO phases.Thus,differences in subsurface process reconstruction are a significant factor in the dissimilarity among ocean reanalyses during warm ENSO events.     

Curvature analysis as a tool for subsidence-related risk zones identification in the city of Tuzla (BiH)

Shallow Miocene salt deposits located below the city of Tuzla (BiH) have been exploited during the last 60 years by means of wells extracting salt-saturated groundwater brines. The massive extraction activities have led to severe subsidence accommodated by collapse and strain localization. Surface topography and geomorphology have been influenced and modified by several faults and fractures.A series of sequential topographical survey data collected during the last 50 years have been used to obtain the total subsidence envelope surface. In order to analyze the salt dissolution-related morphological development of the Tuzla topography and to identify the location of faults and fractures in the deformed area, we applied curvature analysis to the subsidence surface by means of two different analytical methods: the directional, two-dimensional curvature and the analytic Gaussian curvature. The comparison of the curvature maps with those of surface fractures visible in the area shows a good spatial agreement between the directional curvature and fracture intensity, whereas the subsurface normal faults are more evident in the Gaussian curvature maps.