Entropy production and coarse graining of the climate fields in a general circulation model

We extend the analysis of the thermodynamics of the climate system by investigating the role played by processes taking place at various spatial and temporal scales through a procedure of coarse graining. We show that the coarser is the graining of the climatic fields, the lower is the resulting estimate of the material entropy production. In other terms, all the spatial and temporal scales of variability of the thermodynamic fields provide a positive contribution to the material entropy production. This may be interpreted also as that, at all scales, the temperature fields and the heating fields resulting from the convergence of turbulent fluxes have a negative correlation, while the opposite holds between the temperature fields and the radiative heating fields. Moreover, we obtain that the latter correlations are stronger, which confirms that radiation acts as primary driver for the climatic processes, while the material fluxes dampen the resulting fluctuations through dissipative processes. We also show, using specific coarse-graining procedures, how one can separate the various contributions to the material entropy production coming from the dissipation of kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes, without resorting to the full three-dimensional time dependent fields. We find that most of the entropy production is associated to irreversible exchanges occurring along the vertical direction, and that neglecting the horizontal and time variability of the fields has a relatively small impact on the estimate of the material entropy production. The approach presented here seems promising for testing climate models, for assessing the impact of changing their parametrizations and their resolution, as well as for investigating the atmosphere of exoplanets, because it allows for evaluating the error in the estimate of their thermodynamical properties due to the lack of high-resolution data. The findings on the impact of coarse graining on the thermodynamic fields on the estimate of the material entropy production deserve to be explored in a more general context, because they provide a way for understanding the relationship between forced fluctuations and dissipative processes in continuum systems.     

Revisiting the Concentration Observations and Source Apportionment of Atmospheric Ammonia

正While China’s Air Pollution Prevention and Control Action Plan on particulate matter since 2013 has reduced sulfate significantly, aerosol ammonium nitrate remains high in East China. As the high nitrate abundances are strongly linked with ammonia, reducing ammonia emissions is becoming increasingly important to improve the air quality of China. Although satellite data provide evidence of substantial increases in atmospheric ammonia concentrations over major agricultural regions, long-term surface observation of ammonia concentrations are sparse. In addition, there is still no consensus on     

A GIS-based approach for gully erosion susceptibility modelling: a test in Sicily, Italy

The aim of this study is to analyze the susceptibility conditions to gully erosion phenomena in the Magazzolo River basin and to test a method that allows for driving the factors selection. The study area is one of the largest (225 km²) watershed of southern Sicily and it is mostly characterized by gentle slopes carved into clayey and evaporitic sediments, except for the northern sector where carbonatic rocks give rise to steep slopes. In order to obtain a quantitative evaluation of gully erosion susceptibility, statistical relationships between the spatial distributions of gullies affecting the area and a set of twelve environmental variables were analyzed. Stereoscopic analysis of aerial photographs dated 2000, and field surveys carried out in 2006, allowed us to map about a thousand landforms produced by linear water erosion processes, classifiable as ephemeral and permanent gullies. The linear density of the gullies, computed on each of the factors classes, was assumed as the function expressing the susceptibility level of the latter. A 40-m digital elevation model (DEM) prepared from 1:10,000-scale topographic maps was used to compute the values of nine topographic attributes (primary: slope, aspect, plan curvature, profile curvature, general curvature, tangential curvature; secondary: stream power index; topographic wetness index; LS-USLE factor); from available thematic maps and field checks three other physical attributes (lithology, soil texture, land use) were derived. For each of these variables, a 40-m grid layer was generated, reclassifying the topographic variables according to their standard deviation values. In order to evaluate the controlling role of the selected predictive variables, one-variable susceptibility models, based on the spatial relationships between each single factor and gullies, were produced and submitted to a validation procedure. The latter was carried out by evaluating the predictive performance of models trained on one half of the landform archive and tested on the other. Large differences of accuracy were verified by computing geometric indexes of the validation curves (prediction and success rate curves; ROC curves) drawn for each one-variable model; in particular, soil texture, general curvature and aspect demonstrated a weak or a null influence on the spatial distribution of gullies within the studied area, while, on the contrary, tangential curvature, stream power index and plan curvature showed high predictive skills. Hence, predictive models were produced on a multi-variable basis, by variously combining the one-variable models. The validation of the multi-variables models, which generally indicated quite satisfactory results, were used as a sensitivity analysis tool to evaluate differences in the prediction results produced by changing the set of combined physical attributes. The sensitivity analysis pointed out that by increasing the number of combined environmental variables, an improvement of the susceptibility assessment is produced; this is true with the exception of adding to the multi-variables models a variable, as slope aspect, not correlated to the target variable. The addition of this attribute produces effects on the validation curves that are not distinguishable from noise and, as a consequence, the slope aspect was excluded from the final multi-variables model used to draw the gully erosion susceptibility map of the Magazzolo River basin. In conclusion, the research showed that the validation of one-variable models can be used as a tool for selecting factors to be combined to prepare the best performing multi-variables gully erosion susceptibility model.     

Mitochondrial DNA and microsatellites reveal significant divergence in the beachflea Orchestia montagui (Talitridae: Amphipoda)

Talitrids are semiterrestrial crustacean amphipods inhabiting sandy and rocky beaches; they generally show limited active dispersal over long distances. In this study we assessed levels of population genetic structure and variability in the talitrid amphipod Orchestia montagui, a species strictly associated to stranded decaying heaps of the seagrass Posidonia oceanica. The study is based on six populations (153 individuals) and covers five basins of the Mediterranean Sea (Tyrrhenian, Ionian, Adriatic, Western and Eastern basins). Samples were screened for polymorphisms at a fragment of the mitochondrial DNA (mtDNA) coding for the cytochrome oxidase subunit I gene (COI; 571 base pairs) and at eight microsatellite loci. MtDNA revealed a relatively homogeneous haplogroup, which clustered together the populations from the Western, Tyrrhenian and Eastern basins, but not the populations from the Adriatic and Ionian ones; microsatellites detected two clusters, one including the Adriatic and Ionian populations, the second grouping all the others. We found a weak geographic pattern in the genetic structuring of the species, with a lack of isolation by distance at either class of markers. Results are discussed in terms of probability of passive dispersal over long distances through heaps of seagrass.     

Statistical correlation of earthquake and ambient noise spectral ratios

The Horizontal-to-Vertical Spectral Ratio from earthquake (HVSR) and from ambient noise (HVN) recordings realistically indicate the fundamental frequency of soil response but, for the majority of the worldwide examined sites, they do not provide reliable amplification curves as calculated by the earthquake standard Spectral Ratio (SSR). Given the fact that HVSR and especially HVN can be easily obtained, it is challenging to search for a meaningful correlation with SSR amplification functions for the entire frequency band and to use the results for the SSR estimate at a further site where only noise measurements are available. To this aim we used recordings from 75 sites worldwide and we applied a multivariate statistical approach (canonical correlation analysis) to investigate and quantify any correlation among spectral ratios. The canonical correlation between SSR and HVN is then used to estimate the expected SSR at each site by a weighted average of the SSR values measured at the other sites; the weights are properly set to account more for sites with similar behaviour in terms of the canonical correlation results between HVN and SSR. This procedure, repeated for all sites in turn, constitutes the basis of a cross validation. The comparison between the inferred and the original SSR highlights the improvements of site response estimation with respect to the use of ambient noise techniques. The goodness and limitations of the reconstruction procedure are explained by specific geological settings.     

Evaluating fracture patterns within a resurgent caldera: Campi Flegrei, Italy

Understanding deformation of active calderas allows their dynamics to be defined and their hazard mitigated. The Campi Flegrei resurgent caldera (Italy) is one of the most active and hazardous volcanoes in the world, characterized by post-collapse resurgence, eruptions, ground deformation, and seismicity. An original structural analysis provides an overview of the main fracture zones. NW-SE and NE-SW fractures (normal or transtensive faults and extensional fractures) predominate along the rim and within the caldera, suggesting a regional control, both during and after the collapses. While the NE-SW fractures are ubiquitous in the deposits of the last ∼37 ka, NW-SE fractures predominate in the last 4.5 ka, during resurgence. The most recently (<4.5 ka) strained area lies in the caldera center (Solfatara area), where the faults, with an overall ∼ENE-WSW extension direction, appear to be associated with the bending due to resurgence. Solfatara lies immediately to the east of the most uplifted part of the caldera (Pozzuoli area), where domes form and culminate both on the long-term (resurgence, accompanied by volcanic activity) and short-term deformation (1982–1984 bradyseism, accompanied by seismic and hydrothermal activity). Similar volcano-tectonic behavior characterizes the short- and long-term uplifts, and only the intensity of the tectonic and volcanic activity varies, being related to varying amounts of uplift. Seismicity and hydrothermal manifestations occur during the bradyseisms, with moderate uplift, while surface faulting and eruptions occur during resurgence, with higher uplift. The features observed at Campi Flegrei are found at other major calderas, suggesting consistent behavior of large magmatic systems.     

Effects of 4-nonylphenol exposure in mussels (Mytilus galloprovincialis) and crabs (Carcinus aestuarii) with particular emphasis on vitellogenin induction

Since it is often difficult to estimate possible adverse effects due to contamination in selected ecosystems, multi-species biomonitoring may provide more information, taking into account different routes of exposure, ecological roles and metabolic capabilities of animals. In this context, we exposed for 7 days the mussel Mytilus galloprovincialis and the crab Carcinus aestuarii to 4-nonylphenol (NP), a well-known xenoestrogen. In mussels (0–0.2 mg NP l⁻¹), we measured NP bioaccumulation in soft tissues and vitellogenin (Vg)-like protein levels in digestive glands from both males and females by the alkali-labile phosphate assay (ALP). As no reference data were available for crab exposure, the NP 96-h LC₅₀ value was previously determined. Then, in sublethally exposed (0–1.0 mg NP l⁻¹) male crabs, NP bioaccumulation and Vg levels were measured in hemolymph, gonads and digestive gland. Bioaccumulation of NP increased from 43 to 371 μg g⁻¹ d.w. in mussels, and from 3.6 to 37 μg g⁻¹ d.w. in crabs, depending on the NP concentration in water. Dose-dependent Vg-like protein induction was observed in both species, appearing to be related to NP bioaccumulation, although a partial decrease was recorded at the highest concentration tested. A similar trend was observed in both digestive gland and gonad of exposed crabs; Vg increased to a lesser extent, although significantly, in hemolymph. Results demonstrated that NP induces Vg synthesis both in male and female mussels, as well as in male crabs. On the basis of the responsiveness of both species investigated, a multi-species approach is indicated in biomonitoring programmes.     

Is the Linné impact crater morphology influenced by the rheological layering on the Moon's surface? Insights from numerical modeling

Linné is a simple crater, with a diameter of 2.23 km and a depth of 0.52 km, located in northwestern Mare Serenitatis. Recent high‐resolution data acquired by the Lunar Reconnaissance Orbiter Camera revealed that the shape of this impact structure is best described by an inverted truncated‐cone. We perform morphometric measurements, including slope and profile curvature, on the Digital Terrain Model of Linné, finding the possible presence of three subtle topographic steps, at the elevation of +20, ?100, and ?200 m relative to the target surface. The kink at ?100 m might be related to the interface between two different rheological layers. Using the iSALE shock physics code, we numerically model the formation of Linné crater to derive hints on the possible impact conditions and target physical properties. In the initial setup, we adopt a basaltic projectile impacting the Moon with a speed of 18 km s^?1. For the local surface, we consider either one or two layers, in order to test the influence of material properties or composite rheologies on the final crater morphology. The one‐layer model shows that the largest variations in the crater shape take place when either the cohesion or the friction coefficient is varied. In particular, a cohesion of 10 kPa marks the threshold between conical‐ and parabolic‐shaped craters. The two‐layer model shows that the interface between the two layers would be exposed at the observed depth of 100 m when an intermediate value (~200 m) for the upper fractured layer is set. We have also found that the truncated‐cone morphology of Linné might originate from an incomplete collapse of the crater wall, as the breccia lens remains clustered along the crater walls, while the high‐albedo deposit on the crater floor can be interpreted as a very shallow lens of fallout breccia. The modeling analysis allows us to derive important clues on the impactor size (under the assumption of a vertical impact and collision velocity equal to the mean value), and on the approximate, large‐scale preimpact target properties. Observations suggest that these large‐scale material properties likely include some important smaller scale variations, disclosed as subtle morphological steps in the crater walls. Furthermore, the modeling results allow advancing some hypotheses on the geological evolution of the Mare Serenitatis region where Linné crater is located (unit S14). We suggest that unit S14 has a thickness of at least a few hundreds of meters up to about 400 m.