Nature and origin of A-type granites with particular reference to southeastern Australia

In the Lachlan Fold Belt of southeastern Australia, Upper Devonian A-type granite suites were emplaced after the Lower Devonian I-type granites of the Bega Batholith. Individual plutons of two A-type suites are homogeneous and the granites are characterized by late interstitial annite. Chemically they are distinguished from I-type granites with similar SiO₂ contents of the Bega Batholith, by higher abundances of large highly charged cations such as Nb, Ga, Y, and the REE and lower Al, Mg and Ca: high Ga/Al is diagnostic. These A-type suites are metaluminous, but peralkaline and peraluminous A-type granites also occur in Australia and elsewhere. Partial melting of felsic granulite is the preferred genetic model. This source rock is the residue remaining in the lower crust after production of a previous granite. High temperature, vapour-absent melting of the granulitic source generates a low viscosity, relatively anhydrous melt containing F and possibly Cl. The framework structure of this melt is considerably distorted by the presence of these dissolved halides allowing the large highly charged cations to form stable high co-ordination structures. The high concentration of Zr and probably other elements such as the REE in peralkaline or near peralkaline A-type melts is a result of the counter ion effect where excess alkali cations stabilize structures in the melt such as alkali-zircono-silicates. The melt structure determines the trace element composition of the granite. Separation of a fluid phase from an A-type magma results in destabilization of co-ordination complexes and in the formation of rare-metal deposits commonly associated with fluorite. At this stage the role of Cl in metal transport is considered more important than F.     

On model validation for meso/submesoscale currents: Metrics and application to ROMS off Central California

ROMS with horizontal grid spacing of 3.5 km for the region off Central California was compared to RAFOS float observations and satellite altimetry on meso/submesoscales. The approach introduced and used two new metrics for model-data comparison, as well as suggested how to calculate these metrics for different spatio-temporal scales. The first metric consisted of the first two moments of exit time and was used to compare ROMS against RAFOS float observations at mid-depths (between 300 m and 350 m). Exit time is the time a float launched at a point takes to leave a domain for the first time. The second metric was spectral entropy and was used to estimate how well ROMS reproduced variability of the sea surface height (SSH) anomaly field extracted from an AVISO data set (1992–2007) for specified temporal and spatial scales. Calculations showed that ROMS reproduced the mid-depth mesoscale/submesoscale currents next to the coast in a very accurate manner (low-order exit time statistics of floats were reproduced by ROMS with an accuracy better than 95%); but ROMS overestimated the speed of westward drift of floats by as much as 20–30% at distances greater than 350 km from the coastline. ROMS predicted the variability of the mesoscale (100–400 km) SSH anomaly field for temporal scales of 1–12 months with a reasonable accuracy. A wavelet transform modulus maxima technique applied to the spectral entropy of SSH anomaly also demonstrated good agreement between ROMS and satellite altimetry for mesoscales characterized by singular exponents and multi-fractal spectra for 1–12 month time scales.     

El Niño- or La Niña-like climate change?

The potential for the mean climate of the tropical Pacific to shift to more El Niño-like conditions as a result of human induced climate change is subject to a considerable degree of uncertainty. The complexity of the feedback processes, the wide range of responses of different atmosphere–ocean global circulation models (AOGCMs) and difficulties with model simulation of present day El Niño southern oscillation (ENSO), all complicate the picture. By examining the components of the climate-change response that projects onto the model pattern of ENSO variability in 20 AOGCMs submitted to the coupled model inter-comparison project (CMIP), it is shown that large-scale coupled atmosphere–ocean feedbacks associated with the present day ENSO also operate on longer climate-change time scales. By linking the realism of the simulation of present day ENSO variability in the models to their patterns of future mean El Niño-like or La Niña-like climate change, it is found that those models that have the largest ENSO-like climate change also have the poorest simulation of ENSO variability. The most likely scenario (p=0.59) in a model-skill-weighted histogram of CMIP models is for no trend towards either mean El Niño-like or La Niña-like conditions. However, there remains a small probability (p=0.16) for a change to El Niño-like conditions of the order of one standard El Niño per century in the 1% per year CO₂ increase scenario.     

Reliability and importance of structural diversity of climate model ensembles

We investigate the performance of the newest generation multi-model ensemble (MME) from the Coupled Model Intercomparison Project (CMIP5). We compare the ensemble to the previous generation models (CMIP3) as well as several single model ensembles (SMEs), which are constructed by varying components of single models. These SMEs range from ensembles where parameter uncertainties are sampled (perturbed physics ensembles) through to an ensemble where a number of the physical schemes are switched (multi-physics ensemble). We focus on assessing reliability against present-day climatology with rank histograms, but also investigate the effective degrees of freedom (EDoF) of the fields of variables which makes the statistical test of reliability more rigorous, and consider the distances between the observation and ensemble members. We find that the features of the CMIP5 rank histograms, of general reliability on broad scales, are consistent with those of CMIP3, suggesting a similar level of performance for present-day climatology. The spread of MMEs tends towards being “over-dispersed” rather than “under-dispersed”. In general, the SMEs examined tend towards insufficient dispersion and the rank histogram analysis identifies them as being statistically distinguishable from many of the observations. The EDoFs of the MMEs are generally greater than those of SMEs, suggesting that structural changes lead to a characteristically richer range of model behaviours than is obtained with parametric/physical-scheme-switching ensembles. For distance measures, the observations and models ensemble members are similarly spaced from each other for MMEs, whereas for the SMEs, the observations are generally well outside the ensemble. We suggest that multi-model ensembles should represent an important component of uncertainty analysis.     

Impact of ocean model resolution on CCSM climate simulations

The current literature provides compelling evidence suggesting that an eddy-resolving (as opposed to eddy-permitting or eddy-parameterized) ocean component model will significantly impact the simulation of the large-scale climate, although this has not been fully tested to date in multi-decadal global coupled climate simulations. The purpose of this paper is to examine how resolved ocean fronts and eddies impact the simulation of large-scale climate. The model used for this study is the NCAR Community Climate System Model version 3.5 (CCSM3.5)—the forerunner to CCSM4. Two experiments are reported here. The control experiment is a 155-year present-day climate simulation using a 0.5° atmosphere component (zonal resolution 0.625 meridional resolution 0.5°; land surface component at the same resolution) coupled to ocean and sea-ice components with zonal resolution of 1.2° and meridional resolution varying from 0.27° at the equator to 0.54° in the mid-latitudes. The second simulation uses the same atmospheric and land-surface models coupled to eddy-resolving 0.1° ocean and sea-ice component models. The simulations are compared in terms of how the representation of smaller scale features in the time mean ocean circulation and ocean eddies impact the mean and variable climate. In terms of the global mean surface temperature, the enhanced ocean resolution leads to a ubiquitous surface warming with a global mean surface temperature increase of about 0.2?°C relative to the control. The warming is largest in the Arctic and regions of strong ocean fronts and ocean eddy activity (i.e., Southern Ocean, western boundary currents). The Arctic warming is associated with significant losses of sea-ice in the high-resolution simulation. The sea surface temperature gradients in the North Atlantic, in particular, are better resolved in the high-resolution model leading to significantly sharper temperature gradients and associated large-scale shifts in the rainfall. In the extra-tropics, the interannual temperature variability is increased with the resolved eddies, and a notable increases in the amplitude of the El Ni?o and the Southern Oscillation is also detected. Changes in global temperature anomaly teleconnections and local air-sea feedbacks are also documented and show large changes in ocean–atmosphere coupling. In particular, local air-sea feedbacks are significantly modified by the increased ocean resolution. In the high-resolution simulation in the extra-tropics there is compelling evidence of stronger forcing of the atmosphere by SST variability arising from ocean dynamics. This coupling is very weak or absent in the low-resolution model.     

Equations for the Drag Force and Aerodynamic Roughness Length of Urban Areas with Random Building Heights

We use a conceptual model to investigate how randomly varying building heights within a city affect the atmospheric drag forces and the aerodynamic roughness length of the city. The model is based on the assumptions regarding wake spreading and mutual sheltering effects proposed by Raupach (Boundary-Layer Meteorol 60:375?C395, 1992). It is applied both to canopies having uniform building heights and to those having the same building density and mean height, but with variability about the mean. For each simulated urban area, a correction is determined, due to height variability, to the shear stress predicted for the uniform building height case. It is found that u _*/u _*R , where u _* is the friction velocity and u _*R is the friction velocity from the uniform building height case, is expressed well as an algebraic function of ?? and ?? _h /h _m , where ?? is the frontal area index, ?? _h is the standard deviation of the building height, and h _m is the mean building height. The simulations also resulted in a simple algebraic relation for z ₀/z _0R as a function of ?? and ?? _h /h _m , where z ₀ is the aerodynamic roughness length and z _0R is z ₀ found from the original Raupach formulation for a uniform canopy. Model results are in keeping with those of several previous studies.     

Modeling damage and deformation in impact simulations

Abstract— Numerical modeling is a powerful tool for investigating the formation of large impact craters but is one that must be validated with observational evidence. Quantitative analysis of damage and deformation in the target surrounding an impact event provides a promising means of validation for numerical models of terrestrial impact craters, particularly in cases where the final pristine crater morphology is ambiguous or unknown. In this paper, we discuss the aspects of the behavior of brittle materials important for the accurate simulation of damage and deformation surrounding an impact event and the care required to interpret the results. We demonstrate this with an example simulation of an impact into a terrestrial, granite target that produces a 10 km‐diameter transient crater. The results of the simulation are shown in terms of damage (a scalar quantity that reflects the totality of fragmentation) and plastic strain, both total plastic strain (the accumulated amount of permanent shear deformation, regardless of the sense of shear) and net plastic strain (the amount of permanent shear deformation where the sense of shear is accounted for). Damage and plastic strain are both greatest close to the impact site and decline with radial distance. However, the reversal in flow patterns from the downward and outward excavation flow to the inward and upward collapse flow implies that net plastic strains may be significantly lower than total plastic strains. Plastic strain in brittle rocks is very heterogeneous; however, continuum modeling requires that the deformation of the target during an impact event be described in terms of an average strain that applies over a large volume of rock (large compared to the spacing between individual zones of sliding). This paper demonstrates that model predictions of smooth average strain are entirely consistent with an actual strain concentrated along very narrow zones. Furthermore, we suggest that model predictions of total accumulated strain should correlate with observable variations in bulk density and seismic velocity.     

A preliminary assessment of the spatial sources of contemporary suspended sediment in the Ohio River basin,United States,using water quality data from the NASQAN programme in a source tracing procedure

Reliable information on catchment scale suspended sediment sources is required to inform the design of management strategies for helping abate the numerous environmental issues associated with enhanced sediment mobilization and off‐site loadings. Since sediment fingerprinting techniques avoid many of the logistical constraints associated with using more traditional indirect measurement methods at catchment scale, such approaches have been increasingly reported in the international literature and typically use data sets collected specifically for sediment source apportionment purposes. There remains scope for investigating the potential for using geochemical data sets assembled by routine monitoring programmes to fingerprint sediment provenance. In the United States, routine water quality samples are collected as part of the US Geological Survey's revised National Stream Quality Accounting Network programme. Accordingly, the geochemistry data generated from these samples over a 10‐year period (1996–2006) were used as the basis for a fingerprinting exercise to assess the key tributary sub‐catchment spatial sources of contemporary suspended sediment transported by the Ohio River. Uncertainty associated with the spatial source estimates was quantified using a Monte Carlo approach in conjunction with mass balance modelling. Relative frequency weighted means were used as an alternative way of summarizing the spatial source contributions, thereby avoiding the need to use confidence limits. The results should be interpreted in the context of the routine, but infrequent nature, of the suspended sediment samples used to assemble geochemistry as a basis for the sourcing exercise. Nonetheless, the study demonstrates how routine monitoring samples can be used to provide some preliminary information on sediment provenance in large drainage basins. Copyright © 2011 John Wiley & Sons, Ltd.