SEARCH · Search USGS
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.
With the advent of small, low-cost environmental sensors, it is now possible to deploy high-density networks of sensors to measure hyper localized temperature variation. Here, we provide a detailed methodology for constructing a compact version of a previously described custom-fabricated radiation shield for use with inexpensive thermochrons.
Laboratory visible/near-infrared multispectral observations of Mars Exploration Rover Pancam calibration target materials coated with different thicknesses of Mars spectral analog dust were acquired under variable illumination geometries using the Bloomsburg University Goniometer. The data were fit with a two-layer radiative transfer model that combines a Hapke formulation for the dust with measured values of the substrate interpolated using a He-Torrance approach. We first determined the single-scattering albedo, phase function, opposition effect width, and amplitude for the dust using the entire data set (six coating thicknesses, three substrates, four wavelengths, and phase angles 3??-117??). The dust exhibited single-scattering albedo values similar to other Mars analog soils and to Mars Pathfinder dust and a dominantly forward scattering behavior whose scattering lobe became narrower at longer wavelengths. Opacity values for each dust thickness corresponded well to those predicted from the particles sizes of the Mars analog dust. We then restricted the number of substrates, dust thicknesses, and incidence angles input to the model. The results suggest that the dust properties are best characterized when using substrates whose reflectances are brighter and darker than those of the deposited dust and data that span a wide range of dust thicknesses. The model also determined the dust photometric properties relatively well despite limitations placed on the range of incidence angles. The model presented here will help determine the photometric properties of dust deposited on the MER rovers and to track the multiple episodes of dust deposition and erosion that have occurred at both landing sites. Copyright 2006 by the American Geophysical Union.
The American Society of Photogrammetry and Remote Sensing defined remote sensing as the measurement or acquisition of information of some property of an object or phenomenon, by a recording device that is not in physical or intimate contact with the object or phenomenon under study (Colwell et al., 1983). Environmental Systems Research Institute (ESRI) in its geographic information system (GIS) dictionary defines remote sensing as “collecting and interpreting information about the environment and the surface of the earth from a distance, primarily by sensing radiation that is naturally emitted or reflected by the earth’s surface or from the atmosphere, or by sending signals transmitted from a device and reflected back to it (ESRI, 2014).” The usual source of passive remote sensing data is the measurement of reflected or transmitted electromagnetic radiation (EMR) from the sun across the electromagnetic spectrum (EMS); this can also include acoustic or sound energy, gravity, or the magnetic field from or of the objects under consideration. In this context, the simple act of reading this text is considered remote sensing. In this case, the eye acts as a sensor and senses the light reflected from the object to obtain information about the object. It is the same technology used by a handheld camera to take a photograph of a person or a distant scenic view. Active remote sensing, however, involves sending a pulse of energy and then measuring the returned energy through a sensor (e.g., Radio Detection and Ranging [RADAR], Light Detection and Ranging [LiDAR]). Thermal sensors measure emitted energy by different objects. Thus, in general, passive remote sensing involves the measurement of solar energy reflected from the Earth’s surface, while active remote sensing involves synthetic (man-made) energy pulsed at the environment and the return signals are measured and recorded.
A physically based model of variations in δ 18 O and δ 2 H in Pyramid Lake is presented. For inputs, the model uses measurements of liquid water inflows and outflows and their associated isotopic compositions and a set of meteorological data (radiative fluxes, air temperature, relative humidity, and windspeed). The model simulates change of lake volume, thermal and isotopic stratification, evaporation, and the isotopic composition of evaporation. A validation of the model for 1987–1989 and 1991 indicates that it can reproduce measured intra‐ and interannual variations of δ 18 O and δ 2 H. Three applications of the model demonstrate its ability to simulate longer term responses of δ 18 O to change in the hydrologic balance and hydrologic characteristics (opening and closing) of the lake.
Introduction: Optically stimulated luminescence (OSL) is one of a class of measurements known as stimulated phenomena. Such phenomena may be stimulated thermally or optically and the reader is referred to works by Aitken (1998) and Botter-Jensen and others (2003) for more detail. In recent years OSL has become a popular procedure for the determination of environmental radiation doses absorbed by archeological and geological materials in an attempt to date these materials. The first OSL measurements on quartz and feldspar were made using an argon ion-laser (Huntley et al., 1985). However, the development of cheaper stimulation systems based first on filtered lamps and then on light- emitting diodes (LEDs) (Spooner, et al., 1990; Botter-Jensen, and others, 1999) has led to a massive expansion in OSL dating applications. The abstracts in this volume represent presentations from a workshop held in May-June 2006, at the Denver Federal Center, Denver, Colorado, in which OSL methodologies and applications were summarized and integrated to provide a current synthesis of the OSL science being applied throughout North America. The workshop, sponsored by the U.S. Geological Survey Crustal Imaging and Characterization Team and North Dakota State University, was open to all scientists interested in OSL dating techniques and radiation dosimetry. Participants included thirty-six research scientists and students in geology, archaeology, and physics from the U.S. Geological Survey, Los Alamos National Labs, Kentucky Geological Survey, eight universities in the United States, one university in Canada, one university in India, and Riso National Labs of Denmark. The workshop included two keynote speakers: Dr. Ashok Singhvi (Physical Research Laboratory, Ahmedabad, India) spoke on 'Some Unexplored Methodological Aspects and Some New Applications of Luminescence Dating,' while Dr. Jim Feathers (University of Seattle, WA) spoke on OSL Dating of Sediments From Paleoindian Sites in Brazil. The workshop encouraged everyone to interact more to develop a broader perspective on the types of research and the problems encountered when reporting OSL ages. This meeting follows the first North American Luminescence Dating Workshop held in Tulsa, OK, by Oklahoma State University (2001), in Albuquerque, NM, by Los Alamos National Labs (2002), and in Halifax, Nova Scotia, by Dalhousie University (2004, with a name change to New World Luminescence Dating Workshop). These workshops were interspersed with the international meetings on luminescence that were held in Reno, NV, (2002), and Cologne, Germany; (2005).
Radiative effects of anthropogenic changes in atmospheric composition are expected to cause climate changes, in particular an intensification of the global water cycle with a consequent increase in flood risk. But the detection of anthropogenically forced changes in flooding is difficult because of the substantial natural variability; the dependence of streamflow trends on flow regime further complicates the issue. Here we investigate the changes in risk of great floods - that is, floods with discharges exceeding 100-year levels from basins larger than 200,000 km2 - using both streamflow measurements and numerical simulations of the anthropogenic climate change associated with greenhouse gases and direct radiative effects of sulphate aerosols. We find that the frequency of great floods increased substantially during the twentieth century. The recent emergence of a statistically significant positive trend in risk of great floods is consistent with results from the climate model, and the model suggests that the trend will continue.
Gross primary production (GPP) is a key indicator of ecosystem performance, and helps in many decision-making processes related to environment. We used the Eddy covariancelight use efficiency (EC-LUE) model for estimating GPP in the Great Plains, United States in order to evaluate the performance of this model. We developed a novel algorithm for computing the photosynthetically active radiation (PAR) based on net radiation. A strong correlation ( R 2 =0.94, N =24) was found between daily PAR and Landsat-based mid-day instantaneous net radiation. Though the Moderate Resolution Spectroradiometer (MODIS) based instantaneous net radiation was in better agreement ( R 2 =0.98, N =24) with the daily measured PAR, there was no statistical significant difference between Landsat based PAR and MODIS based PAR. The EC-LUE model validation also confirms the need to consider biological attributes (C 3 versus C 4 plants) for potential light use efficiency. A universal potential light use efficiency is unable to capture the spatial variation of GPP. It is necessary to use C 3 versus C 4 based land use/land cover map for using EC-LUE model for estimating spatiotemporal distribution of GPP.
Concentrations of trace elements, radiochemicals, and pesticides in the Bowdoin National Wildlife Refuge lakes generally were not substantially larger than those in the water supplied from Dodson South Canal or in irrigation drainage. Concentrations of arsenic (47 micrograms/L), uranium (43 microg/L), and vanadium (51 microg/L) in Dry Lake Unit, and boron (1,000 microg/L) in Lake Bowdoin were notably larger than at other sites. Zinc concentrations in an irrigation drain (56 microg/L) and two shallow domestic wells (40 and 47 microg/L) were elevated relative to other sites. Concentrations of gross alpha radiation (64 picocuries/L) and gross beta radiation (71 picocuries/L) were elevated in Dry Lake Unit. Pesticides concentrations at all sites were 0.08 microg/L or less. Water use guidelines concentrations for boron, cadmium, uranium, zinc, and gross alpha radiation were slightly exceeded at several sites. In general, trace-constituent concentrations measured in the water do not indicate any potential toxicity problems in Bowdoin National Wildlife Refuge; however, highwater conditions in 1986 probably caused dilution of dissolved constituents compared to recent dry years. Trace element concentrations in bottom sediments of the refuge lakes were generally similar to background concentrations in the soils. The only exception was Dry Lake Unit, which had concentrations of chromium (99 micrograms/g), copper (37 microg/g), nickel (37 microg/g), vanadium (160 microg/g), and zinc (120 microg/g) that were about double the mean background concentrations. The maximum selenium concentration in bottom sediment was 0.6 microg/g. Pesticide concentrations in bottom sediments were less than analytical detection limits at all sites. With few exceptions, concentrations of trace elements and pesticides in biota generally were less than values known to produce harmful effects on growth or reproduction.
The vast majority of surface water resources in the semi-arid western United States start as winter snowpack. Solar radiation is a primary driver of snowmelt, making snowpack water resources especially sensitive to even small increases in concentrations of light absorbing particles such as mineral dust and combustion-related black carbon (BC). Here we show, using fresh snow measurements and snowpack modeling at 51 widely distributed sites in the Rocky Mountain region, that BC dominated impurity-driven radiative forcing in 2018. BC contributed three times more radiative forcing on average than dust, and up to 17 times more at individual locations. Evaluation of 2015 to 2018 archived samples from most of the same sites yielded similar results. These findings, together with long-term observations of atmospheric concentrations and atmospheric model studies, indicate that BC rather than dust has dominated radiative forcing by light absorbing impurities on snow for decades, indicating that mitigation strategies to reduce radiative forcing on headwater snow-water resources would need to focus on reducing winter and spring BC emissions.
This report describes a hydrologic model for the upper Deschutes Basin in central Oregon developed using the U.S. Geological Survey (USGS) integrated Groundwater and Surface-Water Flow model (GSFLOW). The upper Deschutes Basin, which drains much of the eastern side of the Cascade Range in Oregon, is underlain by large areas of permeable volcanic rock. That permeability, in combination with the large annual precipitation at high elevations, results in a substantial regional aquifer system and a stream system that is heavily groundwater dominated. The upper Deschutes Basin is also an area of expanding population and increasing water demand for public supply and agriculture. Surface water was largely developed for agricultural use by the mid-20th century, and is closed to additional appropriations. Consequently, water users look to groundwater to satisfy the growing demand. The well‑documented connection between groundwater and the stream system, and the institutional and legal restrictions on streamflow depletion by wells, resulted in the Oregon Water Resources Department (OWRD) instituting a process whereby additional groundwater pumping can be permitted only if the effects to streams are mitigated, for example, by reducing permitted surface-water diversions. Implementing such a program requires understanding of the spatial and temporal distribution of effects to streams from groundwater pumping. A groundwater model developed in the early 2000s by the USGS and OWRD has been used to provide insights into the distribution of streamflow depletion by wells, but lacks spatial resolution in sensitive headwaters and spring areas. The integrated model developed for this project, based largely on the earlier model, has a much finer grid spacing allowing resolution of sensitive headwater streams and important spring areas, and simulates a more complete set of surface processes as well as runoff and groundwater flow. In addition, the integrated model includes improved representation of subsurface geology and explicitly simulates the effects of hydrologically important fault zones not included in the previous model. The upper Deschutes Basin GSFLOW model was calibrated using an iterative trial and error approach using measured water-level elevations (water levels) from 800 wells, 144 of which have time series of 10 or more measurements. Streamflow was calibrated using data from 21 gage locations. At 14 locations where measured flows are heavily influenced by reservoir operations and irrigation diversions, so called “ naturalized ” flows, with the effects of reservoirs and diversion removed, developed by the Bureau of Reclamation, were used for calibration. Surface energy and moisture processes such as solar radiation, snow accumulation and melting, and evapotranspiration were calibrated using national datasets as well as data from long-term measurement sites in the basin. The calibrated Deschutes GSFLOW model requires daily precipitation, minimum and maximum air temperature data, and monthly data describing groundwater pumping and artificial recharge from leaking irrigation canals (which are a significant source of groundwater recharge). The calibrated model simulates the geographic distribution of hydraulic head over the 5,000 ft range measured in the basin, with a median absolute residual of about 53 ft. Temporal variations in head resulting from climate cycles, pumping, and canal leakage are well simulated over the model area. Simulated daily streamflow matches gaged flows or calculated naturalized flows for streams including the Crooked and Metolius Rivers, and lower parts of the mainstem Deschutes River. Seasonal patterns of runoff are less well fit in some upper basin streams. Annual water balances of streamflow are good over most of the model domain. Model fit and overall capabilities are appropriate for the objectives of the project. The integrated model results confirm findings from other studies and models indicating that most streamflow in the upper Deschutes Basin comes directly from groundwater discharge. The integrated model provides additional insights about the components of streamflow including direct groundwater discharge to streams, interflow, groundwater discharge to the land surface (Dunnian flow), and direct runoff (Hortonian flow). The new model provides improved capability for exploring the timing and distribution of streamflow capture by wells, and the hydrologic response to changes in other external stresses such as canal operation, irrigation, and drought. Because the model uses basic meteorological data as the primary input; and simulates surface energy and moisture balances, groundwater recharge and flow, and all components of streamflow; it is well suited for exploring the hydrologic response to climate change, although no such simulations are included in this report. The model was developed as a tool for future application; however, example simulations are provided in this report. In the example simulations, the model is used to explore the influence of well location and geologic structure on stream capture by pumping wells. Wells were simulated at three locations within a 12-mi area close to known groundwater discharge areas and crossed by a regional fault zone. Simulations indicate that the magnitude and timing of stream capture from pumping is largely controlled by the geographic location of the wells, but that faults can have a large influence on the propagation of pumping stresses.
We document the composition of a natural zircon gemstone sourced from Madagascar, MAD‐559 – a new reference material for calibrating trace element mass fractions in zircon measured by SIMS. The composition of MAD‐559 was quantified by calibration relative to the well‐documented zircon reference material 91500, for which we compiled existing published data (Mg, Al, Y, rare earth elements, Hf, U, Th) and performed new measurements to characterise the mass fraction of less commonly measured elements (Li, Be, B, F, Na, P, K, Ca, Sc, Ti, Fe, Nb). Measurement results of SL13, CZ3 and MAD‐1 zircons and NIST SRM glasses were performed as quality control materials to test measurement bias and repeatability. We show the intermediate precision for most trace element measurement results of MAD‐559 to be between ± 3% and ± 5% RSD based on 139 measurements by SIMS on twenty‐five individual polished zircon chips measured during a 24‐h period, as well as repeat measurements performed over five separate analytical sessions. Trace element mass fractions were also measured by LA‐ICP‐MS in two different laboratories, and major element compositions measured by electron microprobe, to compare with results measured by SIMS. Based on laser Raman and hyperspectral cathodoluminescence spectroscopy, we show MAD‐559 to have high crystal disorder due to radiation damage relative to crystalline zircon (e.g., SL13 and 91500 zircon). Although the high cumulative alpha dose of MAD‐559 zircon makes it a poor reference material for geochronology, the consistency of the trace element mass fraction results measured in multiple sessions and by various measurement methods shows that it is an ideal reference material for microanalytical trace element mass fraction quantification of zircon.
Snowmelt is the principal source for soil moisture, ground-water re-charge, and stream-flow in mountainous regions of the western US, Canada, and other similar regions of the world. Information on the timing, magnitude, and contributing area of melt under variable or changing climate conditions is required for successful water and resource management. A coupled energy and mass-balance model ISNOBAL is used to simulate the development and melting of the seasonal snowcover in several mountain basins in California, Idaho, and Utah. Simulations are done over basins varying from 1 to 2500 km 2 , with simulation periods varying from a few days for the smallest basin, Emerald Lake watershed in California, to multiple snow seasons for the Park City area in Utah. The model is driven by topographically corrected estimates of radiation, temperature, humidity, wind, and precipitation. Simulation results in all basins closely match independently measured snow water equivalent, snow depth, or runoff during both the development and depletion of the snowcover. Spatially distributed estimates of snow deposition and melt allow us to better understand the interaction between topographic structure, climate, and moisture availability in mountain basins of the western US. Application of topographically distributed models such as this will lead to improved water resource and watershed management. Copyright © 1999 John Wiley & Sons, Ltd.
High radium (Ra) concentrations in potable portions of the Cambrian-Ordovician (C-O) aquifer system were investigated using water-quality data and environmental tracers ( 3 H, 3 He trit , SF 6 , 14 C and 4 He rad ) of groundwater age from 80 public-supply wells (PSWs). Groundwater ages were estimated by calibration of tracers to lumped parameter models and ranged from modern (<50 yr) in upgradient, regionally unconfined areas to ancient (>1 Myr) in the most downgradient, confined portions of the potable system. More than 80 and 40 percent of mean groundwater ages were older than 1000 and 50,000 yr, respectively. Anoxic, Fe-reducing conditions and increased mineralization develop with time in the aquifer system and mobilize Ra into solution resulting in the frequent occurrence of combined Ra (Ra c = 226 Ra + 228 Ra) at concentrations exceeding the USEPA MCL of 185 mBq/L (5 pCi/L). The distribution of the three Ra isotopes comprising total Ra (Ra t = 224 Ra + 226 Ra + 228 Ra) differed across the aquifer system. The concentrations of 224 Ra and 228 Ra were strongly correlated and comprised a larger proportion of the Ra t concentration in samples from the regionally unconfined area, where arkosic sandstones provide an enhanced source for progeny from the 232 Th decay series. 226 Ra comprised a larger proportion of the Ra t concentration in samples from downgradient confined regions. Concentrations of Ra t were significantly greater in samples from the regionally confined area of the aquifer system because of the increase in 226 Ra concentrations there as compared to the regionally unconfined area. 226 Ra distribution coefficients decreased substantially with anoxic conditions and increasing ionic strength of groundwater (mineralization), indicating that Ra is mobilized to solution from solid phases of the aquifer as adsorption capacity is diminished. The amount of 226 Ra released from solid phases by alpha-recoil mechanisms and retained in solution increases relative to the amount of Ra sequestered by adsorption processes or co-precipitation with barite as adsorption capacity and the concentration of Ba decreases. Although 226 Ra occurred at concentrations greater than 224 Ra or 228 Ra, the ingestion exposure risk was greater for 228 Ra owing to its greater toxicity. In addition, 224 Ra added substantial alpha-particle radioactivity to potable samples from the C-O aquifer system. Thus, monitoring for Ra isotopes and gross-alpha-activity (GAA) is important in upgradient, regionally unconfined areas as downgradient, and GAA measurements made within 72 h of sample collection would best capture alpha-particle radiation from the short-lived 224 Ra.
Snow and watershed models typically do not account for forest structure and shading; therefore, they display substantial uncertainty when attempting to account for forest change or when comparing hydrological response between forests with varying characteristics. This study collected snow water equivalent (SWE) measurements in a snow-dominated forest in Colorado, the United States, with variable canopy structure. The SWE measurements were integrated with 1 m Lidar derived canopy structure metrics and incoming solar radiation to create empirical SWE offset equations for four canopy structure groupings (forest gaps, south-facing forest edges, north-facing forest edges, and the interior forest) that varied in size compared to an open area. These simple equations indirectly integrate terrain shading and canopy shading and were able to estimate 40 to 70% of SWE variation in a heterogenous forested environment. The equations were then applied to a snow melt model with a 100 m grid size by applying the area-weighted average of SWE offsets from the four canopy structure groupings in each model cell. This tiled model configuration allowed for the model to better represent the subgrid heterogeneity of a forest environment that can be seen through an ensemble or range of potential outputs rather than a singular estimate.
A paper by two U.S. scientists explores the environmental/health hazard posed by abandoned missile launch sites and control facilities (dismantled by 1998 as part of the Nunn-Lugar Cooperative Threat Reduction Program) at an ICBM base in north-central Kazakhstan. It summarizes the findings of Environmental Site Assessment Reports based on a program of water and soil sampling at the sites, with a particular focus on testing for the presence of polychlorinated biphenyls, various toxic metals, and radiation. The study is important in that it documents levels of contamination (and describes abatement measures) at a former Soviet missile base.
The diffuse attenuation coefficient of photosynthetically active radiation (KdPAR) is commonly used to predict light attenuation in aquatic productivity models, but obtaining measurements of PAR to compute KdPAR is difficult. In situ calculations of KdPAR require multiple measurements of PAR through the water column, and these measurements are infeasible for real-time recording. Instead, predictive models using surface-water measurements may be used. Traditional KdPAR models are based on open-ocean habitats and rely on chlorophyll—as a proxy measurement for phytoplankton abundance—as the main predictive parameter. However, elevated suspended sediments and dissolved organic materials may also affect KdPAR values of inland water bodies and estuaries. In this study, we leverage KdPAR calculations derived from in situ light measurements collected along with surface-water-quality parameters across the Sacramento-San Joaquin River Delta in California, USA (the Delta). Sampling occurred between January of 2013 and May of 2014. We also explored regional and seasonal effects, but these did not clearly affect the model. Ultimately, the best-performing model included surface-level turbidity only (R2 = 0.91). The simplicity of the model facilitates use of KdPAR estimates for a variety of purposes throughout the Delta, including euphotic depth calculations, and as inputs to primary-productivity and habitat-suitability models. We demonstrate the model’s usability with two open-sources data sets (one spatially dense, and one temporally dense), and estimate KdPAR, euphotic depth, and primary productivity within the Delta. We provide calculations for each estimation, allowing users to easily adopt these models and apply them to their own data or with open-sourced data, which are abundant.
The University of Wisconsin (UW)-Madison Discovery Farms (Discovery Farms) and UW-Platteville Pioneer Farm (Pioneer Farm) programs were created in 2000 to help Wisconsin farmers meet environmental and economic challenges. As a partner with each program, and in cooperation with the Wisconsin Department of Natural Resources and the Sand County Foundation, the U.S. Geological Survey (USGS) Wisconsin Water Science Center (WWSC) installed, maintained, and operated equipment to collect water-quantity and water-quality data from 25 edge-offield, 6 streamgaging, and 5 subsurface-tile stations at 7 Discovery Farms and Pioneer Farm. The farms are located in the southern half of Wisconsin and represent a variety of landscape settings and crop- and animal-production enterprises common to Wisconsin agriculture. Meteorological stations were established at most farms to measure precipitation, wind speed and direction, air and soil temperature (in profile), relative humidity, solar radiation, and soil moisture (in profile). Data collection began in September 2001 and is continuing through the present (2008). This report describes methods used by USGS WWSC personnel to collect, process, and analyze water-quantity, water-quality, and meteorological data for edge-of-field, streamgaging, subsurface-tile, and meteorological stations at Discovery Farms and Pioneer Farm from September 2001 through October 2007. Information presented includes equipment used; event-monitoring and samplecollection procedures; station maintenance; sample handling and processing procedures; water-quantity, waterquality, and precipitation data analyses; and procedures for determining estimated constituent concentrations for unsampled runoff events.
Aim  Extrapolation of tower CO 2 fluxes will be greatly facilitated if robust relationships between flux components and remotely sensed factors are established. Long-term measurements at five Northern Great Plains locations were used to obtain relationships between CO 2 fluxes and photosynthetically active radiation ( Q ), other on-site factors, and Normalized Difference Vegetation Index ( NDVI ) from the SPOT VEGETATION data set. Location  CO 2 flux data from the following stations and years were analysed: Lethbridge, Alberta 1998–2001; Fort Peck, MT 2000, 2002; Miles City, MT 2000–01; Mandan, ND 1999–2001; and Cheyenne, WY 1997–98. Results  Analyses based on light-response functions allowed partitioning net CO 2 flux ( F ) into gross primary productivity ( P g ) and ecosystem respiration ( R e ). Weekly averages of daytime respiration, γ day , estimated from light responses were closely correlated with weekly averages of measured night-time respiration, γ night ( R 2 0.64 to 0.95). Daytime respiration tended to be higher than night-time respiration, and regressions of γ day on γ night for all sites were different from 1 : 1 relationships. Over 13 site-years, gross primary production varied from 459 to 2491 g CO 2 m −2 year −1 , ecosystem respiration from 996 to 1881 g CO 2 m −2 year −1 , and net ecosystem exchange from −537 (source) to +610 g CO 2 m −2 year −1 (sink). Maximum daily ecological light-use efficiencies, ɛ d , max = P g /Q , were in the range 0.014 to 0.032 mol CO 2 (mol incident quanta) −1 . Main conclusions  Ten-day average P g was significantly more highly correlated with NDVI than 10-day average daytime flux, P d ( R 2 = 0.46 to 0.77 for P g -NDVI and 0.05 to 0.58 for P d -NDVI relationships). Ten-day average R e was also positively correlated with NDVI , with R 2 values from 0.57 to 0.77. Patterns of the relationships of P g and R e with NDVI and other factors indicate possibilities for establishing multivariate functions allowing scaling-up local fluxes to larger areas using GIS data, temporal NDVI, and other factors.