Vertical profile of aerosol single scattering albedo over west coast of India during W_ICARB

Altitude profile of aerosol Single Scattering Albedo (SSA), derived from simultaneous in-situ airborne measurements of the coefficients of aerosol absorption and scattering off the west coast of India over the Arabian Sea (AS), during January 2009 is presented. While both the absorption and scattering coefficients decreased with altitude, their vertical structure differed significantly. Consequently, the derived SSA, with a surface value of 0.94, decreased with altitude, illustrating increasing relative dominance of aerosol absorption at higher altitudes. Altitude profile of SSA, when examined in conjunction with that of hemispheric backscatter fraction, revealed that the continental influence on the aerosol properties was higher at higher altitude, rather than the effect of marine environment. During an east–west transect across the peninsular India at an altitude of ～2500 m (free troposphere), it was found that the aerosol scattering coefficients remained nearly the same over both east and west coasts.     

Aerosol microphysics over a tropical coastal station inferred from the spectral dependence of Angstrom wavelength exponent and inversion of spectral aerosol optical depths

Six years of spectral aerosol optical depth (AOD) measurements have been analyzed from a tropical coastal site, Trivandrum (8.55°N, 76.9°E, 3 m msl) to infer on the seasonal changes in the microphysical properties of columnar aerosols, by examining the derivatives of the Angstrom wavelength exponent (α) in the wavelength domain (α′_λ) as well as in AOD domain (α′_τ) and by retrieving the columnar size distribution by numerical inversion of the AODs. The inference of the changes in the aerosol microphysics drawn from the features of the derivatives α′_λ and α′_τ is consistent with the pattern revealed by the aerosol properties obtained from the columnar size distributions retrieved from the AOD spectra as well as from the surface measurements of mass-size distributions, which are supported by the back-trajectory cluster analysis and the results of chemical species analysis.     

The Himalayan foreland basin crust and upper mantle

Modeling of multimode surface wave group velocity dispersion data sampling the eastern and the western Ganga basins, reveals a three layer crust with an average V_s of 3.7 km s^?1, draped by ～2.5 km foreland sediments. The Moho is at a depth of 43 ± 2 km and 41 ± 2 km beneath the eastern and the western Ganga basins respectively. Crustal V_p/V_s shows a felsic upper and middle crust beneath the eastern Ganga basin (1.70) compared to a more mafic western Ganga basin crust (1.77). Due to higher radiogenic heat production in felsic than mafic rocks, a lateral thermal heterogeneity will be present in the foreland basin crust. This heterogeneity had been previously observed in the north Indian Shield immediately south of the foreland basin and must also continue northward below the Himalaya. The high heat producing felsic crust, underthrust below the Himalayas could be an important cause for melting of midcrustal rocks and emplacement of leucogranites. This is a plausible explanation for abundance of leucogranites in the east-central Himalaya compared to the west. The uppermost mantle V_s is also significantly lower beneath the eastern Ganga basin (4.30 km s^?1) compared to the west (4.44 km s^?1).     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

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.     

Temperature dependence of the elastic moduli of ringwoodite

The elastic moduli of polycrystalline ringwoodite, (Mg_0.91Fe_0.09)₂SiO₄, were measured up to 470 K by means of the resonant sphere technique. The adiabatic bulk (K_S) and shear (μ) moduli were found to be 185.1(2) and 118.22(6) GPa at room temperature, and the average slopes of dK_S/dT and dμ/dT in the temperature range of the study were determined to be −0.0193(9) and −0.0148(3) GPa/K, respectively. Using these results, we estimate seismic wave velocity jumps for a pure olivine mantle model at 520 km depth. We find that the jump for the S-wave velocity is about 1.5 times larger than that for the P-wave velocity at this depth. This suggests that velocity jumps at the 520 km discontinuity are easier to detect using S-waves than P-waves.     

Density and flow structure in the Clyde Sea front

A well-defined front in temperature and salinity separates the stratified Clyde Sea water from the vertically well mixed water of the North Channel. The detailed structure of the front was observed in autumn 1990 by a combination of, repeated crossings of the front using a ship-borne ADCP and a towed undulating CTD system, and the deployment of a fixed mooring system with temperature, salinity and velocity sensors for a period of 12 days. The results show that the front was situated on the Great Plateau near a contour of log₁₀(H/U³₂)=2.7∼3.7 where H is the water depth and U₂ the amplitude of M₂ tidal velocity. The temperature structure in the Clyde Sea was inverted and the Clyde Sea surface temperature was lower than that of the vertically well mixed water in the North Channel. Since the salinity gradient was stronger than the temperature gradient with fresher water on the surface, the density structure was predominantly controlled by salinity. There were indications of warm and saline bottom water upwelling on the mixed side of the front during spring tides. This upwelling disappeared and the salinity and temperature structure at the front was more diffuse during the neap tide period. A jet-like along-front residual current was observed flowing to the northwest in the surface layer with a counter flow to the southeast in the bottom layer. The vertical difference in velocity was about 9 cm s⁻¹ and was approximately consistent with the shear determined from the thermal wind relation. Both cross- and along-front components of the current observed at the mooring station varied in response to the advection of the front, although both components had large variations with periods of less than one day and several days. The front was advected past the mooring system by a mean flow from the North Channel to the inner basin, while oscillating 3–5 km back and forth with the tidal currents. From the velocity at a current meter mooring and CTD data, the front was estimated to have moved up to 20 km during the observational period and the cross frontal velocity was inferred to be 3–4 cm s⁻¹.     

Coupled silicon–oxygen isotope fractionation traces Archaean silicification

Silica alteration zones and cherts are a conspicuous feature of Archaean greenstone belts worldwide and provide evidence of extensive mobilisation of silica in the marine environment of the early Earth. In order to understand the process(es) of silicification we measured the silicon and oxygen isotope composition of sections of variably silicified basalts and overlying bedded cherts from the Theespruit, Hooggenoeg and Kromberg Formations of the Barberton Greenstone Belt, South Africa.The δ³⁰Si and δ¹⁸O values of bulk rock increase with increasing amount of silicification from unsilicified basalts (?0.64‰ < δ³⁰Si < ?0.01‰ and + 8.6‰ < δ¹⁸O < + 11.9‰) to silicified basalts (δ³⁰Si and δ¹⁸O values as high as + 0.81‰ and + 15.6‰, respectively). Cherts generally have positive isotope ratios (+ 0.21‰ < δ³⁰Si < + 1.05‰ and + 10.9 < δ¹⁸O < + 17.1), except two cherts, which have negative δ³⁰Si values, but high δ¹⁸O (up to + 19.5‰).The pronounced positive correlations between δ³⁰Si, δ¹⁸O and SiO₂ imply that the isotope variation is driven by the silicification process which coevally introduced both ¹⁸O and ³⁰Si into the basalts. The oxygen isotope variation in the basalts from about 8.6‰ to 15.6‰ is likely to represent temperature-dependent isotope fractionation during alteration. Our proposed model for the observed silicon isotope variation relies on a temperature-controlled basalt dissolution vs. silica deposition process.     

Latitudinal patterns of export production recorded in surface sediments of the Chilean Patagonian fjords (41–55°S) as a response to water column productivity

The Chilean Patagonian fjords region (41–56°S) is characterized by highly complex geomorphology and hydrographic conditions, and strong seasonal and latitudinal patterns in precipitation, freshwater discharge, glacier coverage, and light regime; all of these directly affect biological production in the water column. In this study, we compiled published and new information on water column properties (primary production, nutrients) and surface sediment characteristics (biogenic opal, organic carbon, molar C/N, bulk sedimentary δ¹³C_org) from the Chilean Patagonian fjords between 41°S and 55°S, describing herein the latitudinal pattern of water column productivity and its imprint in the underlying sediments. Based on information collected at 188 water column and 118 sediment sampling sites, we grouped the Chilean fjords into four main zones: Inner Sea of Chiloé (41° to ～44°S), Northern Patagonia (44° to ～47°S), Central Patagonia (48–51°S), and Southern Patagonia (Magellan Strait region between 52° and 55°S). Primary production in the Chilean Patagonian fjords was the highest in spring–summer, reflecting the seasonal pattern of water column productivity. A clear north–south latitudinal pattern in primary production was observed, with the highest average spring and summer estimates in the Inner Sea of Chiloé (2427 and 5860 mg C m^?2 d^?1) and Northern Patagonia (1667 and 2616 mg C m^?2 d^?1). This pattern was closely related to the higher availability of nutrients, greater solar radiation, and extended photoperiod during the productive season in these two zones. The lowest spring value was found in Caleta Tortel, Central Patagonia (91 mg C m^?2 d^?1), a site heavily influenced by glacier meltwater and river discharge loaded with glacial sediments. Biogenic opal, an important constituent of the Chilean fjord surface sediments (Si_OPAL ～1–13%), reproduced the general north–south pattern of primary production and was directly related to water column silicic acid concentrations. Surface sediments were also rich in organic carbon content and the highest values corresponded to locations far away from glacier influence, sites within fjords, and/or semi-enclosed and protected basins, reflecting both autochthonous (water column productivity) and allochthonous sources (contribution of terrestrial organic matter from fluvial input to the fjords). A gradient was observed from the more oceanic sites to the fjord heads (west–east) in terms of bulk sedimentary δ¹³C_org and C/N ratios; the more depleted (δ¹³C_org ?26‰) and higher C/N (23) values corresponded to areas close to rivers and glaciers. A comparison of the Chilean Patagonian fjords with other fjord systems in the world revealed high variability in primary production for all fjord systems as well as similar surface sediment geochemistry due to the mixing of marine and terrestrial organic carbon.     

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region,Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V_p/V_s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.     

Magnesium isotopes constraints on the origin of Mg-rich olivines from the Allende chondrite: Nebular versus planetary?

High precision Mg isotope measurements by multi-collector ion microprobe show that refractory olivines from the Allende chondrite, either olivines isolated in the matrix (2 samples studied) or olivines in type I chondrules (6 samples studied), have variable δ²⁶Mg* enrichments and deficits (calculated in permil as the ²⁶Mg deviation from the instrumental mass fractionation line) relative to the Earth. Most average δ²⁶Mg* (noted δ²⁶Mg*_av) values (between 10 and 20 analyses per chondrule) are negative but the total range is from ?0.029 (± 0.010) ‰ (2 sigma errors) to + 0.011 (± 0.011) ‰ with an exception of one olivine at + 0.043 (± 0.023) ‰. These variations in δ²⁶Mg*_av reflect the formation of the olivines from reservoirs enriched in various amounts of ²⁶Mg by the decay of short-lived ²⁶Al (T_1/2 = 0.73 Ma). Similarly, 30 analyses of olivines from the Eagle Station pallasite show a δ²⁶Mg*_av value of ?0.033 ± 0.008‰, as negative as some olivines from Allende chondrules and the Solar system initial δ²⁶Mg* value of ?0.038 ± 0.004‰ (defined at the time of formation of type B Ca–Al-rich inclusions – CAIs – when ²⁶Al/²⁷Al = 5.23 × 10^?5, Jacobsen et al., 2008).Because olivines are Al-poor and because their Mg isotopic compositions are not reset during the chondrule forming events, their δ²⁶Mg*_av can be used to calculate model crystallization ages relative to various theoretical Mg isotope growth curves. The two end-member scenarios considered are (i) a “nebular” growth in which the Al/Mg ratio remains chondritic and (ii) a “planetary” growth in which a significant increase of the Al/Mg ratio can be due to, for instance, olivine magmatic fractionation. The low δ²⁶Mg*_av value of olivines from the Eagle Station pallasite demonstrate that metal-silicate differentiation occurred as early as ~ 0. 15_- 0. 23^+ 0. 29 Ma after CAIs in either of the growth scenarios. Similarly the variable δ²⁶Mg*_av values of refractory olivines can be understood if they were formed in planetesimals which started to differentiate as early as the Eagle Station parent body. Accretion of these planetesimals must have been coeval to the formation of CAIs and their disruption could explain why their fragments (Mg-rich olivines) were distributed in the chondrule forming regions of the disk.     

Latitudinal dependence of the seasonal variation of particulate extinction in the UTLS over the Indian longitude sector during volcanically quiescent period based on lidar and SAGE-II observations

The altitude profiles of particulate extinction in the upper troposphere and lower stratosphere (UTLS) obtained from SAGE-II in the latitude region 0–30°N over the Indian longitude sector (70–90°E) are used to study the latitudinal variation of its annual pattern in this region during the volcanically quiescent period of 1998–2003. The SAGE-II data is compared with the lidar measurements from Gadanki (13.5°N, 79.2°E) when the satellite had an overhead occultation pass over a small geographical grid centered at this location. The particulate optical depth (τ_p) in the UT region shows a general decrease with increase in latitude and a pronounced summer–winter contrast with relatively low values during winter and high values during summer. In general, these variations are in accordance with the latitudinal variation of convective available potential energy (CAPE) and thunderstorm activity, which are good representative indices of tropospheric convection. While the particulate extinction (and τ_p) in the 18–21 km (LS₁) region is relatively low in the equatorial region up to 15°N, it shows an increase in the off-equatorial region, beyond 15°N. While the annual variation of τ_p in the LS₁ region is almost insignificant near the equator, it is rather well pronounced in latitude region between 10 and 15°N with relatively high values during winter and low values during summer. Beyond 20°N, this shows a prominent peak during summer. At a higher altitude, the 21–30 km (LS₂) region, the latitude variation of τ_p shows a different pattern with high values near the equator and low values in the off-equatorial region confirming the existence of a stratospheric aerosol reservoir. Low values of τ_p at lower regime (LS₁) near the equator could be due to rapid transport of particulates from the near equatorial region to higher latitudes, while the equatorial high at upper regime (LS₂) could be due to lofting and subsequent accumulation.     

Seismotectonic of Southern Apennines from recent passive seismic experiments

We used data of local earthquakes collected during two recent passive seismic experiments carried out in southern Italy in order to study the seismotectonic setting of the Lucanian Apennine and the surrounding areas. Based on continuous recordings of the temporary stations we extracted over 15,600 waveforms, which were hand-picked along with those recorded by the permanent stations of the Italian national seismic network obtaining a dense, high-quality dataset of P- and S-arrival times. We examined the seismicity occurring in the period 2001–2008 by relocating 566 out of 1047 recorded events with magnitudes M_L ≥ 1.5 and computing 162 fault-plane solutions. Earthquakes were relocated using a minimum one-dimensional velocity model previously obtained for the region and a V_p/V_s ratio of 1.83. Background seismicity is concentrated within the upper crust (between 5 and 20 km of depth) and it is mostly clustered along the Lucanian Apennine chain axis. A significant feature extracted from this study relates to the two E–W trending clusters located in the Potentino and in the Abriola–Pietrapertosa sector (central Lucania region). Hypocentral depths in both clusters are slightly deeper than those observed beneath the Lucanian Apennine. We suggest that these two seismic features are representative of the transition from the inner portion of the chain to the external margin characterized by dextral strike-slip kinematics. In the easternmost part of the study area, below the Bradano foredeep and the Apulia foreland, seismicity is generally deeper and more scattered. The sparse seismicity localized in the Sibari Plain, in the offshore area along the northeastern Calabrian coast and in the Taranto Gulf is also investigated thanks to the new recordings. This seismicity shows hypocenters between 12 and 20 km of depth below the Sibari Plain and is deeper (foci between 10 and 35 km of depth) in the offshore area of the Taranto Gulf. 102 well-constrained fault-plane solutions, showing predominantly normal and strike-slip character with tensional axes (T-axes) generally NE oriented, were selected for the stress tensor analysis. We investigated stress field orientation inverting focal mechanism belonging to the Lucanian Apennine and the Pollino Range, both areas characterized by a more concentrated background seismicity.     

Time-Domain Moment Tensors for shallow (h ≤ 40 km) earthquakes in the broader Aegean Sea for the years 2006 and 2007: The database of the Aristotle University of Thessaloniki

We present a catalog of moment tensor (MT) solutions and moment magnitudes, M_w, for 119 shallow (h ≤ 40 km) earthquakes in Greece and its surrounding lands (34°N–42°N, 19°E–30°E) for the years 2006 and 2007, computed with the 1D Time-Domain Moment Tensor inversion method (TDMT_INV code of Dreger, 2003). Magnitudes range from 3.2 ≤ M_w ≤ 5.7. Green's functions (GF) have been pre-computed to build a library, for a number of velocity profiles applicable to the broader Aegean Sea region, to be used in the inversion of observed broad band waveforms (10–50 s). All MT solutions are the outcome of a long series of tests of different reported source locations and hypocenter depths. Quality factors have been assigned to each MT solution based on the number of stations used in the inversion and the goodness of fit between observed and synthetic waveforms. In general, the focal mechanisms are compatible with previous knowledge on the seismotectonics of the Aegean area. The new data provide evidence for strike-slip faulting along NW–SE trending structures at the lower part of Axios basin, close to the heavily industrialized, and presently subsiding, region of the city of Thessaloniki. Normal faulting along E–W trending planes is observed at the Strimon basin, and in Orfanou Gulf in northern Greece. A sequence of events in the east Aegean Sea close to the coastline with western Anatolia sheds light on an active structure bounding the north coastline of Psara–Chios Islands about 20–25 km in length exhibiting right lateral strike-slip faulting.     

Evolution of the oceanic calcium cycle during the late Mesozoic: Evidence from δ44/40Ca of marine skeletal carbonates

The Ca isotope compositions of 37 late Mesozoic skeletal carbonates, belemnites and brachiopods, from the Tethyan realm were analyzed by thermal (TIMS) and plasma (MC-ICP-MS) ionization mass spectrometry. A poor correlation between δ^44/40Ca and δ¹⁸O values of belemnites suggests only a weak temperature dependency for the Ca isotope composition of belemnites, likely less than 0.02‰/°C. The δ^44/40Ca record of belemnites was therefore used to reconstruct the Ca isotope composition of paleo-seawater (δ^44/40Ca_SW), based on an experimentally determined fractionation factor between seawater Ca and belemnite calcite (α_CC–SW) of ∼ 0.9986. The inferred δ^44/40Ca_SW record, with an average stratigraphic resolution of 1 Ma, shows systematic temporal variation of ∼ 0.5‰ with the Middle/Late Jurassic (∼ 154 Ma) minimum of ∼ 1.4‰ and a subsequent general increase to the Early Cretaceous (∼ 124 Ma) maximum of ∼ 1.9‰. The global nature of the δ^44/40Ca_SW record is supported by identical Ca isotope compositions of coeval (Kimmeridgian) belemnites collected from two distinct paleogeographic regions, the southern (New Zealand) and northern (Germany) margin of the Tethys Ocean. The observed late Mesozoic δ^44/40Ca_SW record was simulated using a simple Ca isotope mass balance model, and the results indicate that the variation in δ^44/40Ca_SW record can be explained by changes in oceanic input fluxes of Ca that were independent of the carbonate ion fluxes, such as the hydrothermal Ca flux or the release of Ca to the oceans via dolomitization of marine carbonates.     

Variations of daytime and nighttime electron temperature and heat flux in the upper ionosphere,topside ionosphere and lower plasmasphere for low and high solar activity

A database of the electron temperature (T_e) comprising of most of the available LEO satellite measurements is used for studying the solar activity variations of T_e. The T_e data are grouped for two levels of solar activity (low LSA and high HSA), five altitude ranges between 350 and 2000 km, and day and night. By fitting a theoretical expression to the T_e values we obtain variation of T_e along magnetic field lines and heat flux for LSA and HSA. We have found that T_e increases with increase in solar activity at low and mid-latitudes during nighttime at all altitudes studied. During daytime the T_e response to solar activity depends on latitude, altitude, and season. This analysis shows existence of anti-correlation between T_e and solar activity at mid-latitudes below 700 km during the equinox and winter day hours. Heat fluxes show small latitudinal dependence for daytime but substantial for nighttime.     

Inferred fault geometry and slip distribution of the 2010 Jiashian,Taiwan, earthquake is consistent with a thick-skinned deformation model

We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 10¹⁸ N-m is equivalent to a M_w = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.     

Characterization of the crust of the Coast Mountains Batholith,British Columbia,from P to S converted seismic waves and petrologic modeling

The late Triassic to early Tertiary Coast Mountains Batholith (CMB) of British Columbia provides an ideal locale to study the processes whereby accreted terranes and subduction-related melts interact to form stable continental crust of intermediate to felsic composition and complementary ultramafic residuals. Seismic measurements, combined with calculated elastic properties of various CMB rock compositions, provide a window into the deep-crustal lithologies that are key to understanding the processes of continental growth and evolution. We use a combination of seismic observations and petrologic modeling to construct hypothetical crustal sections at representative locations across the CMB, then test the viability of these sections via forward modeling with synthetic seismic data. The compositions that make up our petrologic forward models are based on calculations using the free energy minimization program Perple_X to predict mineral assemblages at depth for the bulk compositions of exposed plutonic rocks collected in the study area. Seismic data were collected along two transects in west-central British Columbia: a southern line that crossed the CMB near the town of Bella Coola (near 52° N), and a northern line centered on the towns of Terrace and Kitimat (near 54° N). Along both transects, seismic receiver functions reveal high V_p/V_s ratios near the Insular/Intermontane terrane boundary and crustal thickness increasing from 26 ± 3 km to 34 ± 3 km (at the 1 sigma certainty level) from west to east across the Coast Shear Zone (CSZ). On the southern line, we observe an anomalous region of complex receiver functions and diminished Moho signals beneath the central portion of the CMB. Our petrologic and seismic profiles show that observed seismic data from much of the CMB can be well-matched in terms of crustal thickness and structure, average V_p/V_s, and amplitude of the Moho converted phase, without including ultramafic residual material in the lower crust.     

On the tidal heating of Enceladus

Enceladus, one of Saturn's moons, shows significant volcanic activity identified by the Cassini spacecraft. The aim of the present study is to investigate – with the adaptation of mathematical tools used in geodynamics – the extent of tidal heating due to the mean motion resonance with Dione. For the purpose of calculations a two-layer model of Enceladus was used. The inner part of the model is a “rocky core” with a relative radius 0.55, while the outer part is composed of water ice. The results of model calculations show that the effective tidal heating is not uniformly distributed within Enceladus. It was found for the selected model of Enceladus, that the tidal heating is maximum within the depth interval (25–75) km. Due to the inhomogeneity within Enceladus, 85% of the tidal energy is generated in a volume that contains just 39% of its mass. In time intervals of 3.0 × 10⁸ and 5.3 × 10⁸ years the temperature increase in the relative depth range 0.70 ≤ r/a_E ≤ 0.90 is approximately 270 and 370 K, respectively.     

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.