Scales of columnar jointing in igneous rocks: field measurements and controlling factors

Columnar jointing is a common feature of solidified lavas, sills and dikes, but the factors controlling the characteristic stoutness of columns remain debated, and quantitative field observations are few in number. In this paper, we provide quantitative measurements on sizing of columnar joint sets and our assessment of the principal factors controlling it. We focus on (1) chemistry, as it is the major determinant of the physical (mechanical and thermal) properties of the lava, and (2) geology, as it influences the style of emplacement and lava geometry, setting boundary conditions for the cooling process and the rate of heat loss. In our analysis, we cover lavas with a broad range of chemical compositions (from basanite to phonolite, for six of which we provide new geochemical analyses) and of geological settings. Our field measurements cover 50 columnar jointing sites in three countries. We provide reliable, manually digitized data on the size of individual columns and focus the mathematical analysis on their geometry (23,889 data on side length, of which 17,312 are from full column sections and 3,033 data on cross-sectional area and order of polygonality). The geometrical observations show that the variation in characteristic size of columns between different sites exceeds one order of magnitude (side length ranging from 8 to 338 cm) and that the column-bounding polygons’ average order is less than 6. The network of fractures is found to be longer than required by a minimum-energy hexagonal configuration, indicating a non-equilibrium, geologically quick process. In terms of the development and characteristic sizing of columnar joint sets, our observations suggest that columns are the result of an interplay between the geological setting of emplacement and magma chemistry. When the geological setting constrains the geometry of the emplaced body, it exerts a stronger control on characteristic column stoutness. At unconstrained geometries (e.g. unconfined lava flows), chemistry plays the major role, resulting in stouter columns in felsic lavas and slenderer columns in mafic lavas.     

Quantitative analysis of numerically induced mixing in a coastal model application

In recent years, various attempts have been made to estimate the amount of numerical mixing in numerical ocean models due to discretisation errors of advection schemes. In this study, a high-resolution coastal model using the ocean circulationmodel GETM is applied to the Western Baltic Sea, which is characterised by energetic and episodic inflows of dense bottom waters originating from the Kattegat. The model is equipped with an easy-to-implement diagnostic method for obtaining the numerical mixing which has recently been suggested. In this diagnostic method, the physical mixing is defined as the mean tracer variance decay rate due to turbulent mixing. The numerical mixing due to discretisation errors of tracer advection schemes is defined as the decay rate between the advected square of the tracer variance and the square of the advected tracer, which can be directly compared to the physical variance decay. The source and location of numerical mixing is further investigated by comparing different advection schemes and analysing the amount of numerical mixing in each spatial dimension during the advection time step. The results show that, for the setup used, the numerically and physically induced mixing have the same orders of magnitude but with different vertical and horizontal distributions. As the main mechanism for high numerical mixing, vertical advection of tracers with strong vertical gradients has been identified. The main reason for high numerical mixing is due to bottom-following coordinates when density gradients, especially for regions of steep slopes, are advected normal to isobaths. With the bottom-following coordinates used here, the horizontal gradients are reproduced by a spurious sawtooth-type profile where strong advection through, but not along, the vertical coordinate levels occurs. Additionally, the well known relation between strong tracer gradients and high velocities on the one and high numerical mixing on the other hand is approved quantitatively within this work.     

One step back for a leap forward: toward operational measurements of elements at risk

This study presents the evaluation of 1 year of operational lightning forecasts provided for Europe, using the Weather Research and Forecasting model coupled with a cloud-top height-based lightning parameterization scheme. Three different convective parameterization schemes were employed for parameterizing sub-grid cloud-top heights and consequently driving the lightning scheme. Triggering of the lightning scheme was controlled by means of a model-resolved microphysics-based masking filter, while the formulation for deriving lightning flash rates was also modified, assuming a single “marine” equation instead of the original equations discriminating between continental and marine lightning. Gridded lightning observations were used for evaluating model performance on a dichotomous decision basis. Analysis showed that the lightning scheme is sensitive to the parameterization of convection. In particular, the Kain–Fritsch convective scheme was found to outperform the Grell–Devenyi and Grell–Freitas schemes, showing a statistically significant better performance with respect to lightning prediction. This was most evident during the warm season, while smaller differences among the schemes were recorded during the cold season. Further, for all examined convective schemes, it was found that the application of the masking filter is desirable for improving model performance in terms of lightning forecasting. Last, the reported results revealed that the refinement of the formulation of the lightning parameterization scheme, adhering to a “global” marine equation instead of distinguishing between land and sea lightning, may be necessary in order to obtain reliable lightning forecasts.     

Structure and evolutionary history of the solar system,III

Parts I and II of our analysis of the evolution of the solar system were devoted mainly to the mechanical processes. The present part (Part III) deals primarily with the plasma processes and the hydromagnetic aspects.Much confusion in the cosmogonic field is due to the treatment of the early phases of the evolution of a circumstellar medium by pre-hydromagnetic methods, or by erroneous application of magnetohydrodynamics. In order to reduce the speculative element as far as possible the present analysis tries to connect the cosmogonic processes as directly as possible to laboratory plasma physics and to space phenomena actually observed today (Section 10).Models of the Laplacian type have been made obsolete by magnetohydrodynamics. Furthermore they are in conflict with observations. A new model is suggested (Section 11).A plasma surrounding a rotating central body may attain a state of partial corotation which is determined by the balance between gravitation and the centrifugal force acting on a plasma in a dipole field. Condensation from a partially corotating plasma results in grains orbiting in ellipses withe=1/3 and finally accreting to bodies at 2/3 of the central distance of the point of condensation (Section 12).An application of the theory to the Saturnian rings and to the asteroidal belt shows that the falldown ratio 2/3 (derived from the geometry of a dipole field) is essential for the understanding of their structure. The structure of the groups of planets and satellites is also discussed but only in a preliminary way. The behavior of volatile substances is a major problem which still awaits an appropriate treatment (Section 13).     

Structure and evolutionary history of the solar system,I

1. Introduction and Survey. The method for studying the structure and evolution of the solar system is discussed. It is pointed out that theories that account for the origin of planets alone are basically insufficient. Instead one ought to aim for a general theory for the formation of secondary bodies around a central body, applicable both to planet and satellite formation. A satisfactory theory should not start from assumed properties of the primitive Sun, which is a very speculative subject, but should be based on an analysis of present conditions and a successive reconstruction of the past states.
2. Orbits of Planets and Satellites. As a foundation for the subsequent analysis, the relevant properties of planets and satellites are presented.
3. The Small Bodies. The motion of small bodies is influenced by non-gravitational forces. Collisions (viscosity) are of special importance for the evolution of the orbits. It is pointed out that the focusing property of a gravitational field (which has usually been neglected) leads to the formation of jet streams. The importance of this concept for the understanding of the comet-meteoroid relations and the structure of the asteroidal belt is shown.
4. Resonance Structure. A survey is given of the resonances in the solar system and their possible explanation. It is concluded that in many cases the resonances must already be produced at the times when the bodies formed. It is shown that resonance effects put narrow limits on the post-accretional changes of orbits.
5. Spin and Tides. Tidal effects on planetary spins and satellite orbits are discussed. It is very doubtful if any satellite except the Moon and possibly Triton has had its orbit changed appreciably by tidal effects. The isochronism of planetary and asteroidal spins is discussed, as well as its bearing on the accretional process.
6. Post-accretional Changes in the Solar System. The stability of the solar system and upper limits for changes in orbital and spin data are examined. It is concluded that much of the present dynamic structure has direct relevance to the primordial processes.

     

Does the complexity in temporal precipitation disaggregation matter for a lumped hydrological model?

ABSTRACT

Flood peaks and volumes are essential design variables and can be simulated by precipitation–runoff (P–R) modelling. The high-resolution precipitation time series that are often required for this purpose can be generated by various temporal disaggregation methods. Here, we compare a simple method (M1, one parameter), focusing on the effective precipitation duration for flood simulations, with a multiplicative cascade model (M2, 32/36 parameters). While M2 aims at generating realistic characteristics of precipitation time series, M1 aims only at accurately reproducing flood variables by P–R modelling. Both disaggregation methods were tested on precipitation time series of nine Swiss mesoscale catchments. The generated high-resolution time series served as input for P–R modelling using a lumped HBV model. The results indicate that differences identified in precipitation characteristics of disaggregated time series vanish when introduced into the lumped hydrological model. Moreover, flood peaks were more sensitive than flood volumes to the choice of disaggregation method.     

Mesozoic exhumation history and palaeolandscape of the Iberian Massif in eastern Galicia from apatite fission-track and (U+Th)/He data

Apatite fission-track (AFT) and (U+Th)/He (AHe) data, combined with time–temperature inverse modelling, reveal the cooling and exhumation history of the Iberian Massif in eastern Galicia since the Mesozoic. The continuous cooling at various rates correlates with variation of tectonic boundary conditions in the adjacent continental margins. The data provide constraints on the 10⁷ timescale longevity of a relict paleolandscape. AFT ages range from 68 to 174 Ma with mean track lengths of 10.7 ± 2.6 to 12.6 ± 1.8 μm, and AHe ages range from 73 to 147 Ma. Fastest exhumation (≈0.25 km/Ma) occurred during the Late Jurassic to Early Cretaceous main episode of rifting in the adjacent western and northern margins. Exhumation rates have decreased since then and have been approximately one order of magnitude lower. Across inland Galicia, the AFT data are consistent with Early Cretaceous movement on post-Variscan NE trending faults. This is coeval with an extensional episode offshore. The AHe data in this region indicate less than 1.7 km of denudation in the last 100 Ma. This low exhumation suggests the attainment of a mature landscape during Late Cretaceous post-rift tectonic stability, whose remains are still preserved. The low and steady rate of denudation prevailed across inland Galicia despite minor N–S shortening in the northern margin since ≈45 Ma ago. In north Galicia, rock uplift in response to NW strike-slip faulting since Early Oligocene to Early Miocene has caused insufficient exhumation (<3 km) to remove the Mesozoic cooling signal recorded by the AFT data.