Estimating Curie Point Depth and Heat Flow Map for Northern Red Sea Rift of Egypt and Its Surroundings, from Aeromagnetic Data

In this study, we aim to map the Curie point depth surface for the northern Red Sea rift region and its surroundings based on the spectral analysis of aeromagnetic data. Spectral analysis technique was used to estimate the boundaries (top and bottom) of the magnetized crust. The Curie point depth (CPD) estimates of the Red Sea rift from 112 overlapping blocks vary from 5 to 20 km. The depths obtained for the bottom of the magnetized crust are assumed to correspond to Curie point depths where the magnetic layer loses its magnetization. Intermediate to deep Curie point depth anomalies (10–16 km) were observed in southern and central Sinai and the Gulf of Suez (intermediate heat flow) due to the uplifted basement rocks. The shallowest CPD of 5 km (associated with very high heat flow, ~235 mW m^?2) is located at/around the axial trough of the Red Sea rift region especially at Brothers Island and Conrad Deep due to its association with both the concentration of rifting to the axial depression and the magmatic activity, whereas, beneath the Gulf of Aqaba, three Curie point depth anomalies belonging to three major basins vary from 10 km in the north to about 14 km in the south (with a mean heat flow of about 85 mW m^?2). Moreover, low CPD anomalies (high heat flow) were also observed beneath some localities in the northern part of the Gulf of Suez at Hammam Fraun, at Esna city along River Nile, at west Ras Gharib in the eastern desert and at Safaga along the western shore line of the Red Sea rift. These resulted from deviatoric tensional stresses developing in the lithosphere which contribute to its further extension and may be due to the opening of the Gulf of Suez and/or the Red Sea rift. Furthermore, low CPD (with high heat flow anomaly) was observed in the eastern border of the study area, beneath northern Arabia, due to the quasi-vertical low-velocity anomaly which extends into the lower mantle and may be related to volcanism in northern Arabia. Dense microearthquakes seem to occur in areas where the lateral gradients of the CPD are steep (e.g. entrance of the Gulf of Suez and Brothers Island in the Red Sea). These areas may correspond to the boundaries between high and low thermal regions of the crust. Thus, the variations in the microseismic activity may be closely related to thermal structures of the crust. Indeed, shallow cutoff depths of seismicity can also be found in some geothermal areas (e.g. western area of Safaga city along the Red Sea coastal region and at Esna city along the River Nile). These facts indicate that the changes in the thickness of the seismogenic layer strongly depend on temperature. Generally, the shallow Curie point depth indicates that some regions in our study area are promising regions for further geothermal exploration particularly in some localities along the River Nile, Red Sea and Gulf of Suez coastal regions.     

Fan and playa areas in southern California and adjacent parts of Nevada

In bolsons in the desert regions of southern California and adjacent parts of Nevada, the area underlain by alluvial fan gravels and playa sediments is generally ?1·2 times the area being eroded to produce those sediments. In certain larger basins in the vicinity of Death Valley, however, the depositional area is only about half the size of the erosional area. This reflects the more active tectonic environment in these bolsons. Of the areas underlain by recent sediments in these bolsons, playas make up 2–6 per cent. Smaller playas are found in the Mojave region, and seem to be associated with sedimentary terranes. Conversely, igneous terranes support larger playas. Larger deviations of playa area from these averages are attributable to incorrect identification of bolson boundaries. Fine sediment is either able to pass through the bolson to the next down stream, or is being collected from areas upstream that were not considered to be part of the system. For example, the playa in Death Valley is unusually large. This is in part because the Death Valley playa has been deformed tectonically so parts of it are now eroding, and in part because the outlet of Lake Tecopa was downcut in the geologically recent past, so sediment once trapped there now reaches Death Valley. The size of the playa in Death Valley is still adjusting to these changes.     

Segmentation of alluvial fans in death valley,california: New insights from surface exposure dating and laboratory modelling

Laboratory experiments, recent paleoenvironmental analyses of rock varnish, and surface exposure dating of geomorphic units have led to new insights into the process of entrenchment and segmentation of alluvial fans, and into the history of Quaternary sedimentation in Death Valley. Entrenchment begins at the fanhead. As the trench deepens, its down-slope end migrates down-fan, taking several tens of thousands of years to reach lower parts of the fan. Laboratory experiments suggest, however, that a new segment begins to grow at the toe long before the trench reaches this part of the fan. Furthermore, the initial slope of the segment is not the equilibrium slope. Field evidence supports this model. The tectonic tilting that caused entrenchment and segmentation in Death Valley may have been triggered by loading of the valley with water. Sedimentation on the salt pan in southern Death Valley is not, at present, in equilibrium with that on the fans. Rather, it seems to be adjusting to an increase in the amount of fine material reaching the playa, due in part to breaching of the outlet of Lake Tecopa somewhat after 600 ka BP, and in part to subsidence of different parts of the valley at different rates. Failure to recognize this disequilibrium resulted in errors in earlier estimates of the age of the segmentation events.     

Hydrogeology of desert springs in the Panamint Range,California, USA: Geologic controls on the geochemical kinetics,flowpaths, and mean residence times of springs

Over 180 springs emerge in the Panamint Range near Death Valley National Park, CA, yet, these springs have received very little hydrogeological attention despite their cultural, historical, and ecological importance. Here, we address the following questions: (1) which rock units support groundwater flow to springs in the Panamint Range, (2) what are the geochemical kinetics of these aquifers, and (3) and what are the residence times of these springs? All springs are at least partly supported by recharge in and flow through dolomitic units, namely, the Noonday Dolomite, Kingston Peak Formation, and Johnnie Formation. Thus, the geochemical composition of springs can largely be explained by dedolomitization: the dissolution of dolomite and gypsum with concurrent precipitation of calcite. However, interactions with hydrothermal deposits have likely influenced the geochemical composition of Thorndike Spring, Uppermost Spring, Hanaupah Canyon springs, and Trail Canyon springs. Faults are important controls on spring emergence. Seventeen of twenty-one sampled springs emerge at faults (13 emerge at low-angle detachment faults). On the eastern side of the Panamint Range, springs emerge where low-angle faults intersect nearly vertical Late Proterozoic, Cambrian, and Ordovician sedimentary units. These geologic units are not present on the western side of the Panamint Range. Instead, springs on the west side emerge where low-angle faults intersect Cenozoic breccias and fanglomerates. Mean residence times of springs range from 33 (±30) to 1,829 (±613) years. A total of 11 springs have relatively short mean residence times less than 500 years, whereas seven springs have mean residence times greater than 1,000 years. We infer that the Panamint Range springs are extremely vulnerable to climate change due to their dependence on local recharge, disconnection from regional groundwater flow (Death Valley Regional Flow System - DVRFS), and relatively short mean residence times as compared with springs that are supported by the DVRFS (e.g., springs in Ash Meadows National Wildlife Refuge). In fact, four springs were not flowing during this campaign, yet they were flowing in the 1990s and 2000s.     

Slip Distribution and Seismic Moment of the 2010 and 1960 Chilean Earthquakes Inferred from Tsunami Waveforms and Coastal Geodetic Data

The slip distribution and seismic moment of the 2010 and 1960 Chilean earthquakes were estimated from tsunami and coastal geodetic data. These two earthquakes generated transoceanic tsunamis, and the waveforms were recorded around the Pacific Ocean. In addition, coseismic coastal uplift and subsidence were measured around the source areas. For the 27 February 2010 Maule earthquake, inversion of the tsunami waveforms recorded at nearby coastal tide gauge and Deep Ocean Assessment and Reporting of Tsunamis (DART) stations combined with coastal geodetic data suggest two asperities: a northern one beneath the coast of Constitucion and a southern one around the Arauco Peninsula. The total fault length is approximately 400 km with seismic moment of 1.7 × 10²² Nm (Mw 8.8). The offshore DART tsunami waveforms require fault slips beneath the coasts, but the exact locations are better estimated by coastal geodetic data. The 22 May 1960 earthquake produced very large, ~30 m, slip off Valdivia. Joint inversion of tsunami waveforms, at tide gauge stations in South America, with coastal geodetic and leveling data shows total fault length of ~800 km and seismic moment of 7.2 × 10²² Nm (Mw 9.2). The seismic moment estimated from tsunami or joint inversion is similar to previous estimates from geodetic data, but much smaller than the results from seismic data analysis.     

Geometry and growth of sill complexes: insights using 3D seismic from the North Rockall Trough

Doleritic sill complexes, which are an important component of volcanic continental margins, can be imaged using 3D seismic reflection data. This allows unprecedented access to the complete 3D geometry of the bodies and an opportunity to test classic sill emplacement models. The doleritic sills associated with basaltic volcanism in the North Rockall Trough occur in two forms. Radially symmetrical sill complexes consist of a saucer-like inner sill at the base with an arcuate inclined sheet connecting it to a gently inclined, commonly ragged, outer rim. Bilaterally symmetrical sill complexes are sourced by magma diverted from a magma conduit feeding an overlying volcano. With an elongate, concave upwards, trough-like geometry bilaterally symmetrical sills climb away from the magma source from which they originate. Both sill complex types can appear as isolated bodies but commonly occur in close proximity and consequently merge, producing hybrid sill complexes. Radial sill complexes consist of a series of radiating primary flow units. With dimensions up to 3 km, each primary flow unit rises from the inner saucer and is fed by primary magma tube. Primary flow units contain secondary flow units with dimensions up to 2 km, each being fed by a secondary magma tube branching from the primary magma tube. Secondary flow units in turn are composed of 100-m scale tertiary flow units. A similar branching hierarchy of flow units can also be seen in bilaterally symmetrical sill complexes, with their internal architecture resembling an enlarged version of a primary flow unit from a radial sill complex. This branching flow pattern, as well as the interaction between flow units of varying orders, provides new insights into the origin of the structures commonly seen within sill complexes and the hybrid sill bodies produced by their merger. The data demonstrate that each radially symmetrical sill complex is independently fed from a source located beneath the centre of the inner saucer, grows by climbing from the centre outwards and that peripheral dyking from the upper surface is a common feature. These features suggest a laccolith emplacement style involving peripheral fracturing and dyking during inner saucer growth and thickening. The branching hierarchy of flow units within bilaterally symmetrical sill complexes is broadly similar to that of primary flow units within a radially symmetrical sill complex, suggesting that the general features of the laccolith emplacement model also apply.Editorial responsibility: J. Stix     

基于InSAR形变数据分析台湾纵谷断层南段现今运动特征

本文利用2007—2010年花东纵谷南段区域的InSAR形变数据作为约束, 采用分段断层模型和层状介质模型, 反演中国台湾东部纵谷断层南段滑动速率空间分布, 并据此分析断层运动特征。 研究结果表明, 纵谷断层南段整体以逆冲运动为主, 兼具左旋走滑运动。 纵谷断层南段的滑动速率具有空间非均匀性, 在空间上可以细分为深浅两个极值区, 浅部(0~15 km)最大滑动速率为10 cm/a, 位于深度2.5 km左右; 深部(15~30 km)最大滑动速率为21 cm/a, 位于深度25 km左右。 反演结果与用重复地震估算的深部滑动速率基本吻合。     

Computation of Solar Radiative Fluxes by 1D and 3D Methods Using Cloudy Atmospheres Inferred from A-train Satellite Data

This study used realistic representations of cloudy atmospheres to assess errors in solar flux estimates associated with 1D radiative transfer models. A scene construction algorithm, developed for the EarthCARE mission, was applied to CloudSat, CALIPSO and MODIS satellite data thus producing 3D cloudy atmospheres measuring 61 km wide by 14,000 km long at 1 km grid-spacing. Broadband solar fluxes and radiances were then computed by a Monte Carlo photon transfer model run in both full 3D and 1D independent column approximation modes. Results were averaged into 1,303 (50 km)² domains. For domains with total cloud fractions A _c  < 0.7 top-of-atmosphere (TOA) albedos tend to be largest for 3D transfer with differences increasing with solar zenith angle. Differences are largest for A _c  > 0.7 and characterized by small bias yet large random errors. Regardless of A _c , differences between 3D and 1D transfer rarely exceed ±30 W m^?2 for net TOA and surface fluxes and ±10 W m^?2 for atmospheric absorption. Horizontal fluxes through domain sides depend on A _c with ～20% of cases exceeding ±30 W m^?2; the largest values occur for A _c  > 0.7. Conversely, heating rate differences rarely exceed ±20%. As a cursory test of TOA radiative closure, fluxes produced by the 3D model were averaged up to (20 km)² and compared to values measured by CERES. While relatively little attention was paid to optical properties of ice crystals and surfaces, and aerosols were neglected entirely, ～30% of the differences between 3D model estimates and measurements fall within ±10 W m^?2; this is the target agreement set for EarthCARE. This, coupled with the aforementioned comparison between 3D and 1D transfer, leads to the recommendation that EarthCARE employ a 3D transport model when attempting TOA radiative closure.     

A study on the uncertainties of the centroid depth of the 2013 Lushan earthquake from teleseimic body wave data

Centroid depth of earthquakes is essential for seismic hazard mitigation. But, various studies provided different solutions for the centroid depth of the damaging 2013 Lushan earthquake, thus hindering further studies of the earthquake processes. To resolve its centroid depth and assess the uncertainties, we apply the teleseismic cut and paste method to invert for centroid depth with teleseismic body waves in the epicentral distance of 30°–90°. We performed the inversion for P waves only as well the case of both P and SH waves and found that both cases lead to depth solutions with difference less than 0.5 km. We also investigated the effects on depth inversion from azimuth gap of seismic stations, source duration, and corner frequency of filter. These various tests show that even azimuthal distribution of seismic stations is helpful for accurate depth inversion. It is also found that estimate of centroid depth is sensitive to source duration. Moreover, the depth is biased to larger values when corner frequency of low-pass filter is very low. The uncertainty in the velocity model can also generate some error in the depth estimation (~1.0 km).With all the above factors considered, the centroid depth of Lushan earthquake is proposed to be around 12 km, with uncertainty about 2 km.     

3-D density model of Mt. Etna Volcano (Southern Italy)

A detailed density model of Mt. Etna and its surrounding areas has been evaluated using a 3-D inversion of the gravimetric data acquired in the 1980's. Several high-density and low-density bodies are found, penetrating from shallow depths as far down as 12 km bsl. A positive correlation (in terms of location, extent, density, and velocity) is established between several anomalies of the density model and features identified in previously published seismic tomographies. A prominent high-density body extending down to 7 km bsl is recognized in the southern part of the Valle del Bove, and interpreted as a solidified magmatic intrusion. On the western boundary of this anomaly, a low-density body is interpreted as a bubble and liquid magma mixture. Outside the central area, three other high-density anomalies are imaged and attributed to the earliest phases of volcanic activity in the area. Several interesting low-density anomalies are also identified and correlated with known fault lines and other structural features of the region.     

3D <Emphasis Type="Italic">v</Emphasis><Subscript>P</Subscript> and <Emphasis Type="Italic">v</Emphasis><Subscript>S</Subscript> models of southeastern margin of the Tibetan plateau from joint inversion of body-wave arrival times and surface-wave dispersion data

A new 3D velocity model of the crust and upper mantle in the southeastern (SE) margin of the Tibetan plateau was obtained by joint inversion of body- and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v _S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 M _S 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes (M ≥ 5) in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.     

Crust and upper mantle structures beneath Northeast China from receiver function studies

P-wave and S-wave receiver function analyses have been performed along a profile consisted of 27 broadband seismic stations to image the crustal and upper mantle discontinuities across Northeast China. The results show that the average Moho depth varies from about 37 km beneath the Daxing’anling orogenic belt in the west to about 33 km beneath the Songliao Basin, and to about 35 km beneath the Changbai mountain region in the east. Our results reveal that the Moho is generally flat beneath the Daxing’anling region and a remarkable Moho offset (about 4 km) exists beneath the basin-mountain boundary, the Daxing’anling-Taihang Gravity Line. Beneath the Tanlu faults zone, which seperates the Songliao Basin and Changbai region, the Moho is uplift and the crustal thickness changes rapidly. We interpret this feature as that the Tanlu faults might deeply penetrate into the upper mantle, and facilitate the mantle upwelling along the faults during the Cenozoic era. The average depth of the lithosphere-asthenosphere boundary (LAB) is ~80 km along the profile which is thinner than an average thickness of a continental lithosphere. The LAB shows an arc-like shape in the basin, with the shallowest part approximately beneath the center of the basin. The uplift LAB beneath the basin might be related to the extensive lithospheric stretching in the Mesozoic. In the mantle transition zone, a structurally complicated 660 km discontinuity with a maximum 35 km depression beneath the Changbai region is observed. The 35 km depression is roughly coincident with the location of the stagnant western pacific slab on top of the 660 km discontinuity revealed by the recent P wave tomography.     

Investigation into regional thermal structure of the Thrace Region, NW Turkey, from aeromagnetic and borehole data

The aeromagnetic values over the study region are relatively uniform except for a few anomalies in the northeastern and southwestern areas. Analyses of aeromagnetic data were performed in NW Turkey, in order to have a look into the subsurface regional thermal structure of the region. For this purpose, power spectra, reduced to pole (RTP), and band-pass filtered anomalies were produced using geophysical techniques. Band-pass filtered data were produced from the RTP aeromagnetic anomalies to isolate near surface and undesired deep effects. Based on the aeromagnetic data interpretation, the thickness of the magnetized crust, named the Curie Point Depth (CPD), in the study area lies between 9.7 and 20.3 km. The CPD estimates in the Thrace region of Turkey indicate two shallow CPD (SCPD1 and SCPD2) zones (the Istranca Massif and the Saros Graben area). The deep CPD are located within the Thrace Basin with sediment thickness of about 9 km. The corresponding heat flow map prepared from the averaged thermal conductivities and thermal gradients from the CPD reveals the existence of one low heat flow zone (75 mW/m²) over the center of Thrace Basin, and two high heat flow zones over the Istranca Masif (100–125 mW/m²) in the northern side and Saros Graben (125–135 mW/m²) areas in the southern side of the Thrace Basin.     

Using multivariate statistical analysis of groundwater major cation and trace element concentrations to evaluate groundwater flow in a regional aquifer

Groundwater samples were collected from 11 springs in Ash Meadows National Wildlife Refuge in southern Nevada and seven springs from Death Valley National Park in eastern California. Concentrations of the major cations (Ca, Mg, Na and K) and 45 trace elements were determined in these groundwater samples. The resultant data were subjected to evaluation via the multivariate statistical technique principal components analysis (PCA), to investigate the chemical relationships between the Ash Meadows and Death Valley spring waters, to evaluate whether the results of the PCA support those of previous hydrogeological and isotopic studies and to determine if PCA can be used to help delineate potential groundwater flow patterns based on the chemical compositions of groundwaters. The results of the PCA indicated that groundwaters from the regional Paleozoic carbonate aquifers (all of the Ash Meadows springs and four springs from the Furnace Creek region of Death Valley) exhibited strong statistical associations, whereas other Death Valley groundwaters were chemically different. The results of the PCA support earlier studies, where potentiometric head levels, δ¹⁸O and δD, geological relationships and rare earth element data were used to evaluate groundwater flow, which suggest groundwater flows from Ash Meadows to the Furnace Creek springs in Death Valley. The PCA suggests that Furnace Creek groundwaters are moderately concentrated Ash Meadows groundwater, reflecting longer aquifer residence times for the Furnace Creek groundwaters. Moreover, PCA indicates that groundwater may flow from springs in the region surrounding Scotty's Castle in Death Valley National Park, to a spring discharging on the valley floor. The study indicates that PCA may provide rapid and relatively cost‐effective methods to assess possible groundwater flow regimes in systems that have not been previously investigated. Copyright © 1999 John Wiley & Sons, Ltd.     

Hydrogeology of desert springs in the Panamint Range,California, USA: Identifying the sources and amount of recharge that support spring flow

Despite its location in the rain shadow of the southern Sierra Nevada, the Panamint Range hosts a complex mountain groundwater system supporting numerous springs which have cultural, historical, and ecological importance. The sources of recharge that support these quintessential desert springs remain poorly quantified since very little hydrogeological research has been completed in the Panamint Range. Here we address the following questions: (i) what is the primary source of recharge that supports springs in the Panamint Range (snowmelt or rainfall), (ii) where is the recharge occurring (mountain-block, mountain-front, or mountain-system) and (iii) how much recharge occurs in the Panamint Range? We answer questions (i) and (ii) using stable isotopes measured in spring waters and precipitation, and question (iii) using a chloride mass-balance approach which is compared to a derivation of the Maxey–Eakin equation. Our dataset of the stable isotopic composition (δ¹⁸O and δ²H) of precipitation is short (1.5 years), but analyses on spring water samples indicate that high-elevation snowmelt is the dominant source of recharge for these springs, accounting for 57 (±9) to 79 (±12) percent of recharge. Recharge from rainfall is small but not insignificant. Mountain-block recharge is the dominant recharge mechanism. However, two basin springs emerging along the western mountain-front of the Panamint Range in Panamint Valley appear to be supported by mountain-front and mountain-system recharge, while Tule Spring (a basin spring emerging at the terminus of the bajada on the eastern side of the Panamint Range) appears to be supported by mountain-front recharge. Calculated recharge rates range from 19 mm year⁻¹ (elevations < 1000 mrsl) to 388 mm year⁻¹ (elevations > 1000 mrsl). The average annual recharge is approximately 91 mm year⁻¹ (equivalent to 19.4 percent of total annual precipitation). We infer that the springs in the Panamint Range (and their associated ecosystems) are extremely vulnerable to changes in snow cover associated with climate change. They are heavily dependent on snowmelt recharge from a relatively thin annual snowpack. These findings have important implications for the vulnerability of desert springs worldwide.     

Volcanic systems and calderas in the Vatnajökull region,central Iceland: Constraints on crustal structure from gravity data

A Bouguer anomaly map is presented of southern central Iceland, including the western part of Vatnajökull and adjacent areas. A complete Bouguer reduction for both ice surface and bedrock topography is carried out for the glaciated regions. Parts of the volcanic systems of Vonarskarð-Hágöngur, Bárðarbunga-Veiðivötn, Grímsvötn-Laki, and to a lesser extent Kverkfjöll, show up as distinct features on the gravity map. The large central volcanoes with calderas: Vonarskarð, Bárðarbunga, Kverkfjöll and Grímsvötn, are associated with 15–20 mGal gravity highs caused by high density bodies in the uppermost 5 km of the crust. Each of these bodies is thought to be composed of several hundred km³ of gabbros that have probably accumulated over the lifetime of the volcano. The Skaftárkatlar subglacial geothermal areas are not associated with major anomalous bodies in the upper crust. The central volcanoes of Vonarskarð and Hágöngur belong to the same volcanic system; this also applies to Bárðarbunga and Hamarinn, and Grímsvötn and Þórðarhyrna. None of the smaller of the two volcanoes sharing a system (Hágöngur, Hamarinn and Þórðarhyrna) is associated with distinct gravity anomalies and clear caldera structures have not been identified. However, ridges in the gravity field extend between each pair of central volcanoes, indicating that they are connected by dense dyke swarms. This suggests that when two central volcanoes share the same system, one becomes the main pathway for magma, forming a long-lived crustal magma chamber, a caldera and large volume basic intrusive bodies in the upper crust. Short residence times of magma in the crust beneath these centres favour essentially basaltic volcanism. In the case of the second, auxillary central volcano, magma supply is limited and occurs only sporadically. This setting may lead to longer residence times of magma in the smaller central volcanoes, favouring evolution of the magma and occasional eruption of rhyolites. The eastern margin of the Eastern Volcanic Zone is marked by a NE–SW lineation in the gravity field, probably caused by accumulation of low density, subglacially erupted volcanics within the volcanic zone. This lineation lies 5–10 km to the east of Grímsvötn.     

南海西南次海盆岩石圈结构及其地质意义

南海地区岩石圈资料稀少,阻碍了其形成演化过程的研究.为此,本次研究结合大地热流、空间重力异常、高程、大地水准面和地震数据,在南海西南次海盆反演了两条2.5维岩石圈剖面.本次计算基于三种假设：岩石圈地幔的密度取决于岩石温度;研究区岩石圈处于热稳定状态;研究区处于重力均衡状态.在剖面A-E中,岩石圈底界面从珠江口盆地的105 km迅速抬升到西沙海槽处的50 km,在西沙海槽、西沙-中沙群岛和西南次海盆变化不大,为50~60 km.在剖面F-I中,岩石圈底界面从西沙群岛-中建地块处的88 km向海盆逐渐抬升,在西南次海盆处为46~50 km,到郑和隆起再逐渐变深至64 km.我们比较了西南次海盆岩石圈的冷却模型和热稳定模型,根据冷却模型由水深和热流数据所推断的西南次海盆年龄比实际年龄差很多,说明冷却模型不适用于西南次海盆.通过对比剖面A-E和剖面F-I,说明了剖面A-E经历了更长时间的拉伸,证明南海西南次海盆在形成演化过程中是从北东向南西逐步打开的渐进式扩张.最后,我们综合分析西南次海盆及其大陆边缘的岩石圈结构、减薄陆壳区范围、碳酸盐台地的分布、下地壳韧性流动、流变结构和沉积层特征等多方面资料,认为西南次海盆在形成演化过程中岩石圈地幔首先破裂而地壳后破裂,属于type Ⅱ型非火山型大陆边缘.     

Crustal velocity structure of central Gansu Province from regional seismic waveform inversion using firework algorithm

The firework algorithm (FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude (M _W) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region, inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15, 15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s, respectively.     

Phenocryst-hosted melt inclusions record stalling of magma during ascent in the conduit and upper magma reservoir prior to vulcanian explosions, Soufrière Hills volcano, Montserrat, West Indies

The mechanics of explosive eruptions influence magma ascent pathways. Vulcanian explosions involve a stop–start mechanism that recurs on various timescales, evacuating the uppermost portions of the conduit. During the repose time between explosions, magma rises from depth and refills the conduit and stalls until the overpressure is sufficient to generate another explosion. We have analyzed major elements, Cl, S, H₂O, and CO₂ in plagioclase-hosted melt inclusions, sampled from pumice erupted during four vulcanian events at Soufrière Hills volcano, Montserrat, to determine melt compositions prior to eruption. Using Fourier transform infrared spectroscopy, we measured values up to 6.7 wt.% H₂O and 80 ppm CO₂. Of 42 melt inclusions, 81 % cluster between 2.8 and 5.4 wt.% H₂O (57 to 173 MPa or 2–7 km), suggesting lower conduit to upper magma reservoir conditions. We propose two models to explain the magmatic conditions prior to eruption. In Model 1, melt inclusions were trapped during crystal growth in magma that was stalled in the lower conduit to upper magma reservoir, and during trapping, the magma was undergoing closed-system degassing with up to 1 wt.% free vapor. This model can explain the melt inclusions with higher H₂O contents since these have sampled the upper parts of the magma reservoir. However, the model cannot explain the melt inclusions with lower H₂O because the timescale for plagioclase crystallization and melt inclusion entrapment is longer than the magma residence time in the conduit. In Model 2, melt inclusions were originally trapped at deeper levels of the magma chamber, but then lost hydrogen by diffusion through the plagioclase host during periodic stalling of the magma in the lower conduit system. In this second scenario, which we favor, the melt inclusions record re-equilibration depths within the lower conduit to upper magma reservoir.     

Seismicity and fluid geochemistry at Lassen Volcanic National Park,California: Evidence for two circulation cells in the hydrothermal system

Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270 °C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240 °C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West, Middle, and East seismic clusters) occur at shallow depths (< 6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤ 161 °C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤ 123 °C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N₂-rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ¹³C–CO₂, δ³⁴S–H₂S, and δ¹⁵N–N₂ values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.