SEARCH · Search USGS

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

At least 775 records · Page 43Linked to original sources

USGS2021

Remote sensing-based actual evapotranspiration assessment in a data-scarce area of Brazil: A case study of the Urucuia Aquifer System

The large groundwater reserves of the Urucuia Aquifer System (UAS) enabled agricultural development and economic growth in the western Bahia State, in northeastern Brazil. Over the last several years, concern has grown around the aquifer’s diminishing water levels, and water balance (WB) studies are in demand. Considering the lack of measured actual evapotranspiration (ET a ), a major component of the water cycle, this work uses the Operational Simplified Surface Energy Balance (SSEBop) model to estimate ET a , and compares it to basin-scale estimates from the Soil Moisture Accounting Procedure (SMAP) monthly model and from an annual WB closure method, based on gridded meteorological data and the Gravity Recovery and Climate Experiment (GRACE) product. Additionally, a comparative assessment of different versions of the SSEBop parameterization was performed. Moderate Resolution Imaging Spectroradiometer (MODIS) imagery was used to implement eight different versions of the SSEBop algorithm over the UAS between 2000 and 2013. SSEBop and SMAP ET a yielded similar seasonal patterns, with correlation coefficient (r) up to 0.65, mean difference (MD) of 0.8 mm/month and mean absolute difference (MAD) of 18.5 mm/month. Comparison of SSEBop annual ET a estimates to annual SMAP and WB closure estimates yielded low MD (12.1 and −7.3 mm/year, respectively) and MAD (82.5 and 82.8 mm/year, respectively), but also low r values (0.00 and 0.37, respectively). The comparison of the different SSEBop versions indicated the need to incorporate a calibration step of the aerodynamic heat resistance (r ah ) parameter. SSEBop results were also used for land cover and drought monitoring. Analysis indicates that agriculture, associated with an increasing trend of atmospheric evaporative demand, is responsible for the decrease in groundwater levels and streamflow in the studied time period.

Urucuia Aquifer System↗

USGS2022

Bathymetry retrieval from CubeSat image sequences with short time lags

The rapid expansion of CubeSat constellations could revolutionize the way inland and nearshore coastal waters are monitored from space. This potential stems from the ability of CubeSats to provide daily imagery with global coverage at meter-scale spatial resolution. In this study, we explore the unique opportunity to improve the retrieval of bathymetry offered by CubeSats, specifically those of the PlanetScope constellation. The orbital design of the PlanetScope constellation enables the acquisition of image sequences with short time lags (from seconds to hours). This characteristic allows multiple images to be captured during a short period of steady bathymetric conditions, especially in dynamic environments like rivers. We hypothesize that taking the ensemble mean of a CubeSat image sequence can enhance bathymetry retrieval compared to standard single-image analysis. Along with the existing optimal band ratio analysis (OBRA) algorithm, we also use a new neural network-based depth retrieval (NNDR) technique to infer bathymetry from both individual and time-averaged images. The two methodologies are evaluated using field data from five different river reaches with depths up to 15 m and both top-of-atmosphere (TOA) radiance and bottom-of-atmosphere (BOA) surface reflectance PlanetScope data products. Despite low spectral resolution and concerns about the radiometric quality of CubeSat imagery, accuracy assessment based on in-situ comparisons indicates the potential (0.52 < R 2 < 0.7 for the NNDR method) of PlanetScope imagery to retrieve depths up to ∼ 10 m in clear water conditions. The proposed image averaging consistently improves bathymetry retrieval over single image analysis. The NNDR technique was found to outperform OBRA, illustrating the importance of leveraging all spectral bands through machine learning approaches. TOA data provided more robust bathymetry results than BOA data for the OBRA technique, but the NNDR technique was minimally impacted by the type of data product.

International Journal of Applied Earth Observation↗

USGS2013

Reevaluation of lake trout and lake whitefish bioenergetics models

Using a corrected algorithm for balancing the energy budget, we reevaluated the Wisconsin bioenergetics model for lake trout (Salvelinus namaycush) in the laboratory and for lake whitefish (Coregonus clupeaformis) in the laboratory and in the field. For lake trout, results showed that the bioenergetics model slightly overestimated food consumption by the lake trout when they were fed low and intermediate rations, whereas the model predicted food consumption by lake trout fed ad libitum without any detectable bias. The slight bias in model predictions for lake trout on restricted rations may have been an artifact of the feeding schedule for these fish, and we would therefore recommend application of the Wisconsin lake trout bioenergetics model to lake trout populations in the field without any revisions to the model. Use of the Wisconsin bioenergetics model for coregonids resulted in overestimation of food consumption by lake whitefish both in the laboratory and in the field by between 20 and 30%, on average. This overestimation of food consumption was most likely due to overestimation of respiration rate. We therefore adjusted the respiration component of the bioenergetics model to obtain a good fit to the observed consumption in our laboratory tanks. The adjusted model predicted the consumption in the laboratory and the field without any detectable bias. Until a detailed lake whitefish respiration study can be conducted, we recommend application of our adjusted version of the Wisconsin generalized coregonid bioenergetics model to lake whitefish populations in the field.

Journal of Great Lakes Research↗

USGS2019

A comparison of chlorophyll a values obtained from an autonomous underwater vehicle to satellite-based measures for Lake Michigan

Accurate methods to track changes in lake productivity through time and space are critical to fisheries management. Chlorophyll a is the most widely studied proxy for ecosystem primary production, and has been the topic of many studies. The main sources of chlorophyll a measurements are ship-based measures or multi-spectral satellite data. Autonomous underwater vehicles can survey large spatial extents approaching the scale of satellite data, but with the accuracy of ship-based water sampling methods. We use several statistical measures to compare measures of chlorophyll a collected in Lake Michigan with spatiotemporally matched satellite-derived measures of chlorophyll a from the MODIS Aqua multi-spectral sensor using NASA’s OC3 and the Great Lakes Fit algorithms. Our findings show a near one to one relationship between AUV data and both satellite-derived data sets when the AUV data are coarsened to the resolution of the satellite data. A comparison of satellite-based chlorophyll a to AUV-derived chlorophyll summarized in discrete water depth bins suggested that, based on decreasing coefficients of determination, satellite estimates of chlorophyll accounted for the most variability in chlorophyll a concentrations in the upper 10 m of the water column, even though satellite sensors may detect past this depth.

Lake Michigan↗

USGS2016

Uncertainty analysis of the Operational Simplified Surface Energy Balance (SSEBop) model at multiple flux tower sites

Evapotranspiration (ET) is an important component of the water cycle – ET from the land surface returns approximately 60% of the global precipitation back to the atmosphere. ET also plays an important role in energy transport among the biosphere, atmosphere, and hydrosphere. Current regional to global and daily to annual ET estimation relies mainly on surface energy balance (SEB) ET models or statistical and empirical methods driven by remote sensing data and various climatological databases. These models have uncertainties due to inevitable input errors, poorly defined parameters, and inadequate model structures. The eddy covariance measurements on water, energy, and carbon fluxes at the AmeriFlux tower sites provide an opportunity to assess the ET modeling uncertainties. In this study, we focused on uncertainty analysis of the Operational Simplified Surface Energy Balance (SSEBop) model for ET estimation at multiple AmeriFlux tower sites with diverse land cover characteristics and climatic conditions. The 8-day composite 1-km MODerate resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) was used as input land surface temperature for the SSEBop algorithms. The other input data were taken from the AmeriFlux database. Results of statistical analysis indicated that the SSEBop model performed well in estimating ET with an R2 of 0.86 between estimated ET and eddy covariance measurements at 42 AmeriFlux tower sites during 2001–2007. It was encouraging to see that the best performance was observed for croplands, where R2 was 0.92 with a root mean square error of 13 mm/month. The uncertainties or random errors from input variables and parameters of the SSEBop model led to monthly ET estimates with relative errors less than 20% across multiple flux tower sites distributed across different biomes. This uncertainty of the SSEBop model lies within the error range of other SEB models, suggesting systematic error or bias of the SSEBop model is within the normal range. This finding implies that the simplified parameterization of the SSEBop model did not significantly affect the accuracy of the ET estimate while increasing the ease of model setup for operational applications. The sensitivity analysis indicated that the SSEBop model is most sensitive to input variables, land surface temperature (LST) and reference ET ( ET o ); and parameters, differential temperature ( dT ), and maximum ET scalar ( K max ), particularly during the non-growing season and in dry areas. In summary, the uncertainty assessment verifies that the SSEBop model is a reliable and robust method for large-area ET estimation. The SSEBop model estimates can be further improved by reducing errors in two input variables ( ET o and LST) and two key parameters ( K max and dT ).

Journal of Hydrology↗

USGS2016

An overview of current applications, challenges, and future trends in distributed process-based models in hydrology

Process-based hydrological models have a long history dating back to the 1960s. Criticized by some as over-parameterized, overly complex, and difficult to use, a more nuanced view is that these tools are necessary in many situations and, in a certain class of problems, they are the most appropriate type of hydrological model. This is especially the case in situations where knowledge of flow paths or distributed state variables and/or preservation of physical constraints is important. Examples of this include: spatiotemporal variability of soil moisture, groundwater flow and runoff generation, sediment and contaminant transport, or when feedbacks among various Earth’s system processes or understanding the impacts of climate non-stationarity are of primary concern. These are situations where process-based models excel and other models are unverifiable. This article presents this pragmatic view in the context of existing literature to justify the approach where applicable and necessary. We review how improvements in data availability, computational resources and algorithms have made detailed hydrological simulations a reality. Avenues for the future of process-based hydrological models are presented suggesting their use as virtual laboratories, for design purposes, and with a powerful treatment of uncertainty.

Journal of Hydrology↗

USGS2021

Machine learning predictions of mean ages of shallow well samples in the Great Lakes Basin, USA

The travel time or “age” of groundwater affects catchment responses to hydrologic changes , geochemical reactions, and time lags between management actions and responses at down-gradient streams and wells. Use of atmospheric tracers has facilitated the characterization of groundwater ages, but most wells lack such measurements. This study applied machine learning to predict ages in wells across a large region around the Great Lakes Basin using well, chemistry, and landscape characteristics. For a dataset of age tracers in 961 samples, the travel time from the land surface to the sample location was estimated for each sample using parametric functions. The mean travel times were then modeled using a gradient boosting machine (GBM) algorithm with cross validation tuning of model metaparameters. The GBM approach was able to closely match estimated ages for the training data (RMSE = 0.26 natural-log scale years) and provided a reasonable match to testing data (RMSE = 0.84). Of the variables tested, well characteristics (e.g. depth), land use, hydrologic indicators (e.g. topographic wetness index), and water chemistry (e.g. nitrate, fluoride, and pH), substantially affected the predictions of age. GBM prediction was applied to 14,335 groundwater samples with median sample depth of 5.4 m, indicating for the Great Lakes Basin a broad distribution of ages among wells with a median of 32.9 years. Lag times of decades are likely for these wells to respond to changing solute fluxes near land surface. While depth variables most strongly affected predicted mean ages, chemical constituents exhibited smooth trends with age, consistent with prevailing conceptual models of evolving sources and geochemistry flowpaths. The results provide proof of concept for use of readily available variables of well, landscape, and chemical characteristics to improve groundwater age estimates across large regions.

Great Lakes basin↗

USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗

USGS2005

Seismicity and tilt associated with the 2003 Anatahan eruption sequence

On May 10, 2003, the first historical eruption of Anatahan volcano in the western Pacific Mariana Islands was fortuitously recorded by a broadband seismograph installed on the island only 4 days prior to the eruption. This station, located 7 km WNW of the active crater, together with another broadband seismograph on Sarigan Island 45 km to the north, continued to operate throughout the 2-month period of major eruptive activity in May and June and throughout the majority of the following year. In June 2003, the Saipan Emergency Management Office and the US Geological Survey installed two telemetered high-gain short-period seismic stations to monitor the activity in real-time. The only earthquakes detected in the 4-day period from the initial seismograph installation until 6 h prior to the eruption occurred approximately 20 km to the northeast of the island on May 8. The first volcano-tectonic (VT) event located near the volcano occurred at 01:53 GMT on May 10. The number of events per hour then increased dramatically and a period of about 80 discrete earthquakes per hour commenced at about 06:20 GMT, immediately prior to the estimated eruption time of 07:30 from the Volcanic Ash Advisory Center. A long-period tilt signal recorded on the horizontal components of the broadband seismograph, indicating upward movement of the crater region, also commenced at about 06:20. Inflation continued until 09:30, when the direction of tilt reversed. Deflation continued until 17:50, coinciding with a reduction in the number of VT events. The larger VT events were located with a linearized least-squares location algorithm. Magnitudes of located VT events on May 10 ranged from 2.0 to 3.2, but a period of larger VT events were recorded on May 11, with the largest M 4.2. After about 36 h of intense earthquake activity, the number of discrete VT events declined and was replaced by nearly continuous volcanic tremor for the next 6 weeks. Differing types of very long-period events may suggest complex non-destructive magmatic source mechanisms, persisting dominantly throughout the first 10 days of the initial eruption. Visual reports indicate that a small craggy dome extruded sometime between May 20 and June 5. From analogy with other volcanic dome extrusions, we believe the dome probably extruded just as the tremor amplitude decreased dramatically about May 24. This dome was then destroyed between June 13 and 14. Reduced displacement of the co-eruption tremor is estimated as approximately 40-80 cm2, suggesting an eruption with a Volcanic Explosivity Index of about 3. 2005 Elsevier B.V. All rights reserved.

Northern Mariana Islands↗

USGS2015

An automated SO2 camera system for continuous, real-time monitoring of gas emissions from Kīlauea Volcano's summit Overlook Crater

SO 2 camera systems allow rapid two-dimensional imaging of sulfur dioxide (SO 2 ) emitted from volcanic vents. Here, we describe the development of an SO 2 camera system specifically designed for semi-permanent field installation and continuous use. The integration of innovative but largely “off-the-shelf” components allowed us to assemble a robust and highly customizable instrument capable of continuous, long-term deployment at Kīlauea Volcano's summit Overlook Crater. Recorded imagery is telemetered to the USGS Hawaiian Volcano Observatory (HVO) where a novel automatic retrieval algorithm derives SO 2 column densities and emission rates in real-time. Imagery and corresponding emission rates displayed in the HVO operations center and on the internal observatory website provide HVO staff with useful information for assessing the volcano's current activity. The ever-growing archive of continuous imagery and high-resolution emission rates in combination with continuous data from other monitoring techniques provides insight into shallow volcanic processes occurring at the Overlook Crater. An exemplary dataset from September 2013 is discussed in which a variation in the efficiency of shallow circulation and convection, the processes that transport volatile-rich magma to the surface of the summit lava lake, appears to have caused two distinctly different phases of lake activity and degassing. This first successful deployment of an SO 2 camera for continuous, real-time volcano monitoring shows how this versatile technique might soon be adapted and applied to monitor SO 2 degassing at other volcanoes around the world.

Hawai'i↗

USGS2018

A retrospective look at the February 1993 east rift zone intrusion at Kīlauea volcano, Hawaii

The February 1993 dike intrusion in the East Rift Zone (ERZ) of Kīlauea Volcano, Hawai'i, was recognized from tilt and seismic data, but ground-based geodetic data were too sparse to constrain the characteristics of the intrusion. Analysis of Interferometric Synthetic Aperture Radar (InSAR) from the Japan Aerospace Exploration Agency (JAXA) JERS-1 satellite reveals a maximum of ~30 cm of line-of-sight (LOS) displacement occurring near Makaopuhi Crater in the middle ERZ of Kīlauea. We model this deformation signal as a subvertical dike using a 3D-Mixed Boundary Element Method (3D-MBEM) paired with a nonlinear inversion algorithm to find the best-fit model. The best-fit dike is located just to the west of Makaopuhi Crater striking N50°W, extends to within 100 m of the surface, is ~1.3 km in length by ~4.2 km in width along strike, and has a total volume of ~7.4 × 10 6 m 3 . In addition, a post-intrusion interferogram from JERS-1 spanning 1993–1997 was analyzed. Guided by previous results, our model for the 4-year period consists of opening of the deep rift zones by about 0.5 m at 3–8.5 km depth beneath the Southwest Rift Zone (SWRZ), ERZ and the summit. A sub-horizontal detachment fault is connected to the seaward side of the vertical dike-like source to mimic the décollement known to exist beneath the volcano. We classify the 1993 dike intrusion as a passive intrusion similar to those that occurred in 1997 and 1999. Passive intrusions lack precursory inflation at Kīlauea's summit, and the likely triggering mechanism is persistent deep rift opening combined with seaward motion of the south flank along the basal décollement. Passive intrusions make forecasting and hazard assessment difficult since they are not preceded by inflation nor by large increases in seismicity.

Hawaii↗

USGS2022

Rock alteration mapping in and around fossil shallow intrusions at Mt. Ruapehu New Zealand with laboratory and aerial hyperspectral imaging

Diagnostic absorption features in hyperspectral data can be used to identify a specific mineral or mineral associations. However, it is unknown how accurate hyperspectral mapping can be for identifying alteration mineral compositions at the resolution required to describe structures such as fossil intrusions, or whether it can accurately quantify the alteration present. This study compared petrographic observation with visible, near-infrared (VNIR), and shortwave infrared (SWIR) hyperspectral remote sensing at laboratory- (centimetre-scale) and aerial- (metre-scale) scales to characterise the abundance of surface hydrothermal rock alteration in and around a shallow fossil intrusion on Pinnacle Ridge, Mt. Ruapehu, New Zealand. We classified a high-resolution aerial hyperspectral image to develop a new surface alteration map using Spectral Angle Mapper (SAM) algorithm. The petrographic thin-section and the laboratory and aerial hyperspectral imaging revealed a spectrum of hydrous alteration phases as indicated by the presence of an absorption feature at 2207 nm. Moderate correlation exists between the depth of the absorption feature at 2207 nm and the point counting-derived alteration percent values, indicating reliability of laboratory-based hyperspectral analytical methods. In contrast, aerial hyperspectral data failed to provide any clear correlations to field-mapped alteration using a band-depth approach, and we interpret this due to ‘oversampling’ of surface (supergene) alteration, spectral mixing, and sensor limitations (e.g., bandwidth, signal-to-noise ratio). The hyperspectral image-derived alteration map, created using supervised image classification, can loosely be translated to a geotechnical map where porosity and permeability play a major role in localizing hydrothermal fluid flow and the formation of alteration mineral associations.

Mt. Ruapehu↗

USGS2012

Determination of sediment provenance by unmixing the mineralogy of source-area sediments: The "SedUnMix" program

Along the margins of areas such as Greenland and Baffin Bay, sediment composition reflects a complex mixture of sources associated with the transport of sediment in sea ice, icebergs, melt-water and turbidite plumes. Similar situations arise in many contexts associated with sediment transport and with the mixing of sediments from different source areas. The question is: can contributions from discrete sediment (bedrock) sources be distinguished in a mixed sediment by using mineralogy, and, if so, how accurately? To solve this problem, four end-member source sediments were mixed in various proportions to form eleven artificial mixtures. Two of the end-member sediments are felsic, and the other two have more mafic compositions. End member and mixed sediment mineralogies were measured for the < 2 mm sediment fractions by quantitative X-ray diffraction (qXRD). The proportions of source sediments in the mixtures then were calculated using an Excel macro program named SedUnMix, and the results were evaluated to determine the robustness of the algorithm. The program permits the unmixing of up to six end members, each of which can be represented by up to 5 alternative compositions, so as to better simulate variability within each source region. The results indicate that we can track the relative percentages of the four end members in the mixtures. We recommend, prior to applying the technique to down-core or to other provenance problems, that a suite of known, artificial mixtures of sediments from probable source areas be prepared, scanned, analyzed for quantitative mineralogy, and then analyzed by SedUnMix to check the sensitivity of the method for each specific unmixing problem.

Baffin Bay↗

USGS2023

Modeling surface wave dynamics in upper Delaware Bay with living shorelines

Living shorelines gain increasing attention because they stabilize shorelines and reduce erosion. This study leverages physics-based models and bagged regression tree (BRT) machine learning algorithm to simulate wave dynamics at a living shoreline composed of constructed oyster reefs (CORs) in upper Delaware Bay. The physics-based models consist of coupled Delft3D-FLOW and SWAN in four-level nested domains. The model accuracy converges with increasing mesh resolution. The simulated wave-induced current circulation substantiates the effectiveness of CORs in trapping sediments. The simulated yearly-averaged wave power correlates qualitatively with historical shoreline retreat rates. BRT is adopted to improve the model accuracy, identify key processes responsible for simulation errors in wave height ( H 8 ) and wave period ( T p ), and quantify their importance. In the CORs sheltered area, BRT reveals that simulation errors of wind seas mainly arise from wind forcing, wave breaking and wave triad interactions. Wave breaking is seven times more important than wind forcing for simulating H 8 , while wind forcing and triad interactions are of equal importance for simulating T p . Simulation errors of swells mostly stem from bottom friction and offshore wave boundary conditions. Results from this study can help the assessment and adaptive management of CORs-based living shoreline restoration projects under climate change.

New Jersey↗

USGS2010

Effect of roughness formulation on the performance of a coupled wave, hydrodynamic, and sediment transport model

A variety of algorithms are available for parameterizing the hydrodynamic bottom roughness associated with grain size, saltation, bedforms, and wave–current interaction in coastal ocean models. These parameterizations give rise to spatially and temporally variable bottom-drag coefficients that ostensibly provide better representations of physical processes than uniform and constant coefficients. However, few studies have been performed to determine whether improved representation of these variable bottom roughness components translates into measurable improvements in model skill. We test the hypothesis that improved representation of variable bottom roughness improves performance with respect to near-bed circulation, bottom stresses, or turbulence dissipation. The inner shelf south of Martha’s Vineyard, Massachusetts, is the site of sorted grain-size features which exhibit sharp alongshore variations in grain size and ripple geometry over gentle bathymetric relief; this area provides a suitable testing ground for roughness parameterizations. We first establish the skill of a nested regional model for currents, waves, stresses, and turbulent quantities using a uniform and constant roughness; we then gauge model skill with various parameterization of roughness, which account for the influence of the wave-boundary layer, grain size, saltation, and rippled bedforms. We find that commonly used representations of ripple-induced roughness, when combined with a wave–current interaction routine, do not significantly improve skill for circulation, and significantly decrease skill with respect to stresses and turbulence dissipation. Ripple orientation with respect to dominant currents and ripple shape may be responsible for complicating a straightforward estimate of the roughness contribution from ripples. In addition, sediment-induced stratification may be responsible for lower stresses than predicted by the wave–current interaction model.

Ocean Modelling↗

USGS2023

HyWaves: Hybrid downscaling of multimodal wave spectra to nearshore areas

Long-term and accurate wave hindcast databases are often required in different coastal engineering projects. The assessment of the nearshore wave climate is often accomplished by using downscaling techniques to translate offshore waves to coastal areas. However, dynamical downscaling approaches may incur huge computational cost. Additionally, the common use of bulk parameterizations are often not accurate for multidimensional waves. To overcome these limitations, we present a hybrid downscaling approach that combines mathematical algorithms (statistical downscaling) and numerical modeling (dynamical downscaling) over the individual spectral partitions. Every wave partition is downscaled and aggregated afterward by using principles of wave linear theory. By assuming linearity in the propagation of the wave celerity, the application of the method is limited from offshore to intermediate water depths. In addition, the method proposed uses a technique to simplify the spectral boundary conditions in complex domains. The methodology has been applied and validated in the island states of Samoa, American Samoa, Majuro, and Kwajalein, showing good skill at reproducing the spectral hourly time series of significant wave height, peak period, and peak direction. Moreover, an accurate representation of the observed energy spectrum was achieved. This study provides insight into the numerical approximation of the combined sea-swell states while improving the quality of fast spectral forecasting and early warning systems.

Ocean Modeling↗

USGS2009

TreeMAC: Localized TDMA MAC protocol for real-time high-data-rate sensor networks

Earlier sensor network MAC protocols focus on energy conservation in low-duty cycle applications, while some recent applications involve real-time high-data-rate signals. This motivates us to design an innovative localized TDMA MAC protocol to achieve high throughput and low congestion in data collection sensor networks, besides energy conservation. TreeMAC divides a time cycle into frames and each frame into slots. A parent node determines the children's frame assignment based on their relative bandwidth demand, and each node calculates its own slot assignment based on its hop-count to the sink. This innovative 2-dimensional frame-slot assignment algorithm has the following nice theory properties. First, given any node, at any time slot, there is at most one active sender in its neighborhood (including itself). Second, the packet scheduling with TreeMAC is bufferless, which therefore minimizes the probability of network congestion. Third, the data throughput to the gateway is at least 1/3 of the optimum assuming reliable links. Our experiments on a 24-node testbed show that TreeMAC protocol significantly improves network throughput, fairness, and energy efficiency compared to TinyOS's default CSMA MAC protocol and a recent TDMA MAC protocol Funneling-MAC. Partial results of this paper were published in Song, Huang, Shirazi and Lahusen [W.-Z. Song, R. Huang, B. Shirazi, and R. Lahusen, TreeMAC: Localized TDMA MAC protocol for high-throughput and fairness in sensor networks, in: The 7th Annual IEEE International Conference on Pervasive Computing and Communications, PerCom, March 2009]. Our new contributions include analyses of the performance of TreeMAC from various aspects. We also present more implementation detail and evaluate TreeMAC from other aspects. ?? 2009 Elsevier B.V.

Pervasive and Mobile Computing↗

USGS2008

The imaging performance of the SRC on Mars Express

The Mars Express spacecraft carries the pushbroom scanner high-resolution stereo camera (HRSC) and its added imaging subsystem super resolution channel (SRC). The SRC is equipped with its own optical system and a 1024×1024 framing sensor. SRC produces snapshots with 2.3 m ground pixel size from the nominal spacecraft pericenter height of 250 km, which are typically embedded in the central part of the large HRSC scenes. The salient features of the SRC are its light-weight optics, a reliable CCD detector, and high-speed read-out electronics. The quality and effective visibility of details in the SRC images unfortunately falls short of what has been expected. In cases where thermal balance cannot be reached, artifacts, such as blurring and “ghost features” are observed in the images. In addition, images show large numbers of blemish pixels and are plagued by electronic noise. As a consequence, we have developed various image improving algorithms, which are discussed in this paper. While results are encouraging, further studies of image restoration by dedicated processing appear worthwhile. The SRC has obtained more than 6940 images at the time of writing (1 September 2007), which often show fascinating details in surface morphology. SRC images are highly useful for a variety of applications in planetary geology, for studies of the Mars atmosphere, and for astrometric observations of the Martian satellites. This paper will give a full account of the design philosophy, technical concept, calibration, operation, integration with HRSC, and performance, as well as science accomplishments of the SRC.

Planetary and Space Science↗