Improving a subgrid runoff parameterization scheme for climate models by the use of high resolution data derived from satellite observations

In this study it is shown that the availability of a very high resolution dataset of land surface characteristics leads to the improvement of a surface runoff parameterization scheme. The improved parameterization scheme was developed for application in global and regional climate models and is a further development of the Arno scheme that is widely used in climate models. Here, surface runoff is computed as infiltration excess from a "bucket" type reservoir which takes the subgrid variability of soil saturation within a model gridbox into account. Instead of prescribing a distribution of subgrid scale soil water capacities as in the original Arno scheme, the array of high resolution soil water capacities taken from a global 1 km dataset of land surface parameters is used to obtain individual fractional saturation curves for each model gridbox. From each saturation curve, the three parameters (a shape parameter describing the shape of the subgrid distribution of soil water capacities, subgrid minimum and maximum soil water capacity) required in the modified formulation of the scheme are derived via optimization. As in the original Arno scheme applied in the ECHAM general circulation model and the REMO regional climate model, topography variations will influence the distribution of saturated subgrid areas within a model gridbox. At most gridboxes the net effect of these changes is such that more runoff is produced for high soil water contents and less runoff for low soil water contents. A validation of simulated discharge computed with a simplified land surface scheme applied to reanalysis data of the European Centre for Medium-Range Weather Forecasts and a hydrological discharge model has shown that these changes lead to a more realistic simulation of the annual cycle of discharge for several catchments. In particular this could be shown for the Yangtze Kiang and Amur catchments where adequate input data are available.     

Impact of AMD on water quality in critical watershed in the Hudson River drainage basin: Phillips Mine,Hudson Highlands,New York

A sulfur and trace element enriched U–Th-laced tailings pile at the abandoned Phillips Mine in Garrison, New York, releases acid mine drainage (AMD, generally pH < 3, minimum pH 1.78) into the first-order Copper Mine Brook (CMB) that drains into the Hudson River. The pyrrhotite-rich Phillips Mine is located in the Highlands region, a critical water source for the New York metro area. A conceptual model for derivation/dissolution, sequestration, transport and dilution of contaminants is proposed. The acidic water interacts with the tailings, leaching and dissolving the trace metals. AMD evaporation during dry periods concentrates solid phase trace metals and sulfate, forming melanterite (FeSO₄·7H₂O) on sulfide-rich tailings surfaces. Wet periods dissolve these concentrates/precipitates, releasing stored acidity and trace metals into the CMB. Sediments along CMB are enriched in iron hydroxides which act as sinks for metals, indicating progressive sequestration that correlates with dilution and sharp rise in pH when mine water mixes with tributaries. Seasonal variations in metal concentrations were partly attributable to dissolution of the efflorescent salts with their sorbed metals and additional metals from surging acidic seepage induced by precipitation.     

Conceptual model of recharge to southeastern Badain Jaran Desert groundwater and lakes from environmental tracers

Sources of groundwater recharge to the Badain Jaran Desert in China have been investigated using geochemical and isotopic techniques. Stable isotope compositions (δ¹⁸O and δ²H) of shallow groundwater and surface water from oasis lakes evolve from a starting composition considerably depleted compared to local unsaturated zone moisture, confirming inferences from chloride mass balance that direct infiltration of precipitation is not a volumetrically important source of recharge to the shallow aquifer in the study area. Shallow phreatic and deeper confined groundwater bodies appear unconnected based on chemical composition and radiocarbon activities. Hydrogeologic evidence points toward a bordering mountain range (Yabulai) as a likely recharge zone, which is consistent with tracer results. A mean residence time in the range of 1–2 ka for the desert’s southeastern margin is inferred from radiocarbon. These results reveal that some replenishment to the desert aquifer is occurring but at a rate much lower than previously suggested, which is relevant for water resources planning in this ecologically sensitive area.     

Impacts of deforestation and afforestation in the Mediterranean region as simulated by the MPI atmospheric GCM

For two reasons it is important to study the sensitivity of the global climate to changes in the vegetation cover over land. First, in the real world, changes in the vegetation cover may have regional and global implications. Second, in numerical simulations, the sensitivity of the simulated climate may depend on the specific parameterization schemes employed in the model and on the model's large-scale systematic errors. The Max-Planck-Institute's global general circulation model ECHAM4 has been used to study the sensitivity of the local and global climate during a full annual cycle to deforestation and afforestation in the Mediterranean region. The deforestation represents an extreme desertification scenario for this region. The changes in the afforestation experiment are based on the pattern of the vegetation cover 2000 years before present when the climate in the Mediterranean was more humid. The comparison of the deforestation integration to the control shows a slight cooling at the surface and reduced precipitation during the summer as a result of less evapotranspiration of plants and less evaporation from the assumption of eroded soils. There is no significant signal during the winter season due to the stronger influence of the mid-latitude baroclinic disturbances. In general, the results of the afforestation experiment are opposite to those of the deforestation case. A significant response was found in the vicinity of grid points where the land surface characteristics were modified. The response in the Sahara in the afforestation experiment is in agreement with the results from other general circulation model studies.     

Arsenic enrichment in unconfined sections of the southern Gulf Coast aquifer system, Texas

Groundwater arsenic concentrations exceeding the federal drinking water standard are common in the southern Gulf Coast aquifer system in Texas, including in aerobic, unconfined groundwater which provides much of the municipal and domestic water supplies for the region. The objective of this study was to determine geochemical factors affecting the occurrence and distribution of groundwater As in unconfined portions of the southern Gulf Coast aquifer system through a comparative transect study of groundwater across three major hydrostratigraphic units (the Catahoula Formation, Jasper aquifer and Evangeline aquifer) and analysis of regional water quality data. Results show that As concentrations decrease with increasing distance from the Catahoula Formation, which is consistent with Miocene volcanic ash as the main source of As to groundwater in the region. Arsenic concentrations correlate with V, SiO₂ and K, all of which were released during weathering of volcanic sediments and their degradation products. In all three units, carbonate weathering and active recharge in the unconfined zones result in circum-neutral pH and oxidizing groundwater, which are typically amenable to As immobilization by adsorption of arsenate onto mineral oxides and clays. However, As concentrations exceed 10 μg/L in approximately 30% of wells. Silica that was co-released with As may compete for sorption sites and reduce the capacity for arsenate adsorption.     

Numerical studies of transient planetary circulations in a wind-driven ocean

Summary The time-dependent primitive equations for a shallow homogeneous ocean with a free surface are solved for a bounded basin on the sphere, driven by a steady zonal wind stress and subject to lateral viscous dissipation. These are the vertically integrated equations for a free-surface model, and are integrated to 60 days from an initial state of rest by an explicit centered-difference method with zero-slip lateral boundary conditions. In a series of comparative numerical solutions it is shown that at least a 2-deg resolution is needed to resolve the western boundary currents adequately and to avoid undue distortion of the transient (Rossby waves. The -plane formulation is shown to be an adequate approximation for the mean circulation in the lower and middle latitudes, but noticeably intensifies the transports poleward of about 50 deg and both slows and distorts the transients in the central basin. The influence of the (southern) zonal boundary on the transport solutions is confined to the southernmost gyre, except in the region of the western boundary currents where its influence spreads to the northern edge of the basin by 30 days. The total boundary current transport is shown to be approximately proportional to the zonal width of the basin and independent of the basin's (uniform) depth, while the elevation of the free water surface is inversely proportional to the basin depth, in accordance with linear theory. The introduction of bottom friction has a marked damping effect on the transient Rossby waves, and also reduces the maximum boundary-current transport. The solutions throughout are approximately geostrophic and are only slightly nonlinear.The root-mean-square (rms) transport variability during the period 30 to 60 days is concentrated in the southwest portion of the basin through the reflection of the transient Rossby waves from the western shore and has a maximum corresponding to an rms current variability of about 3 cm sec^–1. The transport variabilities are about 10 percent of the mean zonal transport and more than 100 percent of the mean meridional transport over a considerable region of the western basin (outside the western boundary current regime). Some 99 percent of the total kinetic energy is associated with the zonal mean and standing zonal waves, which are also responsible for the bulk of the meridional transport of zonal angular momentum. Although the transient Rossby waves systematically produce a momentum flux convergence at the latitude of the maximum eastward current, much in the manner of their atmospheric counterparts, this is only a relatively small contribution to the zonal oceanic momentum balance; the bulk of the mean zonal stress is here balanced by a nearly stationary net pressure torque exerted against the meridional boundaries by the wind-raised water. In an ocean without such boundaries the role of the transient circulations may be somewhat more important.