Geology ReportsSearch

SEARCH · Geology Reports

Results for “Algorithms”

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.

At least 325 records · Page 18Linked to original sources

An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks

A highly automated algorithm called vegetation change tracker (VCT) has been developed for reconstructing recent forest disturbance history using Landsat time series stacks (LTSS). This algorithm is based on the spectral–temporal properties of land cover and forest change processes, and requires little or no fine tuning for most forests with closed or near close canopy cover. It was found very efficient, taking 2–3 h on average to analyze an LTSS consisting of 12 or more Landsat images using an average desktop PC. This LTSS-VCT approach has been used to examine disturbance patterns with a biennial temporal interval from 1984 to 2006 for many locations across the conterminous U.S. Accuracy assessment over 6 validation sites revealed that overall accuracies of around 80% were achieved for disturbances mapped at individual year level. Average user's and producer's accuracies of the disturbance classes were around 70% and 60% in 5 of the 6 sites, respectively, suggesting that although forest disturbances were typically rare as compared with no-change classes, on average the VCT detected more than half of those disturbances with relatively low levels of false alarms. Field assessment revealed that VCT was able to detect most stand clearing disturbance events, including harvest, fire, and urban development, while some non-stand clearing events such as thinning and selective logging were also mapped in western U.S. The applicability of the LTSS-VCT approach depends on the availability of a temporally adequate supply of Landsat imagery. To ensure that forest disturbance records can be developed continuously in the future, it is necessary to plan and develop observational capabilities today that will allow continuous acquisition of frequent Landsat or Landsat-like observations.

Remote Sensing of Environment

Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data

There are many types of changes occurring over the Earth's landscapes that can be detected and monitored using Landsat data. Here we focus on monitoring “within-state,” gradual changes in vegetation in contrast with traditional monitoring of “abrupt” land-cover conversions. Gradual changes result from a variety of processes, such as vegetation growth and succession, damage from insects and disease, responses to shifts in climate, and other factors. Despite the prevalence of gradual changes across the landscape, they are largely ignored by the remote sensing community. Gradual changes are best characterized and monitored using time-series analysis, and with the successful launch of Landsat 8 we now have appreciable data continuity that extends the Landsat legacy across the previous 43 years. In this study, we conducted three related analyses: (1) comparison of spectral values acquired by Landsats 7 and 8, separated by eight days, to ensure compatibility for time-series evaluation; (2) tracking of multitemporal signatures for different change processes across Landsat 5, 7, and 8 sensors using anniversary-date imagery; and (3) tracking the same type of processes using all available acquisitions. In this investigation, we found that data representing natural vegetation from Landsats 5, 7, and 8 were comparable and did not indicate a need for major modification prior to use for long-term monitoring. Analyses using anniversary-date imagery can be very effective for assessing long term patterns and trends occurring across the landscape, and are especially good for providing insights regarding trends related to long-term and continuous trends of growth or decline. We found that use of all available data provided a much more comprehensive level of understanding of the trends occurring, providing information about rate, duration, and intra- and inter-annual variability that could not be readily gleaned from the anniversary date analyses. We observed that using all available clear Landsat 5–8 observations with the new Continuous Change Detection and Classification (CCDC) algorithm was very effective for illuminating vegetation trends. There are a number of potential challenges for assessing gradual changes, including atmospheric impacts, algorithm development and visualization of the changes. One of the biggest challenges for studying gradual change will be the lack of appropriate data for validating results and products.

Remote Sensing of Environment

Lessons learned implementing an operational continuous United States national land change monitoring capability: The Land Change Monitoring, Assessment, and Projection (LCMAP) approach

Growing demands for temporally specific information on land surface change are fueling a new generation of maps and statistics that can contribute to understanding geographic and temporal patterns of change across large regions, provide input into a wide range of environmental modeling studies, clarify the drivers of change, and provide more timely information for land managers. To meet these needs, the U.S. Geological Survey has implemented a capability to monitor land surface change called the Land Change Monitoring, Assessment, and Projection (LCMAP) initiative. This paper describes the methodological foundations and lessons learned during development and testing of the LCMAP approach. Testing and evaluation of a suite of 10 annual land cover and land surface change data sets over six diverse study areas across the United States revealed good agreement with other published maps (overall agreement ranged from 73% to 87%) as well as several challenges that needed to be addressed to meet the goals of robust, repeatable, and geographically consistent monitoring results from the Continuous Change Detection and Classification (CCDC) algorithm. First, the high spatial and temporal variability of observational frequency led to differences in the number of changes identified, so CCDC was modified such that change detection is dependent on observational frequency. Second, the CCDC classification methodology was modified to improve its ability to characterize gradual land surface changes. Third, modifications were made to the classification element of CCDC to improve the representativeness of training data, which necessitated replacing the random forest algorithm with a boosted decision tree. Following these modifications, assessment of prototype Version 1 LCMAP results showed improvements in overall agreement (ranging from 85% to 90%).

Remote Sensing of Environment

Need and vision for global medium-resolution Landsat and Sentinel-2 data products

Global changes in climate and land use are threatening natural ecosystems, biodiversity, and the ecosystem services people rely on. This is why it is necessary to track and monitor spatiotemporal change at a level of detail that can inform science, management, and policy development. The current constellation of multiple Landsat and Sentinel-2 satellites collecting imagery at predominantly ≤ "> ≤ 30-m spatial resolution affords an opportunity for the generation of global medium- resolution products every few days. Our goal is to both identify the information needs and provide direction towards the generation of a suite of global, high-level, systematically-generated, medium-resolution products designed for both management and science. Our vision builds on the success of the NASA MODIS/VIIRS product suite, while recognizing the unique strengths of medium-resolution satellite data given their higher spatial resolution and longer time series. We propose a suite of 13 essential products that enable the characterization of the current state and changes in the biosphere, cryosphere, and hydrosphere, and would fill information needs identified by the Committee on Earth Observation Satellites for the Global Climate Observing System and the Global Terrestrial Observing System, by the National Research Council of the US National Academies in the decadal survey, and by others. These products are: land cover, land cover change, burned area, forest loss, vegetation indices, phenology, dynamic habitat indices, albedo, land surface temperature, snow cover, ice extent, surface water extent, and evapotranspiration. Furthermore, we provide a list of desirable products poised for addition to the essential products (e.g., crop type, emissivity, and ice sheet velocity). Lastly, we suggest aspirational products requiring further algorithm development (e.g., forest structure and crop yield). For the identified essential products, algorithms are in place, making it feasible to begin generating products systematically. These products should be accompanied by quality and accuracy assessments undertaken following consensus protocols. Five decades after the first Landsat satellite, and two decades after the MODIS products were first produced, it is time now for readily available, standardized, and consistent high-level products built upon medium-resolution imagery, thereby fulfilling the promise and the vision that inspired the Landsat program since its inception.

Remote Sensing of Environment

Tracking mangrove condition changes using dense Landsat time series

Mangroves in tropical and subtropical coasts are subject to episodic disturbances, notably from severe storms, leading to potential widespread vegetation mortality. The ability of vegetation to recover varies, and with disturbances becoming more frequent and severe, it is vital to track and project vegetation responses to support management and policy decisions. Prior studies have largely focused on binary mangrove mapping (i.e., presence or absence), while tracking conditions and condition change have not received sufficient attention. In this paper, we demonstrate a method based on dense time series Landsat images for continuous monitoring of mangrove conditions, where we track three kinds of post-disturbance mangrove conditions, including disturbed (disturbed, with rebound to the previous state within one growing season), recovering (undergoing natural recovery in longer than one growing season), and declining (showing long-term decline after disturbance). The method starts with disturbance detection using the DEtection and Characterization Of the tiDal wEtland change (DECODE) algorithm, an existing dense time series model designed to detect disturbances in tidal wetlands with adaptation to tidal fluctuations. This algorithm is well suited for the detection of tidal wetland disturbances but does not provide satisfactory post-disturbance monitoring results, due to the substantial variability in post-disturbance Landsat observations. To better monitor post-disturbance conditions, a new time series fitting approach, DECODER (DECODE and Recovery), is proposed for the recovery stage. Additionally, for temporal segments divided by disturbance events, we built a random forest classifier with temporal-spectral variables derived from the time series model to characterize mangrove conditions. Employing this approach in Florida's mangroves, we generated condition maps, such as dieback and recovery, with an overall accuracy of approximately 97.96 ± 0.86- [95 % confidence intervals]. Comparing post-hurricane conditions in Florida revealed that the increased frequency and severity of disturbances are challenging mangrove resilience, potentially diminishing their ability to recover and sustain ecosystem functions.

Florida

The Harmonized Landsat and Sentinel-2 version 2.0 surface reflectance dataset

Frequent multispectral observations of sufficient spatial detail from well-calibrated spaceborne sensors are needed for large-scale terrestrial monitoring. To meet this demand, the NASA Harmonized Landsat and Sentinel-2 (HLS) project was initiated in early 2010s to produce comparable 30-m surface reflectance from the US Landsat 8 Operational Land Imager (OLI) and the European Copernicus Sentinel-2A MultiSpectral Instrument (MSI), and currently from two OLI and two MSI sensors, by applying atmospheric correction to top-of-atmosphere (TOA) reflectance, masking out clouds and cloud shadows, normalizing bi-directional reflectance view angle effects, adjusting for sensor bandpass differences with the OLI as the reference, and providing the harmonized data in a common grid. Several versions of HLS dataset have been produced in the last few years. The recent improvements on almost all the harmonization algorithms had prompted a production of a new HLS dataset, tagged Version 2.0, which was completed in the summer of 2023 and for the first time takes on a global coverage (except for Antarctica). The HLS V2.0 data record starts in April 2013, two months after Landsat 8 launch. For 2022, the first whole year two Landsat and two Sentinel-2 satellites were available, HLS provides a global median of 66 cloud-free observations over land, substantially more than from a single sensor. This paper describes the HLS algorithm improvements and assesses the harmonization efficacy by examining how the reflectance difference between contemporaneous Landsat and Sentinel-2 observations was successively reduced by each harmonization step. The assessment was conducted on 545 pairs of globally distributed same-day Landsat/Sentinel-2 images from 2021 to 2022. Compared to the TOA data, the HLS atmospheric correction slightly increased the reflectance relative difference between Landsat and Sentinel-2 for most of the spectral bands, especially for the two blue bands and the green bands. The subsequent bi-directional reflectance view angle effect normalization effectively reduced the between-sensor reflectance difference present in the atmospherically corrected data for all the spectral bands, and notably to a level below the TOA differences for the red, near-infrared (NIR), and the two shortwave infrared (SWIR) bands. The bandpass adjustment only had a modest effect on reducing the between-sensor reflectance difference. In the final HLS products, the same-day reflectance difference between Landsat and Sentinel-2 was below 4.2% for the red, NIR, and the two SWIR bands, all smaller than the difference in the TOA data. However, the between-sensor differences for the two blue and the green bands remain slightly higher than in TOA data, and this reflects the difficulty in accurately correcting for atmospheric effects in the shorter wavelength visible bands. The data consistency evaluation on a suite of commonly used vegetation indices (VI) calculated from the HLS V2.0 reflectance data showed that the between-sensor VI difference is below 4.5% for most of the indices. These results suggest that the harmonization is robust and the HLS V2.0 data are adequate for quantitative terrestrial applications.

Remote Sensing of Environment

Enhancements to the USGS Landsat Level 2 surface temperature and emissivity product for Collection 3 reprocessing

The Landsat program provides the longest continuous global record of thermal infrared observations of the Earth's surface, underpinning critical applications in climate monitoring, water resources, ecosystem dynamics, urban heat analysis, and natural hazard assessment. The release of a global inventory of Landsat Collection 2 Level 2 surface temperature products by the U.S. Geological Survey (USGS) marked a major milestone in operational provision of Landsat thermal infrared analysis-ready data. Ongoing validations and community uses of Collection 2 have identified opportunities to further improve accuracy, uncertainty characterization, and emissivity correction across diverse atmospheric and surface conditions. In preparation for the planned Landsat Collection 3 reprocessing of the Landsat data record in the late 2020s, the USGS is implementing a coordinated set of enhancements to the Level 2 surface temperature products. These include revised emissivity estimation that leverages external datasets, improved atmospheric characterization and uncertainty propagation, expanded dynamic range for high temperature targets, consideration of split window atmospheric correction algorithm for Landsat 8 and 9, and decoupling of thermal infrared processing from visible to shortwave infrared constraints to enable surface temperature retrievals under low or no solar illumination conditions. These changes are designed to improve product quality and consistency across the Landsat record. Beyond near-term performance gains, the Collection 3 design establishes a scalable processing architecture to accommodate the expanded spectral and radiometric measurement capabilities of the forthcoming Landsat 10 mission. By preserving continuity across the Landsat 4–9 record while enabling future algorithm evolution, Landsat Collection 3 will provide a foundation for long-term, multi-decadal Earth system thermal infrared observations.

Remote Sensing of Environment

Geographic determination of Pinus ponderosa using DART TOFMS, ICP-MS, and LIBS handheld analyzer

Due to legal requirements on international imports, it is important for law enforcement and regulatory agencies to identify the geographical provenance of timber. Current methods for geographic identification utilize data generated by direct analysis in real time time-of-flight mass spectrometry (DART TOFMS), genetics, and isotope-ratio mass spectrometry (IRMS), but identification methods based on genetics and IRMS data require months to years to create usable databases. This study used machine learning algorithms to compare the results of DART TOFMS, inductively coupled plasma mass spectrometry (ICP-MS), and a handheld laser-induced breakdown spectroscopy (LIBS) analyzer for use in geographic identification of five populations of Pinus ponderosa spaced between 14 to 72 km apart. The results of the study showed comparable performances from machine learning algorithms applied to the ICP-MS and LIBS data with accuracy and kappa values over 90% while the DART TOFMS had an accuracy of 76% and a kappa value of 70%. This study demonstrated that data from the LIBS handheld analyzer is a viable and intriguing alternative to ICP-MS and DART TOFMS analyses in generating training databases and further indicates that trace elemental analysis via ICP-MS is a promising method for generating databases used to identify the origin of timber.

California, Oregon

Lineage diversification of fringe-toed lizards (Phrynosomatidae: Uma notata complex) in the Colorado Desert: Delimiting species in the presence of gene flow

Multi-locus nuclear DNA data were used to delimit species of fringe-toed lizards of the Uma notata complex, which are specialized for living in wind-blown sand habitats in the deserts of southwestern North America, and to infer whether Quaternary glacial cycles or Tertiary geological events were important in shaping the historical biogeography of this group. We analyzed ten nuclear loci collected using Sanger sequencing and genome-wide sequence and single-nucleotide polymorphism (SNP) data collected using restriction-associated DNA (RAD) sequencing. A combination of species discovery methods (concatenated phylogenies, parametric and non-parametric clustering algorithms) and species validation approaches (coalescent-based species tree/isolation-with-migration models) were used to delimit species, infer phylogenetic relationships, and to estimate effective population sizes, migration rates, and speciation times. Uma notata , U. inornata , U. cowlesi , and an undescribed species from Mohawk Dunes, Arizona ( U. sp.) were supported as distinct in the concatenated analyses and by clustering algorithms, and all operational taxonomic units were decisively supported as distinct species by ranking hierarchical nested speciation models with Bayes factors based on coalescent-based species tree methods. However, significant unidirectional gene flow (2 NM >1) from U. cowlesi and U. notata into U. rufopunctata was detected under the isolation-with-migration model. Therefore, we conservatively delimit four species-level lineages within this complex ( U. inornata , U. notata , U. cowlesi , and U. sp.), treating U. rufopunctata as a hybrid population ( U. notata x cowlesi ). Both concatenated and coalescent-based estimates of speciation times support the hypotheses that speciation within the complex occurred during the late Pleistocene, and that the geological evolution of the Colorado River delta during this period was an important process shaping the observed phylogeographic patterns.

Molecular Phylogenetics and Evolution

Estimation of hectare-scale soil-moisture characteristics from aquifer-test data

Analysis of a 72-h, constant-rate aquifer test conducted in a coarse-grained and highly permeable, glacial outwash deposit on Cape Cod, Massachusetts revealed that drawdowns measured in 20 piezometers located at various depths below the water table and distances from the pumped well were significantly influenced by effects of drainage from the vadose zone . The influence was greatest in piezometers located close to the water table and diminished with increasing depth. The influence of the vadose zone was evident from a gap, in the intermediate-time zone, between measured drawdowns and drawdowns computed under the assumption that drainage from the vadose zone occurred instantaneously in response to a decline in the elevation of the water table. By means of an analytical model that was designed to account for time-varying drainage, simulated drawdowns could be closely fitted to measured drawdowns regardless of the piezometer locations. Because of the exceptional quality and quantity of the data and the relatively small aquifer heterogeneity, it was possible by inverse modeling to estimate all relevant aquifer parameters and a set of three empirical constants used in the upper-boundary condition to account for the dynamic drainage process. The empirical constants were used to define a one-dimensional (1D) drainage versus time curve that is assumed to be representative of the bulk material overlying the water table. The curve was inverted with a parameter estimation algorithm and a 1D numerical model for variably saturated flow to obtain soil-moisture retention curves and unsaturated hydraulic conductivity relationships defined by the Brooks and Corey equations. Direct analysis of the aquifer-test data using a parameter estimation algorithm and a two-dimensional, axisymmetric numerical model for variably saturated flow yielded similar soil-moisture characteristics. Results suggest that hectare-scale soil-moisture characteristics are different from core-scale predictions and even relatively small amounts of fine-grained material and heterogeneity can dominate the large-scale soil-moisture characteristics and aquifer response.

Journal of Hydrology

Multi-level slug tests in highly permeable formations: 2. Hydraulic conductivity identification, method verification, and field applications

Using the developed theory and modified Springer-Gelhar (SG) model, an identification method is proposed for estimating hydraulic conductivity from multi-level slug tests. The computerized algorithm calculates hydraulic conductivity from both monotonic and oscillatory well responses obtained using a double-packer system. Field verification of the method was performed at a specially designed fully penetrating well of 0.1-m diameter with a 10-m screen in a sand and gravel alluvial aquifer (MSEA site, Shelton, Nebraska). During well installation, disturbed core samples were collected every 0.6 m using a split-spoon sampler. Vertical profiles of hydraulic conductivity were produced on the basis of grain-size analysis of the disturbed core samples. These results closely correlate with the vertical profile of horizontal hydraulic conductivity obtained by interpreting multi-level slug test responses using the modified SG model. The identification method was applied to interpret the response from 474 slug tests in 156 locations at the MSEA site. More than 60% of responses were oscillatory. The method produced a good match to experimental data for both oscillatory and monotonic responses using an automated curve matching procedure. The proposed method allowed us to drastically increase the efficiency of each well used for aquifer characterization and to process massive arrays of field data. Recommendations generalizing this experience to massive application of the proposed method are developed.Using the developed theory and modified Springer-Gelhar (SG) model, an identification method is proposed for estimating hydraulic conductivity from multi-level slug tests. The computerized algorithm calculates hydraulic conductivity from both monotonic and oscillatory well responses obtained using a double-packer system. Field verification of the method was performed at a specially designed fully penetrating well of 0.1-m diameter with a 10-m screen in a sand and gravel alluvial aquifer (MSEA site, Shelton, Nebraska). During well installation, disturbed core samples were collected every 0.6 m using a split-spoon sampler. Vertical profiles of hydraulic conductivity were produced on the basis of grain-size analysis of the disturbed core samples. These results closely correlate with the vertical profile of horizontal hydraulic conductivity obtained by interpreting multi-level slug test responses using the modified SG model. The identification method was applied to interpret the response from 474 slug tests in 156 locations at the MSEA site. More than 60% of responses were oscillatory. The method produced a good match to experimental data for both oscillatory and monotonic responses using an automated curve matching procedure. The proposed method allowed us to drastically increase the efficiency of each well used for aquifer characterization and to process massive arrays of field data. Recommendations generalizing this experience to massive application of the proposed method are developed.

Journal of Hydrology

Classification methods for monitoring Arctic sea ice using OKEAN passive/active two-channel microwave data

This paper presents methods for classifying Arctic sea ice using both passive and active (2-channel) microwave imagery acquired by the Russian OKEAN 01 polar-orbiting satellite series. Methods and results are compared to sea ice classifications derived from nearly coincident Special Sensor Microwave Imager (SSM/I) and Advanced Very High Resolution Radiometer (AVHRR) image data of the Barents, Kara, and Laptev Seas. The Russian OKEAN 01 satellite data were collected over weekly intervals during October 1995 through December 1997. Methods are presented for calibrating, georeferencing and classifying the raw active radar and passive microwave OKEAN 01 data, and for correcting the OKEAN 01 microwave radiometer calibration wedge based on concurrent 37 GHz horizontal polarization SSM/I brightness temperature data. Sea ice type and ice concentration algorithms utilized OKEAN's two-channel radar and passive microwave data in a linear mixture model based on the measured values of brightness temperature and radar backscatter, together with a priori knowledge about the scattering parameters and natural emissivities of basic sea ice types. OKEAN 01 data and algorithms tended to classify lower concentrations of young or first-year sea ice when concentrations were less than 60%, and to produce higher concentrations of multi-year sea ice when concentrations were greater than 40%, when compared to estimates produced from SSM/I data. Overall, total sea ice concentration maps derived independently from OKEAN 01, SSM/I, and AVHRR satellite imagery were all highly correlated, with uniform biases, and mean differences in total ice concentration of less than four percent (sd<15%).

Remote Sensing of Environment

Comments on "Failures in detecting volcanic ash from a satellite-based technique"

The recent paper by Simpson et al. [Remote Sens. Environ. 72 (2000) 191.] on failures to detect volcanic ash using the 'reverse' absorption technique provides a timely reminder of the danger that volcanic ash presents to aviation and the urgent need for some form of effective remote detection. The paper unfortunately suffers from a fundamental flaw in its methodology and numerous errors of fact and interpretation. For the moment, the 'reverse' absorption technique provides the best means for discriminating volcanic ash clouds from meteorological clouds. The purpose of our comment is not to defend any particular algorithm; rather, we point out some problems with Simpson et al.'s analysis and re-state the conditions under which the 'reverse' absorption algorithm is likely to succeed. ?? 2001 Elsevier Science Inc. All rights reserved.

Remote Sensing of Environment

Singularity and Nonnormality in the Classification of Compositional Data

Geologists may want to classify compositional data and express the classification as a map. Regionalized classification is a tool that can be used for this purpose, but it incorporates discriminant analysis, which requires the computation and inversion of a covariance matrix. Covariance matrices of compositional data always will be singular (noninvertible) because of the unit-sum constraint. Fortunately, discriminant analyses can be calculated using a pseudo-inverse of the singular covariance matrix; this is done automatically by some statistical packages such as SAS. Granulometric data from the Darss Sill region of the Baltic Sea is used to explore how the pseudo-inversion procedure influences discriminant analysis results, comparing the algorithm used by SAS to the more conventional Moore-Penrose algorithm. Logratio transforms have been recommended to overcome problems associated with analysis of compositional data, including singularity. A regionalized classification of the Darss Sill data after logratio transformation is different only slightly from one based on raw granulometric data, suggesting that closure problems do not influence severely regionalized classification of compositional data.

Mathematical Geology

Pyroclastic flow hazard at Volcán Citlaltépetl

Volc&aacute;n Citlalt&eacute;petl (Pico de Orizaba) with an elevation of 5,675 m is the highest volcano in North America. Its most recent catastrophic events involved the production of pyroclastic flows that erupted approximately 4,000, 8,500, and 13,000 years ago. The distribution of mapped deposits from these eruptions gives an approximate guide to the extent of products from potential future eruptions. Because the topography of this volcano is constantly changing computer simulations were made on the present topography using three computer algorithms: energy cone, FLOW2D, and FLOW3D. The Heim Coefficient (&mu;), used as a code parameter for frictional sliding in all our algorithms, is the ratio of the assumed drop in elevation (H) divided by the lateral extent of the mapped deposits (L). The viscosity parameter for the FLOW2D and FLOW3D codes was adjusted so that the paths of the flows mimicked those inferred from the mapped deposits. We modeled two categories of pyroclastic flows modeled for the level I and level II events. Level I pyroclastic flows correspond to small but more frequent block-and-ash flows that remain on the main cone. Level II flows correspond to more widespread flows from catastrophic eruptions with an approximate 4,000-year repose period. We developed hazard maps from simulations based on a National Imagery and Mapping Agency ( NIMA ) DTED-1 DEM with a 90 m grid and a vertical accuracy of &plusmn;30 m. Because realistic visualization is an important aid to understanding the risks related to volcanic hazards we present the DEM as modeled by FLOW3D. The model shows that the pyroclastic flows extend for much greater distances to the east of the volcano summit where the topographic relief is nearly 4,300 m. This study was used to plot hazard zones for pyroclastic flows in the official hazard map that was published recently.

Natural Hazards

Extracting low‐resolution river networks from high‐resolution digital elevation models

Including a global river network in the land component of global climate models (GCMs) is necessary in order to provide a more complete representation of the hydrologic cycle. The process of creating these networks is called river network upscaling and consists of lowering the resolution of already available fine networks to make them compatible with GCMs. Fine‐resolution river networks have a level of detail appropriate for analysis at the watershed scale but are too intensive for global hydrologic studies. A river network upscaling algorithm, which processes fine‐resolution digital elevation models to determine the flow directions that best describe the flow patterns in a coarser user‐defined scale, is presented. The objectives of this study were to develop an algorithm that advances the previous work in the field by being applicable at a global scale, allowing for the upscaling to be performed in a projected environment, and generating evenly distributed flow directions.

Water Resources Research

A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

[1] We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green's functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake ( Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth.

Tottori

Interferogram formation in the presence of complex and large deformation

Sierra Negra volcano in Isabela island, Galápagos, erupted from October 22 to October 30 in 2005. During the 8 days of eruption, the center of Sierra Negra's caldera subsided about 5.4 meters. Three hours prior to the onset of the eruption, an earthquake (Mw 5.4) occurred, near the caldera. Because of the large and complex phase gradient due to the huge subsidence and the earthquake, it is difficult to form an interferogram inside the caldera that spans the eruption. The deformation is so large and spatially variable that the approximations used in existing InSAR software (ROI, ROI_PAC, DORIS, GAMMA) cannot properly coregister SAR image pairs spanning the eruption. We have developed here a two-step algorithm that can form intra-caldera interferograms from these data. The first step involves a "rubber-sheeting" SAR image coregistration. In the second step we use range offset estimates to mitigate the steep phase gradient. Using this new algorithm, we retrieve an interferogram with the best coverage to date inside the caldera of Sierra Negra.

Isabela island, Sierra Negra volcano