On determining the noon polar cap boundary from SuperDARN HF radar backscatter characteristics

Previous work has shown that ionospheric HF radar backscatter in the noon sector can be used to locate the footprint of the magnetospheric cusp particle precipitation. This has enabled the radar data to be used as a proxy for the location of the polar cap boundary, and hence measure the flow of plasma across it to derive the reconnection electric field in the ionosphere. This work used only single radar data sets with a field of view limited to 2 h of local time. In this case study using four of the SuperDARN radars, we examine the boundary determined over 6 h of magnetic local time around the noon sector and its relationship to the convection pattern. The variation with longitude of the latitude of the radar scatter with cusp characteristics shows a bay-like feature. It is shown that this feature is shaped by the variation with longitude of the poleward flow component of the ionospheric plasma and may be understood in terms of cusp ion time-of-flight effects. Using this interpretation, we derive the time-of-flight of the cusp ions and find that it is consistent with approximately 1 keV ions injected from a subsolar reconnection site. A method for deriving a more accurate estimate of the location of the open-closed field line boundary from HF radar data is described.     

Solar Microwave Bursts from Electron Populations with a 'Broken’ Energy Spectrum

Usually the gyrosynchrotron emission of microwave bursts from electron populations with a power-law (PL) energy distribution has been considered under the assumption that the spectral index of the distribution is constant over a wide range of energies. Meanwhile, there is strong evidence, in particular from hard X-ray and -ray, but also from cm/mm wavelength radio observations, that in many solar flare events the spectrum of the emitting electrons is characterized by a significant hardening at energies above 100–500 keV. We present some examples of calculated microwave burst spectra at cm/mm wavelengths taking into account the above evidence. It is shown that a break in the energy spectrum of the PL electrons can indeed result in a spectral hardening sometimes observed in microwave bursts at frequencies above 10–30 GHz.     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Rock magnetic properties and palaeomagnetic results of sediments from a stone implement layer in the Bose Basin,Guangxi

Directional samples were taken to study rock magnetism and palaeomagnetic records from the Dong- sheng profile, which is 5 m thick and on the northwest edge of the Bose Basin. Mineralogy and rock magnetism of typical samples indicate that coarse granular titanomagnetite, and fine-grained hematite, superparamagnetic maghemite formed by pedogenesis are in the sediment, which has undergone many transformative processes during different stages of pedogenesis. Parallel samples were taken for thermal demagnetization (TH) (0 to 680℃) and alternating field (AF) demagnetization (0 to 80 mT) respectively. Experimental results of these two kinds of demagnetization illustrate that there are two or more magnetic components in the samples. Intensity of NRM decreases by almost 60% to 90% rapidly when the temperature ranges from 100℃ to 350℃, with a steady magnetic component. It is impossible to analyze the magnetic components at high temperature because those fluctuate widely when the temperature is higher than 400℃. Steady magnetic components from 100℃ to 350℃ indicate that the remanence was mainly carried by fine-grained hematite formed by pedogenesis, reflecting a change in the geomagnetic field while the magnetite was being oxidized into hematite by chemical weathering after deposition. The formative age of the sediments cannot be obtained by magnetic methods in this profile.     

Bizarre geochemistry of komatiites from the Crixás greenstone belt,Brazil

Komatiite lava flows in the Crixás greenstone belt, Goiás, Brazil, have textures and volcanic structures typical of Archean komatiites, but are geochemically most unusual. The flows are porphyritic and massive, or layered with spinifex upper parts and olivine cumulate lower parts. MgO contents range from 18 to 40%. In such lavas, only olivine (and minor chromite) can have crystallized, but neither major nor trace elements fall on olivine control lines. In MgO variation diagrams, CaO and Sr fall on lines with slopes steeper than olivine control lines; SiO₂, FeO, Na₂O, K₂O and Y show little systematic variation; Zr shows a large variation that does not correlate with MgO; and Al₂O₃ decreases markedly with decreasing MgO. The aberrant behaviour is highlighted by the REE (rare earth elements) in spinifex and olivine cumulate layers from three flows: in the spinifex layers, chondrite-normalized REE patterns are hump-shaped with maxima at Nd or Sm ((La/Sm)_N=0.6, (Gd/Yb)_N=1.6–2.3), whereas cumulate zones in the same flows have steadily sloping patterns, with LREE enriched relative to HREE ((La/Sm)_N=1.3, (Gd/Yb)_N=1.4). Neither normal magmatic processes acting during emplacement of the komatiites, nor thermal erosion and wall-rock assimilation can explain these effects, and we speculate that elements commonly thought of as “immobile” (e.g. Al, Zr, REE) migrated during hydrothermal alteration or metamorphism. A Pb-Pb whole rock isochron gave an age of 2,728±140 Ma and selected Sm-Nd analyses an apparent isochron age of 2,825±98 Ma (ɛ_Nd≈0). The Pb-Pb age is believed to be the approximate time of emplacement. Interpretation of the Sm-Nd data is complicated by the evidence of mobility of REE.     

Preliminary study on assessment of nutrient transport in the Taihu Basin based on SWAT modeling

With the Taihu Basin as a study area, using the spatially distributed and mechanism-based SWAT model, preliminary simulations of nutrient transport in the Taihu Basin during the period of 1995:_2002 has been carried out. The topography, soil, meteorology and land use with industrial point pollution discharge, the loss of agricultural fertilizers, urban sewerage, and livestock drainages were all considered in the boundary conditions of the simulations. The model was calibrated and validated against water quality monitoring data from 2001 to 2002. The results show that the annual total productions of nitrogen (TN) and phosphorus (TP) into Lake Taihu are 40000t and 2000t respectively. Nutrient from the Huxi Region is a major resource for Lake Taihu. The non-point source (surface source) pollution is the main form of catchment sources of nutrients into Lake Taihu, occupied TN 53% and TP 56% respectively. TN and TP nutrients from industrial point pollution discharge are 30% and 16%, and sewerage in both forms of point source and non-point source are TN 31% and TP 47%. Both the loss of agricultural fertilizers and livestock drainages from the catchment should be paid more attention as an important nutrient source. The results also show that SWAT is an effective model for the simulation of temporally and spatially nutrient changes and for the assessment of the trends in a catchment scale.     

Empirical ground-motion models for point- and extended-source crustal earthquake scenarios in Europe and the Middle East

This article presents the latest generation of ground-motion models for the prediction of elastic response (pseudo-) spectral accelerations, as well as peak ground acceleration and velocity, derived using pan-European databases. The models present a number of novelties with respect to previous generations of models (Ambraseys et al. in Earthq Eng Struct Dyn 25:371–400, 1996, Bull Earthq Eng 3:1–53, 2005; Bommer et al. in Bull Earthq Eng 1:171–203, 2003; Akkar and Bommer in Seismol Res Lett 81:195–206, 2010), namely: inclusion of a nonlinear site amplification function that is a function of $\text{ V }_\mathrm{S30}$ and reference peak ground acceleration on rock; extension of the magnitude range of applicability of the model down to $\text{ M }_\mathrm{w}$ 4; extension of the distance range of applicability out to 200 km; extension to shorter and longer periods (down to 0.01 s and up to 4 s); and consistent models for both point-source (epicentral, $\text{ R }_\mathrm{epi}$ , and hypocentral distance, $\text{ R }_\mathrm{hyp}$ ) and finite-fault (distance to the surface projection of the rupture, $\text{ R }_\mathrm{JB}$ ) distance metrics. In addition, data from more than 1.5 times as many earthquakes, compared to previous pan-European models, have been used, leading to regressions based on approximately twice as many records in total. The metadata of these records have been carefully compiled and reappraised in recent European projects. These improvements lead to more robust ground-motion prediction equations than have previously been published for shallow (focal depths less than 30 km) crustal earthquakes in Europe and the Middle East. We conclude with suggestions for the application of the equations to seismic hazard assessments in Europe and the Middle East within a logic-tree framework to capture epistemic uncertainty.     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

An Improved Approach for Parameterizing Surface-Layer Turbulent Transfer Coefficients in Numerical Models

Based on classic iterative computation results, new equations to calculate the surface turbulent transfer coefficients are proposed, which allow for large ratios of the momentum and heat roughness lengths. Compared to the Launiainen scheme, our proposed scheme generates results closer to classical iterative computations. Under unstable stratification, the relative error in the Launiainen scheme increases linearly with increasing instability, even exceeding 15%, while the relative error of the present scheme is always less than 8.5%. Under stable stratification, the Launiainen scheme uses two equations, one for 0 < Ri _B ≤ 0.08 and another for 0.08 < Ri _B ≤ 0.2, and does not consider the condition that Ri _B > 0.2, while its relative errors in the region 0 < Ri _B ≤ 0.2 exceed 31 and 24% for momentum and heat transfer coefficients, respectively. In contrast, the present scheme uses only one equation for 0 < Ri _B ≤ 0.2 and another equation for Ri _B > 0.2, and the relative error of the present scheme is always less than 14%.     

Sources of the Late Mesozoic magmatic associations in the northeastern part of the Amurian microcontinent

Based on generalization of available geochronological data, Late Mesozoic magmatic associations in the northeastern part of the Amurian microcontinent are divided into three groups: 142–125, 124–115, and <110 Ma. The age of these associations decreases with approaching the Pacific margin of Asia. In the same direction, they show a change in sources of their parental melts: continental crust (142–125 Ma) → continental crust + PREMA (DM) (124–115 Ma) → continental crust + PREMA (DM) + EMII (<110 Ma). Isotope-geochemical (Sr-Nd) study indicates that intrusive and volcanic rocks of the Late Mesozoic magmatic associations in the northeastern part of the Amurian microcontinent were originated in geodynamic settings that provided access of enriched mantle sources to magma formation. The most probable of these settings are as follows: (1) plate sliding accompanying by the formation of slab window beneath continental margin; (2) passage of the Asian margin over the East Asian mantle hot field in the Late Mesozoic; (3) asthenospheric upwelling due to delamination of the lower crust during closure of the Mongolian-Okhotsk ocean caused by collision between the Amurian microcontinent, Dzhugdzhur-Stanovoy, and Selenga-Stanovoy superterranes in the Central Asian fold belt.