A laboratory model of surface crust formation and disruption on lava flows through non-uniform channels

Crust formation on basaltic lava flows dictates conditions of both flow cooling and emplacement. For this reason, flow histories are dramatically different depending on whether lava is transported through enclosed lava tubes or through open channels. Recent analog experiments in straight uniform channels (Griffiths et al. J Fluid Mech 496:33–62, 2003) have demonstrated that tube flow, dictated by a stationary surface crust, can be distinguished from a mobile crust regime, where a central solid crust is separated from channel walls by crust-free shear zones, by a simple dimensionless parameter ϑ, such that ϑ<25 produces tube flow and ϑ>25 describes the mobile crust regime. ϑ combines a previously determined parameter ψ, which describes the balance between the formation rate of surface solid and the shear strain that disrupts the solid crust, with the effects of thermal convection (described by the Rayleigh number Ra).Here we explore ways in which ϑ can be used to describe the behavior of basaltic lava channels. To do this we have extended the experimental approach to examine the effects of channel irregularities (expansions, contractions, sinuosity, and bottom roughness) on crust formation and disruption. We find that such changes affect local flow behavior and can thus change channel values of ϑ. For example, gradual widening of a channel results in a decrease in flow velocity that causes a decrease in ϑ and may allow a down-flow transition from the mobile crust to the tube regime. In contrast, narrowing of the channel causes an increase in flow velocity (increasing ϑ), thus inhibiting tube formation.We also quantify the fraction of surface covered by crust in the mobile crust regime. In shallow channels, variations in crust width (d _c) with channel width (W) are predicted to follow d _c∼W ^5/3. Analysis of channelized lava flows in Hawaii shows crustal coverage consistent with this theoretical result along gradually widening or narrowing channel reaches. An additional control on crustal coverage in both laboratory and basaltic flows is disruption of surface crust because of flow acceleration through constrictions, around bends, and over breaks in slope. Crustal breakage increases local rates of cooling and may cause local blockage of the channel, if crusts rotate and jam in narrow channel reaches. Together these observations illustrate the importance of both flow conditions and channel geometry on surface crust development and thus, by extension, on rates and mechanisms of flow cooling. Moreover, we note that this type of analysis could be easily extended through combined use of FLIR and LiDAR imaging to measure crustal coverage and channel geometry directly.Editorial responsibility: A. Harris     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

Thermal structure and heat loss at the summit crater of an active lava dome

Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 10⁴  s⁻¹m⁻²), followed by the hot central core (0.1–0.4 × 10⁴ J s⁻¹m⁻²) and cold annulus (0.04–0.1 × 10⁴ J s⁻¹m⁻²). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s⁻¹), but the cold annulus contributed equal power (2.41–7.07 MJ s⁻¹) as the hot central core (2.68–6.92 MJ s⁻¹) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.     

The changing morphology of an open lava channel on Mt. Etna

An open channel lava flow on Mt. Etna (Sicily) was observed during May 30–31, 2001. Data collected using a forward looking infrared (FLIR) thermal camera and a Minolta-Land Cyclops 300 thermal infrared thermometer showed that the bulk volume flux of lava flowing in the channel varied greatly over time. Cyclic changes in the channel's volumetric flow rate occurred over several hours, with cycle durations of 113–190 min, and discharges peaking at 0.7 m³ s⁻¹ and waning to 0.1 m³ s⁻¹. Each cycle was characterized by a relatively short, high-volume flux phase during which a pulse of lava, with a well-defined flow front, would propagate down-channel, followed by a period of waning flow during which volume flux lowered. Pulses involved lava moving at relatively high velocities (up to 0.29 m s⁻¹) and were related to some change in the flow conditions occurring up-channel, possibly at the vent. They implied either a change in the dense rock effusion rate at the source vent and/or cyclic-variation in the vesicle content of the lava changing its bulk volume flux. Pulses would generally overspill the channel to emplace pāhoehoe overflows. During periods of waning flow, velocities fell to 0.05 m s^–1. Blockages forming during such phases caused lava to back up. Occasionally backup resulted in overflows of slow moving ‘a‘ā that would advance a few tens of meters down the levee flank. Compound levees were thus a symptom of unsteady flow, where overflow levees were emplaced as relatively fast moving pāhoehoe sheets during pulses, and as slow-moving ‘a‘ā units during backup. Small, localized fluctuations in channel volume flux also occurred on timescales of minutes. Volumes of lava backed up behind blockages that formed at constrictions in the channel. Blockage collapse and/or enhanced flow under/around the blockage would then feed short-lived, wave-like, down-channel surges. Real fluctuations in channel volume flux, due to pulses and surges, can lead to significant errors in effusion rate calculations. Editorial responsibility: A. Woods     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

Morphological and textural diversity of the Steens Basalt lava flows,Southeastern Oregon,USA: implications for emplacement style and nature of eruptive episodes

This study focuses on Middle Miocene tholeiitic flood basalt lava flows from the Oregon Plateau, northwestern USA (Steens Basalt), and is the first to comprehensively document and evaluate their morphology. Field observations of flows from several sections within and proximal to the main exposures at Steens Mountain have been supplemented with textural and geochemical data, and are used to offer preliminary insights into their emplacement. Compound pahoehoe flows of variable thickness appear to be common throughout the study area, laterally and vertically. These tend to be plagioclase phyric and the morphology and disposition of constituent flow lobes are quite similar to those from other provinces such as Hawaii and the Snake River Plain. Classic a’a flows with brecciated upper and basal crusts are not abundant, but by no means uncommon. Flows with characters different from typical pahoehoe and a’a are also common. Such flows display a range in morphology; flows with preserved upper crusts but brecciated basal crusts, as well as those displaying well-developed flow-top breccias and preserved basal crusts (rubbly pahoehoe) are observed. The Steens Basalt appears to display greater morphological and textural diversity at the outcrop scale than that described for some other flood basalt provinces. The abundant compound pahoehoe flows (often rich in plagioclase phenocrysts) were likely emplaced during slow but sustained eruptive episodes; their constituent lobes show clear evidence for endogenous growth. The relatively aphyric flows with brecciated surfaces (including a’a) hint at higher strain rates and/or higher viscosity, probably caused by higher effusion rates. A couple of sections are characterized by compositionally similar, but morphologically different flows that were possibly part of the same eruption. While differences in pre-eruptive topography could explain this, it is also possible that certain physical parameters changed substantially and abruptly during eruption and that such changes were accompanied by differentiation processes within the plumbing system. It is possible that such observations indicate temporal fluctuations within complex magmatic and eruptive systems, and deserve closer scrutiny.     

Strombolian explosive styles and source conditions: insights from thermal (FLIR) video

Forward Looking Infrared Radiometer (FLIR) cameras offer a unique view of explosive volcanism by providing an image of calibrated temperatures. In this study, 344 eruptive events at Stromboli volcano, Italy, were imaged in 2001–2004 with a FLIR camera operating at up to 30 Hz. The FLIR was effective at revealing both ash plumes and coarse ballistic scoria, and a wide range of eruption styles was recorded. Eruptions at Stromboli can generally be classified into two groups: Type 1 eruptions, which are dominated by coarse ballistic particles, and Type 2 eruptions, which consist of an optically-thick, ash-rich plume, with (Type 2a) or without (Type 2b) large numbers of ballistic particles. Furthermore, Type 2a plumes exhibited gas thrust velocities (>15 m s⁻¹) while Type 2b plumes were limited to buoyant velocities (<15 m s⁻¹) above the crater rim. A given vent would normally maintain a particular gross eruption style (Type 1 vs. 2) for days to weeks, indicating stability of the uppermost conduit on these timescales. Velocities at the crater rim had a range of 3–101 m s⁻¹, with an overall mean value of 24 m s⁻¹. Mean crater rim velocities by eruption style were: Type 1 = 34 m s⁻¹, Type 2a = 31 m s⁻¹, Type 2b = 7 m s⁻¹. Eruption durations had a range of 6–41 s, with a mean of 15 s, similar among eruption styles. The ash in Type 2 eruptions originates from either backfilled material (crater wall slumping or ejecta rollback) or rheological changes in the uppermost magma column. Type 2a and 2b behaviors are shown to be a function of the overpressure of the bursting slug. In general, our imaging data support a broadening of the current paradigm for strombolian behavior, incorporating an uppermost conduit that can be more variable than is commonly considered.     

Distribution of thermal areas on an active lava flow field: Landsat observations of Kilauea,Hawaii, July 1991

A Landsat Thematic Mapper (TM) image acquired on 23 July 1991 recorded widespread activity associated with the Episode 48 of the Pu'u 'O'o-Kupaianaha eruption of Kilauea Volcano, Hawaii. The scene contains a very large number (>3500) of thermally elevated near infrared (0.8–2.35 m) pixels (each 900 m²), which enable the spatial distribution of volcanic activity to be identified. This activity includes a lava lake within Pu'u 'O'o cone, an active lava tube system (7.9 km in length) with skylights between the Kupaianaha lava shield and several ocean entry points, and extensive active surface flows (total area of 1.3 km²) within a much larger area of cooling flows (total16 km²). The production of an average flux density map from the TM data of the flow field, wherein the average flux density is defined in units of Wm^-2, allows for the chronology of emplacement of active and cooling flows to be determined. The flux density map reveals that there were at least three breakouts (>5000 Wm^-2) feeding active flows, but on the day that the data were collected the TM recorded a waning phase of surface activity in this area, based on the relatively large amount of intermediate power-emitting (cooling) flows compared to high power-emitting (active) flows. The production of a comparable flux density map for future eruptions would aid in the assessment of volcanic hazards if the data were available in near-real time.     

Voluminous lava flows at Oldoinyo Lengai in 2006: chronology of events and insights into the shallow magmatic system

The largest natrocarbonatite lava flow eruption ever documented at Oldoinyo Lengai, NW Tanzania, occurred from March 25 to April 5, 2006, in two main phases. It was associated with hornito collapse, rapid extrusion of lava covering a third of the crater and emplacement of a 3-km long compound rubbly pahoehoe to blocky aa-like flow on the W flank. The eruption was followed by rapid enlargement of a pit crater. The erupted natrocarbonatite lava has high silica content (3% SiO₂). The eruption chronology is reconstructed from eyewitness and news media reports and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data, which provide the most reliable evidence to constrain the eruption’s onset and variations in activity. The eruption products were mapped in the field and the total erupted lava volume estimated at 9.2 ± 3.0 × 10⁵ m³. The event chronology and field evidence are consistent with vent construct instability causing magma mixing and rapid extrusion from shallow reservoirs. It provides new insights into and highlights the evolution of the shallow magmatic system at this unique natrocarbonatite volcano.     

Comparison of the response of the Landsat 7 Enhanced Thematic Mapper Plus and the Earth Observing-1 Advanced Land Imager over active lava flows

Using near simultaneously acquired Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Earth Observing-1 Advanced Land Imager (ALI) data we assess the relative radiant responses over active lava flows from the Mt. Etna July/August 2001 flank eruption. By assessing the extent of saturation between the two instruments and using the dual-band method of extracting sub pixel thermal information, we show that the ALI represents an improvement over the ETM+ in the present ability to assess temperatures of hot active lava flows for a number of reasons. (1) The extra spectral channels provided by ALI compliment the current SWIR channels on ETM+ by providing a greater number of paired channel combinations for input into the dual-band method. Thus, dual-band temperature solutions can be determined for a greater range of lava flow types than previously possible using the two paired channel combinations available with the ETM+. (2) The ALI instrument is less susceptible than ETM+ to saturation within the SWIR, especially when using channels 5, 5p and 4p at wavelengths of 1.65, 1.25 and 0.87 μm respectively. (3) The greater radiometric sensitivity of the ALI 12 bit electronics coupled with a significantly higher signal to noise ratio aid in obtaining successful dual-band solutions.     

Volcanic structure of the southern sector of Mt. Etna after the 2001 and 2002 eruptions defined by magnetotelluric measurements

Two magnetotelluric (MT) surveys were carried out on the Mt. Etna volcano after two of the most intense eruptions of the last 30 years which took place in summer 2001 and winter 2002–2003. Surveying was pursued for two main reasons. First, we sought to contribute to the definition of the first-order structure and physico-chemical state (temperature, fluids, melts) of a volcano that has been extensively explored and monitored by means of various geophysical methods, but where only few electrical and electromagnetic surveys have been performed. Secondly, we acquired MT data in the same sites in the two different surveys with the aim of monitoring the possible changes of the first-order structure, since conditions are expected to vary on an active volcano such as Etna, and are supposed to be linked to the eruptive events. Soundings have been acquired in an E-W 10 km-long profile across the southern flank of Mt. Etna, at a distance of almost 6 km south from the Central Crater. The first survey was carried out three months after the 2001 eruption. Inverse models define a pronounced (4 km thickness) low resistivity section at a depth of about 1 km b.s.l. to the west. To the east, a low resistivity section is still present, but appears deeper, thinner and more resistive, and a shallow low resistivity anomaly also exists. The shallow anomaly to the east is tentatively correlated with altered and clayey volcanic units and/or temporary groundwater storage. The deep anomalies are interpreted as being due to melt storage at shallow depths which was not exhausted during the eruption. This would be confirmed by the abundance of lava erupted within one year from the end of the survey. The few good sites retrieved in the second survey, carried out a few weeks after the eruption of 2002–2003, confirm the picture defined in the first survey, and provide a better definition of the bottom of the deep anomaly located in the sedimentary basement.     

Historical measurements of the Earth's magnetic field compared with remanence directions from lava flows in Italy over the last four centuries

Direct measurements of the Earth's magnetic field in Italy since 1640 a.d. have been used to check the remanence directions derived from historically dated volcanic rocks of Etna and Vesuvius. Direct measurements consist of the records of L’Aquila and Pola geomagnetic observatories, the repeat stations of the Italian Magnetic Network and the data base of the Historical Italian Geomagnetic Data Catalogue. All have been relocated to the same reference site (Viterbo — lat. 42.45°N, long. 12.03°E) in order to draw a reference secular variation (SV) curve. The direction of the Earth's field at Viterbo has also been calculated from the historical records (2000-1600) of ref. [Jackson, A., Jonkers, A.R.T., Walker, M.R., 2000. Four centuries of geomagnetic secular variation from historical records. Phil. Trans. R. Soc. London, Ser. A 358, 957-990] database. The remanence directions from Etna show a general agreement with the trend of the SV curve, although their inclination is usually lower than that from the direct measurement. The directions from Vesuvius are more scattered. Large discrepancies occur at both volcanoes and in some cases have been ascribed in the literature to poor geographic information, making it difficult to identify the flows actually emplaced during the eruptions reported in the chronicles. Closer examination shows that the great majority of the best-defined remanence directions (semi-angle of confidence α₉₅ < 2.5°) deviate significantly from the geomagnetic direction measured at the time of the emplacement, the angle between the two directions being larger than the α₉₅ value. The value of 2.5-3.0° can thus be regarded as a conservative evaluation of the error when dealing with dating Etna and Vesuvius lava flows older than 17th century, even when the accuracy attained in remanence measurements is higher. In default of a SV curve for Italy derived from archaeological artefacts, a further error in dating is introduced when reference is made to SV curves of other countries, even if well-established, as these are from regions too far from Italy (>600 km) to confidently relocate magnetic directions.     

Morphologies and emplacement mechanisms of the lava flows of the Faroe Islands Basalt Group,Faroe Islands,NE Atlantic Ocean

The Palaeogene Faroe Islands Basalt Group (FIBG) comprises three eruptive sequences or formations, all emplaced into a subaerial environment during the development of the extensive continental flood basalt province that stretches from East Greenland through the Faroe Islands and into the Faroe-Shetland Basin. The Beinisvørð Formation, having a tabular-classic facies architecture, is composed of a sequence of simple flows each comprising a single sheet lobe. The Beinisvørð Formation is overlain by the distinctly contrasting Malinstindur Formation that has a compound-braided facies architecture. The Enni Formation occurs at the top of the sequence and consists of a mixture of simple and compound flows with tabular-classic and compound-braided facies architectures, respectively. Surface and internal characteristics of the sheet lobes of the Beinisvørð and Enni formations indicate emplacement through inflation, which is more obvious for the tube-fed compound flows of the Malinstindur and Enni formations. The difference between the simple and compound flow sequences of the FIBG is, most likely, linked to the manner in which the lava was supplied during the eruption and the eruptive style of the volcanic system. The sheet lobes were erupted over laterally extensive areas from fissure systems which had a continuous supply of lava, which contrasts with the tube-fed compound flows which were erupted in a gradual, piecemeal manner from point-sourced, low shield volcanoes with limited areal extents.     

Effects of accelerating and decelerating flows in a channel with vegetated banks and gravel bed

Experiments were conducted in an 8 m long,40 cm wide,and 60 cm deep re-circulating fiume with vegetated banks and gravel bed to study the effects of accelerating and decelerating flows on the flow stru...     

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows on infiltrating surfaces with lateral inflow neglected in the momentum equation

Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow, infiltration and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.     

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows with lateral inflow neglected in the momentum equation

Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; and it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.     

Identification of appropriate lags and temporal resolutions for low flow indicators in the River Rhine to forecast low flows with different lead times

The aim of this paper is to assess the relative importance of low flow indicators for the River Rhine and to identify their appropriate temporal lag and resolution. This is done in the context of low flow forecasting with lead times of 14 and 90 days. First, the Rhine basin is subdivided into seven sub‐basins. By considering the dominant processes in the sub‐basins, five low flow indicators were selected: precipitation, potential evapotranspiration, groundwater storage, snow storage and lake storage. Correlation analysis was then carried out to determine the relationship between observed low flows and preselected indicators with varying lags (days) and temporal resolutions (from 1 day to 7 months). The results show that the most important low flow indicators in the Alpine sub‐basins for forecasts with a lead time of 14 days are potential evapotranspiration with a large lag and temporal resolution, and lake levels with a small lag and temporal resolution. In the other sub‐basins groundwater levels with a small lag and temporal resolution are important in addition to potential evapotranspiration with a large lag and temporal resolution. The picture is slightly different for forecasts with a lead time of 90 days. The snow storage in the Alpine sub‐basins and the precipitation in the other sub‐basins also become relevant for low flows. Consequently, the most important low flow indicators in the Alpine sub‐basins for forecasts with a lead time of 90 days are potential evapotranspiration with a large lag and temporal resolution, lake levels with a small lag and temporal resolution and snow storage with a small lag and large temporal resolution. The resultant correlation maps provide appropriate lags and temporal resolutions for indicators to forecast low flows in the River Rhine with different lead times. Copyright © 2012 John Wiley & Sons, Ltd.     

珠江口盆地大地热流特征及其与热岩石圈厚度的关系

沉积盆地现今热流特征是岩石圈构造-热演化过程的综合反映和盆地热史恢复的必要约束条件,其总体变化趋势与热岩石圈厚度密切相关.本文根据新收集的珠江口盆地19口钻井温度数据,新增计算了19个大地热流数据,其中12个数据位于深水区（水深大于300 m）,丰富了该盆地深水区钻井地热数据.结合前人研究成果,绘制了该盆地的大地热流图,并分析了其热流分布特征.在此基础上通过求解一维热传导方程,计算得到36口井位处的热岩石圈厚度,量化了盆地大地热流与热岩石圈厚度间的关系.结果显示,珠江口盆地大地热流值介于24.2~121.0 mW·m-2,平均71.8±13.6 mW·m-2,新增盆地深水区钻井平均热流值高达84.5±4.4 mW·m-2.大地热流分布整体上从陆架区到陆坡区升高,而热岩石圈厚度整体分布趋势与大地热流相反.大地热流与热岩石圈厚度间存在良好的指数相关性.