Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

Century-Scale Change in Water Availability: CO2-Quadrupling Experiment

It has been suggested that, unless a major effort is made, the atmospheric concentration of carbon dioxide may rise above four times the pre-industrial level in a few centuries. Here we use a coupled atmosphere-ocean-land model to explore the response of the global water cycle to such a large increase in carbon dioxide, focusing on river discharge and soil moisture. Our results suggest that water is going to be more plentiful in those regions of the world that are already `water-rich'. However, water stresses will increase significantly in regions and seasons that are already relatively dry. This could pose a very challenging problem for water-resource management around the world. For soil moisture, our results indicate reductions during much of the year in many semi-arid regions of the world, such as the southwestern region of North America, the northeastern region of China, the Mediterranean coast of Europe, and the grasslands of Australia and Africa. In some of these regions, soil moisture values are reduced by almost a factor of two during the dry season. The drying in semi-arid regions is likely to induce the outward expansion of deserts to the surrounding regions. Over extensive regions of both the Eurasian and North American continents in high and middle latitudes, soil moisture decreases in summer but increases in winter, in contrast to the situation in semi-arid regions. For river discharge, our results indicate an average increase of ～ 15% during the next few centuries. The discharges from Arctic rivers such as the Mackenzie and Ob' increase by much larger fractions. In the tropics, the discharges from the Amazonas and Ganga-Brahmaputra also increase considerably. However, the percentage changes in runoff from other tropical and many mid-latitude rivers are smaller.     

Assessment of predictions of the response of the Hualien containment model during forced vibration tests

Blind predictions for the response of the 1/4-scale reinforced concrete Hualien (Taiwan) containment model during forced vibration tests are compared with the observed data. The predictions obtained by the CLASSI approach reflect the experimental conditions prior to and after backfill of the soil surrounding the embedded foundation. The experimental data show a strong and unexpected coupling between the response in the NS and EW directions which is not present in the results for the axisymmetric theoretical models. Also, significant differences can be seen between the experimental responses in the two orthogonal horizontal directions which minimize cross-axis coupling. Although these differences are not accounted for in the theoretical models, the discrepancies between predictions and observations are within the uncertainty of the structural and geotechnical data. The obtained differences between predictions and observations give an excellent measure of the prediction error that can be expected in this type of analysis from uncertainty in the data. A detailed assessment of the initial structural and geotechnical data based on extensive comparisons with the results of previous identification studies is also presented. Finally, comparisons between the observed response and calculations based on revised models for the structure and the soil show that current methods of analysis can account accurately for the observed response.     

Evolution of sand waves in the Messina Strait, Italy

In the present paper the morphodynamics of sand waves in the Messina Strait, Italy, is analysed by comparing data gathered during two different surveys carried out in 1991 and 2001, respectively. In particular, a morphometric analysis on the most recent data and a qualitative analysis of the differences between bottom bedform patterns, are carried out. At locations characterised by greater depths, only minor changes to the planimetric configuration of the field, i.e. crest orientation (which is seen to be orthogonal to the direction of net gross sand transport) and wave length are observed, while differences in wave length and crest direction are more evident in more shallow areas. On the other hand, wave height has significantly increased in the whole field. A possible explanation of such a change, based on the previsions of a theory which relates sand-wave growth and migration to the main components of the tidal ellipse, is provided.Responsible Editor: Jens Kappenberg     

An accurate and robust multigrid algorithm for 2D forward resistivity modelling

We present an adaptation of the full multigrid algorithm in DC resistivity modelling in an effort to increase its accuracy. There is a great difficulty with conventional multigrid solvers in representing the physics of an arbitrary distribution of electrical conductivity on a very coarse grid. In general, conventional rectangular finite‐difference or 5‐point approximations of Poisson's equation cannot represent, at a coarse grid level, the effective anisotropy on a coarse scale which results from fine structure in the model. An exception to this generalization occurs when the principal axes of structural anisotropy are aligned with the coordinate axis. Additional and similarly generated problems arise when a coarse cell is obliged to represent fine structure containing very high conductivity contrasts. We have developed an adaptation of the usual resistive‐network representation of the continuum, which avoids some of these problems, and have compared it with the traditional resistive network currently used. The network adaptation consists of replacing the usual 5‐point Laplacian operator stencil used on the finite‐difference grid with a 9‐point stencil, and the conductivity scalar with a 6‐parameter conductivity parametrization. This parametrization permits representation of arbitrarily orientated anisotropy as well as more complex behaviour related to high conductivity contrasts. The importance of multigrid solvers does not lie in their speed at forward modelling (which is comparable with other methods), but rather in their potential for inverse modelling. Inverse solvers which proceed by refinement of an initially very coarse solution can, in principle, take time only linearly proportional to the number of gridpoints in the final desired model.     

Tectonics and metallogenesis of Proterozoic rocks of the Reading Prong

Detailed geologic mapping, petrography, and major and trace-element analyses of Proterozoic rocks from the Greenwood Lake Quadrangle, New York are compared with chemical analyses and stratigraphic information compiled for the entire Reading Prong. A persistent regional stratigraphy is evident in the mapped area whose geochemistry indicates protoliths consistent with a back-arc marginal basin sequence. The proposed marginal basin may have been floored by an older sialic basement and overlain by a basin-fill sequence consisting of a basal tholeiitic basalt, basic to intermediate volcanic or volcaniclastic rocks and carbonate sediments, a bimodal calc-alkaline volcanic sequence, and finally volcaniclastic, marine, and continental sediments. The presence of high-chlorine biotite and scapolite may indicate circulation of brine fluids or the presence of evaporite layers in the sequence. Abundant, stratabound magnetite deposits with a geologic setting very unlike that of cratonic, Proterozoic banded-iron formations are found throughout the proposed basin sequence. Associated with many of the magnetite deposits is unusual uranium and rare-earth element mineralization. It is proposed here that these deposits formed in an exhalative, volcanogenic, depositional environment within an extensional back-arc marginal basin. Such a tectonic setting is consistent with interpretations of protoliths in other portions of the Reading Prong, the Central Metasedimentary Belt of the Canadian Grenville Province, and recent interpretation of the origin of the Franklin lead-zinc deposits, suggesting a more cohesive evolving arc/back-arc tectonic model for the entire Proterozoic margin of the north-eastern portion of the North American craton.     

Physical disturbance creates bacterial dominance of benthic biological communities in tropical deltaic environments of the Gulf of Papua

Unlike many reactive continental shelf mud deposits in temperate regions, bacteria and microfauna rather than macrofauna typically dominate benthic biomass and activities over large areas of the Gulf of Papua (GoP) deltaic complex, Papua New Guinea. During mid NW monsoon periods (Jan–Feb), macrofaunal densities at Gulf stations were relatively low (), large macroinfauna were absent (upper 25 cm), and small (), surface deposit-feeding polychaetes and tubiculous amphipods were dominant, reflecting a frequently destabilized seabed and high sedimentation/erosion rates. Although frequent physical disturbance generally inhibits development of macrobenthic communities, some regions of the Gulf deposits are periodically colonized and extensively bioturbated during quiescent periods, as shown by preserved biogenic sedimentary structures. Bacterial inventories integrated over the top 20 cm were extremely variable within each sub region of the clinoform complex. A possible bimodal pattern with bathymetric depth and distance offshore may occur: lowest-inventories within the sandy, proximal Fly River delta, an open Gulf inner topset zone (10–20 m) having sites of relatively high inventories, an open Gulf mid-topset region with intermediate values and less extreme variation, and the outer topset—upper foreset zone (40–50 m) where highest values are attained (). Various measures of microbial activity, including measures proportional to the cellular rRNA content and the proportion of dividing cells, indicate extremely productive populations over the upper 1-m of the seabed throughout the Gulf of Papua region. Bacterial biomass (0–20 cm) including data of Alongi et al. (1991, 1992, 1995) varied from a low of in intertidal mud banks to a high of in the topset—foreset zone. Macrofaunal biomass did not exceed in any sampled region, ranging from 0.009±0 to with no obvious correlation with bathymetric depth (1–63 m). Meiofaunal biomass was generally an order of magnitude lower than macrofaunal biomass. Relatively elevated bacterial biomass and high turnover rates are consistent with high measured rates of benthic remineralization, presumably reflecting the rapid response time of bacteria to physical reworking, the associated entrainment of organic substrate, and flushing of metabolites. Solute exchange is also enhanced below the directly mixed surface region, possibly producing ‘far field’ stimulation of microbes in underlying deposits. Physical reworking and reoxidation of sediments between 10 and 50 m water depth maintain suboxic, nonsulfidic conditions in the upper 0.5–1 m despite active microbial communities and high benthic remineralization rates.     

The use of GPR and VES in delineating a contamination plume in a landfill site: a case study in SE Brazil

This paper presents the results of the application of the Ground Penetrating Radar (GPR) method, or Georadar, in outlining a zone of contamination due to solid residues at the waste burial site of Rio Claro in the state of São Paulo, SE Brazil. A total of eight GPR profiles with 50- and 100-MHz antennae were surveyed. Six profiles were located within the landfill site and the remaining two were outside. The main objective of the GPR survey was to evaluate the side extension of contamination. A Vertical Electric Sounding (VES) survey was performed at four points within the site in order to investigate the vertical extent of the contamination plume and to define the bottom of the landfill. Two additional VESs were done outside the landfill with the purpose of determining the top of the ground water table and the geoelectric stratigraphy of the background. From the interpretation of the GPR profiles, it was possible to locate the top of the contamination plume and to infer that it was migrating laterally beyond the limits of the waste disposal site. This was observed along the profile situated close to the highway SP-127, which was about 20 m from the limit of the site. The signature of the contaminant appears as a discontinuous reflector that is believed to be a shallow ground water table. The discontinuity is marked by a shadow zone, which is characteristic of conductive contaminant residues. The contamination did not move far enough to reach a sugar cane plantation located at approximately 100 m from the border of the site. In the regions free from contamination, the ground water table was mapped at approximately 10 m of depth, and it was characterized by a strong and continuous reflector. The radar signal penetrated deep enough and enabled the identification of a second reflector at approximately 14 m deep, interpreted as the contact between the Rio Claro and the Corumbataí formations. The contact is marked by the presence of gravel characterized by ferruginous concretes, which cause the strong amplitude reflection in the GPR profile. Within the landfill site, the quantitative interpretation of the VES results showed the contamination zone. The base of the landfill varies between 11 and 15 m deep. Outside the landfill site, the VES results showed no indication of contamination and allowed the determination of the top of the ground water table and the contact between the Rio Claro and the Corumbataí formations. The results of GPR and VES showed a good agreement and the integrated interpretations were supported by local geology and information from several boreholes, about 17 m depth, on average. The bottom of the landfill reaches a maximum of 14.5 m depth.