Seismology software: state of the practice

We analyzed the state of practice for software development in the seismology domain by comparing 30 software packages on four aspects: product, implementation, design, and process. We found room for improvement in most seismology software packages. The principal areas of concern include a lack of adequate requirements and design specification documents, a lack of test data to assess reliability, a lack of examples to get new users started, and a lack of technological tools to assist with managing the development process. To assist going forward, we provide recommendations for a document-driven development process that includes a problem statement, development plan, requirement specification, verification and validation (V&V) plan, design specification, code, V&V report, and a user manual. We also provide advice on tool use, including issue tracking, version control, code documentation, and testing tools.     

Liquefaction of silts: the Adapazari criteria

Ground failures in the form of liquefaction, loss of bearing capacity and soil softening have been observed during the 1999 Marmara (Turkey) earthquake. Research to understand the failure phenomena has been carried out since the earthquakes. This paper attempts to provide explanations to the liquefaction failure of silts in seismic conditions. Findings from a large amount of data collected in the city of Adapazari on the physical and mechanical properties of soils is presented. A geomorphological map of the city has shown that there are surprising horizontal and vertical variations of the facies due to the activity of rivers in the past. Cases of liquefaction appear to have concentrated in former backswamp areas where silts and sandy silts were deposited by crevasse splays. Properties of the soils in zones of liquefaction and non-liquefaction have been determined down to a reasonable depth by measuring the average size, clay content and liquidity index as well as cone penetration resistances with porewater pressures to discover that there is significant discrepancy among those profiles susceptible to liquefaction and non-liquefying deposits. A set of “Adapazari Criteria” is proposed which is intended to improve over the “Chinese Criteria” and is simple enough to be universally applicable. This classification is similar to the existing criteria but emphasizes more on the clay content in addition to measuring the liquid limit and the liquidity index as well as the average size.     

Study of the air-sea interactions at the mesoscale: the SEMAPHORE experiment

The SEMAPHORE (Structure des Echanges Mer-Atmosphère, Propriétés des Hétérogénéités Océaniques: Recherche Expérimentale) experiment has been conducted from June to November 1993 in the Northeast Atlantic between the Azores and Madeira. It was centered on the study of the mesoscale ocean circulation and air-sea interactions. The experimental investigation was achieved at the mesoscale using moorings, floats, and ship hydrological survey, and at a smaller scale by one dedicated ship, two instrumented aircraft, and surface drifting buoys, for one and a half month in October-November (IOP: intense observing period). Observations from meteorological operational satellites as well as spaceborne microwave sensors were used in complement. The main studies undertaken concern the mesoscale ocean, the upper ocean, the atmospheric boundary layer, and the sea surface, and first results are presented for the various topics. From data analysis and model simulations, the main characteristics of the ocean circulation were deduced, showing the close relationship between the Azores front meander and the occurrence of Mediterranean water lenses (meddies), and the shift between the Azores current frontal signature at the surface and within the thermocline. Using drifting buoys and ship data in the upper ocean, the gap between the scales of the atmospheric forcing and the oceanic variability was made evident. A 2 °C decrease and a 40-m deepening of the mixed layer were measured within the IOP, associated with a heating loss of about 100 W m⁻². This evolution was shown to be strongly connected to the occurrence of storms at the beginning and the end of October. Above the surface, turbulent measurements from ship and aircraft were analyzed across the surface thermal front, showing a 30% difference in heat fluxes between both sides during a 4-day period, and the respective contributions of the wind and the surface temperature were evaluated. The classical momentum flux bulk parameterization was found to fail in low wind and unstable conditions. Finally, the sea surface was investigated using airborne and satellite radars and wave buoys. A wave model, operationally used, was found to get better results compared with radar and wave-buoy measurements, when initialized using an improved wind field, obtained by assimilating satellite and buoy wind data in a meteorological model. A detailed analysis of a 2-day period showed that the swell component, propagating from a far source area, is underestimated in the wave model. A data base has been created, containing all experimental measurements. It will allow us to pursue the interpretation of observations and to test model simulations in the ocean, at the surface and in the atmospheric boundary layer, and to investigate the ocean-atmosphere coupling at the local and mesoscales.     

Special Contribution: Pictorial Series of the Manifestations of the Dynamics of the Earth: 4. Historical Images of Rockbursts' Effects in Mines

Studia Geophysica et Geodaetica -     

On the fast movements of the thermosphere: an interpretation

We examine the oscillations of the meridional neutral wind in the F region as seen by the EISCAT radar. We propose an interpretation in term of eddies (tourbillons) of typical size of a few tens to a few hundreds of kilometers. The observed rotation velocity is a few hundreds of meters per second. We suggest that the tourbillons are a common feature of thermospheric movements. We propose an optical experiment to check the validity of this assumption.     

Microseisms: Its dependence on the movement of the source

Summary The displacements on the free surface of a homogeneous isotropic semi-infinite medium has been obtained assuming a moving normal stress on the surface. It has been studied how the displacements behave with (1) the direction of propagation of the source with respect to the recording station and (2) the velocity of the source.     

Special Contribution: Pictorial Series of the Manifestations of the Dynamics of the Earth: 3. Historical Images of Landslides and Rock Falls

Studia Geophysica et Geodaetica -     

Special Contribution: Pictorial Series of the Manifestations of the Dynamics of the Earth: 6. South Pacific and Antarctica - the Last Explored Regions

Studia Geophysica et Geodaetica -     

Hydrological implications of the Southern Oscillation: variability of the rainfall-runoff relationship

Abstract

A connection between El Niño Southern Oscillation (ENSO) and weather phenomena in eastern Australia has been recognized for several decades. However, little work has been devoted to addressing how this correlation affects hydrological system behaviour within regional-scale catchments. In this study, spatially distributed ENSO effects are evaluated in terms of monthly rainfall, evaporation, streamflow and runoff characteristics for a 1300 km² catchment. The catchment is located in southeastern Australia where previous studies have indicated only modest ENSO influences on rainfall variability. Spatial and temporal analysis indicates that strongest ENSO-induced rainfall variability occurs during summer months. Additionally, the strength of the relationship is variable in space indicating that topographic controls may affect ENSO influences on rainfall totals and intensities. However, analysis of runoff shows substantially magnified ENSO-induced variability in comparison to the induced variability in rainfall. This may be attributable to the nonlinearity of runoff generation. Differences in antecedent moisture storage conditions will exist but may also be enhanced by complementary ENSO influences on daily rainfall intensities and mean monthly evaporation and temperature totals. The degree of the nonlinearity displayed by the hydrometeorological processes presented demonstrates that the significance of ENSO forecasts for surface water resource management should be assessed with direct regard to streamflow generation rather than on the basis of rainfall totals alone.     

Defining the moment of erosion: the principle of thermal consonance timing

Geomorphological process research demands quantitative information on erosion and deposition event timing and magnitude, in relation to fluctuations in the suspected driving forces. This paper establishes a new measurement principle – thermal consonance timing (TCT) – which delivers clearer, more continuous and quantitative information on erosion and deposition event magnitude, timing and frequency, to assist understanding of the controlling mechanisms. TCT is based on monitoring the switch from characteristically strong temperature gradients in sediment, to weaker gradients in air or water, which reveals the moment of erosion. The paper (1) derives the TCT principle from soil micrometeorological theory; (2) illustrates initial concept operationalization for field and laboratory use; (3) presents experimental data for simple soil erosion simulations; and (4) discusses initial application of TCT and perifluvial micrometeorology principles in the delivery of timing solutions for two bank erosion events on the River Wharfe, UK, in relation to the hydrograph. River bank thermal regimes respond, as soil temperature and energy balance theory predicts, with strong horizontal thermal gradients (often >1 K cm^?1 over 6·8 cm). TCT fixed the timing of two erosion events, the first during inundation, the second 19 h after the discharge peak and 13 h after re‐emergence from the flow. This provides rare confirmation of delayed bank retreat, quantifies the time‐lag involved, and suggests mass failure processes rather than fluid entrainment. Erosion events can be virtually instantaneous, implying ‘catastrophic retreat’ rather than ‘progressive entrainment’. Considerable potential exists to employ TCT approaches for: validating process models in several geomorphological contexts; assisting process identification and improving discrimination of competing hypotheses of process dominance through high‐resolution, simultaneous analysis of erosion and deposition events and driving forces; defining shifting erodibility and erosion thresholds; refining dynamic linkages in event‐based sediment budget investigations; and deriving closer approximations to ‘true’ erosion and deposition rates, especially in self‐concealing scour‐and‐fill systems. Copyright © 2005 John Wiley & Sons, Ltd.     

Monitoring of the Arctic region: contribution of the Arkhangelsk seismic network

We evaluate the contribution of the Arkhangelsk seismic network to the recording of seismic events of the Arctic region. The evaluation is performed by means of the general analysis of the network seismic catalog, evaluation of the magnitude sensitivity, and comparative analysis of the Arkhangelsk network catalogs with NORSAR seismological data services and the Geophysical Service of the Russian Academy of Sciences. It is established that the Arkhangelsk network has a rather high sensitivity and records the intra-plate and inter-plate seismic activity in the Arctic Region. Due to the start of the seismic station operation on the Franz Josef Land Archipelago, there is a capability to record earthquakes, data on which are missing in the catalogs of well-known seismological data services.     

Multipole models of the geomagnetic field: Construction of the Nth approximation

The fact that the geomagnetic potential, as well as magnetic induction and its gradient, can be represented with an arbitrary accuracy in the scope of multipole models has been justified programmatically and analytically. The obtained expressions are brief and can be utilized in analytical and computer-based studies using methods of computer algebra and numerical computations. Previously unknown analytical expressions for the components of the first seven multipole tensors, expressed in terms of the Gaussian coefficients, have been obtained for applied problems of space dynamics. An algorithm that makes it possible to construct analytical expressions for the arbitrary-order multipole tensor and magnetic induction vector components in any finite approximation has been developed.     

Shielded by the wind: the influence of the interstellar medium on the environment of Earth

We report on the recent studies on the long-term influence of cosmic rays on the Earth's environment. While on short time-scales solar activity is the driver for atmospheric changes suspected to be due to cosmic rays, for long time-scales the heliosphere, i.e. the circumsolar region occupied by the expanding part of the Sun's atmosphere, has to be considered. The heliosphere is identified as an important shield against interstellar influences and hazards. It has been demonstrated by quantitative modelling that a change of the interstellar medium surrounding the heliosphere as a result of the Sun's quasi-Keplerian motion around the galactic center triggers significant changes of planetary environments caused by enhanced fluxes of neutral atoms as well as by the increased cosmic ray fluxes. We give a compilation of recent space science results of interest to the atmospheric science community.     

Modeling the formation of undulations of the coastline: The role of tides

An idealized model is developed and analyzed to investigate the relevance of tidal motion for the emergence of undulations of a sandy coastline. The model describes feedbacks between tidal and steady flow on the inner shelf, sand transport in the nearshore zone and an irregular coastline. It is demonstrated that an initially straight coastline can become unstable with respect to perturbations with a rhythmic structure in the alongshore direction. The mechanism causing the growth of these perturbations is explained in terms of vorticity concepts. The relative importance of tide-related and wave-driven sediment fluxes in generating undulations of the coastline is investigated for the Dutch coast. Using parameter values that are appropriate for the Dutch coast it is found that tides can render a straight coastline unstable. The model predicts a fastest growing mode (FGM) with a wavelength that is in the order of the observed length of barrier islands. The mode grows on a time scale of 50 yr and it migrates 200 m per year. The wavelength of the FGM decreases with increasing amplitude of the tidal currents. This result is consistent with data of tides, waves and the lengths of barrier islands that are located along the Dutch and German Wadden coast.     

Seismicity of the North of the Russian Plate: Relocation of Recent Earthquakese

The hypocenters of the earthquakes recorded in the north of the Russian Plate from 1982 to 2013 are relocated. The relocation of the hypocenters is based on the common velocity section, common methodology, and the entire set of the initial data and bulletins available from the Russian and foreign seismic stations. The efficiency of the algorithm for calculating the hypocentral parameters and the velocity section is demonstrated by the example of two nonmilitary nuclear explosions in July 18, 1985 and September 6, 1988 in the northern part of the European Russia. For the first time, two earthquakes of July 19, 1982 and October 7, 2012, which have not been previously reported in the catalogs for the north of the Russian plate, are included in the seismic catalog.     

GPR modelling by the Fourier method: improvement of the algorithm

We improve two aspects of the modelling scheme for the simulation of electromagnetic radio waves, based on the Fourier pseudospectral method.
When there are large contrasts in the material properties, use of the standard algorithm (regular grid) causes a series of artefacts, as, for instance, ringing and acausal events. These problems, due to the non-locality of the differential operator, are solved by using the staggered Fourier method (staggered grid).
Realistic radiation patterns can be obtained from simple combinations of magnetic and electric sources. If the directivity pattern of the antenna is known, from either a finite-difference simulation or an analytic evaluation or an experimental characterization, it can then be simulated by a composite-source concept. This effective source is implemented in the modelling algorithm by means of a perturbation technique, which first computes the intensity and directional spectra of the single electromagnetic sources. Their location is optimized to obtain the best fit with a minimum number of sources. The approach is, in principle, valid for the far-field radiation pattern of the antenna.     

Plagioclase ultraphyric basalts in Iceland: the mush of the rift

Glassy, plagioclase ultraphyric basalts from six locations in Iceland have bimodal phenocryst size distributions where microphenocrysts (ol+plg±cpx±mt) are in equilibrium with the matrix glass, but macrophenocrysts (ol+plg±cpx) are too primitive to be so. Matrix glass compositions are similar to those of other rift zone glasses from Iceland, and oxygen isotope variations suggest they interacted with the Icelandic crust. A lack of negative Eu-anomalies in matrix glasses precludes large amounts of plagioclase crystallisation from their parental liquids. Compositions of glass inclusions in plagioclase and olivine macrophenocrysts indicate that parental magma compositions of the macrophenocryst assemblage are similar to those of primitive, Icelandic rift zone glasses. Temperatures for plagioclase macrophenocryst crystallisation obtained from Linkam^® heating stage experiments, and from glass inclusion compositions corrected for post-entrapment crystallisation, give temperatures up to 1260°C, corresponding to crystallisation at middle to deep crustal levels. Temperature differences of less than 100°C between plagioclase-hosted glass inclusions before and after post-entrapment plagioclase crystallisation show that the macrophenocrysts must have been kept at elevated temperatures prior to incorporation in their present host magmas. We suggest that the macrophenocrysts of the plagioclase ultraphyric basalts accumulated in crystal mush bodies underneath the rift zone and were picked up by their present hosts during a rifting event with increased magma supply from the mantle.