Geology Reports⌕ Search

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 991 records · Page 55Linked to original sources
USGS2021

Preliminary assessment of the geometric improvements to the Landsat Collection-2 archive

The U.S. Geological Survey (USGS) has completed processing of the historical Landsat archive to Collection-2 as of December of 2020 and has released it to the public. As part of Collection-2, several geometric changes have been implemented, including changes to the ground control points (GCPs) and elevation datasets. These datasets are used as a geometric reference for all missions. In addition, mission specific improvements were included in Collection-2, such as improvements to the precision correction algorithms and updates of the calibration parameters. This paper discusses a preliminary analysis of a comparison between the Collection-1 and Collection-2 products of the entire Landsat archive. Compared to the Level 1 products in Collection-1, the number of Tier-1 precision- and terrain-corrected (L1TP) Level 1 products in Collection-2 increased by 6.26% across all sensors. Landsat 8 products showed an increase of Tier-1 L1TP products by 5.33%; Landsat 7 products showed an increase of Tier-1 products by 6.94%; and Landsat 5 and Landsat 4 Thematic Mapper (TM) products showed an increase of Tier-1 products by 9.35%. The geometric accuracy of the Tier-1 terrain corrected products also improved by 2 meters or more.

Conference Paper↗
USGS2000

Landsat 7 on-orbit geometric calibration and performance

The Landsat 7 Image Assessment System was developed to characterize and calibrate the radiometric and geometric performance of the Landsat 7 Enhanced Thematic Mapper Pius (ETM+) instrument. Algorithms and software assess the geometric performance of the Landsat 7 spacecraft and ETM+ sensor system and perform geometric calibration by estimating sensor and spacecraft geometric parameters. Following the initial on-orbit calibration, performed during the Landsat 7 on- orbit initialization and verification period, all geometric performance goals were met. Geometric characterization and calibration activities will continue for the life of the Landsat 7 mission.

Conference Paper↗
USGS2002

Absolute irradiance of the Moon for on-orbit calibration

The recognized need for on-orbit calibration of remote sensing imaging instruments drives the ROLO project effort to characterize the Moon for use as an absolute radiance source. For over 5 years the ground-based ROLO telescopes have acquired spatially-resolved lunar images in 23 VNIR (Moon diameter ???500 pixels) and 9 SWIR (???250 pixels) passbands at phase angles within ??90 degrees. A numerical model for lunar irradiance has been developed which fits hundreds of ROLO images in each band, corrected for atmospheric extinction and calibrated to absolute radiance, then integrated to irradiance. The band-coupled extinction algorithm uses absorption spectra of several gases and aerosols derived from MODTRAN to fit time-dependent component abundances to nightly observations of standard stars. The absolute radiance scale is based upon independent telescopic measurements of the star Vega. The fitting process yields uncertainties in lunar relative irradiance over small ranges of phase angle and the full range of lunar libration well under 0.5%. A larger source of uncertainty enters in the absolute solar spectral irradiance, especially in the SWIR, where solar models disagree by up to 6%. Results of ROLO model direct comparisons to spacecraft observations demonstrate the ability of the technique to track sensor responsivity drifts to sub-percent precision. Intercomparisons among instruments provide key insights into both calibration issues and the absolute scale for lunar irradiance.

Conference Paper↗
USGS2004

Digital elevation extraction from multiple MTI data sets

The Digital Elevation Model (DEM) extraction process traditionally uses a stereo pair of aerial photographs that are sequentially captured using an airborne metric camera. Standard DEM extraction techniques have been naturally extended to utilize satellite imagery. However, the particular characteristics of satellite imaging can cause difficulties in the DEM extraction process. The ephemeris of the spacecraft during the collects, with respect to the ground test site, is the most important factor in the elevation extraction process. When the angle of separation between the stereo images is small, the extraction process typically produces measurements with low accuracy. A large angle of separation can cause an excessive number of erroneous points in the output DEM. There is also a possibility of having occluded areas in the images when drastic topographic variation is present, making it impossible to calculate elevation in the blind spots. The use of three or more images registered to the same ground area can potentially reduce these problems and improve the accuracy of the extracted DEM. The pointing capability of the Multispectral Thermal Imager (MTI) allows for multiple collects of the same area to be taken from different perspectives. This functionality of MTI makes it a good candidate for the implementation of DEM extraction using multiple images for improved accuracy. This paper describes a project to evaluate this capability and the algorithms used to extract DEMs from multi-look MTI imagery.

Conference Paper↗
USGS2006

Using self-organizing maps to determine observation threshold limit predictions in highly variant data

A significant data quality challenge for highly variant systems surrounds the limited ability to quantify operationally reasonable limits on the data elements being collected and provide reasonable threshold predictions. In many instances, the number of influences that drive a resulting value or operational range is too large to enable physical sampling for each influencer, or is too complicated to accurately model in an explicit simulation. An alternative method to determine reasonable observation thresholds is to employ an automation algorithm that would emulate a human analyst visually inspecting data for limits. Using the visualization technique of self-organizing maps (SOM) on data having poorly understood relationships, a methodology for determining threshold limits was developed. To illustrate this approach, analysis of environmental influences that drive the abundance of a target indicator species (the pink shrimp, Farfantepenaeus duorarum) provided a real example of applicability. The relationship between salinity and temperature and abundance of F. duorarum is well documented, but the effect of changes in water quality upstream on pink shrimp abundance is not well understood. The highly variant nature surrounding catch of a specific number of organisms in the wild, and the data available from up-stream hydrology measures for salinity and temperature, made this an ideal candidate for the approach to provide a determination about the influence of changes in hydrology on populations of organisms.

Conference Paper↗
USGS2007

Consistency of L4 TM absolute calibration with respect to the L5 TM sensor based on near-simultaneous image acquisition

The Landsat archive provides more than 35 years of uninterrupted multispectral remotely sensed data of Earth observations. Since 1972, Landsat missions have carried different types of sensors, from the Return Beam Vidicon (RBV) camera to the Enhanced Thematic Mapper Plus (ETM+). However, the Thematic Mapper (TM) sensors on Landsat 4 (L4) and Landsat 5 (L5), launched in 1982 and 1984 respectively, are the backbone of an extensive archive. Effective April 2, 2007, the radiometric calibration of L5 TM data processed and distributed by the U.S. Geological Survey (USGS) Center for Earth Resources Observation and Science (EROS) was updated to use an improved lifetime gain model, based on the instrument's detector response to pseudo-invariant desert site data and cross-calibration with the L7 ETM+. However, no modifications were ever made to the radiometric calibration procedure of the Landsat 4 (L4) TM data. The L4 TM radiometric calibration procedure has continued to use the Internal Calibrator (IC) based calibration algorithms and the post calibration dynamic ranges, as previously defined. To evaluate the "current" absolute accuracy of these two sensors, image pairs from the L5 TM and L4 TM sensors were compared. The number of coincident image pairs in the USGS EROS archive is limited, so the scene selection for the cross-calibration studies proved to be a challenge. Additionally, because of the lack of near-simultaneous images available over well-characterized and traditionally used calibration sites, alternate sites that have high reflectance, large dynamic range, high spatial uniformity, high sun elevation, and minimal cloud cover were investigated. The alternate sites were identified in Yuma, Iraq, Egypt, Libya, and Algeria. The cross-calibration approach involved comparing image statistics derived from large common areas observed eight days apart by the two sensors. This paper summarizes the average percent differences in reflectance estimates obtained between the two sensors. The work presented in this paper is a first step in understanding the current performance of L4 TM absolute calibration and potentially serves as a platform to revise and improve the radiometric calibration procedures implemented for the processing of L4 TM data.

Conference Paper↗
USGS2008

MODIS and SeaWIFS on-orbit lunar calibration

The Moon plays an important role in the radiometric stability monitoring of the NASA Earth Observing System's (EOS) remote sensors. The MODIS and SeaWIFS are two of the key instruments for NASA's EOS missions. The MODIS Protoflight Model (PFM) on-board the Terra spacecraft and the MODIS Flight Model 1 (FM1) on-board the Aqua spacecraft were launched on December 18, 1999 and May 4, 2002, respectively. They view the Moon through the Space View (SV) port approximately once a month to monitor the long-term radiometric stability of their Reflective Solar Bands (RSB). SeaWIFS was launched on-board the OrbView-2 spacecraft on August 1, 1997. The SeaWiFS lunar calibrations are obtained once a month at a nominal phase angle of 7??. The lunar irradiance observed by these instruments depends on the viewing geometry. The USGS photometric model of the Moon (the ROLO model) has been developed to provide the geometric corrections for the lunar observations. For MODIS, the lunar view responses with corrections for the viewing geometry are used to track the gain change for its reflective solar bands (RSB). They trend the system response degradation at the Angle Of Incidence (AOI) of sensor's SV port. With both the lunar observation and the on-board Solar Diffuser (SD) calibration, it is shown that the MODIS system response degradation is wavelength, mirror side, and AOI dependent. Time-dependent Response Versus Scan angle (RVS) Look-Up Tables (LUT) are applied in MODIS RSB calibration and lunar observations play a key role in RVS derivation. The corrections provided by the RVS in the Terra and Aqua MODIS data from the 412 nm band are as large as 16% and 13%, respectively. For SeaWIFS lunar calibrations, the spacecraft is pitched across the Moon so that the instrument views the Moon near nadir through the same optical path as it views the Earth. The SeaWiFS system gain changes for its eight bands are calibrated using the geometrically-corrected lunar observations. The radiometric corrections to the SeaWiFS data, after more than ten years on orbit, are 19% at 865 nm, 8% at 765 nm, and 1-3% in the other bands. In this report, the lunar calibration algorithms are reviewed and the RSB gain changes observed by the lunar observations are shown for all three sensors. The lunar observations for the three instruments are compared using the USGS photometric model. The USGS lunar model facilitates the cross calibration of instruments with different spectra bandpasses whose measurements of the Moon differ in time and observing geometry.

Conference Paper↗
USGS2012

Spectroscopic remote sensing for material identification, vegetation characterization, and mapping

Identifying materials by measuring and analyzing their reflectance spectra has been an important procedure in analytical chemistry for decades. Airborne and space-based imaging spectrometers allow materials to be mapped across the landscape. With many existing airborne sensors and new satellite-borne sensors planned for the future, robust methods are needed to fully exploit the information content of hyperspectral remote sensing data. A method of identifying and mapping materials using spectral feature analyses of reflectance data in an expert-system framework called MICA (Material Identification and Characterization Algorithm) is described. MICA is a module of the PRISM (Processing Routines in IDL for Spectroscopic Measurements) software, available to the public from the U.S. Geological Survey (USGS) at http://pubs.usgs.gov/of/2011/1155/. The core concepts of MICA include continuum removal and linear regression to compare key diagnostic absorption features in reference laboratory/field spectra and the spectra being analyzed. The reference spectra, diagnostic features, and threshold constraints are defined within a user-developed MICA command file (MCF). Building on several decades of experience in mineral mapping, a broadly-applicable MCF was developed to detect a set of minerals frequently occurring on the Earth's surface and applied to map minerals in the country-wide coverage of the 2007 Afghanistan HyMap data set. MICA has also been applied to detect sub-pixel oil contamination in marshes impacted by the Deepwater Horizon incident by discriminating the C-H absorption features in oil residues from background vegetation. These two recent examples demonstrate the utility of a spectroscopic approach to remote sensing for identifying and mapping the distributions of materials in imaging spectrometer data.

Conference Paper↗
USGS2025

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗
USGS2005

Pattern-oriented modeling of agent-based complex systems: Lessons from ecology

Agent-based complex systems are dynamic networks of many interacting agents; examples include ecosystems, financial markets, and cities. The search for general principles underlying the internal organization of such systems often uses bottom-up simulation models such as cellular automata and agent-based models. No general framework for designing, testing, and analyzing bottom-up models has yet been established, but recent advances in ecological modeling have come together in a general strategy we call pattern-oriented modeling. This strategy provides a unifying framework for decoding the internal organization of agent-based complex systems and may lead toward unifying algorithmic theories of the relation between adaptive behavior and system complexity.

Science↗
USGS2020

3D fault architecture controls the dynamism of earthquake swarm

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. We leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.

Science↗
USGS2025

An empirical test of maximum depositional age through paired LA-ICP-MS and CA-ID-TIMS detrital zircon analysis with implications for the “Nevadan orogeny” in the Sierra Nevada foothills, California, U.S.A.

The interpretation of maximum depositional age (MDA) from U-Pb detrital zircon data acquired via laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) or secondary ion mass spectrometry (SIMS) is now routine; however, to date only a few studies have presented tests of such MDAs by subsequently analyzing a subset of the same detrital zircon grains with the more accurate and precise chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) method. We first generated LA-ICP-MS dates and MDAs from three turbidite sandstone samples of the Late Jurassic Mariposa Formation (Sierra Nevada foothills, California) containing large proportions of young detrital zircon grains (at the time of sediment deposition). We then removed 4–5 of the apparent youngest grains from the epoxy mount and analyzed them with CA-ID-TIMS. Despite a lack of low Th/U or high U (ppm) characteristics that might indicate U-Pb system disturbance, all LA-ICP-MS dates are younger than corresponding CA-ID-TIMS dates, on the same grain, by up to 11.8 m.y. when considered as point estimates. Only two out of 13 paired analyses overlap within error at 95% confidence. We interpret MDAs from our CA-ID-TIMS dates to be the youngest CA-ID-TIMS dates (YTDs). LA-ICP-MS MDAs based on fewer dates from the young tail of the youngest date distribution perform the worst (i.e., furthest from CA-ID-TIMS MDAs), whereas LA-ICP-MS MDA methods that incorporate more dates from the youngest date distribution perform better; the best performing method in all cases is the maximum likelihood algorithm–minimum. The performance of tested LA-ICP-MS MDA methods is improved by removing statistical outlier dates and by removing visually young “outlier” dates that drift away from the primary distribution of young dates. Our paired LA-ICP-MS/CA-ID-TIMS MDA workflow shows that the accuracy of MDA can be significantly improved by conducting CA-ID-TIMS on as few as four grains from a sample. The combination of our new CA-ID-TIMS MDAs with published CA-ID-TIMS analyses and petrochronology on the nearby Guadalupe igneous complex provides sub-million-year resolution of contemporaneous igneous and sedimentary systems during deformation and clarifies the timing of regional deformation that defines the local “Nevadan orogeny.” Rocks in the upper plate of the Bear Mountains fault zone have a CA-ID-TIMS YTD MDA of 151.71 ± 0.23 Ma, and rocks of the lower plate have a CA-ID-TIMS YTD MDA of 149.92 ± 0.11 Ma. These new MDAs, in conjunction with the observation of fabric-bearing Mariposa Formation xenoliths in the ca. 149.65 ± 0.10 Ma Guadalupe igneous complex, suggest that rocks in the upper plate of the Bear Mountains fault zone represent a slightly older (~2 m.y.) section of Mariposa Formation that was deformed and intruded prior to being juxtaposed against, and further deformed with, a slightly younger lower-plate section of Mariposa Formation in actively deforming, fault-bounded basins. Our observations are not consistent with traditional models that require that sedimentation of the Mariposa Formation ended by ca. 155 Ma. Instead, we interpret our data to be consistent with other evidence for a continuum of deformation in Late Jurassic to Early Cretaceous time and document that the regional “slatey cleavage” observed in the greater Mariposa Formation and used to define the “Nevadan orogeny” in our study area is largely younger than 149.92 ± 0.11 Ma.

California↗
USGS2016

Persistent slip rate discrepancies in the eastern California (USA) shear zone

Understanding fault slip rates in the eastern California shear zone (ECSZ) using GPS geodesy is complicated by potentially overlapping strain signals due to many sub-parallel strike-slip faults and by inconsistencies with geologic slip rates. The role of fault system geometry in describing ECSZ deformation may be investigated with total variation regularization, which algorithmically determines a best-fitting geometry from an initial model with numerous faults, constrained by a western United States GPS velocity field. The initial dense model (1) enables construction of the first geodetically constrained block model to include all ECSZ faults with geologic slip rates, allowing direct geologic-geodetic slip rate comparisons, and (2) permits fault system geometries with many active faults that are analogous to distributed interseismic deformation. Beginning with 58 ECSZ blocks, a model containing 10 ECSZ blocks is most consistent with geologic slip rates, reproducing five of 11 within their reported uncertainties. The model fits GPS observations with a mean residual velocity of 1.5 mm/yr. Persistent geologic-geodetic slip rate discrepancies occur on the Calico and Garlock faults, on which we estimate slip rates of 7.6 mm/yr and <2 mm/yr, respectively, indicating that inconsistencies between geology and geodesy may be concentrated on or near these faults and are not due to pervasive distributed deformation in the region. Discrepancies may in part be due to postseismic relaxation following the A.D. 1992 M w 7.3 Landers and 1999 M w 7.1 Hector Mine earthquakes. Otherwise, resolving geologic-geodetic discrepancies would require as much as 11.4 mm/yr of off-fault deformation within <10 km of the main ECSZ faults, with ∼5 mm/yr concentrated near the Calico fault.

California↗
USGS2011

Surface (sea floor) and near-surface (box cores) sediment mineralogy in Baffin Bay as a key to sediment provenance and ice sheet variations

To better understand the glacial history of the ice sheets surrounding Baffin Bay and to provide information on sediment pathways, samples from 82 seafloor grabs and core tops, and from seven box cores were subjected to quantitative X-ray diffraction weight percent (wt.%) analysis of the <2 mm sediment fraction. The samples were collected between 67°N and 78°N, in water depths of 155 to 2375 m and were retrieved on cruises between A.D. 1964 and 2009. Grain size, magnetic characteristics, and colour reflectance data were also obtained on many of the samples. Twenty-one non-clay and 10 clay mineral species were identified; the average wt.% of the non-clay minerals was 70% and was dominated by quartz, various feldspars, and dolomite, whereas the dominant clay minerals were 1 M illite, biotite, and chlorite. Cluster analysis on principal component scores identified three main mineral groups, which also had strong associations with grain size and sediment magnetic properties. Box cores from the deep central basin (>2000 m) all show an abrupt drop in calcite wt.% (post-5 cal ka BP?) following a major peak in detrital carbonate (mainly dolomite). This dolomite-rich detrital carbonate (DC) event in JR175BC06 is possibly coeval with the Younger Dryas cold event. Four possible glacial-sourced end members were employed in a compositional unmixing algorithm to gain insight into down core changes in sediment provenance at the deep central basin. Estimates of the rates of sediment accumulation in the central basin are only in the range of 2 to 4 cm/cal ka, surprisingly low given the glaciated nature of the surrounding land.

Canadian Journal of Earth Sciences↗
USGS2022

Modeling morphodynamics of coastal response to extreme events: What shape are we in?

This review focuses on recent advances in process-based numerical models of the impact of extreme storms on sandy coasts. Driven by larger-scale models of meteorology and hydrodynamics, these models simulate morphodynamics across the Sallenger storm-impact scale, including swash, collision, overwash, and inundation. Models are becoming both wider (as more processes are added) and deeper (as detailed physics replaces earlier parameterizations). Algorithms for wave-induced flows and sediment transport under shoaling waves are among the recent developments. Community and open-source models have become the norm. Observations of initial conditions (topography, land cover, and sediment characteristics) have become more detailed, and improvements in tropical cyclone and wave models provide forcing (winds, waves, surge, and upland flow) that is better resolved and more accurate, yielding commensurate improvements in model skill. We foresee that future storm-impact models will increasingly resolve individual waves, apply data assimilation, and be used in ensemble modeling modes to predict uncertainties.

Annual Review of Marine Science↗
USGS2004

Accuracy of a pulse-coherent acoustic Doppler profiler in a wave-dominated flow

The accuracy of velocities measured by a pulse-coherent acoustic Doppler profiler (PCADP) in the bottom boundary layer of a wave-dominated inner-shelf environment is evaluated. The downward-looking PCADP measured velocities in eight 10-cm cells at 1 Hz. Velocities measured by the PCADP are compared to those measured by an acoustic Doppler velocimeter for wave orbital velocities up to 95 cm s-1 and currents up to 40 cm s-1. An algorithm for correcting ambiguity errors using the resolution velocities was developed. Instrument bias, measured as the average error in burst mean speed, is -0.4 cm s-1 (standard deviation = 0.8). The accuracy (root-mean-square error) of instantaneous velocities has a mean of 8.6 cm s-1 (standard deviation = 6.5) for eastward velocities (the predominant direction of waves), 6.5 cm s-1 (standard deviation = 4.4) for northward velocities, and 2.4 cm s-1 (standard deviation = 1.6) for vertical velocities. Both burst mean and root-mean-square errors are greater for bursts with ub ??? 50 cm s-1. Profiles of burst mean speeds from the bottom five cells were fit to logarithmic curves: 92% of bursts with mean speed ??? 5 cm s-1 have a correlation coefficient R2 > 0.96. In cells close to the transducer, instantaneous velocities are noisy, burst mean velocities are biased low, and bottom orbital velocities are biased high. With adequate blanking distances for both the profile and resolution velocities, the PCADP provides sufficient accuracy to measure velocities in the bottom boundary layer under moderately energetic inner-shelf conditions.

Journal of Atmospheric and Oceanic Technology↗
USGS2015

Characterizing Congo Basin rainfall and climate using Tropical Rainfall Measuring Mission (TRMM) satellite data and limited rain gauge ground observations

Quantitative understanding of Congo River basin hydrological behavior is poor because of the basin’s limited hydrometeorological observation network. In cases such as the Congo basin where ground data are scarce, satellite-based estimates of rainfall, such as those from the joint NASA/JAXA Tropical Rainfall Measuring Mission (TRMM), can be used to quantify rainfall patterns. This study tests and reports the use of limited rainfall gauge data within the Democratic Republic of Congo (DRC) to recalibrate a TRMM science product (TRMM 3B42, version 6) in characterizing precipitation and climate in the Congo basin. Rainfall estimates from TRMM 3B42, version 6, are compared and adjusted using ground precipitation data from 12 DRC meteorological stations from 1998 to 2007. Adjustment is achieved on a monthly scale by using a regression-tree algorithm. The output is a new, basin-specific estimate of monthly and annual rainfall and climate types across the Congo basin. This new product and the latest version-7 TRMM 3B43 science product are validated by using an independent long-term dataset of historical isohyets. Standard errors of the estimate, root-mean-square errors, and regression coefficients r were slightly and uniformly better with the recalibration from this study when compared with the 3B43 product (mean monthly standard errors of 31 and 40 mm of precipitation and mean r 2 of 0.85 and 0.82, respectively), but the 3B43 product was slightly better in terms of bias estimation (1.02 and 1.00). Despite reasonable doubts that have been expressed in studies of other tropical regions, within the Congo basin the TRMM science product (3B43) performed in a manner that is comparable to the performance of the recalibrated product that is described in this study.

Journal of Applied Meteorology and Climatology↗
USGS2019

A high-resolution 1983-2016 Tmax climate data record based on InfraRed Temperatures and Stations by the Climate Hazard Center

Understanding the dynamics and physics of climate extremes will be a critical challenge for 21st century climate science. Increasing temperatures and saturation vapor pressures may exacerbate heat waves, droughts and precipitation extremes. Yet our ability to monitor temperature variations is limited and declining. Between 1983 and 2016 the number of observations in the CRU Tmax product declined precipitously (5,900→1,000); 1,000 poorly distributed measurements are not enough to resolve regional Tmax variations. Here we show that long (1983-near present) high resolution (0.05°), cloud-screened archives of geostationary satellite Thermal Infrared Red (TIR) observations, when combined with a dense set of ~15,000 station observations, can explain 23, 40, 30, 11% more variance than the CRU over the globe, South America, Africa, India and areas north of 50°N, with even greater levels of improvement over the 2011-2016 period (28, 45, 39, 52, 28%). Described here for the first time, the TIR Tmax algorithm uses sub-daily TIR distributions to screen out cloud contaminated observations, providing accurate (correlation≈0.8) gridded emission Tmax estimates. Blending these gridded fields with ~15,000 station observations provides a seamless, high-resolution source of accurate Tmax estimates that performs well in areas lacking dense in situ observations and even better where in situ observations are available. Cross-validation results indicates that the satellite-only, station-only and combined products all perform accurately (R≈0.8-0.9, mean absolute errors ≈0.8-1.0). Hence, the Climate Hazards center InfraRed Temperature with Stations (CHIRTSmax) data set should provide a valuable resource for climate change studies, climate extreme analyses, and early warning applications.

Journal of Climate↗