Hybrid frequency–time domain models for dynamic response analysis of marine structures

Time-domain models of marine structures based on frequency domain data are usually built upon the Cummins equation. This type of model is a vector integro-differential equation which involves convolution terms. These convolution terms are not convenient for analysis and design of motion control systems. In addition, these models are not efficient with respect to simulation time, and ease of implementation in standard simulation packages. For these reasons, different methods have been proposed in the literature as approximate alternative representations of the convolutions. Because the convolution is a linear operation, different approaches can be followed to obtain an approximately equivalent linear system in the form of either transfer function or state-space models. This process involves the use of system identification, and several options are available depending on how the identification problem is posed. This raises the question whether one method is better than the others. This paper therefore has three objectives. The first objective is to revisit some of the methods for replacing the convolutions, which have been reported in different areas of analysis of marine systems: hydrodynamics, wave energy conversion, and motion control systems. The second objective is to compare the different methods in terms of complexity and performance. For this purpose, a model for the response in the vertical plane of a modern containership is considered. The third objective is to describe the implementation of the resulting model in the standard simulation environment Matlab/Simulink.     

A Coastal Bay Summer Breeze Study,Part 2: High-resolution Numerical Simulation of Sea-breeze Local Influences

We complete the analysis of the data obtained during the experimental campaign around the semi circular bay of Quiberon, France, during two weeks in June 2006 (see Part 1). A reanalysis of numerical simulations performed with the Advanced Regional Prediction System model is presented. Three nested computational domains with increasing horizontal resolution down to 100 m, and a vertical resolution of 10 m at the lowest level, are used to reproduce the local-scale variations of the breeze close to the water surface of the bay. The Weather Research and Forecasting mesoscale model is used to assimilate the meteorological data. Comparisons of the simulations with the experimental data obtained at three sites reveal a good agreement of the flow over the bay and around the Quiberon peninsula during the daytime periods of sea-breeze development and weakening. In conditions of offshore synoptic flow, the simulations demonstrate that the semi-circular shape of the bay induces a corresponding circular shape in the offshore zones of stagnant flow preceding the sea-breeze onset, which move further offshore thereafter. The higher-resolution simulations are successful in reproducing the small-scale impacts of the peninsula and local coasts (breeze deviations, wakes, flow divergences), and in demonstrating the complexity of the breeze fields close to the surface over the bay. Our reanalysis also provides guidance for numerical simulation strategies for analyzing the structure and evolution of the near-surface breeze over a semi-circular bay, and for forecasting important flow details for use in upcoming sailing competitions.     

Mineral equilibria constraints on open‐system melting in metamafic compositions

The recent development of activity–composition relations for mineral and melt phases in high‐grade metamafic rocks allows mineral equilibria tools to be used to further aid our understanding of partial melting and the mineralogical consequences of melt segregation in these rocks. We show that bulk compositional data from natural amphibolites cover a wide compositional range, with particular variability in the content and ratios of Ca, Na and K indicating that low‐grade metasomatic alteration can substantially alter the igneous protolith chemistry and potentially affect the volume and composition of melt generated. Mineral equilibria calculations for five samples that span the compositional variability in our data set indicate that melting occurs primarily via the fluid‐absent breakdown of amphibole+quartz to produce a pressure‐sensitive peritectic assemblage of augite, orthopyroxene and/or garnet. The introduction of orthopyroxene at the onset of the amphibolite‐to‐granulite‐facies transition at lower pressure results in an increased rate of melt production until quartz is typically exhausted, and this is similarly seen for the introduction of garnet at higher pressure. Calculated melt compositions are dependent on the protolith composition, but initial solidus melting and biotite breakdown produce 1–3 mol.% of K‐rich granitic melts. As hornblende melting proceeds, 15–20 vol.% of either more granodioritic‐to‐tonalitic or granodioritic‐to‐trondhjemitic melt is produced. Once quartz is exhausted, intermediate to mafic melt compositions are produced at ultrahigh‐temperature conditions. Quartz‐rich lithologies with high Ca coupled to low Na and K are the most fertile under orogenic conditions, yielding up to 25 mol.% of sub‐alkalic granitic melt by 850°C. Such rocks did not experience significant subsolidus alteration. Altered compositions with low Ca and elevated Na and K are not as fertile, yielding less than 15 mol.% of alkalic granitic melt by 850°C. These melt volumes are enough to be segregated, and can make a contribution to granite magmatism and intracrustal differentiation that should not be overlooked.     

An urban neighborhood temperature and energy study from the CAPITOUL experiment with the Solene model

A methodology is proposed to analyze the radiative and thermal exchanges between a small urban neighborhood and the atmosphere based on the use of the thermoradiative model SOLENE and radiometric measurements to optimize the effective values of its constant parameters. Applied to the center of Toulouse, France, the optimization data are building surface temperatures measured with handheld radiometers and a downward-facing pyranometer during one of the intense observation periods of the CAnopy and Particles Interactions in TOulouse Urban Layer (CAPITOUL) experimental campaign. The quality of the simulations is assessed by comparing, without any other model adjustment, the model outputs for two diurnal cycles (1?day in summer and 1?day in winter) against two independent experimental datasets from fixed permanent radiometers and from sensors measuring the solar, infrared and sensible heat fluxes to the atmosphere at the top of a mast 30?m above the roofs. These simulations allow us to further analyze the separate contributions of the different surface classes, roofs, facades and pavement to these fluxes and to compare them with their counterparts observed over a neighborhood of Marseilles city center during a previous experiment. The partition is remarkably similar in Toulouse and Marseilles: 6–7% for the solar radiation, 73–76% for the infrared radiation and 17–20% for the sensible heat flux. The contribution of roofs to the infrared flux appears proportional to their plan area proportion (built density) but not the contributions to the other two fluxes. The contributions of facades to all three fluxes are roughly proportional to their fraction of the total surface area.     

High-resolution analysis of submarine lobes deposits: Seismic-scale outcrops of the Lauzanier area (SE Alps, France)

The Lauzanier area represents the northernmost extension of the Annot Sandstone series and contains deposits between 650 and 900 m-thick. This basin was active from upper Bartonian or lower Priabonian to early Rupelian. It is composed of two superposed units separated by a major unconformity. The sediment supply is due to channelled flows coming from the south. Flow processes include mass flow to turbidity currents. The size of the particles and the absence of fine-grained sediment suggest a transport over a short distance. The Lower Unit is made of coarse-grained tabular beds interpreted as non-channelled lobe deposits. The Upper Unit is made of massive conglomerates interpreted as the channelled part of lobes. These lobe deposits settle in a tectonically confined basin according to topographic compensation that occurs from bed scale to unit scale. The abrupt progradation between the lower and the upper unit seems related to a major tectonic uplift in the area. This uplift is also suggested by a change in the petrographic nature of the source and an abrupt coarsening of the transported clasts.This field example allows providing high resolution analysis for depositional sedimentary sequences of terminal lobe deposits in a coarse-grained turbidite system. The outcrop analysis shows the lateral evolution of deposits and the system progradation allows a longitudinal analysis of facies evolution by superposing on the same outcrops the channelled lobe system and the non-channelled lobe system. These results of high-resolution outcrop analysis can be extrapolated to results obtained on sedimentary lobes in recent deep-sea turbidite system that are either restricted to cores, or with a lesser resolution (seismic).     

K–Ar dating and δO–δD tracing of illitization within and outside the Shea Creek uranium prospect,Athabasca Basin,Canada

Isotope analyses (K–Ar, δ¹⁸O and δD) were performed on illite from both the sandstone cover and the underlying basement, close to and distant from Shea Creek, an unconformity-type U deposit (Athabasca Basin, Canada); the illite had previously been characterized crystallographically. In the barren areas away from deposit, illite is mainly of the cis-vacant 1M polytype occurring as relatively coarse-grained lath-shaped particles, while it occurs as fine-grained particles of the trans-vacant 1M type next to and in the U mineralized strata. The tectonic-induced hydrothermal system that favored illite crystallization was multi-episodic 1453 ± 2, 1330 ± 20 and probably about 1235 Ma ago. These illite-forming episodes appear to have occurred contemporaneously to those favoring the concentration of the associated U oxides, which were dated independently by the U–Pb method in the Shea Creek deposits and elsewhere in the Athabasca Basin.     

Assessing climate change impacts,sea level rise and storm surge risk in port cities: a case study on Copenhagen

This study illustrates a methodology to assess the economic impacts of climate change at a city scale and benefits of adaptation, taking the case of sea level rise and storm surge risk in the city of Copenhagen, capital of Denmark. The approach is a simplified catastrophe risk assessment, to calculate the direct costs of storm surges under scenarios of sea level rise, coupled to an economic input–output (IO) model. The output is a risk assessment of the direct and indirect economic impacts of storm surge under climate change, including, for example, production and job losses and reconstruction duration, and the benefits of investment in upgraded sea defences. The simplified catastrophe risk assessment entails a statistical analysis of storm surge characteristics, geographical-information analysis of population and asset exposure combined with aggregated vulnerability information. For the city of Copenhagen, it is found that in absence of adaptation, sea level rise would significantly increase flood risks. Results call for the introduction of adaptation in long-term urban planning, as one part of a comprehensive strategy to manage the implications of climate change in the city. Mitigation policies can also aid adaptation by limiting the pace of future sea level rise.     

Practical aspects of frequency-domain identification of dynamic models of marine structures from hydrodynamic data

The motion response of marine structures in waves can be studied using finite-dimensional linear-time-invariant approximating models. These models, obtained using system identification with data computed by hydrodynamic codes, find application in offshore training simulators, hardware-in-the-loop simulators for positioning control testing, and also in initial designs of wave-energy conversion devices. Different proposals have appeared in the literature to address the identification problem in both time and frequency domains, and recent work has highlighted the superiority of the frequency-domain methods. This paper summarises practical frequency-domain estimation algorithms that use constraints on model structure and parameters to refine the search of approximating parametric models. Practical issues associated with the identification are discussed, including the influence of radiation model accuracy in force-to-motion models, which are usually the ultimate modelling objective. The illustration examples in the paper are obtained using a freely available MATLAB toolbox developed by the authors, which implements the estimation algorithms described.     

Late Glacial and Holocene avulsions of the Rio Pastaza Megafan (Ecuador–Peru): frequency and controlling factors

The geomorphological study by means of remote sensing imagery of the Rio Pastaza Megafan (Ecuador and northern Peru) reveals the traces of numerous avulsions. One hundred and eight avulsion sites have been defined. The location of these sites, the available radiocarbon ages as well as historical maps of the seventeenth century, enable us to propose an evolution history of the migration and avulsions of the Rio Pastaza since the Last Glacial Maximum. The first avulsions of the Río Pastaza occurred after the LGM in a zone close to and roughly parallel to the sudandean front, where the developed avulsion gave a distributive pattern to the ancient stream of the Río Pastaza in an area located between the modern Río Morona and Pastaza, where they caused the Rio Pastaza to develop a fan-like distributary pattern. This is interpreted as a response to thrust-related forelimb tilt, progressively shifting eastward the Rio Pastaza and the apex of the megafan. This sequence of events ended with the Great Diversion of the Rio Pastaza towards the modern Rios Corrientes and Tigre. Avulsions occurred in the Tigre-Corrientes Area between 9200 and 8,500 years Cal BP. Afterwards, the Río Pastaza was diverted to its present-day north–south course. This last significant avulsion occurred before AD 1691. In the area located between the modern Río Morona and Pastaza, avulsion frequency—probably overestimated—ranges between 100 and 200 years. In the Ríos Tigre and Corrientes area, avulsion frequency—probably underestimated—ranges from 300 to 400 years. Regional tectonics is likely to have triggered most of the avulsions in the Morona Pastaza area but its influence is restricted to this area. The factors controlling the avulsions in the Tigre-Corrientes area are less clear because the frequently described “hydrologic”-driven avulsion as observed in areas characterized by contrasted hydrologic cycles are inconsistent with the characteristics of the hydrologic cycles of the Rio Pastaza.