The dynamic behaviour of a river-dominated tidal inlet, River Murray, Australia

Australia's largest river, the River Murray, discharges to the southern ocean through a coastal lagoon and river-dominated tidal inlet. Increased water extractions upstream for irrigation have led to significantly reduced flows at the mouth and, as a result, the area is undergoing rapid change, particularly with regard to the rate at which sediment is being transported into the lagoon. Based on detailed and accurate bathymetric surveys it has been possible to estimate that the rate of lagoon in-filling is of the order of 100,000 m³ per year for the period June 2000 to May 2003, although the actual rate shows significant year to year variability. Dredging of the lagoon commenced in 2000 in an attempt to reverse the trend.In an effort to understand the behaviour of the inlet a one-dimensional numerical model of the inlet has been developed. The model extends the original of van de Kreeke by including a dynamic inlet throat area based on predicted river flows and a sediment transport module to predict the resulting net sediment transport. Comparisons with water level data collected on both the ocean and lagoon sides of the mouth have shown that the model is able to predict the attenuation and lag of the tidal signal reasonably well. The sediment transport model was based on predicted sediment concentrations in the surf zone and was found to predict the rate of sediment in-filling to an acceptable level of accuracy. It is envisaged that the model will be a useful management tool, especially since it is possible to manipulate river discharges to the mouth.     

A singular evolutive extended Kalman filter to assimilate ocean color data in a coupled physical–biochemical model of the North Atlantic ocean

Within the European DIADEM project, a data assimilation system for coupled ocean circulation and marine ecosystem models has been implemented for the North Atlantic and the Nordic Seas. One objective of this project is to demonstrate the relevance of sophisticated methods to assimilate satellite data such as altimetry, surface temperature and ocean color, into realistic ocean models. In this paper, the singular evolutive extended Kalman (SEEK) filter, which is an advanced assimilation scheme where three-dimensional, multivariate error statistics are taken into account, is used to assimilate ocean color data into the biological component of the coupled system. The marine ecosystem model, derived from the FDM model [J. Mar. Res. 48 (1990) 591], includes 11 nitrogen and carbon compartments and describes the synthesis of organic matter in the euphotic zone, its consumption by animals of upper trophic levels, and the recycling of detritic material in the deep ocean. The circulation model coupled to the ecosystem is the Miami isopycnic coordinate ocean model (MICOM), which covers the Atlantic and the Arctic Oceans with an enhanced resolution in the North Atlantic basin. The model is forced with realistic ECMWF ocean/atmosphere fluxes, which permits to resolve the seasonal variability of the circulation and mixed layer properties. In the twin assimimation experiments reported here, the predictions of the coupled model are corrected every 10 days using pseudo-measurements of surface phytoplankton as a substitute to chlorophyll concentrations measured from space. The diagnostics of these experiments indicate that the assimilation is feasible with a reduced-order Kalman filter of small rank (of order 10) as long as a sufficiently good identification of the error structure is available. In addition, the control of non-observed quantities such as zooplankton and nitrate concentrations is made possible, owing to the multivariate nature of the analysis scheme. However, a too severe truncation of the error sub-space downgrades the propagation of surface information below the mixed layer. The reduction of the actual state vector to the surface layers is therefore investigated to improve the estimation process in the perspective of sea-viewing wide field-of-view sensor (SeaWiFS) data assimilation experiments.     

Intercomparison of coastal area morphodynamic models

Five different coastal area morphodynamic models have been set up to run on the same offshore breakwater layout and an intercomparison carried out on the hydrodynamic and morphodynamic output produced by each scheme. In addition, the predicted morphodynamics was checked against available laboratory and field data.It is concluded that the models are capable of producing realistic estimates for the dominant morphodynamic features associated with offshore breakwaters. Coupling of the wave, current and sediment transport components of each scheme is shown to yield bathymetry which attains a state of equilibrium, unlike models which are based on the initial transport field only.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

Rare Earth Deposits of North America

Rare earth elements (REE) have been mined in North America since 1885, when placer monazite was produced in the southeast USA. Since the 1960s, however, most North American REE have come from a carbonatite deposit at Mountain Pass, California, and most of the world’s REE came from this source between 1965 and 1995. After 1998, Mountain Pass REE sales declined substantially due to competition from China and to environmental constraints. REE are presently not mined at Mountain Pass, and shipments were made from stockpiles in recent years. Chevron Mining, however, restarted extraction of selected REE at Mountain Pass in 2007. In 1987, Mountain Pass reserves were calculated at 29 Mt of ore with 8.9% rare earth oxide based on a 5% cut‐off grade. Current reserves are in excess of 20 Mt at similar grade. The ore mineral is bastnasite, and the ore has high light REE/heavy REE (LREE/HREE). The carbonatite is a moderately dipping, tabular 1.4‐Ga intrusive body associated with ultrapotassic alkaline plutons of similar age. The chemistry and ultrapotassic alkaline association of the Mountain Pass deposit suggest a different source than that of most other carbonatites. Elsewhere in the western USA, carbonatites have been proposed as possible REE sources. Large but low‐grade LREE resources are in carbonatite in Colorado and Wyoming. Carbonatite complexes in Canada contain only minor REE resources. Other types of hard‐rock REE deposits in the USA include small iron‐REE deposits in Missouri and New York, and vein deposits in Idaho. Phosphorite and fluorite deposits in the USA also contain minor REE resources. The most recently discovered REE deposit in North America is the Hoidas Lake vein deposit, Saskatchewan, a small but incompletely evaluated resource. Neogene North American placer monazite resources, both marine and continental, are small or in environmentally sensitive areas, and thus unlikely to be mined. Paleoplacer deposits also contain minor resources. Possible future uranium mining of Precambrian conglomerates in the Elliott Lake–Blind River district, Canada, could yield by‐product HREE and Y. REE deposits occur in peralkaline syenitic and granitic rocks in several places in North America. These deposits are typically enriched in HREE, Y, and Zr. Some also have associated Be, Nb, and Ta. The largest such deposits are at Thor Lake and Strange Lake in Canada. A eudialyte syenite deposit at Pajarito Mountain in New Mexico is also probably large, but of lower grade. Similar deposits occur at Kipawa Lake and Lackner Lake in Canada. Future uses of some REE commodities are expected to increase, and growth is likely for REE in new technologies. World reserves, however, are probably sufficient to meet international demand for most REE commodities well into the 21st century. Recent experience shows that Chinese producers are capable of large amounts of REE production, keeping prices low. Most refined REE prices are now at approximately 50% of the 1980s price levels, but there has been recent upward price movement for some REE compounds following Chinese restriction of exports. Because of its grade, size, and relatively simple metallurgy, the Mountain Pass deposit remains North America’s best source of LREE. The future of REE production at Mountain Pass is mostly dependent on REE price levels and on domestic REE marketing potential. The development of new REE deposits in North America is unlikely in the near future. Undeveloped deposits with the most potential are probably large, low‐grade deposits in peralkaline igneous rocks. Competition with established Chinese HREE and Y sources and a developing Australian deposit will be a factor.     

Glacio-lacustrine stratigraphy,aquifer characterization and contaminant transport: a case study in South Lake Tahoe,California, USA

A hydrogeologic model that has been used by many researchers and consultants to describe an area of South Lake Tahoe, California, USA impacted by MTBE contamination describes a relatively homogeneous unconfined aquifer comprised of poorly sorted glacial outwash deposits, within which water-supply wells are able to exert significant alteration in natural groundwater flow. A re-examination of the area’s hydrogeology is presented, which supports a layered heterogeneous aquifer system constructed of alternating fine and coarser-grained glacio-lacustrine depositional units. This re-evaluation was accomplished through a review of lithologic logs across an area of approximately 1 km², combined with observations of significant hydraulic head differences and knowledge of the depositional environments controlled by Pleistocene Lake Tahoe high stands. Many of the fine-grained units observed at depths from 6 to 15 m, although relatively thin, are generally continuous and serve as significant barriers to groundwater flow. The vertical migration of contamination across these fine-grained units to deeper groundwaters was facilitated by cross-screened monitoring wells installed as part of site investigation activities. This conclusion highlights the importance of geologic characterization and proper monitoring well construction at contaminated site investigations.

---

Electronic supplementary material   The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Comparison of numerical models of two debris flows in the Cortina d’ Ampezzo area,Dolomites, Italy

The accurate prediction of runout distances, velocities and the knowledge of flow rheology can reduce the casualties and property damage produced by debris flows, providing a means to delineate hazard areas, to estimate hazard intensities for input into risk studies and to provide parameters for the design of protective measures. The application of most of models that describe the propagation and deposition of debris flow requires detailed topography, rheological and hydrological data that are not always available for the debris-flow hazard delineation and estimation. In the Cortina d’Ampezzo area, Eastern Dolomites, Italy, most of the slope instabilities are represented by debris flows; 325 debris-flow prone watersheds have been mapped in the geomorphological hazard map of this area. We compared the results of simulations of two well-documented debris flows in the Cortina d’Ampezzo area, carried on with two different single-phase, non-Newtonian models, the one-dimensional DAN-W and the two-dimensional FLO-2D, to test the possibility to simulate the dynamic behaviour of a debris flow with a model using a limited range of input parameters. FLO-2D model creates a more accurate representation of the hazard area in terms of flooded area, but the results in terms of runout distances and deposits thickness are similar to DAN-W results. Using DAN-W, the most appropriate rheology to describe the debris-flow behaviour is the Voellmy model. When detailed topographical, rheological and hydrological data are not available, DAN-W, which requires less detailed data, is a valuable tool to predict debris-flow hazard. Parameters obtained through back-analysis with both models can be applied to predict hazard in other areas characterized by similar geology, morphology and climate.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.