Evaluating Source Area Contributions from Aircraft Flux Measurements Over Heterogeneous Land Using Large-Eddy Simulation

The estimation of spatial patterns in surface fluxes from aircraft observations poses several challenges in the presence of heterogeneous land cover. In particular, the effects of turbulence on scalar transport and the different behaviour of passive (e.g. water vapour) versus active (e.g. temperature) scalars may lead to large uncertainties in the source area/flux-footprint estimation for sensible (H) and latent (LE) heat-flux fields. This study uses large-eddy simulation (LES) of the land–atmosphere interactions to investigate the atmospheric boundary-layer (ABL) processes that are likely to create differences in airborne-estimated H and LE footprints. We focus on 32~m altitude aircraft flux observations collected over a study site in central Oklahoma during the Southern Great Plains experiment in 1997 (SGP97). Comparison between the aircraft data and traditional model estimates provide evidence of a difference in source area for turbulent sensible and latent heat fluxes. The LES produces reasonable representations of the observed fluxes, and hence provides credible evidence and explanation of the observed differences in the H and LE footprints. Those differences can be quantified by analyzing the change in the sign of the spatial correlation of the H and LE fields provided by the LES model as a function of height. Dry patterns in relatively moist surroundings are able to generate strong, but localized, sensible heating. However, whereas H at the aircraft altitude is still in phase with the surface, LE presents a more complicated connection to the surface as the dry updrafts force a convergence of the surrounding moist air. Both the observational and LES model evidence support the concept that under strongly advective conditions, H and LE measured at the top of the surface layer (≈50 m) can be associated with very different upwind source areas, effectively contradicting surface-layer self-similarity theory for scalars. The results indicate that, under certain environmental conditions, footprint models will need to predict differing source area/footprint contributions between active (H) and passive (LE) scalar fluxes by considering land-surface heterogeneity and ABL dynamics.     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

Site effects observed in alluvial basins: the case of Norcia (Central Italy)

In this work, we investigate the site amplification effects observed in the Norcia plain, Central Italy. Data from 30 selected local earthquakes (2 ≤ Ml ≤ 4.1) recorded by a temporary seismic network composed by 15 stations, are analyzed to determine the spatial variability of site effects. Both the Horizontal-to-Vertical spectral ratio and the Standard Spectral Ratio techniques are applied to estimate the site amplification effects. The results show that most of the sites in the valley are affected by strong amplifications (up to a factor of 20) in the frequency range 0.5–5 Hz. The value of the fundamental frequency of resonance is strictly dependent on the location within the basin and on the sediment thickness. Strong amplifications also affect the vertical components. The time-frequency analysis performed on a station located inside the basin shows the presence of a large spectral amplitudes after the S-wave phase, not observed on a station located on the bedrock, suggesting the presence of locally generated wave trains. Then, in agreement with earlier observations for other alluvial basins in Central Italy, 2D–3D effects play an important role in determining the site amplification effects in Norcia.     

Recent and active tectonics of the external zone of the Northern Apennines (Italy)

We present a comprehensive study of the recent and active tectonics of the external part of the Northern Apennines (Italy) by using morphotectonic, geological–structural, and stratigraphic analysis, compared with the current seismicity of the region. This analysis suggests that the external part of the Northern Apennines is characterised by presence of three major systems of Quaternary compressive structures corresponding to (1) the Apenninic watershed, (2) the Apennines–Po Plain margin (pede-Apenninic thrust front), and (3) the Emilia, Ferrara, and Adriatic Fold systems buried below the Po Plain. Geological data and interpreted seismic sections indicate a roughly N–S Quaternary deformation direction, with rates <2.5 mm/year. The shortening decreased since the Pliocene, when our data indicate compression in a NNW–SSE direction and rates up to 7 mm/year. The trend and kinematics of the structures affecting the Apennines–Po Plain margin and the Po Plain subsoil fit well the pattern of the current seismicity of the area, as well as recent GPS and geodetic levelling data, pointing to a current activity of these thrust systems controlled by an overall compressive stress field. Close to the Apenninic watershed, earthquake focal mechanisms indicate that shallow extension is associated to deep compression. The extensional events may be related to a secondary extensional stress field developing on the hangingwall of the thrust system affecting the Apenninic watershed; alternatively, this thrust system may have been recently deactivated and overprinted by active normal faulting. Deeper compressive events are related to the activity of both a major basement thrust that connects at surface with the pede-Apenninic thrust front and a major Moho structure.     

Compound and hybrid clinothems of the last lowstand Mid-Adriatic Deep: Processes,depositional environments,controls and implications for stratigraphic analysis of prograding systems

Clinoforms with a range of scales are essential elements of prograding continental margins. Different types of clinoforms develop during margin growth, depending on combined changes in relative sea level, sediment supply and oceanographic processes. In studies of continental margin stratigraphy, trajectories of clinoform ‘rollover’ points are often used as proxies for relative sea-level variation and as predictors of the character of deposits beyond the shelf-break. The analysis of clinoform dynamics and rollover trajectory often suffers from the low resolution of geophysical data, the small scale of outcrops with respect to the dimensions of clinoform packages and low chronostratigraphic resolution. Here, through high-resolution seismic reflection data and sediment cores, we show how compound clinoforms were the most common architectural style of margin progradation of the late Pleistocene lowstand in the Adriatic Sea. During compound clinoform development, the shoreline was located landward of the shelf-break. It comprised a wave-dominated delta to the west and a barrier and back-barrier depositional system in the central and eastern area. Storm-enhanced hyperpycnal flows were responsible for the deposition of a sandy lobe in the river mouth, whereas a heterolithic succession formed elsewhere on the shelf. The storm-enhanced hyperpycnal flows built an apron of sand on the slope that interrupted an otherwise homogeneous progradational mudbelt. Locally, the late lowstand compound clinoforms have a flat to falling shelf-break trajectory. However, the main phase of shelf-break bypass and basin deposition coincides with a younger steeply rising shelf-break trajectory. We interpret divergence from standard models, linking shelf-break trajectory to deep-sea sand deposition, as resulting from a great efficiency of oceanographic processes in reworking sediment in the shelf, and from a high sediment supply. The slope foresets had a large progradational attitude during the late lowstand sea-level rise, showing that oceanographic processes can inhibit coastal systems to reach the shelf-edge. In general, our study suggests that where the shoreline does not coincide with the shelf-break, trajectory analysis can lead to inaccurate reconstruction of the depositional history of a margin.     

Modelling of magnetic field measurements at Mercury

Mapping Mercury's internal magnetic field with a magnetometer in closed orbit around the planet will provide valuable information about its internal structure. By measuring magnetic field multipoles of order higher than the dipole we could, in principle, determine some properties, such as size and location, of the internal source. Here we try to quantify these expectations. Using conceptual models, we simulate the actual measurement during the BepiColombo mission, and then we analyze the simulated data in order to estimate the measurement errors due to the limited spatial sampling. We also investigate our ability to locate the field generating current system within the planet. Finally, we address the main limitation of our model, due to the presence of time-varying external magnetospheric currents.     

Magma-induced strain localization in centrifuge models of transfer zones

Scaled centrifuge experiments have been used to investigate the dynamic relations between deformation and magma distribution in rift-related transfer zones. The physical models were built using suitable analogue materials, such as sand to represent the brittle upper crust, various kinds of silicone mixtures to simulate the lower crust and upper mantle and glycerol to reproduce magma. Models simulated the development of transfer zones across pre-existing glycerol reservoirs placed at the base of the analogue continental crust. In plan view, different geometries, dimensions and positions of subcrustal reservoirs were reproduced in three different sets of experiments; to compare results, models were also performed without magma-simulating glycerol.Set 1 experiments, incorporating a narrow rectangular glycerol reservoir, show that the low-viscosity material is able to localise deformation into the overlying crust, giving rise to discrete transfer zones. This concentrated surface deformation corresponds at depth to major magma accumulation. Set 2 experiments, with an initial wide squared glycerol reservoir, show instead that deformation is distributed across the whole model surface, corresponding at depth to relatively minor magma accumulation. Set 3 experiments explored various positions of a small squared reservoir that invariably localised faulting in the overlying analogue brittle crust at the onset of model deformation.The overall model behaviour suggests that magma distribution at depth can effectively control the strain distribution in the overlying crust and the deformative pattern of transfer zones. Strain distribution, in turn, may control magma emplacement as localized deformation would favour major accumulation of magma at transfer zones. Coupled to a strong thermal weakening of the country rocks, this process may ultimately lead to a positive feedback interaction between magma and deformation.