CrystalMom: a new model for the evolution of crystal size distributions in magmas with the quadrature-based method of moments

Nucleation and growth of crystals, and the resulting crystal size distribution, play a fundamental role in controlling the physical properties of magmas and consequently the dynamics of the eruptions. In the past decades, laboratory experiments demonstrated that size and shape of crystals strongly control the physical properties of magma and lava. Additionally, natural and experimental samples are usually characterized in terms of their crystal size distribution to link it with physical processes that are not directly observable, such as cooling or decompression mechanisms. In this paper, we present CrystalMoM, a new predictive model, based on the quadrature-based method of moments, developed for studying the kinetic of crystallization in volcanic systems. The quadrature-based method of moments, well established in the field of chemical engineering, represents a mesoscale modelling approach that rigorously simulates the space–time evolution of a distribution of particles, by considering its moments. The method is applied here, for the first time, for studying the equilibrium/disequilibrium crystallization in magma, modelling the temporal evolution of the moments of a crystal size distribution. The model, verified against numerical and experimental data, represents a valuable tool to infer the cooling and decompression rates from the crystal size distribution observed in natural samples.     

Seasonal dynamics of water and nutrient fluxes in an agricultural peatland

The American cranberry (Vaccinium macrocarpon Ait.) is an important part of the cultural heritage and economy of Southeastern Massachusetts, yet water quality concerns and wetland protection laws challenge its commercial production. Here, we report inputs and outputs of water, nitrogen (N), and phosphorus (P) for a 2.12‐ha cranberry bed over a 2‐year period from 2013 to 2015. Water‐budget analysis indicated that precipitation contributed 40%, floodwater 37%, irrigation 15%, and groundwater 8% of water inputs to the cranberry bed. Minor annual variation in surface water discharge (~90 mm·year^?1 or 3%) contrasted with large decreases in net (= outputs ? inputs) nutrient export, from 16.2 to 9.1 kg N·ha^?1·year^?1 for total (dissolved + suspended particulate) nitrogen (TN) and from 3.34 to 1.47 kg P·ha^?1·year^?1 for total phosphorus (TP) between Years 1 and 2. Annual variation in net TN and TP export was tied to decreases in spring and summer nutrient export and controlled by the combined effects of fertilizer management, soil biogeochemistry, and hydrology. The relatively high spring TN export in Year 1 was associated with coincident increases in soil temperature and rainfall. A second factor was the timing of fertilizer application, which occurred 1 day prior to a major summer storm (i.e., third largest daily rainfall since 1926) and was responsible for up to 15% and 9% of the Year 1 TN and TP export, respectively. Nutrient budgets, which balanced water and fertilizer inputs with water, fruit, and vegetative outputs, were consistent with the burial of 21.6 kg N·ha^?1·year^?1 and 7.27 kg P·ha^?1·year^?1. Field measurements indicated that burial would increase TN and TP in the shallow (0–5 cm) rooting zone by 14% and 6%, respectively, which seemed plausible based on the relatively young age of the bed (4–5 years) and new root growth patterns in Vaccinium plants.     

P wave velocity structure of upper and middle crust beneath the Three Gorges reservoir dam and adjacent regions

The three-dimensional P wave velocity structures of upper and middle crust beneath the Three Gorges reservoir dam and adjacent regions were reconstructed by means of the seismic tomography, using the absolute P-arrivals from the events observed by the Three Gorges seismic network. Our preliminary study of the deep structure shows: (1) the lithology of the Huangling anticline, integral and highly solidified, exhibits a uniformity. Down to a depth of 20 km, the Huangling anticline still shows high-velocity anomalies. Note that the characteristics of the velocity images with the different depths are not the same; (2) the depth of the sedimentary rocks beneath the entire Zigui basin is not constant, which is less than 6 km; (3) the top boundary of the crystalline basement varies from 0 to 7 km, whereas the bottom one between 14 and 16 km; (4) the partially molten substances of the upper mantle are located in the Badong-Zigui and its adjacent regions as well as nearby Peishi area. The development and evolution of the current tectonics are associated with the cooling and solidifying process of the molten substances within the upper mantle; (5) the different velocity structures of the deep crust in both east and west side of the Jiuwanxi fault are to be investigated in order to understand the origin of gravitational gradient zone in the Three Gorge areas; and (6) the velocities within the upper crust would be affected by reservoir water permeability. Supported by National Natural Science Foundation of China (Grant No. 40574039), National Key Technology R& D Program (Grant No. 2008BAC38B02)     

Warm frontal structure in association with a rapidly deepening extratropical cyclone

Abstract

The detailed characteristics of a CASPII warm frontal passage are presented in this article. This storm, Intensive Operating Period (IOP) 13 (February 26–27, 1992), was observed in detail with an array of diverse instruments. It has the advantage over earlier freezing precipitation studies of having simultaneous, in situ and remote sensing measurements by aircraft and ground‐based Doppler radar.

The associated precipitation was in the form of banded structures parallel to the front. Within these bands were embedded precipitation cores, some parallel to the band, some perpendicular. The warm front itself was characterized by major perturbations in its kinematic and thermodynamic features. The cores oriented parallel to the front were the result of embedded convection generated, at least in part, by the irregularities in the frontal surface.

The cores oriented perpendicular to the front were closely associated with the 0°C isotherm on the underside of the frontal inversion. Precipitation phase changes played a significant role in the occurrence of wide near 0°C regions, both vertically and horizontally. These regions had a profound influence on the observed precipitation types and led to complex precipitation‐thermodynamic‐dynamic interactions. Instabilities produced by these interactions are seen in wave‐like features observed by the Doppler radar in these regions, both parallel and perpendicular to the frontal zone.     

Hydrogeochemical investigations and groundwater provinces of the Friuli Venezia Giulia Plain aquifers,northeastern Italy

Water resources are a key factor, particularly for the planning of the sustainable regional development of agriculture, as well as for socio-economic development in general. A hydrochemical investigation was conducted in the Friuli Venezia Giulia aquifer systems to identify groundwater evolution, recharge and extent of pollution. Temperature, pH, electric conductivity, total dissolved solids, alkalinity, total hardness, SAR, Ca²⁺, Na⁺, K⁺, Mg²⁺, Cl⁻, SO₄ ²⁻, NO₃ ⁻, HCO₃ ⁻ _, water quality and type, saturation indexes and the environmental stable isotope δ¹⁸O were determined in 149 sampling stations. The pattern of geochemical and oxygen stable isotope variations suggests that the sub-surface groundwater (from phreatic and shallow confined aquifers) is being recharged by modern precipitations and local river infiltrations. Four hydrogeological provinces have been recognised and mapped in the Friuli Venezia Giulia Plain having similar geochemical signatures. These provinces have different degrees of vulnerability to contamination. The deep confined groundwater samples are significantly less impacted by surface activities; and it appears that these important water resources have very low recharge rates and would, therefore, be severely impacted by overabstraction.     

Artificial Neural Network (ANN) Based Inversion of Benthic Substrate Bottom Type and Bathymetry in Optically Shallow Waters - Initial Model Results

Ocean-colour remote sensing in optically shallow waters is influenced by contribution from the water column depth as well as by the substrate type. Therefore, it is required to include the contribution from the water column and substrate bottom type for bathymetry estimation. In this report we demonstrate the use of Artificial Neural Network (ANN) based approach to spectrally distinguish various benthic bottom types and estimate depth of substrate bottom simultaneously in optically shallow waters. We have used in-water radiative transfer simulation modeling to generate simulated top-of-the-water column reflectance the four major benthic bottom types viz. sea grass, coral sand, green algae and red algae using Hydrolight simulation model. The simulated remote sensing reflectance, for the four benthic bottom types having benthic bottom depth up to 30 m were generated for moderately clear waters. A multi-layer perceptron (MLP) type neural network was trained using the simulated data. ANN based approach was used for classification of the benthic bottom type and simultaneous inversion of bathymetry. Simulated data was inverted to yield benthic bottom type classification with an accuracy of ~98% for the four benthic substrate types and the substrate depth were estimated with an error of 0% for sea grass, 1% for coral sand and 1–3% for green and red algae up to 25 m, whereas for substrate bottom deeper than 25 m depth the classification errors increased by 2–5% for three substrate bottom types except sea grass bottom type. The initial results are promising which needs validation using the in-situ measured remote sensing reflectance spectra for implementing further on satellite data.     

Organized Roll Circulation and Transport of Mineral Aerosols in the Atmospheric Boundary Layer

Izvestiya, Atmospheric and Oceanic Physics - An investigation into mesoscale roll circulation and its transport characteristics in the atmospheric boundary layer (ABL) is carried out. The case...     

Distribution and transport of organic matter along a filament-eddy system in the Canaries – NW Africa coastal transition zone region

The distribution and transport of chlorophyll a (Chla), particulate (POC) and dissolved (DOC) organic carbon, and the respiratory ETS activity of the microplankton community were studied along a filament-eddy system located in the transition zone between the NW Africa upwelling and Canary Islands waters. Two independent filaments (F1 and F2) stemming from the coastal jet, between Cape Juby and Cape Bojador, merged about 100 km offshore, turning southward and onshore forced by the circulation of a recurrent oceanic cyclonic eddy. In general, the coastal upwelling waters presented higher Chla, but lower POC, DOC and ETS activity than filament waters. However, differences in organic carbon distribution and respiratory activity were observed among stations from the two filaments. The bio-chemical fields were strongly influenced by a complex sub-mesoscale hydrography resulting from the interaction of cyclonic and anticyclonic island eddies with the filaments. The combined F1 + F2 filament system transported 97.1 kg s⁻¹ of excess (non-refractory) total organic carbon (e-TOC), a value comparable to other published estimates from upwelling filaments in the NE Atlantic. About 90% of e-TOC was exported as DOC, since eddy re-circulation precluded the offshore transport of POC. Assuming that the calculated transport of e-TOC is representative of the annual average, the yearly offshore transport (3.1 x 10⁹ kg C) would represent about 25% of the upwelling primary production of the region of study.     

A tangential wind profile for simulating strong tropical cyclones with MM5

A new tangential wind profile for simulating strong tropical cyclones is put forward and planted into the NCARAFWA tropical cyclone bogussing scheme in MM5. The scheme for the new profile can make full use of the information from routine typhoon reports, including not only the maximum wind, but also the additional information of the wind speeds of 25.7 and 15.4 ms^-1 and their corresponding radii, which are usually provided for strong cyclones. Thus, the new profile can be used to describe the outer structure of cyclones more accurately than by using the earlier scheme of MM5 in which only the maximum wind speed is considered. Numerical experimental forecasts of two strong tropical cyclones are performed to examine the new profile. Results show that by using the new profile the prediction of both cyclones‘ intensity can be obviously improved, but the effects on the track prediction of the two cyclones are different. It seems that the new profile might be more suitable for strong cyclones with shifted tracks. However, the conclusion is drawn from only two typhoon cases, so more cases are needed to evaluate the new profile.     

Shadow mechanism and the opposition effect of brightness of atmosphereless celestial bodies

We consider the Irvine-Yanovistkii modification of the shadow model developed by Hapke for the opposition effect of brightness. The relation between the single scattering albedo ω and the transparency coefficient of particles κ is suggested to be used in the form κ = (1 ? ω) ⁿ , which allows the number of unknowns in the model to be reduced to two parameters (the packing density of particles g and ω) and the single-scattering phase function χ(α). The analysis of spectrophotometric measurements of the moon and Mars showed that the data on the observed opposition effect and the changes in the color index with the phase angle α well agree if the values of n = 0.25 and g = 0.4 (the moon) and 0.6 (Mars) are assumed in calculations. When being applied to asteroids of several types, this method also yielded a satisfactory agreement. For the E-type asteroids, the sets of parameters are [g = 0.6, ω = 0.6, A _g = 0.21, and q = 0.83] or [g = 0.3, ω = 0.4, A _g = 0.15, and q = 0.71] under the Martian single-scattering phase function; for the M-type asteroids, it is [g = 0.4, ω ≤ 0.1, A _g ≤ 0.075, and q ≤ 0.42] under the lunar single-scattering phase function; for the S-type asteroids, it is [g = 0.4, ω = 0.4, A _g = 0.28, and q = 0.49] under the lunar single-scattering phase function; and for the C-type asteroids, it is [g = 0.6, ω ≤ 0.1, A _g ≤ 0.075, and q = 0.43] under the modified lunar single-scattering phase function. The polarization measurements fulfilled by Gehrels et al. (1964) for the bright feature on the lunar surface, Copernicus (L = -20°08′, φ = +10°11′), at a phase angle α = 1.6° revealed the deviations in the position of the polarization plane from that typical for the negative branch. They were 22° and 12° in the G and I filters, respectively. At the same time, the deviation was within the error (±3°) in the U filter and for the dark feature Plato (L = -10°32′, φ = +51°25′), which can be caused by the coherent mechanism of the formation of the polarization peak.