Distribution and morphology of sinkholes triggered by flooding following Tropical Storm Alberto at Albany, Georgia, USA

Flooding of the Flint River in July 1994 triggered the collapse of at least 312 sinkholes in the karstic Dougherty Plain at Albany, Georgia. We examined the distribution and morphology of these new sinkholes to evaluate the mode of formation, to characterize early stages of the evolution of sinkhole form, and to estimate the lowering of the surface associated with the development of new sinkholes.Eighty-eight percent of sinkholes occur inside the limits of flooding, especially in areas of sandy overburden, and they often follow joint-controlled linear trends. Sinkhole dimensions are log-normally distributed with median values of circumference = 5.7 m, length = 1.8 m, width = 1.6 m, and depth = 0.7 m; asymmetry (L:W) = 1.2. Cross-sectional forms range from narrow cylinders to large bowls, with many sinkholes having undercut sides.Flooding triggered the formation of sinkholes by saturating and liquefying overburden, which caused soil arches to collapse and flow into cavities in bedrock. The prevalence of sinkholes near the periphery of flooding suggests that drainage and loss of buoyant support as flood waters subsided may also have contributed to failure. A volume ratio index is used to quantify the three dimensional geometric form of sinkholes. Initially, small cylindrical shafts open over a bedrock joint, followed by progressive slumping that leads to widening and increases in volume to a final bowl form. Estimates of the aggregate volume of overburden transported underground in flooded areas range from 7,990 to 11,130 m3. Averaged over flooded areas, this accounts for 0.26 to 0.37 mm/km² lowering of the surface. Based on a 500 year recurrence interval for the flood event, values for lowering of the surface range from 0.52 to 0.74 mm per 1,000 years. These values are an order of magnitude less than estimates of carbonate dissolution and suggest that transport of overburden underground is limited by triggering events.     

Improved Representation of Land-surface Heterogeneity in a Non-hydrostatic Numerical Weather Prediction Model

This study focuses on the relevance of accurate surface parameters, in particular soil moisture, and of parameterizations for heterogeneous land surfaces, for the prediction of sensible and latent heat fluxes by a mesoscale weather forecast model with horizontal grid resolution of 7 km. The analysis is based on model integrations for a 30-day period, which are compared both to flux measurements obtained from the LITFASS-2003 field experiment and to high-resolution-model (1-km grid spacing) results. At first, the relevance of improved parameter sets and input data compared to usual operational practice for an accurate prediction of near-surface fluxes is shown and discussed. It is demonstrated that an observation-based land-surface assimilation scheme leads to an improved soil moisture analysis, which is shown to be essential for the realistic simulation of surface fluxes. Secondly, the implementation of two efficient parameterization strategies for subgrid-scale variability of the surface, the mosaic and the tile approach, is presented. Using these methods, the simulations are in better agreement with measurements than simulations with simple aggregation methods that use effective surface parameters. Integrations with the mosaic approach reproduce high resolution simulations very well and more accurately than simulations with the tile method. Finally, the high resolution simulations are analyzed to justify and discuss the approximations underlying both methods.     

A furrow runs through it: An example of sustainable traditional irrigation in Western Cape Province,South Africa

Water and water related matters are important issues in water scarce countries like the Republic of South Africa. This is especially true for the agricultural sector that relies heavily on the availability of water for irrigation to ensure a secure food supply. In South Africa, present day irrigation projects rely on government subsidies and are usually large undertakings beyond the reach of individuals or small communities. The sustainability of such large irrigation projects is still a matter of debate. A 100 years ago a different set of rules applied. Irrigation could only be achieved by hand-dug furrows that supplied water in relatively small quantities to small farming communities. This paper investigated such a small community, the Leeukloof farming community in the Western Cape Province of the Republic of South Africa. For the past 93 years this community has maintained a furrow that supplies water not only to farms bordering the Perdeberg River, which supplies the furrow with water, but also to farms not bordering the river. The authors try to reconstruct the changes in agricultural practices brought about by the furrow, methods of water extraction during each phase, as well as the development and role of the furrow system. The perception by the community of the necessity of this method of water extraction from the river system is also investigated. This revised version was published online in July 2006 with corrections to the Cover Date.     

Comparison of Positron Emission Tomography and X-ray radiography for studies of physical processes in sandstone

The migration of moisture in cores of porous homogeneous sandstone of Ledian age (Belgian Eocene, Tertiary) is monitored with PET (Positron Emission Tomography) and micro-focus X-ray radiography. In the case of PET, a nuclear medical imaging technique, ⁵⁵Co-EDTA (Ethylenediamine Tetraacetic acid) and R-¹⁸F were used as water-soluble tracers. The X-ray projection method has evolved from the better-known medical technique and allows a fast and accurate determination of the two-dimensional transient moisture content profiles. Results indicate that both techniques can deliver important information concerning physical processes in situ.     

Cycling the representer algorithm for variational data assimilation with a nonlinear reduced gravity ocean model

The representer method was used by [Ngodock, H.E., Jacobs, G.A., Chen, M., 2006. The representer method, the ensemble Kalman filter and the ensemble Kalman smoother: a comparison study using a nonlinear reduced gravity ocean model. Ocean Modelling 12, 378–400] in a comparison study with the ensemble Kalman filter and smoother involving a 1.5 nonlinear reduced gravity idealized ocean model simulating the Loop Current (LC) and the Loop Current eddies (LCE) in the Gulf of Mexico. It was reported that the representer method was more accurate than its ensemble counterparts, yet it had difficulties fitting the data in the last month of the 4-month assimilation window when the data density was significantly decreased. The authors attributed this failure to increased advective nonlinearities in the presence of an eddy shedding causing the tangent linear model (TLM) to become inaccurate. In a separate study [Ngodock, H.E., Smith, S.R., Jacobs, G.A., 2007. Cycling the representer algorithm for variational data assimilation with the Lorenz attractor. Monthly Weather Review 135 (2), 373–386] applied the cycling representer algorithm to the Lorenz attractor and demonstrated that the cycling solution was able to accurately fit the data within each cycle and beyond the range of accuracy of the TLM, once adjustments were made in the early cycles, thus overcoming the difficulties of the non-cycling solution. The cycling algorithm is used here in assimilation experiments with the nonlinear reduced gravity model. It is shown that the cycling solution overcomes the difficulties encountered by the non-cycling solution due to a limited time range of accuracy of the TLM. Thus, for variational assimilation applications where the TLM accuracy is limited in time, the cycling representer becomes a very powerful and attractive alternative, given that its computational cost is significantly lower than that of the non-cycling algorithm.     

A high-resolution radio study of neutral gas in the starburst galaxy NGC 520

We present subarcsec angular resolution observations of the neutral gas in the nearby starburst galaxy NGC 520. The central kpc region of NGC 520 contains an area of significantly enhanced star formation. The radio continuum structure of this region resolves into ∼10 continuum components. By comparing the flux densities of the brightest of these components at 1.4 GHz with published 15-GHz data we infer that these components detected at 1.4 and 1.6 GHz are related to the starburst and are most likely to be collections of several supernova remnants within the beam. None of these components is consistent with emission from an active galactic nuclei. Both neutral hydrogen (H  i ) and hydroxyl (OH) absorption lines are observed against the continuum emission, along with a weak OH maser feature probably related to the star formation activity in this galaxy. Strong H  i absorption  ( N _H∼ 10²² atoms cm⁻²)  traces a velocity gradient of 0.5 km s⁻¹ pc⁻¹ across the central kpc of NGC 520. The H  i absorption velocity structure is consistent with the velocity gradients observed in both the OH absorption and in CO emission observations. The neutral gas velocity structure observed within the central kpc of NGC 520 is attributed to a kpc-scale ring or disc. It is also noted that the velocity gradients observed for these neutral gas components appear to differ with the velocity gradients observed from optical ionized emission lines. This apparent disagreement is discussed and attributed to the extinction of the optical emission from the actual centre of this source hence implying that optical ionized emission lines are only detected from regions with significantly different radii to those sampled by the observations presented here.     

Monthly Variations of Water Masses in the Yellow and East China Seas, November 6, 1998

The monthly water mass variations in the Yellow Sea and the East China Sea are investigated using over 40 years of historical temperature and salinity observations via a cluster analysis that incorporates geographical distance and depth separation in addition to the temperature and salinity. Results delineate monthly variations in the major water masses and provide some insight into formation mechanisms and intermixing. The major water masses include the Kuroshio-East China Sea water (KE), the Yellow Sea surface water (YSS) and bottom cold water (YSB), mixed water (MW), and coastal water (CW). The distribution of the KE water mass reveals the intrusion pattern into the area west of Cheju. A separate mixed water type appears between the KE water mass and the Yellow Sea water masses during winter. The formation mechanism of the YSB appears to be the surface cooling and active mixing in winter. In the East China Sea, during summer, surface water is differentiated from the subsurface water while there is no differentiation during winter. In the Yellow Sea, a three layer system exists in the summer and fall (May–November) while a two layer system exists during the rest of the year. A fresh water mass generated by Yangtze River discharge (YD) is present over the northern East China Sea and the southern Yellow Sea during summer. This revised version was published online in August 2006 with corrections to the Cover Date.     

Quantum-thermodynamic treatment of intrinsic anharmonicity; Wallace’s theorem revisited

Wallace (in Thermodynamics of crystals, 1972) developed a theorem, rooted in rigid lattice dynamics, which incorporates intrinsic anharmonic effects in solids. The practical application of this theorem in mineral physics is computationally involved and this is the main reason for the theorem not getting the attention it deserves. Because intrinsic anharmonicity is an important issue at the extreme conditions in planetary mantles, we derived a method which removes the computational obstacles in applying this theorem. We extended the theorem to incorporate details of the phonon spectrum and tested our algorithm on forsterite (Mg₂SiO₄). Using a least squares inversion technique applied to all available experimental data, we show that it results in an accurate representation of thermodynamic properties and sound wave velocities of Mg₂SiO₄ in its complete pressure–temperature stability range. We also show that the accuracy of our results is not significantly affected by the use of a different equation of state.