Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes simulation

Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and physical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length z_0m and excess resistance κB^-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB^-1 at the two stations. Results indicate that z_0m has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of z_0m in the land surface model and that obtained from this study. κB^-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of z_0m and κB^-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.     

Influence of aquatic plants on the hydrogen isotope composition of sedimentary long-chain n-alkanes in the Lake Qinghai region,Qinghai-Tibet Plateau

The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.     

Experimental study on the shear behavior of the interface between cushion materials and the concrete raft

Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

Long-Term Trends in the Critical Frequency of the E-layer

A search for trends k(foE) in the critical frequency of the ionospheric E layer at Juliusruh and Slough stations is performed by the method often used by the authors to analyze trends in the F2-layer parameters. It is found that k(foE) could differ in both magnitude and even sign within different time intervals. However, the k(foE) trends have been stably negative over the last two decades for both stations and all months of the year. The k(foE) values averaged over a year are ?0.012 and ?0.005 MHz per year for Juliusruh and Slough stations, respectively. The method used in the recent paper by La?tovi?ka et al. (2016) to determine foE trends is analyzed, and it is shown that the difference in linear approximation of the dependence of the observed foE values on F10.7 within different time intervals could be interpreted not as the presence of a different foE dependence on the F10.7 index within these intervals but as the presence within them of foE trends that change the slope of the linear approximation.     

Application of the dynamic calibration method to international monitoring system stations in Central Asia using natural seismicity data

The dynamic calibration method (DCM), using natural seismicity data and initially elaborated in [Kedrov, 2001; Kedrov et al., 2001; Kedrov and Kedrov, 2003], is applied to International Monitoring System (IMS) stations in Central Asia. The algorithm of the method is refined and a program is designed for calibrating diagnostic parameters (discriminants) that characterize a seismic source on the source-station traces. The DCM calibration of stations in relation to the region under study is performed by the choice of attenuation coefficients that adapt the diagnostic parameters to the conditions in a reference region. In this method, the stable Eurasia region is used as the latter. The calibration used numerical data samples taken from the archive of the International Data Centre (IDC) for the IMS stations MKAR, BVAR, EIL, ASF, and CMAR. In this paper, we used discriminants in the spectral and time domains that have the form

$D_i = X_i - a_m m_b - b_\Delta \log \Delta $

and are independent of the magnitude m _b and the epicentral distance Δ; these discriminants were elaborated in [Kedrov et al., 1990; Kedrov and Lyuke, 1999] on the basis of a method used for identification of events at regional distances in Eurasia. Prerequisites of the DCM are the assumptions that the coefficient a _m is regionindependent and the coefficient b _Δ depends only on the geotectonic characteristics of the medium and does not depend on the source type. Thus, b _Δ can be evaluated only from a sample of earthquakes in the region studied; it is used for adapting the discriminants D(X _i ) in the region studied to the reference region. The algorithm is constructed in such a way that corrected values of D(X i) are calculated from the found values of the calibration coefficients b _Δ, after which natural events in the region under study are selected by filtering. Empirical estimates of the filtering efficiency as a function of a station vary in a range of 95–100%. The DCM was independently tested using records obtained at the IRIS (Incorporated Research Institutions for Seismology) stations BRVK and MAKZ from explosions detonated in India on May 11, 1998, and Pakistan on May 28, 1998; these stations are similar in location and recording instrumentation characteristics to the IMS stations BVAR and MKAR. This test resulted in correct recognition of the source type and thereby directly confirmed the validity of the proposed calibration method of stations with the use of natural seismicity data. It is shown that the calibration coefficients b _Δ for traces similar in the conditions of signal propagation (e.g., the traces from Iran to the stations EIL and ASF) are comparable for nearly all diagnostic parameters. We arrive at the conclusion that the method of dynamic calibration of stations using natural seismicity data in a region where no explosions were detonated can be significant for a rapid and inexpensive calibration of IMS stations. The DCM can also be used for recognition of industrial chemical explosions that are sometimes erroneously classified in regional catalogs as earthquakes.

     

Characterization of the negative carbon isotope shift in segment C2, its global implications as a harbinger of OAE1a

Lower Cretaceous C-isotope records show intermittent negative/positive spikes, and consistent patterns of coeval chemostratigraphic curves thus document shifts that signal simultaneous responses of temporal changes in the global carbon reservoir. The standard pattern registered by the δ ¹³C_org and δ ¹³C_carb in Lower Aptian sediments includes distinct isotope segments C1 to C8 (Menegatti et al., 1998). In the El Pui section, Organyà Basin, Spain, C-isotope segment C2 is the longest interval preceding segments C3–C6 associated with oceanic anoxic event 1a (OAE 1a), and reveals a distinct negative shift of ~1.8‰ to ~2.23‰ defining the C-isotope pattern within that interval. Total inorganic carbon (TIC), total organic carbon (TOC), δ ¹³C_org, microfacies, n-alkanes show no difference before, during, or after the negative inflection. The biomarkers indicate that organic matter (OM) mainly originates from algal/microbial sources because short-chain length homologues (≤nC₁₉) dominate. nC₂₀ through nC₂₅ indicate some contribution from aquatic vegetation, but little from higher plants (>nC₂₅), as also suggested by the terrestrial/aquatic ratio of n-alkanes or (TAR) = [(nC₂₇+nC₂₉+nC₃₁)/(nC₁₅+nC₁₇+nC₁₉)] (averages 0.085). We suggest that conjoint pulses of contemporaneous LIPs (Ontong Java) and massive explosive volcanism in northeast Asia, the Songliao Basin (SB-V), best conform to plausible causes of the negative intra-C2 carbon isotopic excursion (CIE) at that time. Because of its apparent common occurrence the intra-C2 inflection could be a useful marker harbinger to the more pronounced CIE C3, the hallmark of OAE1a.     

The generalized kinetic equation as a model for the nonlinear transfer in third-generation wave models

An alternative model for the nonlinear interaction term S_nl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of S_nl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of S_nl, and to compare it to the existing models for S_nl in WAVEWATCH; the DIA and WRT. Some differences between the different models for S_nl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of S_nl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.     

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.     

Long-term trends in the ratio of the daytime and nighttime values of <Emphasis Type="Italic">foF</Emphasis>2

The consideration of the relation between the daytime and nighttime values of the critical frequency F2, foF2 of the ionospheric F2 layer, started in the previous publication of the authors, is continued. The main regularities in variations in the correlation coefficient R(foF2) characterizing this relation are confirmed using larger statistical material (more ionospheric stations and longer observational series). Long-term trends in the R(foF2) value are found: at all stations the negative value of R(foF2) increases with time after 1980.     

Variation of <Emphasis Type="Italic">Q</Emphasis> value before and after the 1999 Xiuyan,Liaoning Province, <Emphasis Type="Italic">M</Emphasis>=5.4 earthquake on the basis of analysis on attenuation dispersion of P waves

A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q _m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q _m started to appear from Nov. 9, 1999. The events giving the deduced high Q _m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q _m versus distance observation obviously deviates from the normal trend of Q _m linearly increasing with distance. (2) The average Q _m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q _m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q _m after the M=5.1 event is 204.     

Long-period irregular pulsations under the conditions of a quiet magnetosphere

Simultaneous observations of high-latitude long-period irregular pulsations at frequencies of 2.0–6.0 mHz (ipcl) and magnetic field disturbances in the solar wind plasma at low geomagnetic activity (Kp ~ 0) have been studied. The 1-s data on the magnetic field registration at Godhavn (GDH) high-latitude observatory and the 1-min data on the solar wind plasma and IMF parameters for 2011–2013 were used in an analysis. Ipcl (irregular pulsations continuous, long), which were observed against a background of the IMF Bz reorientation from northward to southward, have been analyzed. In this case other solar wind plasma and IMF parameters, such as velocity V, density n, solar wind dynamic pressure P = ρV² (ρ is plasma density), and strength magnitude B, were relatively stable. The effect of the IMF Bz variation rate on the ipcl spectral composition and intensity has been studied. It was established that the ipcl spectral density reaches its maximum (~10–20 min) after IMF Bz sign reversal in a predominant number of cases. It was detected that the ipcl average frequency (f) is linearly related to the IMF Bz variation rate (ΔBz/Δt). It was shown that the dependence of f on ΔBz/Δt is controlled by the α = arctan(By/Bx) angle value responsible for the MHD discontinuity type at the front boundary of magnetosphere. The results made it possible to assume that the formation of the observed ipcl spectrum, which is related to the IMF Bz reorientation, is caused by solar wind plasma turbulence, which promotes the development of current sheet instability and surface wave amplification at the magnetopause.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

Spatial inhomogeneity of variations in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer critical frequencies under the conditions of winter 1995

The daily samples of the hourly measurements of the foF2 critical frequency, obtained on January 5–21, 1995, at the midlatitude and high-latitude automated ionospheric stations (geographic latitude higher than 60°), are considered. The {fo} sets are transformed into the {δfoF2} sets of relative variations, for which asymmetry (A) and excess (E) are calculated. The selected stations are grouped into 20 pairs of automated ionospheric stations (AISs) located at distances of 200–10 000 km from one another. Sign estimates of the cross-correlation coefficients between the sets of 16 A and E values for different pairs of stations are used. Two types of structures of the statistical invariant spatial distribution are established: the structures with a scale of about 300 km, invariant with respect to latitude, and with a scale of about 6000 km (for only high latitudes).     

Geomagnetically Induced Currents and Space Weather: Pi3 Pulsations and Extreme Values of Time Derivatives of the Geomagnetic Field’s Horizontal Components

The problem of estimating the time derivatives of the horizontal components of the geomagnetic field and forecasting the probability of the occurrence of perturbations that exceed a given threshold level (the over-threshold perturbations) arises in the applications concerned with the geomagnetically induced currents (GICs). In this work, we consider the temporal and spatial structure of the Pi3 pulsations with quasi-periods of 10² to 10³ s during which the auroral and subauroral stations of the IMAGE network record over-threshold values in the derivatives of the meridional (along the longitudinal circle) B_X component and latitudinal (along the latitudinal circle) B_Y component. The extreme |dB_X/dt| values mainly develop against the background of the Pi3 pulsations with a complex frequency content, whereas the extreme |dB_Y/dt| values appear when the buildup (decay) phases of the bay-like disturbance associated with the evolution of a substorm coincide with the respective phases of the field of pulsations. The conditions under which the derivatives |dB_X/dt| and |dB_Y/dt| reach their over-threshold values are studied for subauroral latitudes by the technique of superposed epoch analysis. The extreme values of the derivatives most frequently occur during the main phase of moderate magnetic storms or beyond the storm—during high substorm activity under the conditions of a negative vertical component of the interplanetary magnetic field. The probability of the occurrence of over-threshold values increases at high amplitudes of the Pi3 pulsations and depends on their spectral content. The problem of analyzing and forecasting the over-threshold |dB_Y/dt| perturbations is complicated by the fact that the scale of the perturbations is small along the lines of latitude and large along the meridians. This can result in GIC excitation in the North–South oriented electric power lines by the geomagnetic perturbations localized within a narrow band in longitude which can be missed during the measurements.     

G: Fracture energy,friction and dissipation in earthquakes

Recent estimates of fracture energy G ^′ in earthquakes show a power-law dependence with slip u which can be summarized as G ^′ ∝ u ^a where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to G _f : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G ^′ and G _f should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s²). The resulting fracture energy G _f is similar to the seismological estimates, with G _f and G ^′ being comparable over most of the slip range. However, G _f appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G ^′ appears to increase further and surpasses G _f at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.     

Gravity Anomalies of Arbitrary 3D Polyhedral Bodies with Horizontal and Vertical Mass Contrasts

During the last 15 years, more attention has been paid to derive analytic formulae for the gravitational potential and field of polyhedral mass bodies with complicated polynomial density contrasts, because such formulae can be more suitable to approximate the true mass density variations of the earth (e.g., sedimentary basins and bedrock topography) than methods that use finer volume discretization and constant density contrasts. In this study, we derive analytic formulae for gravity anomalies of arbitrary polyhedral bodies with complicated polynomial density contrasts in 3D space. The anomalous mass density is allowed to vary in both horizontal and vertical directions in a polynomial form of \(\lambda =ax^m+by^n+cz^t\), where m, n, t are nonnegative integers and a, b, c are coefficients of mass density. First, the singular volume integrals of the gravity anomalies are transformed to regular or weakly singular surface integrals over each polygon of the polyhedral body. Then, in terms of the derived singularity-free analytic formulae of these surface integrals, singularity-free analytic formulae for gravity anomalies of arbitrary polyhedral bodies with horizontal and vertical polynomial density contrasts are obtained. For an arbitrary polyhedron, we successfully derived analytic formulae of the gravity potential and the gravity field in the case of \(m\le 1\), \(n\le 1\), \(t\le 1\), and an analytic formula of the gravity potential in the case of \(m=n=t=2\). For a rectangular prism, we derive an analytic formula of the gravity potential for \(m\le 3\), \(n\le 3\) and \(t\le 3\) and closed forms of the gravity field are presented for \(m\le 1\), \(n\le 1\) and \(t\le 4\). Besides generalizing previously published closed-form solutions for cases of constant and linear mass density contrasts to higher polynomial order, to our best knowledge, this is the first time that closed-form solutions are presented for the gravitational potential of a general polyhedral body with quadratic density contrast in all spatial directions and for the vertical gravitational field of a prismatic body with quartic density contrast along the vertical direction. To verify our new analytic formulae, a prismatic model with depth-dependent polynomial density contrast and a polyhedral body in the form of a triangular prism with constant contrast are tested. Excellent agreements between results of published analytic formulae and our results are achieved. Our new analytic formulae are useful tools to compute gravity anomalies of complicated mass density contrasts in the earth, when the observation sites are close to the surface or within mass bodies.     

Anomalous solar-diurnal variations in cosmic rays related to crossing of the IMF sector boundaries by the Earth and the problem of earthquakes

Based on a comparison of the cases of a decrease in the ratio of A _n/A _μ (where A _n and A _μ are the amplitudes of the diurnal variations of the neutron and hard cosmic ray components) to the instants of the Earth crossing the neutral IMF, it has been indicated that the process of such crossing is most effective for stimulating large destructive earthquakes with a magnitude of M ≥ 6. The 11-year period in the cyclicity of the occurrence probability of the above earthquakes has been revealed.     

Seismic performance evaluation of traditional timber <Emphasis Type="Italic">Hımış</Emphasis> frames: capacity spectrum method based assessment

Timber constructions have been widely suggested to be seismically resistant based on post-disaster reconnaissance studies. This observation has, however, remained to a large extent anecdotal due to the lack of experimental work supporting it, especially for certain timber architectural forms, including traditional timber frame “h?m??” structures. To fill this gap, the authors carried out an extensive full-scale testing scheme using frames of various geometrical configurations, tested under reverse-cyclic lateral loading with/without infill (brick and adobe) or cladding (ba?dadi and ?amdolma) (Aktas et al. in Earthq Spectra 30(4):1711–1732, 2014a, b). The tests concluded that h?m?? frames had high energy dissipation capabilities due mostly to nailed connections. Infill/cladding significantly helped improve stiffness and lateral load strength of the frames, and timber type did not seem to make a remarkable impact on the overall behaviour. The current paper, on the other hand, uses test data to calculate capacity/demand ratios based on capacity spectrum method and Eurocode 8 to elaborate more on the performance of “h?m??” structures under seismic loading. The obtained results are discussed to draw important conclusions with regards to how frame geometry and infill/cladding techniques affect the overall performance.