Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Families of periodic orbits in Hill’s problem with solar radiation pressure: application to Hayabusa 2

This work studies periodic solutions applicable, as an extended phase, to the JAXA asteroid rendezvous mission Hayabusa 2 when it is close to target asteroid 1999 JU3. The motion of a spacecraft close to a small asteroid can be approximated with the equations of Hill’s problem modified to account for the strong solar radiation pressure. The identification of families of periodic solutions in such systems is just starting and the field is largely unexplored. We find several periodic orbits using a grid search, then apply numerical continuation and bifurcation theory to a subset of these to explore the changes in the orbit families when the orbital energy is varied. This analysis gives information on their stability and bifurcations. We then compare the various families on the basis of the restrictions and requirements of the specific mission considered, such as the pointing of the solar panels and instruments. We also use information about their resilience against parameter errors and their ground tracks to identify one particularly promising type of solution.     

A model for oxygen and sulfur isotope fractionation in sulfate during bacterial sulfate reduction processes

We present a model of bacterial sulfate reduction that includes equations describing the fractionation relationship between the sulfur and the oxygen isotope composition of residual sulfate (δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}) and the amount of residual sulfate. The model is based exclusively on oxygen isotope exchange between cell-internal sulfur compounds and ambient water as the dominating mechanism controlling oxygen isotope fractionation processes. We show that our model explains δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns observed from natural environments and from laboratory experiments, whereas other models, favoring kinetic isotope fractionation processes as dominant process, fail to explain many (but not all) observed δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns. Moreover, we show that a “typical” δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} slope does not exist. We postulate that measurements of δ³⁴S_{SO4_residual} and δ¹⁸O_{SO4_residual} can be used as a tool to determine cell-specific sulfate reduction rates, oxygen isotope exchange rates, and equilibrium oxygen isotope exchange factors. Data from culture experiments are used to determine the range of sulfur isotope fractionation factors in which a simplified set of equations can be used. Numerical examples demonstrate the application of the equations. We postulate that, during denitrification, the oxygen isotope effects in residual nitrate are also the result of oxygen isotope exchange with ambient water. Consequently, the equations for the relationship between δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}, and the amount of residual sulfate could be modified and used to calculate the fractionation-relationship between δ¹⁵N_{NO3_residual}, δ¹⁸O_{NO3_residual}, and the amount of residual nitrate during denitrification.     

Paleogeography and structure of the central Mediterranean: Sicily and its offshore area

The geology of the mainland and offshore of Sicily is illustrated by a few geologic sections and seismic profiles across the late Cenozoic orogenic belt of central and western Sicily and across the Sardinia Channel and Sicily Straits. This belt is the result of several tectonic events. Deformation involved mainly the sedimentary cover of the old African continental margin characterized by a broad basinal domain, flanked along its external (southern) margin by a shallow-water carbonate platform attached to Africa in the Triassic. Compressional deformation started in the more internal basinal rock assemblages overlying a thinned crust. The most important structural characteristic of the early phase of thrusting is the duplex pile forming the bulk of the chain in Central Western Sicily. The structure consists of a basal allochthon, made up of Permian to Middle Triassic layers, an intermediate duplex wedge, composed of competent Mesozoic carbonates, and a roof complex, including Upper Mesozoic-Lower Tertiary less competent rocks. Large-scale clockwise rotation of the thrusts accompanied transpressional movements in the hinterland during the Pliocene. Right oblique reverse faults modified the previous tectonic contacts between the allochthons in the hinterland zones. Contemporaneous south-directed imbrications affected the southern external areas, progressively incorporating foreland and piggyback basins. Development of the Gela Thrust System appears to be linked to the transpressional event; its accretion is also related to contemporaneous underthrusting at deeper levels of Mesozoic carbonate substratum. The older buried thrust sheets were pushed up to the surface breaching the deformed Tertiary cover of the Gela TS. Northwards in the belt post-Messinian normal growth faults opened half graben whose sedimentary fill underwent structural inversion. Alternation of extension and compression tectonics characterizes the Sicilian continental margin in the last million years.     

The effect of Congo River freshwater discharge on Eastern Equatorial Atlantic climate variability

The surface ocean explains a considerable part of the inter-annual Tropical Atlantic variability. The present work makes use of observational datasets to investigate the effect of freshwater flow on sea surface salinity (SSS) and temperature (SST) in the Gulf of Guinea. In particular, the Congo River discharges a huge amount of freshwater into the ocean, affecting SSS in the Eastern Equatorial Atlantic (EEA) and stratifying the surface layers. The hypothesis is that an excess of river runoff emphasize stratification, influencing the ocean temperature. In fact, our findings show that SSTs in the Gulf of Guinea are warmer in summers following an anomalously high Congo spring discharge. Vice versa, when the river discharges low freshwater, a cold anomaly appears in the Gulf. The response of SST is not linear: temperature anomalies are considerable and long-lasting in the event of large freshwater flow, while in dry years they are less remarkable, although still significant. An excess of freshwater seems able to form a barrier layer, which inhibits vertical mixing and the entrainment of the cold thermocline water into the surface. Other processes may contribute to SST variability, among which the net input of atmospheric freshwater falling over EEA. Likewise the case of continental runoff from Congo River, warm anomalies occur after anomalously rainy seasons and low temperatures follow dry seasons, confirming the effect of freshwater on SST. However, the two sources of freshwater anomaly are not in phase, so that it is possible to split between atypical SST following continental freshwater anomalies and rainfall anomalies. Also, variations in air-sea fluxes can produce heating and cooling of the Gulf of Guinea. Nevertheless, atypical SSTs cannot be ascribed to fluxes, since the temperature variation induced by them is not sufficient to explain the SST anomalies appearing in the Gulf after anomalous peak discharges. The interaction processes between river runoff, sea surface salinity and temperature play an effective role in the interannual variability in the EEA region. Our results add a new source of variability in the area, which was often neglected by previous studies.     

The June 4th 1999 severe weather episode in San Quirino, Italy: a tornado event?

On the morning of June 4th 1999, a severe weather event took place in San Quirino, a small village of Friuli-Venezia Giulia in the northeast of Italy. This village is located near the piedmont of the Alps, 40 km west from Udine and 60 km north from Venice.Around 0900 UTC (1100 local time), a thunderstorm with an intense hail fall affected the area of San Quirino. A few minutes later (around 0920 UTC, source: a farmer), a funnel cloud from a cumulonimbus touched the ground, producing damages to houses, trees and sheds. The damaged area was quite narrow (about 300 m) and short (less than 10 km). No injuries to people were reported.In spite of the smallness of the area interested by the phenomenon, this storm is studied here starting from the synoptic scale, moving to the mesoscale and finishing with the storm scale, trying to underline its characteristics. These analyses, especially those coming from the Doppler radar images, bring us to the conclusion that the San Quirino episode was produced by a supercell storm.     

Large wood transport modelling by a coupled Eulerian–Lagrangian approach

The paper discusses a model which predicts the trajectory of floating rigid bodies and may be applied to compute the motion of woody “debris” mobilized during floods. The model couples a Discrete Element (DE) Lagrangian approach for the calculation of motion of rigid bodies with the Eulerian solution of the shallow water equations (SWE), in order to simulate the transport of a cylinder in a two-dimensional stream. It differs from existing models since it is based on a dynamic approach, adapting the Basset–Boussinesq–Oseen equation. In a first step, forces are computed from flow and log velocities; then, the equations of dynamics are solved to model the planar roto-translation of the body. Model results and physical reliability are clearly affected by the values of the drag and side coefficients, especially since logs, modelled as cylinders, are able to change their orientation towards the flow. Experimental studies to evaluate drag and side coefficients can be found in the literature for a submerged cylinder, with various orientations. To extend such results to the case of a floating log, the authors performed a series of laboratory tests on partially submerged cylinders, implementing the outcomes in the proposed DE-SWE model. The coupled model is validated against existing laboratory data concerning spheres and wooden cylinder transport.