Spatial and temporal effects of grazing management and rainfall on the vertebrate fauna of a tropical savanna

Grazing by domestic livestock is one of the most widespread uses of the rangelands of Australia. There is limited information on the effects of grazing by domestic livestock on the vertebrate fauna of Australia and the establishment of a long-term grazing experiment in north-eastern Queensland at Wambiana provided an opportunity to attempt an examination of the changes in vertebrate fauna as a consequence of the manipulation of stocking rates. The aim was to identify what the relative effects of vegetation type, stocking rate and other landscape-scale environmental factors were on the patterns recorded. Sixteen 1-ha sites were established within three replicated treatments (moderate, heavy and variable stocking rates). The sites were sampled in the wet and dry seasons in 1999-2000 (T-0) and again in 2003-04 (T-1). All paddocks of the treatments were burnt in 1999. Average annual rainfall declined markedly between the two sampling periods, which made interpretation of the data difficult. A total of 127 species of vertebrate fauna comprising five amphibian, 83 bird, 27 reptile and 12 mammal species were recorded. There was strong separation in faunal composition from T-0 to T-1 although changes in mean compositional dissimilarity between the grazing stocking rate treatments were less well defined. There was a relative change in abundance of 24 bird, four mammal and five reptile species from T-0 to T-1. The generalised linear modelling identified that, in the T-1 data, there was significant variation in the abundance of 16 species explained by the grazing and vegetation factors. This study demonstrated that vertebrate fauna assemblage did change and that these changes were attributable to the interplay between the stocking rates, the vegetation types on the sites surveyed, the burning of the experimental paddocks and the decrease in rainfall over the course of the two surveys. It is recommended that the experiment is sampled again but that the focus should be on a rapid survey of abundant taxa (i.e. birds and reptiles) to allow an increase in the frequency of sampling and replication of the data. This would help to articulate more clearly the trajectory of vertebrate change due to the relative effects of stocking rates compared with wider landscape environmental changes. Given the increasing focus on pastoral development in northern Australia, any opportunity to incorporate the collection of data on biodiversity into grazing manipulation experiments should be taken for the assessment of the effects of land management on faunal species.

Stochastic modelling of monthly rainfall

A transformation is suggested which can transform a non-Gaussian monthly hydrological time series into a Gaussian one. The suggested approach is verified with data of ten Indian rainfall time series. Incidentally, it is observed that once the deterministic trends are removed, the transformation leads to an uncorrelated process for monthly rainfall. The procedure for normalization is general enough in that it should be also applicable to river discharges. This is verified to a limited extent by considering data of two Indian river discharges.

Coherent rainfall zones of the Indian region

We have delineated rainfall zones for the Indian region that are coherent with respect to the variations of the summer monsoon rainfall. Within each zone, the time series of the summer monsoon rainfall at every pair of stations are significantly positively correlated, and the mean interseries correlation for each zone is high. The interseries correlation data set is analysed in order to delineate the rainfall zones, using an objective method specifically developed for the purpose. Each of the zonal averages are shown to be representative of the zone as a whole. We suggest that this regionalization is appropriate for study of the variation of the summer monsoon rainfall over the Indian region on interannual and larger scales.

Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Regionalization of precipitation in data sparse areas using large scale atmospheric variables - A fuzzy clustering approach

Delineation of homogeneous precipitation regions (regionalization) is necessary for investigating frequency and spatial distribution of meteorological droughts. The conventional methods of regionalization use statistics of precipitation as attributes to establish homogeneous regions. Therefore they cannot be used to form regions in ungauged areas, and they may not be useful to form meaningful regions in areas having sparse rain gauge density. Further, validation of the regions for homogeneity in precipitation is not possible, since the use of the precipitation statistics to form regions and subsequently to test the regional homogeneity is not appropriate. To alleviate this problem, an approach based on fuzzy cluster analysis is presented. It allows delineation of homogeneous precipitation regions in data sparse areas using large scale atmospheric variables (LSAV), which influence precipitation in the study area, as attributes. The LSAV, location parameters (latitude, longitude and altitude) and seasonality of precipitation are suggested as features for regionalization. The approach allows independent validation of the identified regions for homogeneity using statistics computed from the observed precipitation. Further it has the ability to form regions even in ungauged areas, owing to the use of attributes that can be reliably estimated even when no at-site precipitation data are available. The approach was applied to delineate homogeneous annual rainfall regions in India, and its effectiveness is illustrated by comparing the results with those obtained using rainfall statistics, regionalization based on hard cluster analysis, and meteorological sub-divisions in India. (C) 2011 Elsevier B.V. All rights reserved.

Assessing future rainfall projections using multiple GCMs and a multi-site stochastic downscaling model

Impact of global warming on daily rainfall is examined using atmospheric variables from five General Circulation Models (GCMs) and a stochastic downscaling model. Daily rainfall at eleven raingauges over Malaprabha catchment of India and National Center for Environmental Prediction (NCEP) reanalysis data at grid points over the catchment for a continuous time period 1971-2000 (current climate) are used to calibrate the downscaling model. The downscaled rainfall simulations obtained using GCM atmospheric variables corresponding to the IPCC-SRES (Intergovernmental Panel for Climate Change - Special Report on Emission Scenarios) A2 emission scenario for the same period are used to validate the results. Following this, future downscaled rainfall projections are constructed and examined for two 20 year time slices viz. 2055 (i.e. 2046-2065) and 2090 (i.e. 2081-2100). The model results show reasonable skill in simulating the rainfall over the study region for the current climate. The downscaled rainfall projections indicate no significant changes in the rainfall regime in this catchment in the future. More specifically, 2% decrease by 2055 and 5% decrease by 2090 in monsoon (HAS) rainfall compared to the current climate (1971-2000) under global warming conditions are noticed. Also, pre-monsoon (JFMAM) and post-monsoon (OND) rainfall is projected to increase respectively, by 2% in 2055 and 6% in 2090 and, 2% in 2055 and 12% in 2090, over the region. On annual basis slight decreases of 1% and 2% are noted for 2055 and 2090, respectively.

Large scale features and assessment of spatial scale correspondence between TMPA and IMD rainfall datasets over Indian landmass

Daily rainfall datasets of 10 years (1998-2007) of Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) version 6 and India Meteorological Department (IMD) gridded rain gauge have been compared over the Indian landmass, both in large and small spatial scales. On the larger spatial scale, the pattern correlation between the two datasets on daily scales during individual years of the study period is ranging from 0.4 to 0.7. The correlation improved significantly (similar to 0.9) when the study was confined to specific wet and dry spells each of about 5-8 days. Wavelet analysis of intraseasonal oscillations (ISO) of the southwest monsoon rainfall show the percentage contribution of the major two modes (30-50 days and 10-20 days), to be ranging respectively between similar to 30-40% and 5-10% for the various years. Analysis of inter-annual variability shows the satellite data to be underestimating seasonal rainfall by similar to 110 mm during southwest monsoon and overestimating by similar to 150 mm during northeast monsoon season. At high spatio-temporal scales, viz., 1 degrees x1 degrees grid, TMPA data do not correspond to ground truth. We have proposed here a new analysis procedure to assess the minimum spatial scale at which the two datasets are compatible with each other. This has been done by studying the contribution to total seasonal rainfall from different rainfall rate windows (at 1 mm intervals) on different spatial scales (at daily time scale). The compatibility spatial scale is seen to be beyond 5 degrees x5 degrees average spatial scale over the Indian landmass. This will help to decide the usability of TMPA products, if averaged at appropriate spatial scales, for specific process studies, e.g., cloud scale, meso scale or synoptic scale.

Copula-Based Modeling of TMI Brightness Temperature With Rainfall Type

Overland rain retrieval using spaceborne microwave radiometer offers a myriad of complications as land presents itself as a radiometrically warm and highly variable background. Hence, land rainfall algorithms of the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) have traditionally incorporated empirical relations of microwave brightness temperature (Tb) with rain rate, rather than relying on physically based radiative transfer modeling of rainfall (as implemented in the TMI ocean algorithm). In this paper, sensitivity analysis is conducted using the Spearman rank correlation coefficient as benchmark, to estimate the best combination of TMI low-frequency channels that are highly sensitive to the near surface rainfall rate from the TRMM Precipitation Radar (PR). Results indicate that the TMI channel combinations not only contain information about rainfall wherein liquid water drops are the dominant hydrometeors but also aid in surface noise reduction over a predominantly vegetative land surface background. Furthermore, the variations of rainfall signature in these channel combinations are not understood properly due to their inherent uncertainties and highly nonlinear relationship with rainfall. Copula theory is a powerful tool to characterize the dependence between complex hydrological variables as well as aid in uncertainty modeling by ensemble generation. Hence, this paper proposes a regional model using Archimedean copulas, to study the dependence of TMI channel combinations with respect to precipitation, over the land regions of Mahanadi basin, India, using version 7 orbital data from the passive and active sensors on board TRMM, namely, TMI and PR. Studies conducted for different rainfall regimes over the study area show the suitability of Clayton and Gumbel copulas for modeling convective and stratiform rainfall types for the majority of the intraseasonal months. Furthermore, large ensembles of TMI Tb (from the most sensitive TMI channel combination) were generated conditional on various quantiles (25th, 50th, 75th, and 95th) of the convective and the stratiform rainfall. Comparatively greater ambiguity was observed to model extreme values of the convective rain type. Finally, the efficiency of the proposed model was tested by comparing the results with traditionally employed linear and quadratic models. Results reveal the superior performance of the proposed copula-based technique.

Numerical Study on Rainfall Infiltration in Rock-Soil Blocks

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

Numerical study on rainfall infiltration in rock-soil slop

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

Comments on hydrographic and topographic LIDAR acquisition and merging with multibeam sounding data acquired in the Olympic Coast National Marine Sanctuary

In April 2005, a SHOALS 1000T LIDAR system was used as an efficient alternative for safely acquiring data to describe the existing conditions of nearshore bathymetry and the intertidal zone over an approximately 40.7 km2 (11.8 nm2) portion of hazardous coastline within the Olympic Coast National Marine Sanctuary (OCNMS). Data were logged from 1,593 km (860 nm) of track lines in just over 21 hours of flight time. Several islands and offshore rocks were also surveyed, and over 24,000 geo-referenced digital still photos were captured to assist with data cleaning and QA/QC. The 1 kHz bathymetry laser obtained a maximum water depth of 22.2 meters. Floating kelp beds, breaking surf lines and turbid water were all challenges to the survey. Although sea state was favorable for this time of the year, recent heavy rainfall and a persistent low-lying layer of fog reduced acquisition productivity. The existence of a completed VDatum model covering this same geographic region permitted the LIDAR data to be vertically transformed and merged with existing shallow water multibeam data and referenced to the mean lower low water (MLLW) tidal datum. Analysis of a multibeam bathymetry-LIDAR difference surface containing over 44,000 samples indicated surface deviations from –24.3 to 8.48 meters, with a mean difference of –0.967 meters, and standard deviation of 1.762 meters. Errors in data cleaning and false detections due to interference from surf, kelp, and turbidity likely account for the larger surface separations, while the remaining general surface difference trend could partially be attributed to a more dense data set, and shoal-biased cleaning, binning and gridding associated with the multibeam data for maintaining conservative least depths important for charting dangers to navigation. (PDF contains 27 pages.)

An ecological characterization of the marine resources of Vieques, Puerto Rico. Part I: Historical data synthesis

From the 1940s until 2003, portions of the island of Vieques, a municipality within the Commonwealth of Puerto Rico, were used by the US Navy as a base and training facility, resulting in development and zoning history that differ in comparison to other Caribbean islands. The majority of former Navy lands are now under the jurisdiction of the Department of the Interior’s Fish and Wildlife Service as a National Wildlife Refuge, while a smaller percentage of land was transferred to the Vieques municipality and the Puerto Rico Conservation Trust. An analysis of the distribution and status of the marine resources is timely in light of the recent land transfer, increases in development and tourism, and potential changes in marine zoning around the island. To meet this need, NOAA’s Biogeography Branch, in cooperation with the Office of Response and Restoration and other local and regional partners, conducted Part I of an ecological characterization to integrate historical data and research into a synthesis report. The overall objective of this report is to provide resource managers and residents a comprehensive characterization of the marine resources of Vieques to support research, monitoring, and management. For example, knowledge of the spatial distribution of physical features, habitats, and biological communities is necessary to make an informed decision of the establishment and placement of a marine protected area (MPA). The report is divided into chapters based on the physical environment (e.g., climate, geology, bathymetry), habitat types (e.g., reefs and hardbottom, seagrasses, mangroves) and major faunal groups (e.g. fish, turtles, birds). Each section includes five subsections: an overview, description of the relevant literature, methods of analysis, information on the distribution, status and trends of the particular resource, and a discussion of ecological linkages with other components of the Vieques marine ecosystem and surrounding environment. The physical environment of Vieques is similar to other islands within the Greater Antilles chain, with some distinctions. The warm, tropical climate of Vieques, mediated by the northeasterly trade winds, is characterized by a dry season (December-April) and a rainy season (May-November), the latter of which is characterized by the occasional passage of tropical cyclones. Compared to mainland Puerto Rico, Vieques is characterized by lower elevation, less annual precipitation, and higher average temperatures. The amount of annual precipitation also varies spatially within Vieques, with the western portion of the island receiving higher amounts of rainfall than further east. While the North Equatorial Current dominates the circulation pattern in the Greater Antilles region, small scale current patterns specific to Vieques are not as well characterized. These physical processes are important factors mitigating the distribution and composition of marine benthic habitats around Vieques. In general, the topography of Vieques is characterized by rolling hills. Mt. Pirata, the tallest point at 301 m, is located near the southwest coast. In the absence of island wide sedimentation measurements, information on land cover, slope, precipitation, and soil type were used to estimate relative erosion potential and sediment delivery for each watershed. While slope and precipitation amount are the primary driving factors controlling runoff, land use practices such as urban development, military activity, road construction, and agriculture can increase the delivery of pollution and sediments to coastal waters. Due to the recent land transfer, increased development and tourism is expected, which may result in changes in the input of sediments to the coastal environment.

Case studies of two contrasting rainfall episodes in Hawaii

This case study, utilizing surface and upper-air data, has attempted to shed light on the mechanisms that exerted control on two contrasting rainfall episodes in Hawaii [in the dry winter of 1981 and wet winter of 1982].

Influence of climatography on rainfall thresholds for initiation of debris-flows on the California coast [abstract]

EXTRACT (SEE PDF FOR FULL ABSTRACT): Since 1986, the U.S. Geological Survey and National Weather Service have operated a warning system for debris flows triggered by severe rainstorms in the San Francisco Bay region. The NWS tracks storm systems as they approach the region, forecasts precipitation, and observes rainfall with a network of radiotelemetered rain gauges (ALERT). The USGS also monitors ALERT data and compares the observed and forecast rainfall to thresholds for debris-flow initiation. Both groups jointly assess debris-flow hazards and issue public advisories when rainfall conditions reach or approach critical levels.

Coupling model uncertainty for coupled rainfall/runoff and surface water quality models in river problems

Coupled hydrology and water quality models are an important tool today, used in the understanding and management of surface water and watershed areas. Such problems are generally subject to substantial uncertainty in parameters, process understanding, and data. Component models, drawing on different data, concepts, and structures, are affected differently by each of these uncertain elements. This paper proposes a framework wherein the response of component models to their respective uncertain elements can be quantified and assessed, using a hydrological model and water quality model as two exemplars. The resulting assessments can be used to identify model coupling strategies that permit more appropriate use and calibration of individual models, and a better overall coupled model response. One key finding was that an approximate balance of water quality and hydrological model responses can be obtained using both the QUAL2E and Mike11 water quality models. The balance point, however, does not support a particularly narrow surface response (or stringent calibration criteria) with respect to the water quality calibration data, at least in the case examined here. Additionally, it is clear from the results presented that the structural source of uncertainty is at least as significant as parameter-based uncertainties in areal models. © 2012 John Wiley & Sons, Ltd.