PREDICTION OF MEAN ANNUAL FLOWS IN NORTH AND CENTRAL FLORIDA1

ABSTRACT: A frequency analysis approach for the prediction of flow characteristics at ungaged locations is applied to a region of high annual precipitation and low topography in north and central Florida. Stationary time series of annual flows are fitted with the lognormal distribution and estimated parameters of the distribution are fitted by third order trend surfaces. These explain 65 and 74 percent of the observed variances in the mean and standard deviation, respectively. Predictions of parameters are then made for several locations previously unused in the study and they are used to estimate the return periods of various flows from the lognormal distribution. Application of the Kolmogorov-Smirnov goodness-of-fit test suggests that only one of the five test stations can be considered significantly different from the observed data, confirming the applicability of this technique.     

Mass conservation and atmospheric dynamics in the Regional Atmospheric Modeling System (RAMS)

This paper examines the spatio-temporal patterns of atmospheric carbon dioxide transport predicted by the Regional Atmospheric Modeling System (RAMS). Forty-eight hour simulations over northern New England incorporating a simple representation of the diurnal summertime surface carbon dioxide forcing arising from biological activity indicate that, in its native formulation, RAMS exhibits a significant degree of mass non-conservation. Domain-wide rates of non-physical mass gain and mass loss are as large as three percent per day which translates into approximately eleven parts per million per day for carbon dioxide — enough to rapidly dilute the signature of carbon dioxide fluxes arising from biological activity. Analysis shows that this is due to the approximation used by RAMS to compute the Exner function. Substitution of the exact, physically complete equation improves mass conservation by two orders of magnitude. In addition to greatly improving mass conservation, use of the complete Exner function equation has a substantial impact on the spatial pattern of carbon dioxide predicted by the model, yielding predictions differing from a conventional RAMS simulation by as much as forty parts per million. Such differences have important implications both for comparisons of modeled atmospheric carbon dioxide concentrations to observations and for carbon dioxide inversion studies, which use estimates of atmospheric transport of carbon dioxide in conjunction with measurements of atmospheric carbon dioxide concentrations to infer the spatio-temporal distribution of surface carbon dioxide fluxes. Furthermore, use of the complete Exner function equation affects the vertical velocity and water mixing ratio fields, causing significant changes in accumulated precipitation over the region.     

An Approach for the Aggregation of Aerodynamic Surface Parameters in Calculating the Turbulent Fluxes over Heterogeneous Surfaces in Atmospheric Models

Experimental evidence indicates that there is a significant departure of the wind profile above the underlying surface consisting of patches of solid and liquid parts, and plant communities with different morphological from that predicted by the logarithmic relationship, which gives the values larger than those observed. This situation can seriously affect the transfer of momentum, heat and water vapor from the surface fluxes into the atmosphere.The object of this paper is to generalize the calculation of the exchange of momentum between the atmosphere and a very heterogeneous surface, find a general equation for the wind speed profile in a roughness sublayer under neutral conditions, and, then, derive aggregated roughness length and displacement height over the grid cell. The suggested expression for the wind profile is compared with some earlier approaches, using a common parameterization of aerodynamic parameters over the grid cell, and the observations obtained at an experimental site in Philadelphia, PA.     

GIS and geostatistics: Essential partners for spatial analysis

Initially, geographical information systems (GIS) concentrated on two issues: automated map making, and facilitating the comparison of data on thematic maps. The first required high quality graphics, vector data models and powerful data bases, the second is based on grid cells that can be manipulated by suites of mathematical operators collectively termed map algebra. Both kinds of GIS are widely available and are taught in many universities and technical colleges. After more than 20 years of development, most standard GIS provide both kinds of functionality and good quality graphic display, but until recently they have not included the methods of statistics and geostatistics as tools for spatial analysis. Recently, standard statistical packages have been linked to GIS for both exploratory data analysis and statistical analysis and hypothesis testing. Standard statistical packages include methods for the analysis of random samples of cases or objects that are not necessarily co-located in space—if the results of statistical analysis display a spatial pattern then that is because the underlying data also share that pattern. Geostatistics addresses the need to make predictions of sampled attributes (i.e., maps) at unsampled locations from sparse, often expensive data. To make up for lack of hard data geostatistics has concentrated on the development of powerful methods based on stochastic theory. Though there have been recent moves to incorporate ancillary data in geostatistical analyses, insufficient attention has been paid to using modern methods of data display for the visualization of results. GIS can serve geostatistics by aiding geo-registration of data, facilitating spatial exploratory data analysis, providing a spatial context for interpolation and conditional simulation, as well as providing easy-to-use and effective tools for data display and visualization. The value of geostatistics for GIS lies in the provision of reliable interpolation methods with known errors, methods of upscaling and generalization, and for supplying multiple realizations of spatial patterns that can be used in environmental modeling. These stochastic methods are improving understanding of how errors in models of spatial processes accrue from errors in data or incompleteness in the structure of the models. New developments in GIS, based on ideas taken from map algebra, cellular automata and image analysis are providing high level programming languages for modeling dynamic processes such as erosion or the development of alluvial fans and deltas. Research has demonstrated that these models need stochastic inputs to yield realistic results. Non-stochastic tools such as fuzzy subsets have been shown to be useful for spatial analysis when probabilistic approaches are inappropriate or impossible. The conclusion is that in spite of differences in history and approach, the linkage of GIS, statistics and geostatistics provides a powerful, and complementary suite of tools for spatial analysis in the agricultural, earth and environmental sciences.     

A method for assessing generalised nitrogen leaching estimates for agricultural land

The purpose of this study was to develop a method for assessing generalised N leaching estimates from large areas of agricultural land. The system developed was based on calculating a number of N leaching estimates for different typical cropping situations. The estimates were normalised with respect to varying weather conditions and crop production. The different cropping situations were described by setting up a matrix consisting of crucial factors influencing leaching such as soils, crops and climate. Nitrogen leaching was then estimated for a number of combinations of these factors. Calculations were made for three different regions where all the major crops were cultivated on soils with seven different textures and four different organic-N classes and two fertilisation regimes. The three regions are representative of climates and agricultural practices in some of the major agricultural areas in Sweden. The model used was the SOILN model. Leaching of nitrogen from the root zone showed large variations. The range was from 1 to 50 kg ha⁻¹ for different soils and crops when only fertiliser N was applied. Leaching varied both due to different climates and differences in cultivation practices between the regions. Leaching decreased in a south-north gradient. Leaching increased as a result of greater mineralisation when the organic matter content in the soils was increased, leaching was less from soils with a high clay content and was very small for the heavy clay soil. This revised version was published online in July 2006 with corrections to the Cover Date.     

The role of life histories and trophic interactions in population recovery

Factors affecting population recovery from depletion are at the focus of wildlife management. Particularly, it has been debated how life‐history characteristics might affect population recovery ability and productivity. Many exploited fish stocks have shown temporal changes towards earlier maturation and reduced adult body size, potentially owing to evolutionary responses to fishing. Whereas such life‐history changes have been widely documented, their potential role on stock's ability to recover from exploitation often remains ignored by traditional fisheries management. We used a marine ecosystem model parameterized for Southeastern Australian ecosystem to explore how changes towards “faster” life histories might affect population per capita growth rate r. We show that for most species changes towards earlier maturation during fishing have a negative effect (3–40% decrease) on r during the recovery phase. Faster juvenile growth and earlier maturation were beneficial early in life, but smaller adult body sizes reduced the lifetime reproductive output and increased adult natural mortality. However, both at intra‐ and inter‐specific level natural mortality and trophic position of the species were as important in determining r as species longevity and age of maturation, suggesting that r cannot be predicted from life‐history traits alone. Our study highlights that factors affecting population recovery ability and productivity should be explored in a multi‐species context, where both age‐specific fecundity and survival schedules are addressed simultaneously. It also suggests that contemporary life‐history changes in harvested species are unlikely to increase their resilience and recovery ability.     

Strategic Rat Control for Restoring Populations of Native Species in Forest Fragments

Forest fragments have biodiversity value that may be enhanced through management such as control of non‐native predators. However, such efforts may be ineffective, and research is needed to ensure that predator control is done strategically. We used Bayesian hierarchical modeling to estimate fragment‐specific effects of experimental rat control on a native species targeted for recovery in a New Zealand pastoral landscape. The experiment was a modified BACI (before‐after‐control‐impact) design conducted over 6 years in 19 forest fragments with low‐density subpopulations of North Island Robins (Petroica longipes). The aim was to identify individual fragments that not only showed clear benefits of rat control, but also would have a high probability of subpopulation growth even if they were the only fragment managed. We collected data on fecundity, adult and juvenile survival, and juvenile emigration, and modeled the data in an integrated framework to estimate the expected annual growth rate (λ) of each subpopulation with and without rat control. Without emigration, subpopulation growth was estimated as marginal (λ = 0.95–1.05) or negative (λ = 0.74–0.90) without rat control, but it was estimated as positive in all fragments (λ = 1.4–2.1) if rats were controlled. This reflected a 150% average increase in fecundity and 45% average increase in adult female survival. The probability of a juvenile remaining in its natal fragment was 0.37 on average, but varied with fragment connectivity. With juvenile emigration added, 6 fragments were estimated to have a high (>0.8) probability of being self‐sustaining (λ > 1) with rat control. The key factors affecting subpopulation growth rates under rat control were low connectivity and stock fencing because these factors were associated with lower juvenile emigration and higher fecundity, respectively. However, there was also substantial random variation in adult survival among fragments, illustrating the importance of hierarchical modeling for fragmentation studies. Control Estratégico de Ratas para Restaurar Poblaciones de Especies Nativas en Fragmentos de Bosque     

Social values and biodiversity conservation in a dynamic world

Understanding what shape values (which ultimately shape human behavior) will help improve the effectiveness of conservation solutions that depend on public support. To contribute to this understanding, we investigated the influence of societal‐level changes, such as modernization, on values in a multilevel framework. We collected survey responses (n = 4183) to questionnaires mailed to a random selection of households within each county in Washington (U.S.A.) (response rate 32%). We used multilevel modeling to determine the relationship between modernization (e.g., county‐level urbanization, wealth, and education) and wildlife value orientations (values that shape thought about wildlife) while controlling for individual‐level sociodemographics. We then explored how values influence conservation support at different levels (e.g., individual and county) and how values explain conservation support in a case study of public responses to wolf (Canis lupis) recovery. We found positive associations between county‐level examples of modernization and mutualism (a wildlife value orientation that prioritizes the perceived needs of wildlife) independent of a respondent's sociodemographics, and negative associations between modernization and domination (a wildlife value orientation that prioritizes human needs). Our results suggest that context has an additive impact on one's values; certain locations exhibited domination values, whereas others exhibited a mix of value types. This finding is important because actions that restrict human interests to promote biodiversity were negatively associated with domination and positively associated with mutualism. In the wolf case study, mutualism was strongly correlated with less social conflict over wolf recovery in many, but not all, counties (e.g., Pearson's r correlation = 0.59 in one county and a nonsignificant correlation in another). Our findings suggest that modernization operates on values within a state with implications for biodiversity, but other factors in addition to values must be investigated to fully understand what leads to proconservation behavior.     

A spatially explicit estimate of the prewhaling abundance of the endangered North Atlantic right whale

The North Atlantic right whale (NARW) (Eubalaena glacialis) is one of the world's most threatened whales. It came close to extinction after nearly a millennium of exploitation and currently persists as a population of only approximately 500 individuals. Setting appropriate conservation targets for this species requires an understanding of its historical population size, as a baseline for measuring levels of depletion and progress toward recovery. This is made difficult by the scarcity of records over this species’ long whaling history. We sought to estimate the preexploitation population size of the North Atlantic right whale and understand how this species was distributed across its range. We used a spatially explicit data set on historical catches of North Pacific right whales (NPRWs) (Eubalaena japonica) to model the relationship between right whale relative density and the environment during the summer feeding season. Assuming the 2 right whale species select similar environments, we projected this model to the North Atlantic to predict how the relative abundance of NARWs varied across their range. We calibrated these relative abundances with estimates of the NPRW total prewhaling population size to obtain high and low estimates for the overall NARW population size prior to exploitation. The model predicted 9,075–21,328 right whales in the North Atlantic. The current NARW population is thus <6% of the historical North Atlantic carrying capacity and has enormous potential for recovery. According to the model, in June–September NARWs concentrated in 2 main feeding areas: east of the Grand Banks of Newfoundland and in the Norwegian Sea. These 2 areas may become important in the future as feeding grounds and may already be used more regularly by this endangered species than is thought.     

IDENTIFICATION OF CHANGES IN STREAMFLOW CHARACTERISTICS1

ABSTRACT: A set of procedures for identifying changes in selected streamflow characteristics at sites having long‐term continuous streamflow records is illustrated by using streamflow data from the Waccamaw River at Freeland, North Carolina for the 55‐year period of 1940–1994. Data were evaluated and compared to streamflow in the adjacent Lumber River Basin to determine if changes in streamflow characteristics in the Waccamaw River were localized and possibly the result of some human activity, or consistent with regional variations. Following 1963, droughts in the Waccamaw Basin seem to have been less severe than in the Lumber Basin, and the annual one‐, seven‐, and 30‐day low flows exhibited a slightly increasing trend in the Waccamaw River. Mean daily flows in the Waccamaw River at the 90 percent exceedance level (low flows) during 1985–194, a relatively dry period, were very nearly equal to flows at the same exceedance level for 1970–1979, which represents the 10‐year period between 1940 and 1994 with the highest flows. Prior to the 1980s, flows per unit drainage area in the Waccamaw Basin were generally less than those in the Lumber Basin, but after 1980, the opposite was true. The ratio of base flow to runoff in the Waccamaw River may have changed relative to that in the Lumber River in the late 1970s. There was greater variability in Waccamaw River streamflow than in Lumber River flow, and flow variability in the Waccamaw River may have increased slightly during 1985–1994.