METHODOLOGIES FOR CALIBRATION AND PREDICTIVE ANALYSIS OF A WATERSHED MODEL1

ABSTRACT: The use of a fitted parameter watershed model to address water quantity and quality management issues requires that it be calibrated under a wide range of hydrologic conditions. However, rarely does model calibration result in a unique parameter set. Parameter nonuniqueness can lead to predictive nonuniqueness. The extent of model predictive uncertainty should be investigated if management decisions are to be based on model projections. Using models built for four neighboring watersheds in the Neuse River Basin of North Carolina, the application of the automated parameter optimization software PEST in conjunction with the Hydrologic Simulation Program Fortran (HSPF) is demonstrated. Parameter nonuniqueness is illustrated, and a method is presented for calculating many different sets of parameters, all of which acceptably calibrate a watershed model. A regularization methodology is discussed in which models for similar watersheds can be calibrated simultaneously. Using this method, parameter differences between watershed models can be minimized while maintaining fit between model outputs and field observations. In recognition of the fact that parameter nonuniqueness and predictive uncertainty are inherent to the modeling process, PEST's nonlinear predictive analysis functionality is then used to explore the extent of model predictive uncertainty.     

Spatiotemporal variability of wet atmospheric nitrogen deposition to the Neuse River Estuary, North Carolina

Excessive nitrogen (N) loading to N-sensitive waters such as the Neuse River estuary (North Carolina) has been shown to promote changes in microbial and algal community composition and function (harmful algal blooms), hypoxia and anoxia, and fish kills. Previous studies have estimated that wet atmospheric deposition of nitrogen (WAD-N), as deposition of dissolved inorganic nitrogen (DIN: NO3-, NH3/NH4+) and dissolved organic nitrogen, may contribute at least 15% of the total externally supplied or "new" N flux to the coastal waters of North Carolina. In a 3-yr study from June 1996 to June 1999, we calculated the weekly wet deposition of inorganic and organic N at eleven sites on a northwest-southeast transect in the watershed. The annual mean total (wet DIN + wet organics) WAD-N flux for the Neuse River watershed was calculated to be 956 mg N/m2/yr (15026 Mg N/yr). Seasonally, the spring (March-May) and summer (June-August) months contain the highest total weekly N deposition; this pattern appears to be driven by N concentration in precipitation. There is also spatial variability in WAD-N deposition; in general, the upper portion of the watershed receives the lowest annual deposition and the middle portion of the watershed receives the highest deposition. Based on a range of watershed N retention and in-stream riverine processing values, we estimate that this flux contributes approximately 24% of the total "new" N flux to the estuary.     

Respiration and lipid content of the Arctic copepod<Emphasis Type="Italic"> Calanus hyperboreus</Emphasis> overwintering 1 m above the seafloor at 2,300 m water depth in the Fram Strait

In August 2000 high concentrations of the dominant herbivorous copepod Calanus hyperboreus were detected in the Arctic Fram Strait, west of Spitsbergen, 1 m above the seafloor at 2,290 m water depth. Individuals from that layer were sampled by a hyper-benthic net attached to the frame of an epi-benthic sledge. For comparison, the vertical distribution of C. hyperboreus in the water column was studied simultaneously by a multiple opening/closing net haul from 2,250 m depth to the surface. Maximum abundance was found close to the surface with 6.6 and 10.0 ind. m^?3 at 0–50 m and 50–100 m depth, respectively. However, the major fraction of the population (>40%) occurred between 1,000 and 1,500 m depth. In the deepest layer (2,000–2,250 m) abundance measured 2.2 ind. m^?3 and was twice as high as between 100 and 1,000 m depth. In comparison to individuals from surface waters, copepods from the hyper-benthic layer were torpid and did not react to mechanical stimuli. Stage CV copepodids and females from the deep sample contained 4–10% less lipid and showed significantly reduced respiration rates of 0.24 and 0.26 ml O₂ h^?1 g^?1 dry mass (DM) as compared to surface samples (0.49 and 0.43 ml O₂ h^?1 g^?1 DM). All these observations indicate that the hyper-benthic part of the population had already started a dormant overwintering phase at great depth. Based on the lipid deposits and energy demands, the potential maximum duration of the non-feeding dormant phase was estimated at 76–110 days for females and at 98–137 days for CV copepodids, depending on what indispensable minimum lipid content was assumed. In any case, the estimated times could not meet the necessary requirements for a starvation period of >6 months until the next phytoplankton bloom in the following spring. The ecological implications of these results are discussed with respect to the life cycle and eco-physiological adaptations of C. hyperboreus to its high-Arctic habitat.     

Global sensitivity analysis of a process-based model for ammonia emissions from manure storage and treatment structures

The development of process-based models to estimate ammonia emissions from animal feeding operations (AFOSs) is sought to replace costly and time-consuming direct measurements. Critical to process-based model development is conducting sensitivity analysis to determine the input parameters and their interactions that contribute most to the variance of the model output. Global and relative sensitivity analyses were applied to a process-based model for predicting ammonia emissions from the surface of anaerobic lagoons for treating and storing manure. The objectives were to compare global sensitivity analysis (GSA) to relative (local) sensitivity analysis (RSA) on a process-based model for ammonia emissions. Based on the first-order coefficient, both GSA and RSA showed the model input parameters in order of importance in process model for ammonia emissions from lagoon surfaces were: (i) pH, (ii) lagoon liquid temperature, (iii) wind speed above the lagoon surface, and (iv) the concentration of ammoniacal nitrogen in the lagoon. The GSA revealed that interactions between model parameters accounted for over two-thirds of the model variance, a result that cannot be achieved using traditional RSA. Also, the GSA showed that parameter interactions involving liquid pH had more impact on the model output variance than the single parameters: (i) temperature, (ii) wind speed, or (iii) total ammoniacal nitrogen. This study demonstrates that GSA provides a more complete analysis of model input parameters and their interactions on the model output compared to RSA. A comprehensive tutorial regarding the application of GSA to a process model is presented.     

Using long-term data for a reintroduced population to empirically estimate future consequences of inbreeding

Inbreeding depression is an important long-term threat to reintroduced populations. However, the strength of inbreeding depression is difficult to estimate in wild populations because pedigree data are inevitably incomplete and because good data are needed on survival and reproduction. Predicting future population consequences is especially difficult because this also requires projecting future inbreeding levels and their impacts on long-term population dynamics, which are subject to many uncertainties. We illustrate how such projections can be derived through Bayesian state-space modeling methods based on a 26-year data set for North Island Robins (Petroica longipes) reintroduced to Tiritiri Matangi Island in 1992. We used pedigree data to model increases in the average inbreeding level (F ) over time based on kinship of possible breeding pairs and to estimate empirically N_e/N (effective/census population size). We used multiple imputation to model the unknown components of inbreeding coefficients, which allowed us to estimate effects of inbreeding on survival for all 1458 birds in the data set while modeling density dependence and environmental stochasticity. This modeling indicated that inbreeding reduced juvenile survival (1.83 lethal equivalents [SE 0.81]) and may have reduced subsequent adult survival (0.44 lethal equivalents [0.81]) but had no apparent effect on numbers of fledglings produced. Average inbreeding level increased to 0.10 (SE 0.001) as the population grew from 33 (0.3) to 160 (6) individuals over the 25 years, giving a ratio of 0.56 (0.01). Based on a model that also incorporated habitat regeneration, the population was projected to reach a maximum of 331–1144 birds (median 726) in 2130, then to begin a slow decline. Without inbreeding, the population would be expected stabilize at 887–1465 birds (median 1131). Such analysis, therefore, makes it possible to empirically derive the information needed for rational decisions about inbreeding management while accounting for multiple sources of uncertainty.     

Empirical Trials of Plant Field Guides

We designed 3 image‐based field guides to tropical forest plant species in Ghana, Grenada, and Cameroon and tested them with 1095 local residents and 20 botanists in the United Kingdom. We compared users’ identification accuracy with different image formats, including drawings, specimen photos, living plant photos, and paintings. We compared users’ accuracy with the guides to their accuracy with only their prior knowledge of the flora. We asked respondents to score each format for usability, beauty, and how much they would pay for it. Prior knowledge of plant names was generally low (<22%). With a few exceptions, identification accuracy did not differ significantly among image formats. In Cameroon, users identifying sterile Cola species achieved 46–56% accuracy across formats; identification was most accurate with living plant photos. Botanists in the United Kingdom accurately identified 82–93% of the same Cameroonian species; identification was most accurate with specimens. In Grenada, users accurately identified 74–82% of plants; drawings yielded significantly less accurate identifications than paintings and photos of living plants. In Ghana, users accurately identified 85% of plants. Digital color photos of living plants ranked high for beauty, usability, and what users would pay. Black and white drawings ranked low. Our results show the potential and limitations of the use of field guides and nonspecialists to identify plants, for example, in conservation applications. We recommend authors of plant field guides use the cheapest or easiest illustration format because image type had limited bearing on accuracy; match the type of illustration to the most likely use of the guide for slight improvements in accuracy; avoid black and white formats unless the audience is experienced at interpreting illustrations or keeping costs low is imperative; discourage false‐positive identifications, which were common; and encourage users to ask an expert or use a herbarium for groups that are difficult to identify. Pruebas Empíricas de Guías de Campo de Plantas Hawthorne, Cable & Marshall     

Two Hundred Years of Local Avian Extinctions in Eastern Amazonia

Local, regional, and global extinctions caused by habitat loss, degradation, and fragmentation have been widely reported for the tropics. The patterns and drivers of this loss of species are now increasingly well known in Amazonia, but there remains a significant gap in understanding of long‐term trends in species persistence and extinction in anthropogenic landscapes. Such a historical perspective is critical for understanding the status and trends of extant biodiversity as well as for identifying priorities to halt further losses. Using extensive historical data sets of specimen records and results of contemporary surveys, we searched for evidence of local extinctions of a terra firma rainforest avifauna over 200 years in a 2500 km² eastern Amazonian region around the Brazilian city of Belém. This region has the longest history of ornithological fieldwork in the entire Amazon basin and lies in the highly threatened Belém Centre of Endemism. We also compared our historically inferred extinction events with extensive data on species occurrences in a sample of catchments in a nearby municipality (Paragominas) that encompass a gradient of past forest loss. We found evidence for the possible extinction of 47 species (14% of the regional species pool) that were unreported from 1980 to 2013 (80% last recorded between 1900 and 1980). Seventeen species appear on the International Union for Conservation of Nature Red List, and many of these are large‐bodied. The species lost from the region immediately around Belém are similar to those which are currently restricted to well‐forested catchments in Paragominas. Although we anticipate the future rediscovery or recolonization of some species inferred to be extinct by our calculations, we also expect that there are likely to be additional local extinctions, not reported here, given the ongoing loss and degradation of remaining areas of native vegetation across eastern Amazonia. Doscientos Años de Extinciones Locales de Aves en la Amazonia Oriental     

Defining and Evaluating the Umbrella Species Concept for Conserving and Restoring Landscape Connectivity

Conserving or restoring landscape connectivity between patches of breeding habitat is a common strategy to protect threatened species from habitat fragmentation. By managing connectivity for some species, usually charismatic vertebrates, it is often assumed that these species will serve as conservation umbrellas for other species. We tested this assumption by developing a quantitative method to measure overlap in dispersal habitat of 3 threatened species—a bird (the umbrella), a butterfly, and a frog—inhabiting the same fragmented landscape. Dispersal habitat was determined with Circuitscape, which was parameterized with movement data collected for each species. Despite differences in natural history and breeding habitat, we found substantial overlap in the spatial distributions of areas important for dispersal of this suite of taxa. However, the intuitive umbrella species (the bird) did not have the highest overlap with other species in terms of the areas that supported connectivity. Nevertheless, we contend that when there are no irreconcilable differences between the dispersal habitats of species that cohabitate on the landscape, managing for umbrella species can help conserve or restore connectivity simultaneously for multiple threatened species with different habitat requirements. Definición y Evaluación del Concepto de Especie Paraguas para Conservar y Restaurar la Conectividad de Paisajes