Stability and dynamics of the continental tectosphere

Continental cratons overlie thick, high-viscosity, thermal and chemical boundary layers, where the chemical boundary layers are less dense than they would be due to thermal effects alone, perhaps because they are depleted in basaltic constituents. If the continental tectosphere is the same age as the overlying Archaean crust, then the continental tectosphere must be able to survive for several billion years without undergoing a convective instability, despite being both cold and thick. Since platforms and shields correlate only weakly with Earth's gravity and geoid anomalies, acceptable models of the continental tectosphere must also satisfy this gravity constraint. We investigate the long-term stability of the continental tectosphere by carrying out a number of numerical convection experiments within a two-dimensional Cartesian domain. We initiate our experiments with a tectosphere (thermal and chemical boundary layers) immersed in a region of uniform composition, temperature, and viscosity, and consider the effects on the stability of the tectosphere of (1) activation energy (used to define the temperature dependence of viscosity), (2) compositional buoyancy, and (3) linear or non-linear rheology. The large lateral thermal gradients required to match oceanic and tectosphere structures initiate the dominant instability, a “drip” which develops at the side of the tectosphere and moves to beneath its center. High activation energies and high background viscosities restrict the amount and rate of entrainment. Compositional buoyancy does not significantly change the flow pattern. Rather, compositional buoyancy slows the destruction process somewhat and reduces the stress within the tectosphere. With a non-Newtonian rheology, this reduction in stress helps to stiffen the tectosphere. In these experiments, dynamical systems that adequately model the present ocean-continent structures have activation energy E*≥180 kJ mole⁻¹ — a value about one third the estimate of activation energy for olivine, E*≈520 kJ mole⁻¹. Although for E*≈520 kJ mole⁻¹, compositional buoyancy is not required for the tectosphere to survive, the joint application of longevity and gravity constraints allows us to reject all models not containing compositional buoyancy, and to predict that the ratio of compositional to thermal buoyancy within the continental tectosphere is approximately unity.     

Chromian spinel as a petrogenetic indicator in abyssal and alpine-type peridotites and spatially associated lavas

The composition of chromian spinel in alpine-type peridotites has a large reciprocal range of Cr and Al, with increasing Cr# (Cr/(Cr+Al)) reflecting increasing degrees of partial melting in the mantle. Using spinel compositions, alpine-type peridotites can be divided into three groups. Type I peridotites and associated volcanic rocks contain spinels with Cr#<0.60; Type III peridotites and associated volcanics contain spinels with Cr#>0.60, and Type II peridotites and volcanics are a transitional group and contain spinels spanning the full range of spinel compositions in Type I and Type II peridotites. Spinels in abyssal peridotites lie entirely within the Type I spinel field, making ophiolites with Type I alpine-type peridotites the most likely candidates for sections of ocean lithosphere formed at a midocean ridge. The only modern analogs for Type III peridotites and associated volcanic rocks are found in arc-related volcanic and intrusive rocks, continental intrusive assemblages, and oceanic plateau basalts. We infer a sub-volcanic arc petrogenesis for most Type III alpine-type peridotites. Type II alpine-type peridotites apparently reflect composite origins, such as the formation of an island-arc on ocean crust, resulting in large variations in the degree and provenance of melting over relatively short distances. The essential difference between Type I and Type III peridotites appears to be the presence or absence of diopside in the residue at the end of melting.Based on an examination of co-existing rock and spinel compositions in lavas, it appears that spinel is a sensitive indicator of melt composition and pressure of crystallization. The close similarity of spinel composition fields in genetically related basalts, dunites and peridotites at localities in the oceans and in ophiolite complexes indicates that its composition reflects the degree of melting in the mantle source region. Accordingly, we infer from the restricted range of spinel compositions in abyssal basalts that the degree of mantle melting beneath mid-ocean ridges is generally limited to that found in Type I alpine-type peridotites. It is apparent, therefore, that the phase boundary OL-EN-DI-SP +meltOL-EN-SP+melt has limited the degree of melting of the mantle beneath mid-ocean ridges. This was clearly not the case for many alpine-type peridotites, implying very different melting conditions in the mantle, probably involving the presence of water.     

Andesite and dacite genesis via contrasting processes: the geology and geochemistry of El Valle Volcano,Panama

The easternmost stratovolcano along the Central American arc is El Valle volcano, Panama. Several andesitic and dacitic lava flows, which range in age 5–10 Ma, are termed the old group. After a long period of quiescence (approximately 3.4 Ma), volcanic activity resumed approximately 1.55 Ma with the emplacement of dacitic domes and the deposition of dacitic pyroclastic flows 0.9–0.2 Ma. These are referred to as the young group. All of the samples analyzed are calc-alkaline andesites and dacites. The mineralogy of the two groups is distinct; two pyroxenes occur in the old-group rocks but are commonly absent in the young group. In contrast, amphibole has been found only in the young-group samples. Several disequilibrium features have been observed in the minerals (e.g., oscillatory zoning within clinopyroxenes). These disequilibrium textures appear to be more prevalent among the old- as compared with the young-group samples and are most likely the result of magma-mixing, assimilation, and/or polybaric crystallization. Mass-balance fractionation models for major and trace elements were successful in relating samples from the old group but failed to show a relationship among the young-group rocks or between the old- and young-group volcanics. We believe that the old-group volcanics were derived through differentiation processes from basaltic magmas generated within the mantlewedge. The young group, however, does not appear to be related to more primitive magmas by differentiation. The young-group samples cannot be related by fractionation including realistic amounts of amphibole. Distinctive geochemical features of the young group, including La/Yb ratios〉15, Yb〈1, Sr/Y〉150, and Y〈6, suggest that these rocks were derived from the partial melting of the subducted lithosphere. These characteristics can be explained by the partial melting of a source with residual garnet and amphibole. Dacitic material with the geochemical characteristics of subducted-lithosphere melting is generated apparently only where relatively hot crust is subducted, based on recent work. The young dacite-genesis at El Valle volcano is related to the subduction of relatively hot lithosphere.     

Magmatic-to-hydrothermal crystallization in the W–Sn mineralized Mole Granite (NSW, Australia): Part I: Crystallization of zircon and REE-phosphates over three million years—a geochemical and U–Pb geochronological study

Zircon, monazite and xenotime crystallized over a temperature interval of several hundred degrees at the magmatic to hydrothermal transition of the Sn and W mineralized Mole Granite. Magmatic zircon and monazite, thought to have crystallized from hydrous silicate melt, were dated by conventional U–Pb techniques at an age of 247.6 ± 0.4 and 247.7 ± 0.5 Ma, respectively. Xenotime occurring in hydrothermal quartz is found to be significantly younger at 246.2 ± 0.5 Ma and is interpreted to represent hydrothermal growth. From associated fluid inclusions it is concluded that it precipitated from a hydrothermal brine ≤ 600 °C, which is below the accepted closure temperature for U–Pb in this mineral. These data are compatible with a two-stage crystallization process: precipitation of zircon and monazite as magmatic liquidus phases in deep crustal magma followed by complete crystallization and intimately associated Sn–W mineralization after intrusion of the shallow, sill-like body of the Mole Granite. Later hydrothermal formation of monazite in a biotite–fluorite–topaz reaction rim around a mineralized vein was dated at 244.4 ± 1.4 Ma, which distinctly postdates the Mole Granite and is possibly related to a younger hidden intrusion and its hydrothermal fluid system.

Obtaining precise age data for magmatic and hydrothermal minerals of the Mole Granite is hampered by uncertainties introduced by different corrections required for multiple highly radiogenic minerals crystallising from evolved hydrous granites, including ²³⁰Th disequilibrium due to Th/U fractionation during monazite and possibly xenotime crystallization, variable Th/U ratios of the fluids from which xenotime was precipitating, elevated contents of common lead, and post-crystallization lead loss in zircon, enhanced by the fluid-saturated environment. The data imply that monazite can also survive as a liquidus phase in protracted magmatic systems over periods of 10⁶ years. The outlined model is in agreement with prominent chemical core-rim variation of the zircon.     

Analysis of the surface temperature and wind forecast errors of the NCAR-AirDat operational CONUS 4-km WRF forecasting system

Investigating the characteristics of model-forecast errors using various statistical and object-oriented methods is necessary for providing useful guidance to end-users and model developers as well. To this end, the random and systematic errors (i.e., biases) of the 2-m temperature and 10-m wind predictions of the NCAR-AirDat weather research and forecasting (WRF)-based real-time four-dimensional data assimilation (RTFDDA) and forecasting system are analyzed. This system has been running operationally over a contiguous United States (CONUS) domain at a 4-km grid spacing with four forecast cycles daily from June 2009 to September 2010. In the result an exceptionally useful forecast dataset was generated and used for studying the error properties of the model forecasts, in terms of both a longer time period and a broader coverage of geographic regions than previously studied. Spatiotemporal characteristics of the errors are investigated based on the 24-h forecasts between June 2009 and April 2010, and the 72-h forecasts between May and September 2010. It was found that the biases of both wind and temperature forecasts vary greatly seasonally and diurnally, with dependency on the forecast length, station elevation, geographical location, and meteorological conditions. The temperature showed systematic cold biases during the daytime at all station elevations and warm biases during the nighttime above 1,000 m above sea level (ASL), while below 600 m ASL cold biases occurred during the nighttime. The forecasts of surface wind speed exhibited strong positive biases during the nighttime, while the negative biases were observed in the spring and summer afternoons. The surface wind speed was mostly over-predicted except for the stations located between 1,000 and 2,100 m ASL, for which negative biases were identified for most forecast cycles. The highest wind-speed errors were found over the high terrain and near sea-level stations. The wind-direction errors were relatively large at the high-terrain elevation in the Rocky and Appalachian mountain ranges and the western coastal areas and the error structure exhibited notable diurnal variability.     

Summer temperature variability across four urban neighborhoods in Knoxville,Tennessee, USA

The urban heat island (UHI) is a well-documented effect of urbanization on local climate, identified by higher temperatures compared to surrounding areas, especially at night and during the warm season. The details of a UHI are city-specific, and microclimates may even exist within a given city. Thus, investigating the spatiotemporal variability of a city’s UHI is an ongoing and critical research need. We deploy ten weather stations across Knoxville, Tennessee, to analyze the city’s UHI and its differential impacts across urban neighborhoods: two each in four neighborhoods, one in more dense tree cover and one in less dense tree cover, and one each in downtown Knoxville and Ijams Nature Center that serve as control locations. Three months of temperature data (beginning 2 July 2014) are analyzed using paired-sample t tests and a three-way analysis of variance. Major findings include the following: (1) Within a given neighborhood, tree cover helps negate daytime heat (resulting in up to 1.19 ^°C lower maximum temperature), but does not have as large of an influence on minimum temperature; (2) largest temperature differences between neighborhoods occur during the day (0.38–1.16 ^°C difference), but larger differences between neighborhoods and the downtown control occur at night (1.04–1.88 ^°C difference); (3) presiding weather (i.e., air mass type) has a significant, consistent impact on the temperature in a given city, and lacks the differential impacts found at a larger-scale in previous studies; (4) distance from city center does not impact temperature as much as land use factors. This is a preliminary step towards informing local planning with a scientific understanding of how mitigation strategies may help minimize the UHI and reduce the effects of extreme weather on public health and well-being.     

Estimating Sea-Level Allowances for Atlantic Canada using the Fifth Assessment Report of the IPCC

Abstract

Sea-level allowances at 22 tide-gauge sites along the east coast of Canada are determined based on projections of regional sea-level rise for the Representative Concentration Pathway 8.5 (RCP8.5) from the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5) and on the statistics of historical tides and storm surges (storm tides). The allowances, which may be used for coastal infrastructure planning, increase with time during the twenty-first century through a combination of mean sea-level rise and the increased uncertainty of future projections with time. The allowances show significant spatial variation, mainly a consequence of strong regionally varying relative sea-level change as a result of glacial isostatic adjustment (GIA). A methodology is described for replacement of the GIA component of the AR5 projection with global positioning system (GPS) measurements of vertical crustal motion; this significantly decreases allowances in regions where the uncertainty of the GIA models is large. For RCP8.5 with GPS data incorporated and for the 1995–2100 period, the sea-level allowances range from about 0.5?m along the north shore of the Gulf of St. Lawrence to more than 1?m along the coast of Nova Scotia and southern Newfoundland.     

On the inclusion of ground-based gravity measurements to the calibration process of a global rainfall-discharge reservoir model: case of the Durzon karst system (Larzac, southern France)

This work examines the relevance of the inclusion of ground-based gravity data in the calibration process of a global rainfall-discharge reservoir model. The analysis is performed for the Durzon karst system (Larzac, France). The first part of the study focuses on the hydrological interpretation of the ground-based gravity measurements. The second part of the study investigates further the information content of the gravity data with respect to water storage dynamics modelling. The gravity-derived information is found unable to either reduce equifinality of the single-objective, discharge-based model calibration process or enhance model performance through assimilation.     

Oxygen isotopic composition of Central Sierra Nevada precipitation, I. Identification of ice-phase water capture regions in winter storms

Measurements of the stable isotopic ratios of oxygen and simultaneously observed ice crystal structure in freshly fallen snow, have been used to estimate the weighted mean altitudes of ice-phase precipitation formation in winter clouds over the Central Sierra Nevada.Observations of dominant, diffusionally grown ice-crystal habits were used to estimate relatively narrow ranges of temperatures of initial formation of the precipitation particles using the Nakaya (1954) ice-crystal classification techniques. The mean oxygen isotope ratio ¹⁸O/¹⁶O for each snow sample, together with local upper air soundings, were used to estimate ranges of temperature-altitude within the clouds where the precipitating particles had captured their ice-phase water. For this initial study, snow samples were collected each five (5) to ten (10) minutes during three snowfall periods on 27 January, 26 February and 28 February, 1983.For the 27 January case, the ice formation mechanism was predominantly vapor deposition and hence the narrow range of temperatures determined by the ice-crystal habits was directly relatable to the mean δ¹⁸O value.The 26 and 28 February cases were more complicated because the solid-phase precipitation was formed by vapor deposition and by freezing of drops and droplets (which occurs without isotopic fractionation). In these cases, the oxygen isotopic composition of the snow reaching the ground was representative of solid phase precipitation which had formed at warmer temperatures than those corresponding to the primary ice-crystal habits alone. There was no apparent relationship between δ¹⁸O values and surface temperature at the sampling site.If relationships between ¹⁸O/¹⁶O and temperature-altitude are established for this geograpic region for winter snowfall conditions, this crystal habit—isotopic composition technique can provide knowledge about the regions of the clouds in which the ice-phase precipitation is forming over the Central Sierra Nevada.