Petrology of the Proterozoic Potato River Layered Intrusion, Northern Wisconsin, USA

The Potato River intrusion is a Keweenawan (1100 Ma) mafic plutonemplaced in Keweenawan volcanics and earlier Proterozoic metasedimentaryrocks along the southeastern flank of the Lake Superior syncline.It comprises the following lithostratigraphic zones: a thinto absent Border zone of altered olivine gabbro; a Lower zoneof olivine gabbro; a Picritic zone of picrite and troctolite;a Middle zone of olivine gabbro and leucogabbro; an Upper zoneof quartz leucogabbro and ferrogabbro; and a Roof zone of granophyricand granitic rocks. Fractional crystallization is evident fromcompositional changes in the rocks and cumulus minerals withstratigraphic height. Elements concentrated in the cumulus mineralsolivine and plagioclase (Mg, Fe²⁺, Al, Ca, Ni, Co, Cr, Sr) decreasewith height; elements concentrated in the trapped liquid (Na,K, La, Y, Zr, Nb, Rb, Ba) increase with height; and other elements(Ti, Fe³⁺, P, Ga, V, Sc, Cu, Zn) show complicated behavior relatedto the appearance of additional cumulus phases such as clinopyroxene,Fe-Ti oxides, and apatite. Lower zone rocks contain some sulfide,probably from sulfur derived from the country rock, and theUpper zone has sulfides probably precipitated from an immisciblesulfide liquid. The sulfide-bearing rocks have similaritiesto those of other intrusions, such as Bushveld, Stillwater,and Skaergaard. The picritic and troctolitic rocks of the Picritic zone indicatethat the intrusion was open to additional injections of maficmagma. Roof zone granophyric rocks are residual liquids intrudedalong the upper margin of the intrusion during regional tilting,but Roof zone granitic rocks are probably melted country rock.An attempt is made to estimate by reverse stratigraphic summationthe compositional path of the magma that solidified above thePicritic zone. The first compositions are highly aluminous,which suggests that the upper part of the intrusion has beenenriched in plagioclase by convection-aided crystal sorting.A complementary unit of mafic rocks is not exposed, but it couldbe present down dip. Some of the later compositions are similarto typical Keweenawan high-Al tholeiites. The magma did notundergo extreme iron enrichment, probably because of oxygenfugacity buffering.     

Changes in American urban floodplain occupancy since 1958: the experiences of nine cities

In 1958, researchers from the University of Chicago documented increases in the number of structures in the floodplains of 17 American cities, a phenomenon attributed to the prevalence of flood control structures. Because federal policies have shifted to managing floodplains, this paper updates the 1958 study for nine of the cities. Several trends are apparent. Development pressures determine flood-plain encroachment; floodplain management regulations have been implemented where it was easiest to do so. Further, structural measures continue to dominate. Thus, while many local officials are becoming more aware of advantages of flood-plain regulation, implementation and enforcement are inconsistent and uneven.     

Reconstruction of palaeosalinity using carbon isotopes and benthic associations: a comparison

Carbon and oxygen isotopes were determined on 40 recrystallized shells of Late Jurassic bivalves from the Lusitanian Basin of Portugal. In contrast with the oxygen isotopes, which exhibited considerable diagenetic distortion, the carbon isotopes are thought to preserve a record of the salinity of the Jurassic marginal marine seas in which these bivalves lived. The reconstructed palaeosalinities range from 35%o (euhaline) to 5% (oligohaline). Comparing these values with the palaeosalinity reconstructed from a palaeoecological analysis of 17 stratigraphic levels within the basin, the independently derived values agree in most cases. Strongly differing values are explained as being due to biotic factors and to diagenetic distortion of the isotopic signal; they are less likely to be due to smallscale time-averaging or insufficient microstratigraphic sampling. On the whole, the carbon isotope analyses are thought to produce reasonable palaeosalinity values, although data from infaunal, originally aragonitic bivalves appear to be less reliable than those from epifaunal bivalves with a predominantly or exclusively calcitic shell. As diagenetic alteration of the carbon isotope signal is, however, unpredictable and biotic effects on the isotopic composition are insufficiently known, palaeosalinity reconstructions based on stable isotope data should be supported by palaeoecological data.     

Improved Petrological Modal Analyses from X-ray Powder Diffraction Data by use of the Rietveld Method I. Selected Igneous, Volcanic, and Metamorphic Rocks

Mineral modes have been determined for specimens of eight rocktypes from CuK X-ray powder diffraction data using the Rietveldmethod. The samples include two granites, a granodiorite, adamellite,gabbro, basalt, trachyte, and two granulite-facies metamorphicrocks. Up to eight individual mineral components have been measuredin each sample (no glassy phases were observed), with a detectionlimit of {small tilde}1 wt.%, depending on the mineral assemblage.Marked variations in grain size (i.e., granite vs. trachyte)provide no difficulties for the X-ray method. The X-ray resultscompare very favourably with (1) optical modes determined forthe medium–coarse-grained samples by point counting, (2)normative calculations obtained using locally enhanced catanormand mesonorm software, and (3) corresponding Rietveld modesdetermined, for two samples, from neutron powder data. Wheredifferences occur, these are discussed in relation to the limitationsof each of the methods. The improved accuracy of the X-ray method is due primarily tothe incorporation of the full diffraction profile in the Rietveldanalysis calculations, and the elimination of preferred orientationby collecting the data from samples packed in glass capillaries(i.e., Debye–Scherrer mode). The good agreement of theX-ray and neutron modes shows that the usual problems encounteredwith microabsorption, extinction, and sampling are of littleconcern in these rocks. The results highlight one of the majoradvantages provided by Rietveld modal analysis over the moretraditional ‘reference intensity’ X-ray methods,namely, that the crystal chemistry (and thus the calibrationconstants) of the individual phases can be adjusted dynamicallyduring each individual analysis. This not only provides moreaccurate phase abundances, but also gives important supplementaryinformation about the mineralogy of the major components.     

A Model of Magmatic Crystallization

A computer model simulating fractional crystallization at oneatmosphere pressure incorporates nine broadly-defined minerals—magnetite,olivine, hypersthene, augite, quartz, plagioclase, orthoclase,leucite, and nepheline. The crystallization temperature of eachmineral is considered to be a smooth function of the compositionof the magmatic liquid. These mineral temperature equationsare obtained by multiple linear regression analysis of informationfrom published silicate systems and rock melting experiments.The nine equations are solved for any primary liquid, withinthe broad range of common magma types, to select the crystallizingmineral or minerals. Partition ratios from published experimentsand analyses of lavas and phenocrysts permit calculation ofthe composition of the crystallizing mineral assemblage. Subtractionof a small amount of that composition from the primary liquidyields a new liquid, which may be recycled to yield a sequenceof liquids during fractional crystallization. The crystallizationmodel handles assemblages of co-precipitating minerals, andcan trace progressive saturation in new minerals, substitutionof a new mineral for an old mineral, and cessation of crystallizationof a mineral. The sequences of minerals and liquids derivedfrom a broad set of primary liquids are geologically realistic,so the model is useful in predicting phenocrysts in volcanicrocks and events during crystallization of shallow intrusions.     

Metadolerite age relationships in the Dalradian of northwest Donegal,Ireland and their orogenic significance

The structural relationships of several intrusive sheets belonging to the extensive metadolerite suite of the Appin (lower) Dalradian rocks of the Creeslough area, northwest Donegal, Ireland have been examined. It is concluded that, relative to the regional chronology of deformation of the enclosing metasediments, sheet emplacement has occurred: pre-D2; syn-D2; between D2 and D3; and syn-D3. The height of activity came in the post-D2, pre-D3 interval and coincided with the peak of regional metamorphism in the area. These relationships in Donegal are compared with those of the pretectonic basic rocks of the southwest Highlands and models are discussed.     

Island shadow effects and the wave climate of the Western Tuamotu Archipelago (French Polynesia) inferred from altimetry and numerical model data

To implement a numerical model of atoll lagoon circulation, we characterized first the significant wave height (Hs) regime of the Western Tuamotu Archipelago and the local attenuation due to the protection offered by large atolls in the south Tuamotu. Altimetry satellite data and a WAVEWATCH III two-way nested wave model at 5km resolution from 2000 to 2010 were used. Correlation between altimetry and model was high (0.88) over the period. According to the wave model, the archipelago inner seas experienced attenuated Hs year-long with a yearly average Hs around 1.3m vs a minimum of 1.6m elsewhere. The island shadow effect is especially significant in the austral winter. In contrast with southern atolls, Western Tuamotu experienced only few days per year of Hs larger than 2.5m generated by very high Hs southern swell, transient western local storms, strong easterly winds, and during the passage of distant hurricanes.     

Numerical simulation of undersea cable dynamics

A fully three-dimensional code has been written to compute the motion of a towed cable. The code is based on a robust and stable finite difference approximation to the differential equations derived from basic dynamics. A 3500-ft (1.07 km) cable pulled at 18.5 knots (34.3 km hr⁻¹) through a circular turn of 700 yd (0.64 km) radius has been computed in about half of the real time of the maneuver. The computed displacements are close to the measured ones; the changes in depth are within 2%.