Temporal and spatial variations of seismicity scaling behavior in Southern México

R/S analysis is used in this work to investigate the fractal correlations in terms of the Hurst exponent for the 1998–2011 seismicity data in Southern Mexico. This region is the most seismically active area in Mexico, where epicenters for severe earthquakes (e.g., September 19, 1985, M_w = 8.1) causing extensive damage in highly populated areas have been located. By only considering the seismic events that meet the Gutenberg–Ritcher law completeness requirement (b = 0.97, M_GR = 3.6), we found time clustering for scales of about 100 and 135 events. In both cases, a cyclic behavior with dominant spectral components at about one cycle per year is revealed. It is argued that such a one-year cycle could be related to tidal effects in the Pacific coast. Interestingly, it is also found that high-magnitude events (M_w ≥ 6.0) are more likely to occur under increased interevent correlations with Hurst exponent values H > 0.65. This suggests that major earthquakes can occur when the tectonic stress accumulates in preferential directions. In contrast, the high-magnitude seismic risk is reduced when stresses are uniformly distributed in the tectonic shell. Such cointegration between correlations (i.e., Hurst exponent) and macroseismicity is confirmed for spatial variations of the Hurst exponent. In this way, we found that, using the Hurst exponent standpoint, the former presumed Michoacan and the Guerrero seismic gaps are the riskiest seismic zones. To test this empirical finding, two Southern Mexico local regions with large earthquakes were considered. These are the Atoyac de Alvarez, Guerrero (M_w = 6.3), and Union Hidalgo, Oaxaca (M_w = 6.6), events. In addition, we used the Loma Prieta, California, earthquake (October 17, 1989, M_w = 6.9) to show that the high-magnitude earthquakes in the San Andreas Fault region can also be linked to the increments of determinism (quantified in terms of the Hurst exponent) displayed by the stochastic dynamics of the interevent period time series. The results revealed that the analysis of seismic activity by means of R/S analysis could provide further insights in the advent of major earthquakes.     

Age estimation of the deposition of the cretaceous Upper Sandstones and a later remagnetization event,from Hirna area,SE Ethiopia

Outcrops of the Cretaceous Upper sandstone formation some 375 km to the East of Addis Ababa on the motor Highway to Harar was paleomagnetically investigated. About seventy core samples were collected at various stratigraphic levels from 250–300 meters thick sedimentary formation. After standard sample preparations in the laboratory the resulting specimens were subjected to routine paleomagnetic demagnetization protocol. In the first steps of demagnetizations process the recent and viscous magnetizations were removed by heating until a temperature of level of 300 °C. Further demagnetization of the samples resulted in the isolation of the final magnetization with stable line segments that is directed towards the origin, which is interpreted as Characteristic Remanent Magnetization (ChRM). Rock – magnetic experiments have identified goethite (αFeOOH), hematite (αFe₂O₃), detritial hematite, and magnetite as the magnetic mineral phases carrying the remanence. The ChRM identified resulted in an average value of (D_s = 0.5°, I_s = ?0.7°, α₉₅ = 4.3°, N = 34) for the red sandstones while an average value of (D_s = 335.8°, I_s = ?31.8°, α₉₅ = 4.7°, N = 14) for the limestone intercalations. The former ChRM in the red sandstone is determined to be secondary while the latter ChRM is known to be primary. Comparison of these directional results and their pole equivalents with the African plate Apparent Polar Wander Path curve established by Besse and Courtillot (2003) give ages of between 115–130 Million years for the limestone intercalation and ages of 30 million years for red sandstone unit. These are interpreted respectively as estimates of the age of deposition and a later remagnetization respectively.     

Exact solutions to radially symmetric two-phase flow for an arbitrary diffusivity

Over the past decade there have been a variety of exact solutions developed for one-dimensional two-phase flow, however when higher dimensions are considered there is a distinct scarcity of solutions. In this paper we consider the problem of radially symmetric two-phase flow, into an infinite medium of uniform initial saturation, subject to a constant flux V from a line source at the origin. We show that in the absence of gravity and when the two-phase diffusivity D is related to the fraction flow function f by βD = V df/dθ, where θ is the water content and β is a constant of proportionality, a new class of exact solutions can be found. In particular, when β = 2, we show that the solution is given by a simple quadrature for arbitrary D, and is fully integrable for specific functional forms of D. It has been shown by Weeks et al. [Weeks SW, Sander GC, Parlange J-Y. n-Dimensional first integral and similarity solutions for two-phase flow. ANZIAM J 2003;44:365–80] that when D obeys the above relation, a saturated zone does not grow around the line of injection, consequently we find for β = 1, the flow equation maps to one-dimensional single-phase flow under a saturated boundary condition. Consequently solutions developed for one-dimensional single-phase flow (exact or approximate) apply to radially symmetric two-phase flow. Solutions for β = 1 or 2 can be derived for either a wetting fluid displacing a non-wetting fluid, or a non-wetting fluid displacing a wetting fluid, however for arbitrary β numerical methods are required.     

Whistlers observed at low-latitude ground-based VLF facility in Fiji

The propagation features of nighttime whistlers to low-latitude station, Suva (−18.2°, 178.3°, geomag. lat. −22.1°, geomag. long. 253.5°, L=1.15), Fiji, from preliminary observations made during the period from September 2003–2005, are reported. The observations of ELF–VLF signals commenced in September 2003 using the VLF set-up of World Wide Lightning Location Network at our station. The whistlers were observed during the severe magnetic storm of 20–22 November 2003 and moderate magnetic storm of 17–19 July 2005. A whistler with dispersion D=12.7 s^1/2 occurred on 22 November at 00:11 h LT. On 20 July at 01:00 h LT, a short whistler with dispersion D=20.9 s^1/2 and two whistler events having two-component whistlers with D=15.8, 16.7 s^1/2 and 16.7, 17.3 s^1/2 were observed. Non-ducted pro-longitudinal mode of the whistler propagation supported by negative latitudinal electron density gradients in the ionosphere that are enhanced by magnetic storms, seems most likely mode of propagation for the whistlers with dispersion of 12.7–17.3 s^1/2 to this low-latitude station.     

Patterns of significant seismic quiescence on the Mexican Pacific coast

Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, M_w = 7.8), Petatlán (March 1979, M_w = 7.6), Michoacán (September 1985, M_w = 8.0, and M_w = 7.6) and Colima (October 1995, M_w = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.     

The equivalent dose of different grain size quartz fractions from lakeshore sediments in the arid region of north China

Equivalent dose (D_e) values were measured by using medium aliquots of different grain size quartz fractions of five lakeshore sediments from the arid region of north China. There are two different relationships between D_e values and grain sizes of these five samples. The first relationship is that the D_e values obtained from various grain sizes are in agreement within 1 delta errors. The second relationship is that D_e values are similar to each other for fractions between 125 and 300 μm, while the D_e value of the 63–90 μm fraction is 40～55% smaller than others. For example, the D_e values obtained for sample #3 are 20.15 ± 1.19 Gy, 19.80 ± 0.83 Gy and 20.93 ± 1.06 Gy for fractions of 90–125, 125–150 and 250–300 μm respectively, but are 10.79 ± 0.84 Gy for the 63–90 μm fraction. The second relationship can't be interpreted by previous studies of both dosimetry and heterogeneous bleaching. It is deduced for sample #2, #3 and #6 that fine particles (<90 μm) intruded after the dominant sedimentation. Comparison of OSL ages from different grain size fractions of sample #2 with a radiocarbon age from the same lithologic layer supports that fractions coarser than 125 μm yield more reliable burial ages, while the fraction finer than 90 μm yields underestimated ages for some lakeshore sediments from this arid region.     

Correlations between electrical and elastic properties of solid–liquid composites with interfacial energy-controlled equilibrium microstructures

Electrical conductivity and seismic velocity are studied for plausible pore geometries in the Earth's interior for reliable quantitative analysis of experimental data such as seismic tomography and magnetotelluric explorations. Electrical conductivity of a two-phase system with equilibrium, interfacial energy-controlled phase geometry is calculated for the dihedral angles θ = 40°–100° that are typical for rock–aqueous fluid and θ = 20°–60° for rock–melt systems of lower crust and upper mantle for the case of tetrakaidecahedral grains. Electrical conductivity vs. seismic velocity correlations are acquired by combining of the simulated electrical conductivities with the seismic velocity calculated with the help of equilibrium geometry model Takei [Takei, Y., Effect of pore geometry on V_P/V_S: From equilibrium geometry to crack. J. Geophys. Res. 107 (2002): 10.1029/2001JB000522.] for the same pore geometries. The results show that electrical conductivity gradually decreases reaching zero when seismic velocities reach seismic velocities of intact rock for rock–melt systems, while for rock–aqueous fluid systems with θ ≥ 60° conductivity drops to zero at velocities up to 10% smaller. This can explain the seeming discrepancy of the low seismic velocity region, attributed to the high fluid fraction, and the low electrical conductivity of the same region, which is sometimes faced at collocated electromagnetic and seismic experiments.     

Systematic deflection and offset of the Yangtze River drainage system along the strike-slip Ganzi-Yushu-Xianshuihe Fault Zone,Tibetan Plateau

The Ganzi-Yushu-Xianshuihe Fault Zone (GYXFZ) is a typical active strike-slip fault that has triggered many large historic earthquakes, including the 2010 M_w 6.9 Yushu earthquake in the central Tibetan Plateau. This fault zone extends for ca. 800 km from the central Tibetan Plateau to its southeastern margin and varies in trend from WNW-ESE in the northwestern segment of the fault zone to NNW-SSE in the southeastern segment, having the geometry of an arc projecting northeastwards. In this study, we present evidence for the systematical sinistral deflection and/or offset of the Yangtze River and its branch stream channels and valleys along the GYXFZ. Topographic analysis of three-dimensional (3D) perspective images constructed using digital elevation model (DEM) data, 0.5 m-resolution WorldView and GeoEye images, and 15 m-resolution Landsat-Enhanced Thematic Mapper (ETM+) images, together with analysis of geological structures, reveals the following: (i) the main river channels and valleys of the Yangtze River drainage system show systematic sinistral deflections and/or offsets along the GYXFZ; (ii) various amounts of sinistral offset have accumulated on the tributary stream channels, valleys, and gullies of the Yangtze River along the fault, with a linear relation, D = aL, between the upstream length L from the deflected point and the offset amount D with a certain coefficient a; (iii) the maximum amount of sinistral offset is up to ca. 60 km, which was accumulated in the past 13–5 Ma; and (iv) the long-term average strike-slip rate is ca. 4.6–12 mm/year. Geological and geomorphic evidence, combined with geophysical data, demonstrates that the GYXFZ is currently active as one of the major seismogenic faults in the Tibetan Plateau, dominated by left-lateral strike-slip motion. Our findings supply important evidence for the tectonic evolution of strike-slip faults in the Tibetan Plateau since the Eurasia-India continental collision.     

Comparison of absorbing boundary conditions: A correlation of reflection coefficients and determination of the optimum grid interval for migration

Modeling and migration couple is one of the most important steps in seismic data processing and interpretation. Absorbing boundary conditions used in the modeling were studied with the wave-equation by different authors. In this study, reflection coefficient analyses of recent solutions are compared to each other for the different incident angles of seismic waves to the modeling boundaries. According to the reflection coefficients correlation, the A3 condition is the most appropriate solution which greatly reduces artificial reflections from the boundaries. However, multi-transmitting Formula is better for relatively high angles between 32–58° with the usage of a special parameter. On the contrary, this formula is not an appropriate condition for angles lower than 32°, although it allows setting the boundary at any preferred angle. Considering that most of the boundaries are set in low angles, A3 solution is still most preferential condition. In this study, it is also aimed to find out the optimum grid intervals for minimizing the ill-posedness arose from the combination of the 45° finite difference migration equation and the B₂ absorbing boundary condition for migration. Appropriate values are determined as ωΔx = 0.2 and ωΔz = 0.4 or neighbouring coarser values. It is also concluded that finer mesh spacing can increase the ill-posedness, in contrast to existence of some fine grid size values providing well-posedness. In addition, ill-posedness is obviously standard after ωΔx = 0.6 for all values of ωΔz.     

Spatial distribution and ecology of the Recent Ostracoda from Tangra Yumco and adjacent waters on the southern Tibetan Plateau: A key to palaeoenvironmental reconstruction

We elucidate the ecology of Recent Ostracoda from a deep brackish lake, Tangra Yumco (30°45′—31°22′N and 86°23′—86°49′E, 4595 m a.s.l.) and adjacent waters on the southern Tibetan Plateau. Ostracod associations (living and empty valves) in sixty-six sediment samples collected from diverse aquatic habitats (lakes, estuary-like water and lagoon-like water waters, rivers, ponds and springs) were quantitatively assessed.Eleven Recent Ostracoda were found (nine living and two as empty valves only). Cluster analysis established two significant (p < 0.05) habitat specific associations; (i) Leucocytherella sinensis, Limnocythere inopinata, Leucocythere? dorsotuberosa, Fabaeformiscandona gyirongensis and Candona xizangensis are lacustrine fauna. (ii) Tonnacypris gyirongensis, Candona candida, Ilyocypris sp., Heterocypris incongruens and Heterocypris salina are temporary water species.Ostracod distribution and abundance are significantly (p < 0.05) correlated to physico-chemical variables. The first two axes of a canonical correspondence analysis (CCA) explain 30.9% of the variation in the species abundance data. Conductivity and habitat types are the most influential ecological factors explaining the presence and abundance of ostracods. Spearman correlation analysis reveals that: (i) Two species, L.? dorsotuberosa (r = 0.25) and L. inopinata (r = 0.36) have a significant positive correlation with conductivity while one species, T. gyirongensis (r = −0.68) displays a significant negative correlation with conductivity. Limnocythere inopinata correlates significantly positive (r = 0.37) with alkalinity. Fabaeformiscandona gyirongensis correlates significantly positive (r = 0.28) with water depth.Key indicator living assemblages are: (i) L. sinensis dominates Ca-depleted brackish waters although ubiquitously distributed; (ii) L.? dorsotuberosa dwells in fresh to brackish waters; (iii) L. inopinata predominates in mesohaline to polyhaline waters; (iv) F. gyirongensis inhabits exclusively brackish-lacustrine deeper waters; (v) C. candida populates freshwaters; (vi) T. gyirongensis and Ilyocypris sp. are restricted to shallow temporary waters; (vii) H. incongruens occurs in ponds.Water depth indicators are F. gyirongensis and L.? dorsotuberosa, useful in ostracod assemblages for palaeo-water depth reconstruction.Our results expand the knowledge of the ecological significance of Recent Tibetan Ostracoda ecology. This is a new insight on habitat chacteristics of both living assemblages and sub-Recent associations of ostracods in mountain aquatic ecosystems. The new modern ostracod dataset can be used for the quantitative reconstruction of past environmental variables (e.g., conductivity) and types of water environment. The key indicator ostracods are relevant in palaeolimnological and climate research on the Tibetan Plateau.     

Coda-Q and its lapse time dependence analysis in the interaction zone of the Dinarides,the Alps and the Pannonian basin

The single backscattering model was used to estimate total attenuation of coda waves (Q_c) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q₀ and n, parameters of the frequency dependent coda-Q using the relation Q_c = Q₀fⁿ. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Q_c were distributed into classes according to their lapse time, t_L. For t_L = 20–50 s we estimated Q₀ = 45–184 and n = 0.49–0.94, and for t_L = 60–100 s we obtained Q₀ = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q₀, and vice versa. The rates of change of both Q₀ and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.     

A similarity approach to determine response times to steady-state saturation in landscapes

Relationships describing response times for landscape saturation and subsurface flow for idealised hillslopes after a change in water balance are derived in terms of similarity parameters given by their topographic, soil and climatic attributes. The study was carried out under a set of specific assumptions.The work quantitatively describes effects of each of these attributes on travel time. For example, divergent hillslopes can exhibit travel times that are double those of convergent hillslopes, and concave slopes tend to have lower travel times than planar or convex slopes. It is shown that the hillslope travel time T is dependent on hillslope length L, hydraulic conductivity K and slope S such that T = f(L/KS, B, CR, Δq/smd) and if saturation does not occur in the lower reaches of a hillslope, this can be simplified to T = f(L/KS, B). Here, Δq is change in net external flux and smd is soil moisture deficit; CR and B are the convergence ratio and profile factor respectively. The relative importance of these attributes on response time is discussed.It is shown that subsurface flow hydrographs of simple hillslopes of different scale, derived numerically, can be collapsed into a single curve by normalising them by means of a form of response time, the time constant.Results from the analytical derivation were compared with numerical analysis with good agreement. The theory was applied to a natural catchment using observed data from a bushfire event. Comparison of analytically calculated response time with observed response time for the event showed a large discrepancy. The reason is discussed.     

Influence of nitrate and dissolved organic carbon loading on the interaction of Microcystis aeruginosa and heterotrophic bacteria from hyper-eutrophic lake (Taihu Lake,China)

To investigate the interaction between cyanobacteria and heterotrophic bacteria under gradients of glucose and nitrate, a cyanobacterial strain of Microcystis aeruginosa was grown in microcosms with and without a freshwater bacterial mixture, which was collected from Lake Taihu. Concentrations of glucose (1350, 975, 600, 300, 150, and 37.5 μmol C L⁻¹) and nitrate (150, 300, and 9000 μmol N L⁻¹) were used in a range of combinations giving 9 different treatments of glucose:nitrate. In the microcosm without the bacterial mixture, M. aeruginosa abundance gradually increased with days in all treatments. However, M. aeruginosa had much lower density in some treatments with the bacterial mixture. The difference in M. aeruginosa growth could be explained by competition with bacteria in the cultures in which these were added. The abundance of M. aeruginosa and bacteria when grown together was nearly equal and the number of the bacterial species was highest in the treatment with 300 μmol C L⁻¹ and 150 μmol N L⁻¹. Our results suggest that at this glucose:nitrate ratio M. aeruginosa and the bacterial mixture maintain a balance, and bacteria maintain diversity. In conclusion, we propose that dissolved organic carbon and nitrate availability fundamentally affects the structure as well as stoichiometry of pelagic associations.     

Empirical relationship between leaf wax n-alkane δD and altitude in the Wuyi,Shennongjia and Tianshan Mountains,China: Implications for paleoaltimetry

Estimating past elevation not only provides evidence for vertical movements of the Earth's lithosphere, but also increases our understanding of interactions between tectonics, relief and climate in geological history. Development of biomarker hydrogen isotope-based paleoaltimetry techniques that can be applied to a wide range of sample types is therefore of continuing importance. Here we present leaf wax-derived n-alkane δD (δD_wax) values along three soil altitudinal transects, at different latitudes, in the Wuyi, Shennongjia and Tianshan Mountains in China, to investigate δD_wax gradients and the apparent fractionation between leaf wax and precipitation (ε_wax-p).We find that soil δD_wax track altitudinal variations of precipitation δD along the three transects that span variable environment conditions and vertical vegetation spectra. An empirical δD_wax-altitude relation is therefore established in which the average δD_wax lapse rate of ? 2.27 ± 0.38‰/100 m is suitable for predicting relative paleoelevation change (relative uplift). The application of this empirical gradient is restricted to phases in the mountain uplift stage when the atmospheric circulation had not distinctly changed and to when the climate was not arid. An empirical δD_wax–latitude–altitude formula is also calculated: δD_wax = 3.483LAT ? 0.0227ALT ? 261.5, which gives the preliminary spatial distribution pattern of δD_wax in modern China.Mean value of ε_wax-p in the extreme humid Wuyi Mountains is quite negative (? 154‰), compared to the humid Shennongjia (? 129‰) and the arid (but with abundant summer precipitation) Tianshan Mountains (? 130‰), which suggests aridity or water availability in the growing season is the primary factor controlling soil/sediment ε_wax-p. Along the Tianshan transects, values of ε_wax-p are speculated to be constant with altitude; while along the Wuyi and Shennongjia transects, ε_wax-p are also constant at the low-mid altitudes, but become slightly more negative at high altitudes which could be attributed to overestimates of precipitation δD or the vegetation shift to grass/conifer.Additionally, a reversal of altitude effect in the vertical variation of δD_wax was found in the alpine zone of the Tianshan Mountains, which might be caused by atmospheric circulation change with altitude. This implies that the paleo-circulation pattern and its changes should also be evaluated when stable isotope-based paleoaltimetry is applied.     

Dynamic analysis of offshore wind turbine in clay considering soil–monopile–tower interaction

A comprehensive study is performed on the dynamic behavior of offshore wind turbine (OWT) structure supported on monopile foundation in clay. The system is modeled using a beam on nonlinear Winkler foundation model. Soil resistance is modeled using American Petroleum Institute based cyclic p–y and t–z curves. Dynamic analysis is carried out in time domain using finite element method considering wind and wave loads. Several parameters, such as soil–monopile–tower interaction, rotor and wave frequencies, wind and wave loading parameters, and length, diameter and thickness of monopile affecting the dynamic characteristics of OWT system and the responses are investigated. The study shows soil–monopile–tower interaction increases response of tower and monopile. Soil nonlinearity increases the system response at higher wind speed. Rotor frequency is found to have dominant role than blade passing frequency and wave frequency. Magnitude of wave load is important for design rather than resonance from wave frequency.     

Seismic microzonation studies for the city of Ragusa (Italy)

The seismic history of the city of Ragusa (Italy), the geotechnical characterisation of the subsoil and the site response analysis should be correctly evaluated for the definition of the Seismic Geotechnical Hazard of the city of Ragusa, through geo-settled seismic microzoning maps. Basing on the seismic history of the city of Ragusa, the following earthquake scenarios have been considered: the “Val di Noto” earthquake of January 11, 1693 (with intensity X–XI on MCS scale, magnitude M_W=7.41 and epicentral distance of about 53 km); the “Etna” earthquake of February 20, 1818 (with intensity IX on MCS scale, magnitude M_W=6.23 and epicentral distance of about 64 km); the Vizzini earthquake of April 13, 1895 (with intensity I=VII–VIII on MCS scale, magnitude M_W=5.86 and epicentral distance of about 26 km); the “Modica” earthquake of January 23, 1980 (with intensity I=V–VI on MCS scale, magnitude M_W=4.58 and epicentral distance of about 10 km); the “Sicilian” earthquake of December 13, 1990 (with intensity I=VII on MCS scale, magnitude M_W=5.64 and epicentral distance of about 50 km). Geotechnical characterisation has been performed by in situ and laboratory tests, with the definition of shear wave velocity profiles in the upper 30 m of soil. Soil response analyses have been evaluated for about 120 borings location by some non-linear 1-D models. Finally the seismic microzonation of the city of Ragusa has been obtained in terms of maps with different peak ground acceleration at the surface; shaking maps for the central area of the city of Ragusa were generated via GIS for the earthquake scenarios.     

On the comparison of tidal gravity parameters with tidal models in central Europe

Two accurately calibrated superconducting gravimeters (SGs) provide high quality tidal gravity records in three central European stations: C025 in Vienna and at Conrad observatory (A) and OSG050 in Pecný (CZ). To correct the tidal gravity factors from ocean loading effects we compared the load vectors from different ocean tides models (OTMs) computed with different software: OLFG/OLMP by the Free Ocean Tides Loading Provider (FLP), ICET and NLOADF. Even with the recent OTMs the mass conservation is critical but the methods used to correct the mass imbalance agree within 0.1 nm/s². Although the different software agrees, FLP probably provides more accurate computations as this software has been optimised. For our final computation we used the mean load vector computed by FLP for 8 OTMs (CSR4, NAO99, GOT00, TPX07, FES04, DTU10, EOT11a and HAMTIDE). The corrected tidal factors of the 3 stations agree better than 0.04% in amplitude and 0.02° in phase. Considering the weighted mean of the three stations we get for O1 δ_c = 1.1535 ± 0.0001, for K1 δ_c = 1.1352 ± 0.0003 and for M2 δ_c = 1.1621 ± 0.0003. These values confirm previous ones obtained with 16 European stations. The theoretical body tides model DDW99/NH provides the best agreement for M2 (1.1620) and MATH01/NH for O1 (1.1540) and K1 (1.1350). The largest discrepancy is for O1 (0.05%). The corrected phase α_c does not differ significantly from zero except for K1 and S2. The calibrations of the two SG's are consistent within 0.025% and agree with Strasbourg results within 0.05%.     

Cave air ventilation and CO2 outgassing by radon-222 modeling: How fast do caves breathe?

In general, the rate and timing of calcite precipitation is in part affected by variations in cave air CO₂ concentrations. Knowledge of cave ventilation processes is required to quantify the effect variations in CO₂ concentrations have on speleothem deposition rates and thus paleoclimate records. In this study we use radon-222 (²²²Rn) as a proxy of ventilation to estimate CO₂ outgassing from the cave to the atmosphere, which can be used to infer relative speleothem deposition rates. Hollow Ridge Cave, a wild cave preserve in Marianna, Florida, is instrumented inside and out with multiple micro-meteorological sensor stations that record continuous physical and air chemistry time-series data. Our time series datasets indicate diurnal and seasonal variations in cave air ²²²Rn and CO₂ concentrations, punctuated by events that provide clues to ventilation and drip water degassing mechanisms. Average cave air ²²²Rn and CO₂ concentrations vary seasonally between winter (²²²Rn = 50 dpm L^? 1, where 1 dpm L^? 1 = 60 Bq m^? 3; CO₂ = 360 ppmv) and summer (²²²Rn = 1400 dpm L^? 1; CO₂ = 3900 ppmv). Large amplitude diurnal variations are observed during late summer and autumn (²²²Rn = 6 to 581 dpm L^? 1; CO₂ = 360 to 2500 ppmv).We employ a simple first-order ²²²Rn mass balance model to estimate cave air exchange rates with the outside atmosphere. Ventilation occurs via density driven flow and by winds across the entrances which create a ‘venturi’ effect. The most rapid ventilation occurs 25 m inside the cave near the entrance: 45 h^? 1 (1.33 min turnover time). Farther inside (175 m) exchange is slower and maximum ventilation rates are 3 h^? 1 (22 min turnover time). We estimate net CO₂ flux from the epikarst to the cave atmosphere using a CO₂ mass balance model tuned with the ²²²Rn model. Net CO₂ flux from the epikarst is highest in summer (72 mmol m^? 2 day^? 1) and lowest in late autumn and winter (12 mmol m^? 2 day^? 1). Modeled ventilation and net CO₂ fluxes are used to estimate net CO₂ outgassing from the cave to the atmosphere. Average net CO₂ outgassing is positive (net loss from the cave) and is highest in late summer and early autumn (about 4 mol h^? 1) and lowest in winter (about 0.5 mol h^? 1). Modeling of ventilation, net CO₂ flux from the epikarst, and CO₂ outgassing to the atmosphere from cave monitoring time-series can help better constrain paleoclimatic interpretations of speleothem geochemical records.     

Ground-based Pc5 ULF wave power: Solar wind speed and MLT dependence

Using over 20 years of ground-based magnetometer data from the CANOPUS/CARISMA magnetometer array, we present a statistical characterisation of Pc5 ultra-low frequency (ULF) power in the 2–10 mHz band as a function of magnetic local time (MLT), L-shell, and solar wind speed. We examine the power across L-shells between 4.2 and 7.9, using data from the PINA, ISLL, GILL and FCHU stations, and demonstrate that there is a significant MLT dependence in both the H- and D-component median 2–10 mHz power during both fast (>500 km/s) and slow (<500 km/s) solar wind speeds. The H-component power consistently dominates over D-component power at all MLTs and during both fast and slow solar wind. At the higher-L stations (L>5.4), there are strong MLT power peaks in the morning and midnight local time sectors; the morning sector dominating midnight during fast solar wind events. At lower L-shells, there is no evidence of the midnight peak and the 2–10 mHz power is more symmetric with respect to MLT except during the fastest solar wind speeds. There is little evidence in the ground-based power of a localised MLT peak in ULF power at dusk, except at the lowest L-shell station, predominantly in the H-component. The median 2–10 mHz power increases with an approximate power law dependence on solar wind speed, at all local times across the L-shell domain studied in both components. The H-component power peaks at the latitude of the GILL station, with significantly lower power at both higher and lower L-shells. Conversely, the D-component power increases monotonically. We believe that this is evidence for 2–10 mHz power accumulating at auroral latitudes in field line resonances. Finally, we discuss how such ULF wave power characterisation might be used to derive empirical radiation belt radial diffusion coefficients based on, and driven by, the solar wind speed dependence of ULF wave power.     

Earthquake scaling,fault segmentation,and structural maturity

Slip and length measurements on earthquakes suggest large stress drop variability. We analyze an extended set of slip-length measurements for large earthquakes (M ≥ 6) to seek for the possible origin(s) of this apparent variability. We propose that such variability arises from earthquakes breaking a variable number of major fault segments. That number depends on the strength of the inter-segment zones, which itself depends on the structural maturity of the faults. We propose new D_max–L parameterizations based on that idea of multiple segment-ruptures. In such parameterizations, each broken segment roughly scales as a crack, while the total multi-segment rupture does not. Stress drop on individual segments is roughly constant, only varying between 3.5 to 9 MPa. The slight variation that is still observed depends on fault structural maturity; more mature faults have lower stress drops than immature ones. The new D_max–L functions that we propose reduce uncertainties with respect to available relationships. They thus provide a more solid basis to estimate seismic hazard by integrating fault properties revealed by geological studies.