USGS Science⌕ Search

SEARCH · USGS Science

Results for “Special Study”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37Linked to original sources

Groundwater-quality and select quality-control data from the National Water-Quality Assessment Project, January through December 2015, and previously unpublished data from 2013 to 2014

Groundwater-quality data were collected from 502 wells as part of the National Water-Quality Assessment Project of the U.S. Geological Survey National Water-Quality Program and are included in this report. Most of the wells (500) were sampled from January through December 2015, and 2 of them were sampled in 2013. The data were collected from five types of well networks: principal aquifer study networks, which are used to assess the quality of groundwater used for public water supply; land-use study networks, which are used to assess land-use effects on shallow groundwater quality; major aquifer study networks, which are used to assess the quality of groundwater used for domestic supply; enhanced trends networks, which are used to evaluate the time scales during which groundwater quality changes; and vertical flow-path study networks, which are used to evaluate changes in groundwater quality from shallow to deeper depths. Groundwater samples were analyzed for a large number of water-quality indicators and constituents, including major ions, nutrients, trace elements, volatile organic compounds, pesticides, radionuclides, and some constituents of special interest (arsenic speciation, chromium [VI], and perchlorate). These groundwater-quality data, along with data from quality-control samples, are tabulated in this report and in an associated data release. Some data from environmental samples collected in 2013 and quality-control samples collected in 2014 also are included in the associated data release; these data are associated with networks described in this report and have not been published previously.

Data Series↗

Factors affecting plant diversity during post-fire recovery and succession of mediterranean-climate shrublands in California, USA

Plant community diversity, measured as species richness, is typically highest in the early post-fire years in California shrublands. However, this generalization is overly simplistic and the present study demonstrates that diversity is determined by a complex of temporal and spatial effects. Ninety sites distributed across southern California were studied for 5 years after a series of fires. Characteristics of the disturbance event, in this case fire severity, can alter post-fire diversity, both decreasing and increasing diversity, depending on life form. Spatial variability in resource availability is an important factor explaining patterns of diversity, and there is a complex interaction between landscape features and life form. Temporal variability in resource availability affects diversity, and the diversity peak in the immediate post-fire year (or two) appears to be driven by factors different from subsequent diversity peaks. Early post-fire diversity is influenced by life-history specialization, illustrated by species that spend the bulk of their life cycle as a dormant seed bank, which is then triggered to germinate by fire. Resource fluctuations, precipitation in particular, may be associated with subsequent post-fire diversity peaks. These later peaks in diversity comprise a flora that is compositionally different from the immediate post-fire flora, and their presence may be due to mass effects from population expansion of local populations in adjacent burned areas. ?? 2005 Blackwell Publishing Ltd.

Diversity and Distributions↗

What can seismology say about hotspots?

Seismological methods offer the highest-resolution views of the structure of the mantle. Since deep mantle plumes were proposed to explain melting anomalies (“hotspots”), increasingly powerful seismologic studies have sought to detect them, but so far without definitive success. This paper summarizes the relevant seismological methods and results for Earth scientists who are not seismologists.

Special Papers of the Geological Society of Americ↗

Location-only and use-availability data: analysis methods converge

This Special Feature arose from a session on a topic of the same name that took place during The Wildlife Society meeting in Kona, Hawaii, from 5 to 10 November, 2011. The purpose of that session and this Special Feature is to compare methods for predictive modelling of species geographical distributions and the modelling of habitat (resource) selection by animals. The predictive modelling of species geographical distributions and the modelling of habitat selection based on the environmental conditions at sites where animals are known to occur are essentially the same problem. Presence-only and used-available data both consist of a sample of locations with known presence of a species or an individual. A separate sample of locations from a study area, with unknown presence (pseudo-absence), is also assumed to exist. The probability or relative probability of presence of a species or individual is modelled and estimated across a certain time implicitly defined by the sampling mechanism, for example, by the time period during which museum specimens or radiotelemetry data were collected. A number of modelling methods have appeared in the literature over the last couple of decades. Many of these methods were made feasible by the availability of geographical information systems (GIS), global positioning system (GPS) radiotelemetry and public online data access initiatives (e.g. global biodiversity information facility). The papers in this Special Feature are intended to present the state of the methodological art in their subject area, with particular attention paid to contrasting the advantages and disadvantages of alternative methods of analysis for data.

Journal of Animal Ecology↗

Ecological islands: Conserving biodiversity hotspots in a changing climate

For decades, botanists have recognized that rare plants are clustered into ecological “islands”: small and isolated habitat patches produced by landscape features such as sinkholes and bedrock outcrops. Insular ecosystems often provide unusually stressful microhabitats for plant growth (eg because of thin soils, high temperatures, extreme pH, or limited nutrients) to which rare species are specially adapted. Climate-driven shifts to these stressors may undermine the competitive advantage of stress-adapted species, allowing them to be displaced by competitors, or may overwhelm their coping strategies altogether. Special features of insular ecosystems—such as extreme habitat fragmentation and association with unusual landscape features—may also affect their climate sensitivity and adaptive capacity. To predict and manage these changes, a simple conceptual framework is presented based on a synthesis of over 300 site-level studies. Using this framework, conservation efforts can leverage existing ecological knowledge to anticipate changes in particular microhabitats and design targeted strategies for conserving rare species.

Frontiers in Ecology and the Environment↗

Benthic boundary layer processes in the Lower Florida Keys

This special issue of Geo-Marine Letters, "Benthic Boundary Layer Processes in the Lower Florida Keys," includes 12 papers that present preliminary results from the Key West Campaign. The Dry Tortugas and Marquesas Keys test sites were selected by a group of 115 scientists and technicians to study benthic boundary layer processes in a carbonate environment controlled by bioturbation and biogeochemical processes. Major activities included remote sediment classification; high-frequency acoustic scattering experiments; sediment sampling for radiological, geotechnical, biological, biogeochemical, physical, and geoacoustic studies; and hydrodynamic studies using an instrumented tetrapod. All these data are being used to improve our understanding of the effects of environmental processes on sediment structure and behavior.

Florida↗

Assessing the impacts of future climate conditions on the effectiveness of winter cover crops in reducing nitrate loads into the Chesapeake Bay Watershed using SWAT model

Winter cover crops (WCCs) have been widely implemented in the Coastal Plain of the Chesapeake Bay watershed (CBW) due to their high effectiveness at reducing nitrate loads. However, future climate conditions (FCCs) are expected to exacerbate water quality degradation in the CBW by increasing nitrate loads from agriculture. Accordingly, the question remains whether WCCs are sufficient to mitigate increased nutrient loads caused by FCCs. In this study, we assessed the impacts of FCCs on WCC nitrate reduction efficiency on the Coastal Plain of the CBW using Soil and Water Assessment Tool (SWAT) model. Three FCC scenarios (2085 – 2098) were prepared using General Circulation Models (GCMs), considering three Intergovernmnental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) greenhouse gas emission scenarios. We also developed six representative WCC implementation scenarios based on the most commonly used planting dates and species of WCCs in this region. Simulation results showed that WCC biomass increased by ~ 58 % under FCC scenarios, due to climate conditions conducive to the WCC growth. Prior to implementing WCCs, annual nitrate loads increased by ~ 43 % under FCC scenarios compared to the baseline scenario (2001 – 2014). When WCCs were planted, annual nitrate loads were substantially reduced by ~ 48 % and WCC nitrate reduction efficiency water ~ 5 % higher under FCC scenarios relative to the baseline. The increase rate of WCC nitrate reduction efficiency varied by FCC scenarios and WCC planting methods. As CO2 concentration was higher and winters were warmer under FCC scenarios, WCCs had greater biomass and therefore showed higher nitrate reduction efficiency. In response to FCC scenarios, the performance of less effective WCC practices (e.g., barley, wheat, and late planting) under the baseline indicated ~ 14 % higher increase rate of nitrate reduction efficiency compared to ones with better effectiveness under the baseline (e.g., rye and early planting), due to warmer temperatures. According to simulation results, WCCs were effective to mitigate nitrate loads accelerated by FCCs and therefore the role of WCCs in mitigating nitrate loads is even more important in the given FCCs.

Chesapeake Bay↗

Theory, modelling and calibration of passive samplers used in water monitoring: Chapter 7

This chapter discusses contaminant uptake by a passive sampling device (PSD) that consists of a central sorption phase, surrounded by a membrane. A variety of models has been used over the past few years to better understand the kinetics of contaminant transfer to passive samplers. These models are essential for understanding how the amounts of absorbed contaminants relate to ambient concentrations, as well as for the design and evaluation of calibration experiments. Models differ in the number of phases and simplifying assumptions that are taken into consideration, such as the existence of (pseudo-) steady-state conditions, the presence or absence of linear concentration gradients within the membrane phase, the way in which transport within the WBL is modeled and whether or not the aqueous concentration is constant during the sampler exposure. The chapter introduces the basic concepts and models used in the literature on passive samplers for the special case of triolein-containing semipermeable membrane devices (SPMDs). These can easily be extended to samplers with more or with less sorption phases. It also discusses the transport of chemicals through the various phases constituting PSDs. the implications of these models for designing and evaluating calibration studies have been discussed.

Book chapter↗

Stable isotope evolution and paleolimnology of ancient Lake Creede

The lacustrine carbonate and travertine (tufa) deposits of ancient Lake Creede preserve a remarkable record of the isotopic evolution of the lake. That record indicates that the δ18O of the lake water, and by analogy its salinity, evolved through evaporation. Limited and less reliable data on hydrous minerals and fluid inclusions in early diagenetic carbonates indicate that the δD of the lake waters also evolved through evaporation. The isotope data place restrictions on models of the physical limnology of the lake and its evolution. The closed-basin Lake Creede formed shortly after collapse of the 26.9 Ma Creede caldera. Throughout most of its history it occupied the northern three quarters of the moat between the resurgent dome and wall of the caldera. The Creede Formation was deposited in the basin, dominantly as lacustrine sediments. Travertine mounds interfinger with Creede Formation sediments along the inner and outer margins of the lake basin. An estimated one-half of the original thickness of the Creede Formation has been lost mainly to erosion although scattered remnants of the upper portion remain on the caldera walls. Two diamond core holes (CCM-1 and CCM-2) sampled the uneroded portion of the Creede Formation as part of the U.S. Continental Drilling Program. Volcaniclastic material, including tuff units deposited directly into the lake and ash washed in from the watershed, compose the main lithologies of the Creede Formation. These volcaniclastic strata were produced by episodic ring-fracture volcanism. Lacustrine carbonates make up about 15% of the section sampled by drill core. They occur as 1 mm to 2 cm low-Mg calcite laminae alternating with siliciclastic laminae in scattered intervals throughout the preserved section. The carbonate laminae are accumulations of 5–20 µm crystallites (microsparites) and brine shrimp fecal pellets (peloids) composed mainly of microsparite particles. Low-Mg calcite also occurs as an early diagenetic replacement of gypsum or ikaite (CaCO3 ·6H2O) crystals grown displacively in the muds and silts near the water-sediment interface (rice grains). Other studies indicate that aragonite was the original CaCO3 precipitate forming the microsparite and peloidal laminae and that it converted to calcite during burial diagenesis. Samples from CCM-2 and nearby outcrop do not appear to have undergone significant isotope exchange during recrystallization. Samples from CCM-1 and nearby outcrop, however, appear to have undergone extensive oxygen isotope exchange with meteoric water-dominated fluids possibly during a local 17.6 Ma hydrothermal event. The δ18O-δ13C data set produced by microsampling of individual carbonate lamellae and rice grains is exceptional in several aspects and provides important clues concerning the evolution of limnologic structure of the lake and its chemical and isotopic composition. Travertine and ikaite pseudomorphs in travertine deposits extend the record an additional 330 m above the collar of CCM-2. The δ18O values on CCM-2 samples range from 10.4‰ to 37.3‰ and δ13C values range from –10.8‰ to 9.6‰. The data fall into two distinct groups, a covariant group and an invariant group. The covariant group shows a strong negative covariance and a large range of δ18O and δ13C values. The negative covariance is opposite that normally reported for lacustrine carbonates. The large range of δ18O and δ13C values requires that the carbonates precipitated from waters have a large range of temperature and carbon and oxygen isotopic composition. The invariant group has a narrow range of large δ18O values (35‰ ± 2‰) and a wide range of δ13C values (–10.8‰ to 9.6‰), indicating precipitation from waters with a narrow range of temperature and δ18O but a wide range in δ13C of aqueous carbon. The ranges of isotope values for microsparite and peloid samples are virtually identical; two-thirds are in the covariant group. By contrast, the values for almost all rice grain samples are in the invariant group. The range in δ18O for all samples reflects precipitation from waters having varying proportions of deep, cold evaporated lake water and shallow, warmer meteoric water. The range for δ13C reflects varying proportions of organic carbon and carbon of volcanic or atmospheric origin, probably dominantly volcanic, in the aqueous carbon. Changes in the detailed carbon-oxygen isotope systematics with stratigraphic position define three periods of isotopic evolution of Lake Creede. Period I is represented by the lowest ~200 m of Creede Formation core in CCM-2. Analyses of individual microsparite and peloidal carbonate laminae within single thin sections of samples from that interval are tightly grouped. The data set as a whole shows a negative covariance. Rice grains are not found in this interval. Period II is represented by the succeeding 120 m of core in CCM-2. In that interval, δ13C-δ18O values for individual microsparite and peloidal carbonate laminae within single thin sections show strong negative covariance, and the set of values for the entire interval also shows strong negative covariance. Rice grains occur near the top of the interval. Period III is represented by the upper 225 m of CCM-2 core. In this interval, rice grains are abundant and δ13C-δ18O values for microsparite and peloidal laminae as well as rice grains fall in the invariant group. During Period I the lake was well mixed and the oxygen isotopic composition of the lake in the productive zone was only slightly influenced by short-term (e.g., annual) variations in the water budget of the lake. In Period II the lake was stratified, possibly with annual overturn. The productive zone included the mixolimnion and the isotopic composition of the microsparites and peloids reflected mixtures of shallow surface (meteoric) water containing volcanic or atmospheric CO2 (epilimnion) and cold underlying waters, the oxygen isotopic compositions of which had evolved through evaporation and were dominated by CO2 produced by the oxidation of organic matter (hypolimnion). During Period III the lake remained stratified. The productive zone was in the hypolimnion, probably due to a thinning of the epilimnion resulting from an increase in the surface area of the lake or a decrease in input waters reflecting a climate change. An upsection increase in values of δ18O for the heaviest samples during Periods I and II indicates evaporative concentration of 18O and, by analogy, salinity in the hypolimnion. The δD-δ18O evolution of the lake is inferred on theoretical evaporation trends, comparison to Mono Lake, and measurement of the δD in fluid inclusions in a calcite pseudomorph after ikaite. The δD-δ18O composition of the lake water followed a curved path that eventually hooked over at a nearly constant δ18O value for the lake of 2‰ ± 2‰ Travertine (tufa) mounds formed along the inner and outer margins of the lake in a zone of mixing of warm, volcanic CO2-bearing, meteoric waters and lake water. Ikaite crystals formed on the mounds from unmixed saline lake water, probably below the thermocline. As the position of the thermocline deepened, likely following the spring runoff, the ikaite was replaced by calcite and the resulting “pearls” were covered with travertine deposited from mixed meteoric and lake waters. The upsection increase in δ18O values of the carbonates, the long period of invariance of large δ18OH2O values, the presence of brine shrimp fecal pellets, and the inferred hooked δD-δ18O path are consistent with evidence from other studies that Lake Creede obtained significant salinity rather early in its history and certainly by the time the lake became permanently stratified.

Colorado↗

Determination of reaeration coefficients for Ohio streams

The hydrocarbon-gas tracer technique was used to determine reaeration coefficients on 30 reaches of Ohio streams. The studies were done from September 1979 through August 1982 to determine the reaeration coefficients for the individual reaches and to develop general equation that could be used to estimate the coefficients. Multiple linear regression was used to determine relationships among the reaeration coefficients and physical stream characteristics. Four special equation based on the general equation were developed from various combinations of discharge, slope, width, depth, measured velocity, and estimated velocity. The standard errors of estimate for these equations ranged from 37 to 47 percent. The variables that resulted in the lowest standard error of estimate were discharge, slope, width, depth, and measured velocity. The most significant variables were depth and velocity.

Water-Resources Investigations Report↗

System designed for issuing landslide alerts in the San Francisco Bay area

A system for forecasting landslides during major storms has been developed for the San Francisco Bay area by the U.S Geological Survey and was successfully tested during heavy storms in the bay area during February 1986. Based on the forecasts provided by the USGS, the National Weather Service (NWS) included landslide warnings in its regular weather forecasts or in special weather statements transmitted to local radio and television stations and other news media. USGS scientists said the landslide forecasting and warning system for the San Francisco Bay area can be used as a prototype in developing similar systems for other parts of the Nation susceptible to landsliding. Studies show damage from landslides in the United States averages an estimated $1.5 billion per year.

California↗

Reconnaissance Seismic Refraction Studies at Calico Hills, Wahmonie, and Yucca Mountain, Southwest Nevada Test Site, Nye County, Nevada

Reconnaissance refraction surveys consisting off a total of 5 spreads were conducted in the Calico Hills, Wahmonie and Yucca Mountain areas, southwestern Nevada Test Site (NTS). Data from Calico Hills and Wahmonie are generally high in quality; data from Yucca Mountain are for the most part low in quality. At Calico Hills and Wahmonie, special attention was focused on the possible occurrence of a major intrusive body at depth. At Calico Hills this occurrence is supported by an inferred dome-shaped velocity interface. possibly associated with the roof of an altered phase of argillite. However, if an intrusive body is present, its top must be buried deeper than 3 km or it must be so pervasively altered that its velocity is similar to that of the calcareous argillite encountered at the bottom of drill hole DE 25a-3. At Wahmonie, the seismic data suggest the occurrence of a massive lenticular unit within 60 m of the ground surface, probably consisting of argillite but possibly consisting of intensively altered intrusive rock. At Yucca Mountain, preliminary interpretations of the most reliable data suggest the occurrence of a major, steeply inclined velocity interface 500 m from the southwest end of the Yucca C spread. This interface may represent a major fault or erosional feature separating the Topopah Spring and Tiva Canyon Members with Paintbrush Tuff at depth. This interface is 800 m east of a previously mapped fault. On the basis of poor-quality data obtained at Yucca Mountain, the subsurface velocity distribution appears to be complex. For example, one spread near drill hole UE25 a-I suggests not only a much thicker section of Tiva but also that this material is down thrown in the valley. This may suggest faulting with throws exceeding 100 meters or an equivalent erosional feature.

Open-File Report↗

Wetlands under global change

Wetlands are among the ecosystem types most threatened by global change, including both climate change and other anthropogenic factors such as sea level rise, urban development, deforestation, agricultural land use, drainage, levees, tidal flow restrictions, pollution, eutrophication, and fires. Wetlands not only store disproportionate amounts of carbon compared to other terrestrial ecosystems, but they lie at the terrestrial-aquatic interface crucial to understanding landscape and global scale biogeochemical cycles. In this chapter, we focus on the major global change factors affecting wetlands and the responses of different wetland types to those global change factors. Special attention is given to direct responses to increasing atmospheric carbon dioxide levels. Because of their hydrological connections and placement at the terrestrial-aquatic interface, the conservation of wetlands involves accounting for uncertainties related to interacting stressors. While the past decades have seen many important experimental and observational studies of wetland responses to global change factors, large uncertainties remain, especially within tropical regions where even the basic extent of wetland ecosystems is not well documented.

Book chapter↗

The Driftwood anticline

Recent interest in the anticlines near the headwaters of the Utukok River has prompted the issuance of two special reports. This, the first of the two reports, deals with the Driftwood anticlines, the second will consider the Awuna anticline and the western part of the Carbon Creek anticline. At present, only trimetrogon photography is available in this area. This limits the accuracy and completeness of the photographic studies. Vertical photos were flown this summer but they will probably not be available for study before the November meeting. At such time as they are retrieved a more thorough photo investigation will be made.

Alaska↗

Palaeoseismology of the Vilariça segment of the Manteigas-Bragança fault in northeastern Portugal

The Manteigas-Bragança fault is a major, 250-km-long, NNE-striking, sinistral strike-slip structure in northern Portugal. This fault has no historical seismicity for large earthquakes, although it may have generated moderate (M5+) earthquakes in 1751 and 1858. Evidence of continued left horizontal displacement is shown by the presence of Cenozoic pull-apart basins as well as late Quaternary stream deflections. To investigate its recent slip history, a number of trenches were excavated at three sites along the Vilariça segment, north and south of the Douro River. At one site at Vale Meão winery, the occurrence of at least two and probably three events in the past 14.5 ka was determined, suggesting an average return period of about 5–7 ka. All three events appear to have occurred as a cluster in the interval between 14.5 and 11 ka, or shortly thereafter, suggesting a return period of less than 2 ka between events within the cluster. In the same area, a small offset rill suggests 2–2.5 m of slip in the most recent event and about 6.1 m after incision below a c . 16 ka alluvial fill event along the Douro River. At another site along the Vilariça River alluvial plain, NE of the Vale Meão site, several trenches were excavated in late Pleistocene and Holocene alluvium, and exposed the fault displacing channel deposits dated to between 18 and 23 ka. In a succession of closely spaced parallel cuts and trenches, the channel riser was traced into and across the fault to resolve c . 6.5 m of displacement after 18 ka and c . 9 m of slip after c . 23 ka. These observations yield a slip rate of 0.3–0.5 mm/a, which is consistent with earlier estimates. Combining the information on timing at Vale Meão winery and displacement at Vilariça argues for earthquakes in the M7+ range, with coseismic displacements of 2–3 m. This demonstrates that there are potential seismic sources in Portugal that are not associated with the 1755 Lisbon earthquake or the Tagus Valley, and, although rare, large events on the Vilariça fault could be quite destructive for the region. This work provides an analogue for the study of active faulting in intracontinental settings and supports the view that earthquakes within intracontinental settings tend to cluster in time. In addition, this study highlights the usefulness and application of multiple field, remote sensing and geochronological techniques for seismic hazard mitigation.

Manteigas - Braganca fault↗

Preliminary report on the geology and ground-water supply of the Newark, New Jersey, area

In the Newark area, ground water is used chiefly for industrial cooling, air-conditioning, general processing, and for sanitary purposes. A small amount is used in the manufacture of beverages. Total ground-water pumpage in Newark is estimated at not less than 20,000,000 gallons daily. The Newark area is underlain by formations of Recent, Pleistocene and Triassic age, and the geology and hydrologic properties of these formations are discussed. Attention is called to the important influence of a buried valley in the rock floor beneath the Newark area on the yield of wells located within it. Data on the fluctuation of the water levels and the variation in pumpage are presented, and their significance discussed. The results of a pumping test made during the investigation were inconclusive. The beneficial results of artificially recharging the aquifers in one part of the area are described. The intrusion of salt water into certain parts of the ground-water body is described and graphically portrayed by a map showing the chloride concentration of the ground water in various parts of the City. Insofar as available data permit, the chemical quality of the ground water is discussed and records are given of the ground-water temperatures in various parts of the City. There has been marked lowering of the water table in the eastern part of the area, accompanied by salt water intrusion, indicating that the safe yield of the formations in this part of Newark has probably been exceeded. It is recommended that the study of the ground-water resources of this area be continued, and that artificial recharging of the aquifers be increased over as wide an area as possible.

New Jersey↗

Adaptation of Selenastrum capricornutum (Chlorophyceae) to copper

Selenastrum capricornutum Printz, growing in a chemically defined medium, was used as a model for studying adaptation of algae to a toxic metal (copper) ion. Cells exhibited lag‐phase adaptation to 0.8 μM total Cu (10 −12 M free ion concentration) after 20 generations of Cu exposure. Selenastrum adapted to the same concentration when Cu was gradually introduced over an 8‐h period using a specially designed apparatus that provided a transient increase in exposure concentration. Cu adaptation was not attributable to media conditioning by algal exudates. Duration of lag phase was a more sensitive index of copper toxicity to Selenastrum that was growth rate or stationary‐phase cell density under the experimental conditions used. Chemical speciation of the Cu dosing solution influenced the duration of lag phase even when media formulations were identical after dosing. Selenastrum initially exposed to Cu in a CuCl 2 injection solution exhibited a lag phase of 3.9 d, but this was reduced to 1.5 d when a CuEDTA solution was used to achieve the same total Cu and EDTA concentrations. Physical and chemical processes that accelerated the rate of increase in cupric ion concentration generally increased the duration of lag phase.

Environmental Toxicology and Chemistry↗

Global modeling of land water and energy balances. Part I: The land dynamics (LaD) model

A simple model of large-scale land (continental) water and energy balances is presented. The model is an extension of an earlier scheme with a record of successful application in climate modeling. The most important changes from the original model include 1) introduction of non-water-stressed stomatal control of transpiration, in order to correct a tendency toward excessive evaporation: 2) conversion from globally constant parameters (with the exception of vegetation-dependent snow-free surface albedo) to more complete vegetation and soil dependence of all parameters, in order to provide more realistic representation of geographic variations in water and energy balances and to enable model-based investigations of land-cover change; 3) introduction of soil sensible heat storage and transport, in order to move toward realistic diurnal-cycle modeling; 4) a groundwater (saturated-zone) storage reservoir, in order to provide more realistic temporal variability of runoff; and 5) a rudimentary runoff-routing scheme for delivery of runoff to the ocean, in order to provide realistic freshwater forcing of the ocean general circulation model component of a global climate model. The new model is tested with forcing from the International Satellite Land Surface Climatology Project Initiative I global dataset and a recently produced observation-based water-balance dataset for major river basins of the world. Model performance is evaluated by comparing computed and observed runoff ratios from many major river basins of the world. Special attention is given to distinguishing between two components of the apparent runoff ratio error: the part due to intrinsic model error and the part due to errors in the assumed precipitation forcing. The pattern of discrepancies between modeled and observed runoff ratios is consistent with results from a companion study of precipitation estimation errors. The new model is tuned by adjustment of a globally constant scale factor for non-water-stressed stomatal resistance. After tuning, significant overestimation of runoff is found in environments where an overall arid climate includes a brief but intense wet season. It is shown that this error may be explained by the neglect of upward soil water diffusion from below the root zone during the dry season. With the exception of such basins, and in the absence of precipitation errors. It is estimated that annual runoff ratios simulated by the model would have a root-mean-square error of about 0.05. The new model matches observations better than its predecessor, which has a negative runoff bias and greater scatter.

Journal of Hydrometeorology↗