Multiple longshore sand bars in the upper Chesapeake Bay

In the upper Chesapeake Bay (Maryland, U.S.A.) field surveys were conducted at 18 multiple longshore sand bar sites. The multiple bar systems were found in water depths less than approximately 2 m (mean sea level), and exhibited mild bottom slopes of 0·0052 or less. The number of bars composing each system ranged from four to 17 and the spacing between the crests typically increased in the offshore direction, ranging from 12 to 70 m. Bar height also typically increased with distance offshore and ranged from 0·03 to 0·61 m. A grain size analysis of crest and trough sediment did not reveal any significant differences and the sediment was categorized as ‘fine sand’. A review of the literature data indicated that the Chesapeake Bay multiple bars possessed similar characteristics to those found in Gelding Bay (Baltic Sea); similarities in fetch, wave height and tidal range between the two bays may account for this finding. The surf-scaling parameter indicated that the multiple bar systems were extremely dissipative with regard to wave energy, and wave height appeared to be an important factor in controlling bar spacing and bar height. A multiple wave break point hypothesis was discussed as a possible mechanism for the formation of Chesapeake Bay multiple longshore bars, and limited observational evidence appeared to support such a mechanism.     

Population growth and the development of a central place system

This paper describes the spatial and functional evolution of a central place system as market conditions change with population growth. Utilizing a partial equilibrium optimization model, we examine the spatial response of two economic sectors to increases in market populations resulting from natural increase and migration. Response in both sectors is conditioned by threshold demand, with factor prices also affecting one of the sectors. As the central place system evolves it exhibits spatial and functional characteristics that are initially consistent with a Löschian landscape, then a Christallerian landscape at higher populations, while at even larger populations Krugman’s landscape emerges.     

Glacial flood pulse effects on benthic fauna in equatorial high‐Andean streams

Equatorial glacier‐fed streams present unique hydraulic patterns when compared to glacier‐fed observed in temperate regions as the main variability in discharge occurs on a daily basis. To assess how benthic fauna respond to these specific hydraulic conditions, we investigated the relationships between flow regime, hydraulic conditions (boundary Reynolds number, Re*), and macroinvertebrate communities (taxon richness and abundance) in a tropical glacier‐fed stream located in the high Ecuadorian Andes (> 4000 m). Both physical and biotic variables were measured under four discharge conditions (base‐flow and glacial flood pulses of various intensities), at 30 random points, in two sites whose hydraulic conditions were representative to those found in other streams of the study catchment. While daily glacial flood pulses significantly increased hydraulic stress in the benthic habitats (appearance of Re* > 2000), low stress areas still persisted even during extreme flood events (Re* < 500). In contrast to previous research in temperate glacier‐fed streams, taxon richness and abundance were not significantly affected by changes in hydraulic conditions induced by daily glacial flood pulses. However, we found that a few rare taxa, in particular rare ones, preferentially occurred in highly stressed hydraulic habitats. Monte‐Carlo simulations of benthic communities under glacial flood reduction scenarios predicted that taxon richness would be significantly reduced by the loss of high hydraulic stress habitats following glacier shrinking. This pioneer study on the relationship between hydraulic conditions and benthic diversity in an equatorial glacial stream evidenced unknown effects of climate change on singular yet endangered aquatic systems. Copyright © 2013 John Wiley & Sons, Ltd.     

改进粒子滤波器的室内地磁匹配精度研究

针对粒子滤波算法在室内地磁分布的连续性变化中会发生滤波发散,加之算法中重采样步骤使该现象恶化,直接导致室内定位过程中目标丢失的问题。该文受Hausdorff距离度量法中点集匹配思想的启发,提出了一种利用初定位误差作为最大粒子匹配距离的约束思想,通过缩小匹配空间的区域范围的方法来解决传统算法的发散问题。最后实验表明,此方法解决了粒子滤波器在地磁匹配中出现的滤波发散现象,同时使得定位结果更加快速收敛。     

Influence of discrete step size on wind setup component of storm surge

Open coast storm surge water levels consist of wind setup due to wind shear at the water surface; a wave setup component caused by wind induced waves transferring momentum to the water column; an atmospheric pressure head component due to the atmospheric pressure deficit over the spatial extent of the storm system; a Coriolis forced setup or setdown component due to the effects of the rotation of the earth acting on the wind driven alongshore current at the coast; a possible seiche component due to resonance effects initiated by moving wind system, and, if astronomical tides are present, an astronomical tide component (although the tide is typically considered to be a forced astronomical event and not really a direct part of the external wind-driven meteorological component of storm surge). Typically the most important component of a storm surge is the wind setup component, especially on the U.S. East Coast and the Gulf of Mexico shorelines. In many approaches to storm surge modeling, a constant depth approximation is invoked over a limited step size in the computational domain. The use of a constant depth approximation has received little attention in the literature although can be very important to the resulting magnitude of the computed storm surge. The importance of discrete step size to the wind setup storm surge component is considered herein with a simple case computation of the wind setup component on a linear slope offshore profile. The present study findings show that the constant depth approximation to wind setup storm surge estimation is biased on the low side (except in extremely shallow water depths) and can provide large errors if discrete step size is not sufficiently resolved. Guidance has been provided on the error that one might encounter for various step sizes on different slopes.     

Propagation of Leaky MHD Waves at Discontinuities with Tilted Magnetic Field

We investigate the characteristics of magneto-acoustic surface waves propagating at a single density interface, in the presence of an inclined magnetic field. For linear wave propagation, the dispersion relation is obtained and analytical solutions are derived for small inclination angle. The inclination of the field renders the frequency of the waves complex, where the imaginary part describes wave attenuation, due to lateral energy leakage.     

Generation of mid-ocean ridge basalts at pressures from 1 to 7 GPa

We propose a model for the generation of average MORBs based on phase relations in the CaO-MgO-Al₂O₃-SiO₂-CO₂ system at pressures from 3 to 7 GPa and in the CaO-MgO-Al₂O₃-SiO₂-Na₂O-FeO (CMASNF) system at pressures from ∼0.9 to 1.5 GPa. The MELT seismic tomography (Forsyth et al., 2000) across the East Pacific Rise shows the largest amount of melt centered at ∼30-km depth and lesser amounts at greater depths. An average mantle adiabat with a model-system potential temperature (T_p) of 1310°C is used that is consistent with this result. In the mantle, additional minor components would lower solidus temperatures ∼50°C, which would lower T_p of the adiabat for average MORBs to ∼1260°C. The model involves generation of carbonatitic melts and melts that are transitional between carbonatite and kimberlite at very small melt fractions (<0.2%) in the low-velocity zone at pressures of ∼2.6 to 7 GPa in the CMAS-CO₂ system, roughly the pressure range of the PREM low-velocity zone. These small-volume, low-viscosity melts are mixed with much larger volumes of basaltic melt generated at the plagioclase-spinel lherzolite transition in the pressure range of ∼0.9 to 1.5 GPa.In this model, solidus phase relations in the pressure range of the plagioclase-spinel lherzolite transition strongly, but not totally, control the major-element characteristics of MORBs. Although the plagioclase-spinel lherzolite transition suppresses isentropic decompression melting in the CMAS system, this effect does not occur in the topologically different and petrologically more realistic CMASNF system. On the basis of the absence of plagioclase from most abyssal peridotites, which are the presumed residues of MORB generation, we calculate melt productivity during polybaric fractional melting in the plagioclase-spinel lherzolite transition interval at exhaustion of plagioclase in the residue. In the CMASN system, these calculations indicate that the total melt productivity is ∼24%, which is adequate to produce the oceanic crust. The residual mineral proportions from this calculation closely match those of average abyssal peridotites.Melts generated in the plagioclase-spinel lherzolite transition are compositionally distinct from all MORB glasses, but do not have a significant fractional crystallization trend controlled by olivine alone. They reach the composition field of erupted MORBs mainly by crystallization of both plagioclase and olivine, with initial crystallization of either one of these phases rapidly joined by the other. This is consistent with phenocryst assemblages and experimental studies of the most primitive MORBs, which do not show an olivine-controlled fractionation trend. The model is most robust for the eastern Pacific, where an adiabat with a T_p of ∼1260°C is supported by the MELT seismic data and where the global inverse correlation of (FeO)₈ with (Na₂O)₈ is weak. Average MORBs worldwide also are well modeled. A heterogeneous mantle consisting of peridotite of varying degrees of major-element depletion combined with phase-equilibrium controls in the plagioclase-spinel lherzolite transition interval would produce the form of the global correlations at a constant T_p, which suggests a modest range of T_p along ridges. Phase-composition data for the CMASNF system are presently not adequate for quantitative calculation of (FeO)₈-(Na₂O)₈-(CaO/Al₂O₃)₈ systematics in terms of this model. The near absence of basalts in the central portion of the Gakkel Ridge suggests a lower bound for T_p along ridges of ∼1240°C, a potential temperature just low enough to miss the solidus for basalt production at ∼0.9 GPa. An upper bound for T_p is poorly constrained, but the complete absence of picritic glasses in Iceland and the global ridge system suggests an upper bound of ∼1400°C. In contrast to some previous models for MORB generation that emphasize large potential temperature variations in a relatively homogeneous peridotitic mantle, our model emphasizes modest potential temperature variations in a peridotitic mantle that shows varying degrees of heterogeneity. Calculations indicate that melt productivity changes from 0 to 24% for a change in T_p from 1240 to 1260°C, effectively producing a rapid increase to full crustal thickness or decrease to none as ridges appear and disappear.     

Origin of the Oceanic Lithosphere

In a global examination of mid-ocean ridge basalt (MORB) glasscompositions, we find that Na₈–Fe₈–depth variationsdo not support modeling of MORBs as aggregates of melt compositionsgenerated over a large range of temperature and pressure. However,the Na₈–Fe₈ variations are consistent with the compositionalsystematics of solidus melts in the plagioclase–spinellherzolite transition in the CaO–MgO–Al₂O₃–SiO₂–Na₂O–FeO(CMASNF) system. For natural compositions, the P–T rangefor melt extraction is estimated to be 1·2–1·5GPa and 1250–1280°C. This P–T range is a closematch with the maximum P–T conditions for explosive pressure-releasevaporization of carbonate-bearing melts. It is proposed thatfracturing of the lithosphere induces explosive formation andescape of CO₂ vapor. This provides the vehicle for extractionof MORBs at a relatively uniform T and P. The upper portionof the CO₂-bearing and slightly melted seismic low-velocityzone flows toward the ridge, rises at the ridge axis to lower-lithospheredepths, melts much more extensively during this rise, and releasesMORB melts to the surface driven by explosively escaping CO₂vapor. The residue and overlying crust produced by this meltingthen migrate away from the ridge axis as new oceanic lithosphere.The entire process of oceanic lithosphere creation involvesonly the upper 140 km. When lithospheric stresses shift fractureformation to other localities, escape of CO₂ ceases, the vehiclefor transporting melt to the surface disappears, and ridgesdie. Inverse correlations of Na₈ vs Fe₈ for MORB glasses areexplained by mantle heterogeneity, and positive variations superimposedon the inverse variations are consistent with progressive extractionof melts from short, ascending melting columns. The uniformlylow temperatures of MORB extraction are not consistent withthe existence of hot plumes on or close to ocean ridges. Inthis modeling, the southern Atlantic mantle from Bouvet to about26°N is relatively homogeneous, whereas the Atlantic mantlenorth of about 26°N shows significant long-range heterogeneity.The mantle between the Charlie Gibbs and Jan Mayen fracturezones is strongly enriched in FeO/MgO, perhaps by a trappedfragment of basaltic crust. Iceland is explained as the productof this enrichment, not a hot plume. The East Pacific Rise,Galapagos Ridge, Gorda Ridge, and Juan de Fuca Ridge samplemantle that is heterogeneous over short distances. The mantlebeneath the Red Sea is enriched in FeO/MgO relative to the mantlebeneath the northern Indian Ocean.     

Automated forecasting of volcanic ash dispersion utilizing Virtual Globes

There are over 100 active volcanoes in the North Pacific (NOPAC) region, most of which are located in sparsely populated areas. Dispersion models play an important role in forecasting the movement of volcanic ash clouds by complementing both remote sensing data and visual observations from the ground and aircraft. Puff is a three-dimensional dispersion model, primarily designed for forecasting volcanic ash dispersion, used by the Alaska Volcano Observatory and other agencies. Since early 2007, the model is in an automated mode to predict the movement of airborne volcanic ash at multiple elevated alert status volcanoes worldwide to provide immediate information when an eruption occurs. Twelve of the predictions are within the NOPAC region, nine more within the southern section of the Pacific ring of fire and the others are in Europe and the Caribbean. Model forecasts are made for initial ash plumes ranging from 4 to 20 km altitude above sea level and for a 24-h forecast period. This information is made available via the Puff model website. Model results can be displayed in Virtual Globes for three-dimensional visualization. Here, we show operational Puff predictions in two and three-dimensions in Google Earth^®, both as iso-surfaces and particles, and study past eruptions to illustrate the capabilities that the Virtual Globes can provide. In addition, we show the opportunity that Google Maps^® provides in displaying Puff operational predictions via an application programming web interface and how radiosonde data (vertical soundings) and numerical weather prediction vertical profiles can be displayed in Virtual Globes for assisting in estimating ash cloud heights.