High-Performance Parallel Solver for Integral Equations of Electromagnetics Based on Galerkin Method

A new parallel solver for the volumetric integral equations (IE) of electrodynamics is presented. The solver is based on the Galerkin method, which ensures convergent numerical solution. The main features include: (i) memory usage eight times lower compared with analogous IE-based algorithms, without additional restrictions on the background media; (ii) accurate and stable method to compute matrix coefficients corresponding to the IE; and (iii) high degree of parallelism. The solver’s computational efficiency is demonstrated on a problem of magnetotelluric sounding of media with large conductivity contrast, revealing good agreement with results obtained using the second-order finite-element method. Due to the effective approach to parallelization and distributed data storage, the program exhibits perfect scalability on different hardware platforms.     

General patterns of circulation, sediment fluxes and ecology of the Palamós (La Fonera) submarine canyon, northwestern Mediterranean

Currents, particle fluxes and ecology were studied in the Palamós submarine canyon (also known as the Fonera canyon), located in the northwestern Mediterranean. Seven mooring arrays equipped with current meters and sediment traps were deployed along the main canyon axis, on the canyon walls and on the adjacent slope. Additionally, local and regional hydrographic cruises were carried out. Current data showed that mean near surface and mid-depth currents were oriented along the mean flow direction (NE–SW), although at 400 and 1200 m depth within the canyon current reversals were significant, indicating a more closed circulation inside the canyon. Mean near-bottom currents were constrained by the local bathymetry, especially at the canyon head. The most significant frequency at all levels was the inertial frequency. A second frequency of about three days, attributed to a topographic wave, was observed at all depths, suggesting that this wave was probably not trapped near the bottom. The current field observed during the most complete survey revealed a meandering pattern with cyclonic vorticity just upstream from and within the canyon. The associated vertical velocity ranged between 10 and 20 m/day and was constrained to the upper 300 m. This latter feature, together with other computations, suggests that during this survey the meander was not induced by the canyon but by some kind of instability of the mean flow.In the canyon, suspended sediment concentration, downward particle fluxes, chlorophyll and particulate C and N were significantly higher up-canyon from about 1200 m depth than offshore, defining, along with the different hydrodynamics, two canyon domains: one from the canyon head to about 1200 m depth more affected by the canyon confinement and the other deeper than 1200 m depth more controlled by the mean flow and the shelf-slope front. The higher near-bottom downward total mass fluxes were recorded in the canyon axis at 1200 m depth along with sharp turbidity increases and are related to sediment gravity flows. During the deployment period, the increase in downward particle fluxes occurred by mid-November, when a severe storm took place. On the canyon walls at 1200 m depth, suspended sediment concentrations, downward particle fluxes, chlorophyll and particulate C and N were higher on the southern wall than on the northern wall inversely to the current’s energy. This could be caused by an upward water supply on the southern canyon wall and/or the mean flow interacting with the canyon bathymetry. In the swimmers collected by the sediment traps, the dominant species was an elasipod holothurian, which has not been recorded in other canyons or elsewhere in the Mediterranean, indicating particular speciation.     

Picoplankton community structure on the Atlantic Meridional Transect: a comparison between seasons

Samples collected from 10 depths at 25 stations in September–October 1996 and 12 depths at 28 stations in April–May 1997 on an Atlantic Meridional Transect between the British Isles and the Falkland Islands were analysed by flow cytometry to determine the numbers and biomass of four categories of picoplankton: Prochlorococcus spp, Synechococcus spp, picoeukaryotic phytoplankton and heterotrophic bacteria. The composition of the picoplankton communities confirmed earlier findings (Zubkov, Sleigh, Tarran, Burkill & Leakey, 1998) about distinctive regions along the transect and indicated that the stations should be grouped into five provinces: northern temperate, northern Atlantic gyre, equatorial, southern Atlantic gyre and southern temperate, with an intrusion of upwelling water off the coast of Mauritania between the northern Atlantic gyre and equatorial waters. Prochlorococcus was the most numerous phototrophic organism in waters of both northern and southern gyres and in the equatorial region, at concentrations in excess of 0.1×10⁶ml⁻¹; it also dominated plant biomass in the gyres, but the biomass of the larger picoeukaryotic algae equalled that of Prochlorococcus in the equatorial region; higher standing stocks of both Prochlorococcus and picoeukaryotes were present in spring than in autumn in waters of both gyres. In temperate waters at both ends of the transect the numbers and biomass of picoeukaryotes and, more locally, of Synechococcus increased, and the Synechococcus, particularly, were more numerous in spring than in autumn. There was a pronounced southward shift of the main populations of both Synechococcus and Prochlorococcus in April–May in comparison to those of September–October, associated with seasonal changes in solar radiation, the abundance of Prochlorococcus dropping sharply near the 17°C contour, while Synechococcus was still present at temperatures below 10°C. Picoeukaryotes were more tolerant of low temperatures and lower light levels, often being more abundant in samples from greater depths, where they contributed to the deep chlorophyll maximum. Heterotrophic bacterial numbers and biomass tended to be highest in those samples where phototrophic biomass was greatest, with peaks in temperate and equatorial waters, which were shifted southwards in April–May compared with September–October.     

Polarization characteristics of Jensch coelostat

The polarization characteristics of the Jensch coelostat, used at the Sayan observatory, are presented. Results of calculations of the instrumental polarization (IP) are compared with the observations. It is shown that structural features of the Jensch coelostat lead to a very important dependence of the Stokes's parameters, that characterize the IP, on solar coordinates.     

远震数据显示贝加尔裂谷带北穆伊斯克(Severomuysk)段地壳和上地幔的精细结构

通过贝加尔裂谷系统北穆伊斯克(Severomuysk)段的密集地震台站线性网络获得的远距离强震记录,并利用P波接收函数技术,揭示了地壳和上地幔顶部的复杂分层块状结构。横波速度的分布表明构成北穆伊斯克地壳的地块具有不同性质。这些地块的西部聚散和地壳下部的分层证实了该地区隆起的堆积-碰撞起源。位于西伯利亚克拉通变薄的倾斜边缘上的北穆伊斯克段解释了该地区地壳碰撞效应的强度。2015年地震的震中深度与Muyakan 凹陷地壳上部的明显速度差异存在令人信服的相关性。     

Appearance of the nonlinearity from the nonlocality in diffusion through multiscale fractured porous media

We shall consider diffusion or single-phase flow in a multiscale porous medium which represents an infinite set of self-similar double-porosity media. At each scale, the medium consists of a highly permeable network of connected channels and low-permeable blocks. The characteristic scale of heterogeneity is ε at the highest level of hierarchy, wherein ε is a small parameter. The ratio between the channel and block permeability at each scale is ε ². The process analyzed is described using a diffusion equation with an oscillating multiscale diffusion parameter. The macroscale behavior is of interest. The transition to the macroscale is performed by means of the two-scale homogenization procedure. One step of averaging at each level of hierarchy leads to the appearance of the memory terms in the averaged equation. The successive averaging steps lead to progressive memory accumulation, so at each step of averaging, the macroscale model changes its type, and even the result of the second step is unknown a priori. The objective was to determine the macroscopic limit model for the infinite number of scales. By the method of induction, we obtained the macroscale model for an arbitrary number of scales and its limit for the infinite hierarchy. The limit model represents the system of two equations with memory terms. The kernel of the memory operator is the solution of a nonlinear integro-differential equation. Its solution is obtained through Laplace transform.     

Geoarchaeological investigations at Studenoe,an Upper Paleolithic site in the Transbaikal Region,Russia

In 1998, geoarchaeological research was carried out at the Upper Paleolithic site of Studenoe in the Transbaikal Region of Russia. The site is divided into three loci situated on two terraces overlooking the Chikoi River. Alluvial sediments beneath the oldest terrace (T₂) consist of two depositional units. Radiocarbon ages indicate that aggradation of the T₂ fill began before 18,000 yr B.P. Alluvial sediments beneath the lower terrace (T₁) range in age from 13,000 to 10,000 yr B.P. and are divided into three depositional units. Both terraces are overlain by Holocene colluvium. Archaeological materials at Studenoe 1/1, 1/2, and 2 include dwellings, hearths, and thousands of bone and stone artifacts assigned to the late Upper Paleolithic through the Bronze Age. Evidence of microblade technology is present in all components of the site. Material from recent excavations of Paleolithic levels in the T₂ fill at Studenoe includes mobiliary art, bone needles, and a large dwelling with four hearth features. © 2003 Wiley Periodicals, Inc.     

How rough sea affects marine seismic data and deghosting procedures

Most seismic processing algorithms generally consider the sea surface as a flat reflector. However, acquisition of marine seismic data often takes place in weather conditions where this approximation is inaccurate. The distortion in the seismic wavelet introduced by the rough sea may influence (for example) deghosting results, as deghosting operators are typically recursive and sensitive to the changes in the seismic signal. In this paper, we study the effect of sea surface roughness on conventional (5–160 Hz) and ultra‐high‐resolution (200–3500 Hz) single‐component towed‐streamer data. To this end, we numerically simulate reflections from a rough sea surface using the Kirchhoff approximation. Our modelling demonstrates that for conventional seismic frequency band sea roughness can distort results of standard one‐dimensional and two‐dimensional deterministic deghosting. To mitigate this effect, we introduce regularisation and optimisation based on the minimum‐energy criterion and show that this improves the processing output significantly. Analysis of ultra‐high‐resolution field data in conjunction with modelling shows that even relatively calm sea state (i.e., 15 cm wave height) introduces significant changes in the seismic signal for ultra‐high‐frequency band. These changes in amplitude and arrival time may degrade the results of deghosting. Using the field dataset, we show how the minimum‐energy optimisation of deghosting parameters improves the processing result.     

Parameters of Tsunami Source Versus Earthquake Magnitude

This study employs an empirical-analytical approach in combination with Monte-Carlo method to establish relationships between earthquake moment magnitude and upper limits of tsunami source parameters: double-amplitude of vertical bottom deformation, displaced water volume, and potential energy of initial elevation. The approach is based on the Okada solution for a finite rectangular fault and empirical scaling laws for earthquake sources. Results are compared to empirical and theoretical relationships published previously. Parameters of some recent tsunami sources, those for which USGS provides slip distribution data, are considered in light of the established relationships.     

Permafrost distribution and temperature in the Elkon Horst,Russia

The Elkon Horst is a geological structure that consists of heterogeneous strata with highly variable geocryological and temperature conditions. Gaining accurate knowledge of permafrost distribution patterns within this structure is of both scientific and practical importance. In mountainous terrain, the ground thermal regime is controlled by both surface and subsurface conditions. Surface conditions include snow cover characteristics, the presence or absence of vegetation, vegetation density, etc.. In contrast, subsurface conditions involve rock lithology or petrography, density, quantity and depth of fissures, groundwater, etc.. This article examines ground thermal regimes in various geomorphological settings based on temperature measurement data from geotechnical boreholes. The occurrence and extent of permafrost were evaluated for the entire horst area using direct and indirect methods. The maximum permafrost thickness measured in the Elkon Horst is 330 m, and the estimated maximum is 450 m at higher elevations. Thermophysical properties were determined for the major rock types, and the geothermal heat flux was estimated for the study area. The thermal conductivities were found to vary from 1.47 to 4.20 W/(m·K), and the dry bulk densities to range between 2,236 kg/m~3 and 3,235 kg/m~3. The average geothermal heat flux was estimated to be 44 mW/m~2.