The physical volcanology of the 1600 eruption of Huaynaputina, southern Peru

Volcán Huaynaputina is a group of four vents located at 16°36'S, 70°51'W in southern Peru that produced one of the largest eruptions of historical times when ~11 km³ of magma was erupted during the period 19 February to 6 March 1600. The main eruptive vents are located at 4200 m within an erosion-modified amphitheater of a significantly older stratovolcano. The eruption proceeded in three stages. Stage I was an ~20-h sustained plinian eruption on 19-20 February that produced an extensive dacite pumice fall deposit (magma volume ~2.6 km³). Throughout medial-distal and distal parts of the dispersal area, a fine-grained plinian ashfall unit overlies the pumice fall deposit. This very widespread ash (magma volume ~6.2 km³) has been recognized in Antarctic ice cores. A short period of quiescence allowed local erosion of the uppermost stage-I deposits and was followed by renewed but intermittent explosive activity between 22 and 26 February (stage II). This activity resulted in intercalated pyroclastic flow and pumice fall deposits (~1 km³). The flow deposits are valley confined, whereas associated co-ignimbrite ash fall is found overlying the plinian ash deposit. Following another period of quiescence, vulcanian-type explosions of stage III commenced on 28 February and produced crudely bedded ash, lapilli, and bombs of dense dacite (~1 km³). Activity ceased on 6 March. Compositions erupted are predominantly high-K dacites with a phenocryst assemblage of plagioclase>hornblende>biotite>Fe-Ti oxides-apatite. Major elements are broadly similar in all three stages, but there are a few important differences. Stage-I pumice has less evolved glass compositions (~73% SiO₂), lower crystal contents (17-20%), lower density (1.0-1.3 g/cm³), and phase equilibria suggest higher temperature and volatile contents. Stage-II and stage-III juvenile clasts have more evolved glass (~76% SiO₂) compositions, higher crystal contents (25-35%), higher densities (up to 2.2 g/cm³), and lower temperature and volatile contents. All juvenile clasts show mineralogical evidence for thermal disequilibrium. Inflections on a plot of log thickness vs area^1/2 for the fall deposits suggest that the pumice fall and the plinian ash fall were dispersed under different conditions and may have been derived from different parts of the eruption column system. The ash appears to have been dispersed mainly from the uppermost parts of the umbrella cloud by upper-level winds, whereas the pumice fall may have been derived from the lower parts of the umbrella cloud and vertical part of the eruption column and transported by a lower-altitude wind field. Thickness half distances and clast half distances for the pumice fall deposit suggests a column neutral buoyancy height of 24-32 km and a total column height of 34-46 km. The estimated mass discharge rate for the ~20-h-long stage-I eruption is 2.4᎒⁸ kg/s and the volumetric discharge rate is ~3.6᎒⁵ m³/s. The pumice fall deposit has a dispersal index (Hildreth and Drake 1992) of 4.4, and its index of fragmentation is at least 89%, reflecting the dominant volume of fines produced. Of the 11 km³ total volume of dacite magma erupted in 1600, approximately 85% was evacuated during stage 1. The three main vents range in size from ~70 to ~400 m. Alignment of these vents and a late-stage dyke parallel to the NNW-SSE trend defined by older volcanics suggest that the eruption initiated along a fissure that developed along pre-existing weaknesses. During stage I this fissure evolved into a large flared vent, vent 2, with a diameter of approximately 400 m. This vent was active throughout stage II, at the end of which a dome was emplaced within it. During stage III this dome was eviscerated forming the youngest vent in the group, vent 3. A minor extra-amphitheater vent was produced during the final event of the eruptive sequence. Recharge may have induced magma to rise away from a deep zone of magma generation and storage. Subsequently, vesiculation in the rising magma batch, possibly enhanced by interaction with an ancient hydrothermal system, triggered and fueled the sustained Plinian eruption of stage I. A lower volatile content in the stage-II and stage-III magma led to transitional column behavior and pyroclastic flow generation in stage II. Continued magma uprise led to emplacement of a dome which was subsequently destroyed during stage III. No caldera collapse occurred because no shallow magma chamber developed beneath this volcano.     

Onshore sandbar migration at Tairua Beach (New Zealand): Numerical simulations and field measurements

We observed the onshore migration (3.5 m/day) of a nearshore sandbar at Tairua Beach, New Zealand during 4 days of low-energy wave conditions. The morphological observations, together with concurrent measurements of waves and suspended sediment concentrations, were used to test a coupled, wave-averaged, cross-shore model. Because of the coarse bed material and the relatively low-energy conditions, the contribution of the suspended transport to the total transport was predicted and observed to be negligible. The model predicted the bar to move onshore because of the feedback between near-bed wave skewness, bedload, and the sandbar under weakly to non-breaking conditions at high tide. The predicted bathymetric evolution contrasts, however, with the observations that the bar migrated onshore predominantly at low tide. Also, the model flattened the bar, while in the observations the sandbar retained its steep landward-facing flank. A comparison between available observations and numerical simulations suggests that onshore propagating surf zone bores in very shallow water (< 0.25 m) may have been responsible for most of the observed bar behaviour. These processes are missing from the applied model and, given that the observed conditions can be considered typical of very shallow sandbars, highlight a priority for further field study and model development. The possibility that the excess water transported by the bores across the bar was channelled alongshore to near-by rip-channels further implies that traditional cross-shore measures to judge the applicability of a cross-shore morphodynamic model may be misleading.     

Topodata: Brazilian full coverage refinement of SRTM data

This work presents the selection of a set of geostatistical coefficients suitable for a unified SRTM data refinement from 3″ to 1″ through kriging over the entire Brazilian territory. This selection aimed at data potential for geomorphometric derivations, given by the preservation of detailed geometric characteristics of the resulting digital elevation models (DEM), which are sensitive to refining procedures. The development contained a long-term experimentation stage, when data refinement through kriging was locally developed to support distinct regional projects, followed by a unified selection stage, where the acquired experience was applied to select a single and unified interpolation scheme. In this stage, the selected geostatistical models with promising performances were tested for unified refinement on 40 Brazilian areas with distinct geological settings. Tested areas encompass reliefs varying from mountainous to plain. The effects of data preparation were observed on the perception of patterns (texture and roughness), as well as of singularities (edges, peaks, thalwegs etc.). Results were evaluated mainly through the examination of shaded reliefs, transects and perspectives observed in different scales. Terrains with low slopes and small amplitudes had their DEM promptly affected by the refining methods, as opposed to mountainous terrains. The evaluation, unambiguously confirmed by all consulted interpreters, converged into a refining model with outstanding performance in all tested conditions.     

Accessibility for High Asia： Comparative Perspectives on Northern Pakistan＇s Traffic Infrastructure and Linkages with Its Neighbours in the Hindukush-Karakoram- Himalaya

Starting with a discussion of development concepts which were applied in practice and which followed the developmentalist paradigm the expansion of traffic infrastructure in colonial and post-colonial periods is presented for the High Asian mountain rim. Selective railways and roads are the major feature of this development, which aimed first on serving the convenience of hill station visitors and followed strategic considerations later on. This bias between regional planning and implementation remains a characteristic feature. At the same time traffic infrastructure without asphalt roads is important for the mountain areas, thus breaking up the strong correlation between development and asphalt roads.     

An application of spatial filtering techniques for land cover mapping using TM‐images

Land cover classification using satellite imagery is commonly based on spectral information in the individual pixels. The information in neighbouring pixels is ignored. Spatial filtering techniques using information present in neighbouring pixels may however, contribute significantly to an improvement of the classification. In this study different methods of spatial filtering are applied to a part of a TM‐scene of Kenya to assess their relative reliability. The study area is characterized by extended, relatively homogeneous areas of eucalyptus forests and tea estates and by fragmentated areas of agricultural land use. Spectral information was combined with the results of different spatial filtering methods and then classified. The spatial filtering techniques applied were texture calculation by means of variance, “median minus original” filtering and fractal dimension computations using several sizes of templates. The obtained classification accuracy of several image combinations is compared using the percentage correctly classified and using an overall accuracy measure: the Kappa coefficient. It is concluded that in this case the spatial filtering techniques only slightly improve the classification. From the applied filtering methods texture calculation by means of variance yielded the best results.     

The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

Recent sediment transport and accumulation on the NW Iberian margin

Five transects across the NW Iberian margin were studied in the framework of the EU-funded Ocean Margin EXchange II (OMEX II) project, to determine and establish recent sediment and organic carbon transport and accumulation processes and fluxes.On the Galician shelf and shelf edge, resuspension of sediments resulting in well-developed bottom nepheloid layers was observed at all stations, but transport of suspended sediment appears largely confined to the shelf. On the continental slope, only very dilute bottom nepheloid layers were present, and intermediate nepheloid layers were only occasionally seen. This suggests that cross-slope transfer of particles is limited by the prevailing northerly directed shelf and slope currents.Optical backscatter and ADCP current measurements by the BOBO lander, deployed at 2152 m depth on the Galician slope, indicated that particles in the bottom boundary layer were kept in suspension by tidal currents with highest speeds between 15–25 cm s⁻¹. Net currents during the recording period August 6th–September 10th 1998, were initially directed along-slope toward the NNW, but later turned off-slope toward the SW.The separation of the water masses on the slope from the sediment-laden shelf water by the along-slope current regime is reflected in the recent sedimentary deposits of the Galician shelf and slope. Apart from compositional differences, shelf deposits differ from those on the slope by their higher flux of excess ²¹⁰Pb (0.57–5.37 dpm cm⁻²y⁻¹ versus 0.11–3.00 dpm cm⁻²y⁻¹), a much higher sediment accumulation rate (315.6–2295.9 g m⁻²y⁻¹ versus 10.9–124.7 g m⁻²y⁻¹) and organic carbon burial rate (1.01–34.30 g m⁻²y⁻¹ versus 0.01–0.69 g m⁻²y⁻¹).In contrast to the observations on the Galician margin, pronounced nepheloid layers occurred in the Nazaré Canyon, which extended to considerably greater water depths. This indicates that significantly greater transport of fine-grained particles in both the INL and the BNL was occurring within the canyon, as reflected in the exceptionally high ²¹⁰Pb excess flux (up to 34.09 dpm cm⁻²y⁻¹), mass accumulation rates (maximum 9623.1 g m⁻²y⁻¹) and carbon burial fluxes (up to 180.91 g m⁻²y⁻¹) in the sediment. However, radioisotope fluxes in the lower canyon were only slightly higher than at comparable depths on the Galician margin. This suggests that transport and rapid accumulation is focused on the upper and middle part of the canyon, from where it is episodically released to the deep sea. Compared to the Galician margin, the Nazaré Canyon may be considered as an important organic carbon depocenter on short time-scales, and a major conduit for particulate matter transport to the deep sea on >100 y time-scales.     

西江流域的有机碳侵蚀通量

在西江下游的马口水文站对径流进行了4个季节的有机碳采样分析。研究表明,西江径流有机碳的断面构成在各个季节均较为一致;季节性变化表现为,有机碳和悬浮物含量随流量的增加而增加。随着水体悬浮物含量的增加,悬浮物的有机碳含量呈对数趋势降低。西江流域的有机碳侵蚀通量为10.18×106gC/km2·yr.,是全球外流域有机碳侵蚀通量的2～3倍,其中以颗粒有机碳的侵蚀通量为主,达到8.30×106gC/km2·yr.,溶解有机碳的侵蚀通量为1.88×106gC/km2·yr.。反映了流域内较强的机械剥蚀过程,这与西江流域典型的季风气候、较大的地形高差,以及农业耕作历史长久、土地利用强度较大等因素有关。     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).