Incorporating space,time, and magnitude measures in a network characterization of earthquake events

We investigate the structural properties of a spatio-temporal network of earthquake events that incorporates magnitude information between the connected events. The network creates temporally directed links from an origin event towards a later event if it breaks the record closest distance from the origin among all the events in the catalog so far. Additionally, the links are conditionally classified based on the magnitude difference between connected events: “up” (“down”) connections point from a weaker (stronger) to a stronger (weaker) event. Using earthquake records from the Philippines from 1973 to 2012 and southern California from 1982 to 2012, we observe that the out-degree distributions show slight deviations from the corresponding Poisson distribution of the same mean. The space and time separations of connected earthquakes both show power-law regimes, suggesting spatio-temporal (self-)organization. More importantly, the conditional distributions of “up” and “down” connections in space, time, and network structure point to a higher likelihood of a stronger event triggering a nearby weaker event for the first few connections, as in the case of aftershocks. The results are captured by a sandpile-based model where a small but finite probability of preferentially targeting the most susceptible grid site is introduced. Our analysis, coupled with the discrete model analog, provides a quantitative picture of the spatio-temporal and magnitude organization of seismicity beyond just the successive events. The technique may be extended to further characterize similar long-period earthquake records to yield a more complete picture of the underlying processes involved in seismicity.     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

Dispersion features of transmitted channel waves and inversion of coal seam thickness

In-seam seismic survey currently is a hot geophysical exploration technology used for the prediction of coal seam thickness in China. Many studies have investigated the relationship between the group velocity of channel wave at certain frequency and the actual thickness of exposed coal beds. But these results are based on statistics and not universally applicable to predict the thickness of coal seams. In this study, we first theoretically analyzed the relationship between the depth and energy distribution of multi-order Love-type channel waves and found that when the channel wave wavelength is smaller than the thickness of the coal seam, the energy is more concentrated, while when the wavelength is greater than the thickness, the energy reduces linearly. We then utilized the numerical simulation technology to obtain the signal of the simulated Love-type channel wave, analyzed its frequency dispersion, and calculated the theoretical dispersion curves. The results showed that the dispersion characteristics of the channel wave are closely related to the thickness of coal seam, and the shear wave velocity of the coal seam and its surrounding rocks. In addition, we for the first time realized the joint inversion of multi-order Love-type channel waves based on the genetic algorithm and inversely calculated the velocities of shear wave in both coal seam and its surrounding rocks and the thickness of the coal seam. In addition, we found the group velocity dispersion curve of the single-channel transmitted channel wave using the time–frequency analysis and obtained the phase velocity dispersion curve based on the mathematical relationship between the group and phase velocities. Moreover, we employed the phase velocity dispersion curve to complete the inversion of the above method and obtain the predicted coal seam thickness. By comparing the geological sketch of the coal mining face, we found that the predicted coal seam thickness is in good agreement with the actual thickness. Overall, adopting the channel wave inversion method that creatively uses the complete dispersion curve can obtain the shear wave velocities of the coal and its surrounding rocks, and analyzing the depth of the abruptly changed shear wave velocity can accurately obtain the thickness of the coal seam. Therefore, our study proved that this inversion method is feasible to be used in both simulation experiments and actual detection.     

Preliminary comparison of trace metals in coastal aerosols between Qingdao and Liverpool

Aerosol samples were collected at two coastal suburban stations, Qingdao （China） in 1995-1996 and Liverpool （U.K.） in 1995, respectively. The samples were analyzed to determine the concentrations of trace metals （Cr, Zn, Cu, Co, Ni, Pb, V, and Cd） as well as A1, Fe and Mn. Data were examined to understand the difference of trace metals in aerosols between coastal zones downwind the developing area （near the Yellow Sea） and developed region （near the Irish Sea）. The results show that most elements at Qingdao have levels 4-5 times higher than those at Liverpool, particularly for the crust-dominated elements （e.g. Al, Fe and Mn）. Moreover, the aerosol composition at Qingdao is higher in spring than in summer, underlying the influence of westerlies and local emissions in combination, whereas seasonal change of aerosol composition is not significant at Liverpool. The enrichment factors for the crustal source elements （EFcrust） at Liverpool are much higher than those at Qingdao. The contributions from the pollutant source （Rp） for some trace metals like Cu, Pb, Zn and Cd are 〉90% at Qingdao and Liverpool, suggesting overwhelming anthropogenic contributions to these metals. The contributions from crustal source （Re） for trace metals tend to increase with higher aerosol levels and Al concentration at Qingdao, indicating a good correlation between the crust-dominated component and the air mass. At Liverpool, the Rc values for trace metals are positively correlated with Al concentrations instead of with aerosol mass, suggesting that Al in aerosols represents the crustal component even though the aerosols come from different sources.     

Joint probability density function of the stochastic responses of nonlinear structures

The joint probability density fimction （PDF） of different structural responses is a very important topic in the stochastic response analysis of nonlinear structures. In this paper, the probability density evolution method, which is successfully developed to capture the instantaneous PDF of an arbitrary single response of interest, is extended to evaluate the joint PDF of any two responses. A two-dimensional partial differential equation in terms of the joint PDF is established. The strategy of selecting representative points via the number theoretical method and sieved by a hyper-ellipsoid is outlined. A two-dimensional difference scheme is developed. The free vibration of an SDOF system is examined to verify the proposed method, and a flame structure exhibiting hysteresis subjected to stochastic ground motion is investigated. It is pointed out that the correlation of different responses results from the fact that randomness of different responses comes from the same set of basic random parameters involved. In other words, the essence of the probabilistic correlation is a physical correlation.     

∼3,850 Ma tonalites in the Nuuk region,Greenland: geochemistry and their reworking within an Eoarchaean gneiss complex

The Eoarchaean (>3,600 Ma) Itsaq Gneiss Complex of southern West Greenland is dominated by polyphase orthogneisses with a complex Archaean tectonothermal history. Some of the orthogneisses have c. 3,850 Ma zircons, and they vary from rare single phase metatonalites to more common complexly banded migmatites. This is due to heterogeneous strain, in situ anatexis and granitic veining superimposed during younger tectonothermal events. In the single-phase tonalites with c. 3,850 Ma zircon, oscillatory-zoned prismatic zircon is all 3,850 Ma old, but shows patchy ancient loss of radiogenic Pb. SHRIMP spot analyses and laser ablation ICP-MS depth profiling show that thin (usually < 10 μm) younger (3,660–3,590 Ma and Neoarchaean) shells of lower Th/U metamorphic zircon are present on these 3,850 Ma zircons. Several samples with this simple zircon population occur on islands near Akilia. In contrast, migmatites usually contain more complex zircon populations, with often more than one generation of igneous zircon present. Additional zircon dating of banded gneisses across the Complex shows that samples with c. 3,850 Ma igneous zircon are not just a phenomenon restricted to Akilia and adjacent islands. For example, migmatites from Itilleq (c. 65 km from Akilia) contain variable amounts of oscillatory-zoned 3,850 Ma and 3,650 Ma zircon, interpreted, respectively, as the rock age and the time of crustal melting under Eoarchaean metamorphism. With only 110–140 ppm Zr in the tonalites and likely magmatic temperatures of >850°C, zircon solubility–melt composition relationships show that they were only one-third saturated in zircon. Any zircon entrained in the precursor magmas would thus have been highly soluble. Combined with the cathodoluminesence imaging, this demonstrates that the c. 3,850 Ma oscillatory zoned zircon crystallised out of the melt and hence gives a magmatic age. Thus the rare well-preserved tonalites and palaeosome in migmatites testify that c. 3,850 Ma quartzo–feldspathic rocks are a widespread (but probably minor) component in the Itsaq Gneiss Complex. C. 3,850 Ma zircon with negative Eu anomalies (showing growth in felsic systems) also occurs as detrital grains in rare c. 3,800 Ma metaquartzites and as inherited grains in some 3,660 Ma granites (sensu stricto). These demonstrate that still more c. 3,850 Ma rocks were present, but were recycled into Eoarchaean sediments and crustally derived granites. The major and trace element characteristics (e.g. LREE enrichment, HREE depletion, low MgO) of the best-preserved c. 3,850 Ma rocks are typical of Archaean TTG suites, and thus argue for crust formation processes involving important contributions from melting of hydrated mafic crust to the earliest Archaean. Five c. 3,850 Ma tonalites were selected as the best preserved on the basis of field criteria and zircon petrology. Four of these samples have overlapping initial ɛ_Nd (3,850 Ma) values from +2.9 to +3.6± 0.5, with the fourth having a slightly lower value of +0.6. These data provide additional evidence for a markedly LREE-depleted early terrestrial mantle reservoir. The role of c. 3,850 Ma crust should be considered in interpreting isotope signatures of the younger (3,800–3,600 Ma) rocks of the Itsaq Gneiss Complex. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Environmental problems associated with the development and operation of a lined and unlined landfill site: a case study demonstrating two landfill sites in Patra,Greece

The present work demonstrates the environmental problems associated with the development and operation of a lined and unlined landfill site, both used for municipal solid waste landfilling, that are situated next to the city of Patra, Greece. Findings from a detailed site investigation as well calculations on hydrologic evaluation of landfill’s performance and measurements on leachate composition proved that the lined landfill site is more secure and environmental friendly compared with the unlined landfill site. Even though, until today, there is no evidence from near boreholes for severe contamination problems generated by any of the two landfills, in the forthcoming future several environmental problems are expected to occur from the unlined site. In addition the prevailing hydrogeotechnical conditions indicated that the unlined site is a potential source of contamination; hence extra remedial measures should be received by the local authorities to prevent severe contamination in soil and groundwater.     

Equatorial electrojet in the south Atlantic anomaly region

Features of the equatorial electrojet are studied at Sao Luiz (2.6°S, 44.2°W, inclination −0.25°) in eastern Brazil and Sikasso (11.3°N, 5.7°W, inclination 0.1°) in the western African sector. The stations are situated on either side of the lowest magnetic field intensity in the region of rapid changes in the declination. The daily variations of ΔX at the two stations are almost similar with the peak around noon with maximum values during equinoxes and minimum values during J-solstices. Daily variations of ΔY differ with the maximum deviation of about −35 nT around noon at Sao Luiz and much smaller value of about −10 nT around 14 h LT for Sikasso. The direction of the H vector varies from 15°W of north at 08 h to more than 30°W of north at 17 h for Sao Luiz and from 14°E of north to 25°W of north at 18 h for Sikasso. The plot of the deviations in ΔX and ΔY at different hours for the two stations shows the points along narrow ellipses with major axis aligned along 22°W of north for Sao Luiz and along 3°W of north for Sikasso as compared to declination of 20°W for Sao Luiz and 6°W for Sikasso. The deviations in ΔX at the two stations are fairly well correlated.     

Remote effects of mixed layer development on ocean acidification in the subsurface layers of the North Pacific

Using the outputs of projections under the highest emission scenario of the representative concentration pathways performed by Earth system models (ESMs), we evaluate the ocean acidification rates of subsurface layers of the western North Pacific, where the strongest sink of atmospheric CO₂ is found in the mid-latitudes. The low potential vorticity water mass called the North Pacific Subtropical Mode Water (STMW) shows large dissolved inorganic carbon (DIC) concentration increase, and is advected southwestward, so that, in the sea to the south of Japan, DIC concentration increases and ocean acidification occurs faster than in adjacent regions. In the STMW of the Izu-Ogasawara region, the ocean acidification occurs with a pH decrease of ~0.004 year^?1 , a much higher rate than the previously estimated global average (0.0023 year^?1), so that the pH decreases by 0.3–0.4 during the twenty-first century and the saturation state of calcite (Ω_Ca) decreases from ~4.8 down to ~2.4. We find that the ESMs with a deeper mixed layer in the Kuroshio Extension region show a larger increase in DIC concentration within the Izu-Ogasawara region and within the Ryukyu Islands region. Comparing model results with the mixed layer depth obtained from the Argo dataset, we estimate that DIC concentration at a depth of ~200 m increases by 1.4–1.6 μmol kg^?1 year^?1 in the Izu-Ogasawara region and by 1.1–1.4 μmol kg^?1 year^?1 in the Ryukyu Islands region toward the end of this century.     

CME Velocity and Acceleration Error Estimates Using the Bootstrap Method

The bootstrap method is used to determine errors of basic attributes of coronal mass ejections (CMEs) visually identified in images obtained by the Solar and Heliospheric Observatory (SOHO) mission’s Large Angle and Spectrometric Coronagraph (LASCO) instruments. The basic parameters of CMEs are stored, among others, in a database known as the SOHO/LASCO CME catalog and are widely employed for many research studies. The basic attributes of CMEs (e.g. velocity and acceleration) are obtained from manually generated height-time plots. The subjective nature of manual measurements introduces random errors that are difficult to quantify. In many studies the impact of such measurement errors is overlooked. In this study we present a new possibility to estimate measurements errors in the basic attributes of CMEs. This approach is a computer-intensive method because it requires repeating the original data analysis procedure several times using replicate datasets. This is also commonly called the bootstrap method in the literature. We show that the bootstrap approach can be used to estimate the errors of the basic attributes of CMEs having moderately large numbers of height-time measurements. The velocity errors are in the vast majority small and depend mostly on the number of height-time points measured for a particular event. In the case of acceleration, the errors are significant, and for more than half of all CMEs, they are larger than the acceleration itself.