Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Names”

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 775 records · Page 43Linked to original sources

Climate warming shifts carbon allocation from stemwood to roots in calcium-depleted spruce forests

Increased greening of northern forests, measured by the Normalized Difference Vegetation Index (NDVI), has been presented as evidence that a warmer climate has increased both net primary productivity (NPP) and the carbon sink in boreal forests. However, higher production and greener canopies may accompany changes in carbon allocation that favor foliage or fine roots over less decomposable woody biomass. Furthermore, tree core data throughout mid- and northern latitudes have revealed a divergence problem (DP), a weakening in tree ring responses to warming over the past half century that is receiving increasing attention, but remains poorly understood. Often, the same sites exhibit trend inconsistency phenomenon (TIP), namely positive, or no trends in growing season NDVI where negative trends in tree ring indexes are observed. Here we studied growth of two Norway spruce (Picea abies) stands in western Russia that exhibited both the DP and TIP but were subject to soil acidification and calcium depletion of differing timing and severity. Our results link the decline in radial growth starting in 1980 to a shift in carbon allocation from wood to roots driven by a combination of two factors: (a) soil acidification that depleted calcium and impaired root function and (b) earlier onset of the growing season that further taxed the root system. The latter change in phenology appears to act as a trigger at both sites to push trees into nutrient limitation as the demand for Ca increased with the longer growing season, thereby causing the shift in carbon allocation.

Global Biogeochemical Cycles↗

What is the Anthropocene?

Since Paul Crutzen and Eugene Stoermer introduced the word “Anthropocene” in 2000, scientists and nonscientists alike have used the word to highlight the concept that we are now living in a time when the global environment, at some level, is shaped by humankind rather than vice versa. Humans have significantly altered Earth’s land surface, oceans, rivers, atmosphere, flora, and fauna. By its emphasis on the here and now and on what humans have done and can do in the future, the word “Anthropocene” has served as a call to action for environmental sustainability and responsibility [ Crutzen and Stoermer , 2000; Waters et al. , 2014; Ruddiman et al. , 2015]. So far, however, the term “Anthropocene” has not been integrated into the official Geologic Time Scale, which geologists use to divide the past into named blocks based on the rock record. In 2016 or thereabouts, the International Commission on Stratigraphy—the scientific body that maintains the official Geologic Time Scale—will consider a proposal to formalize a definition of this term. It’s a decision that has both semantic and scientific implications and may have legal implications as well.

Eos, Earth and Space Science News↗

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

Overview of spirit microscopic imager results

This paper provides an overview of Mars Exploration Rover Spirit Microscopic Imager (MI) operations and the calibration, processing, and analysis of MI data. The focus of this overview is on the last five Earth years (2005–2010) of Spirit 's mission in Gusev crater, supplementing the previous overview of the first 450 sols of the Spirit MI investigation. Updates to radiometric calibration using in‐flight data and improvements in high‐level processing are summarized. Released data products are described, and a table of MI observations, including target/feature names and associated data sets, is appended. The MI observed natural and disturbed exposures of rocks and soils as well as magnets and other rover hardware. These hand‐lens‐scale observations have provided key constraints on interpretations of the formation and geologic history of features, rocks, and soils examined by Spirit . MI images complement observations by other Spirit instruments, and together show that impact and volcanic processes have dominated the origin and evolution of the rocks in Gusev crater, with aqueous activity indicated by the presence of silica‐rich rocks and sulfate‐rich soils. The textures of some of the silica‐rich rocks are similar to terrestrial hot spring deposits, and observations of subsurface cemented layers indicate recent aqueous mobilization of sulfates in places. Wind action has recently modified soils and abraded many of the rocks imaged by the MI, as observed at other Mars landing sites.

Journal of Geophysical Research E: Planets↗

Environmental controls, emergent scaling, and predictions of greenhouse gas (GHG) fluxes in coastal salt marshes

Coastal salt marshes play an important role in mitigating global warming by removing atmospheric carbon at a high rate. We investigated the environmental controls and emergent scaling of major greenhouse gas (GHG) fluxes such as carbon dioxide (CO 2 ) and methane (CH 4 ) in coastal salt marshes by conducting data analytics and empirical modeling. The underlying hypothesis is that the salt marsh GHG fluxes follow emergent scaling relationships with their environmental drivers, leading to parsimonious predictive models. CO 2 and CH 4 fluxes, photosynthetically active radiation (PAR), air and soil temperatures, well water level, soil moisture, and porewater pH and salinity were measured during May–October 2013 from four marshes in Waquoit Bay and adjacent estuaries, MA, USA. The salt marshes exhibited high CO 2 uptake and low CH 4 emission, which did not significantly vary with the nitrogen loading gradient (5–126 kg · ha −1 · year −1 ) among the salt marshes. Soil temperature was the strongest driver of both fluxes, representing 2 and 4–5 times higher influence than PAR and salinity, respectively. Well water level, soil moisture, and pH did not have a predictive control on the GHG fluxes, although both fluxes were significantly higher during high tides than low tides. The results were leveraged to develop emergent power law‐based parsimonious scaling models to accurately predict the salt marsh GHG fluxes from PAR, soil temperature, and salinity (Nash‐Sutcliffe Efficiency = 0.80–0.91). The scaling models are available as a user‐friendly Excel spreadsheet named Coastal Wetland GHG Model to explore scenarios of GHG fluxes in tidal marshes under a changing climate and environment.

Massachusetts↗

A lacustrine paleoenvironment recorded at Vera Rubin ridge, Gale crater: Overview of the sedimentology and stratigraphy observed by the Mars Science Laboratory Curiosity rover

For ~500 Martian solar days (sols), the Mars Science Laboratory team explored Vera Rubin ridge (VRR), a topographic feature on the northwest slope of Aeolis Mons. Here we review the sedimentary facies and stratigraphy observed during sols 1,800–2,300, covering more than 100 m of stratigraphic thickness. Curiosity's traverse includes two transects across the ridge, which enables investigation of lateral variability over a distance of ~300 m. Three informally named stratigraphic members of the Murray formation are described: Blunts Point, Pettegrove Point, and Jura, with the latter two exposed on VRR. The Blunts Point member, exposed just below the ridge, is characterized by a recessive, fine-grained facies that exhibits extensive planar lamination and is crosscut by abundant curvi-planar veins. The Pettegrove Point member is more resistant, fine-grained, thinly planar laminated, and contains a higher abundance of diagenetic concretions. Conformable above the Pettegrove Point member is the Jura member, which is also fine-grained and parallel stratified, but is marked by a distinct step in topography, which coincides with localized meter-scale inclined strata, a thinly and thickly laminated facies, and occasional crystal molds. All members record low-energy lacustrine deposition, consistent with prior observations of the Murray formation. Uncommon outcrops of low-angle stratification suggest possible subaqueous currents, and steeply inclined beds may be the result of slumping. Collectively, the rocks exposed at VRR provide additional evidence for a long-lived lacustrine environment (in excess of 10 6 years via comparison to terrestrial records of sedimentation), which extends our understanding of the duration of habitable conditions in Gale crater.

Journal of Geophysical Research: Planets↗

Sea‐level rise will drive divergent sediment transport patterns on fore reefs and reef flats, potentially causing erosion on atoll islands

Atoll reef islands primarily consist of unconsolidated sediment, and their ocean‐facing shorelines are maintained by sediment produced and transported across their reefs. Changes in incident waves can alter cross‐shore sediment exchange and, thus, affect the sediment budget and morphology of atoll reef islands. Here we investigate the influence of sea level rise and projected wave climate change on wave characteristics and cross‐shore sediment transport across an atoll reef at Kwajalein Island, Republic of the Marshall Islands. Using a phase‐resolving model, we quantify the influence on sediment transport of quantities not well captured by wave‐averaged models, namely, wave asymmetry and skewness and flow acceleration. Model results suggest that for current reef geometry, sea level, and wave climate, potential bedload transport is directed onshore, decreases from the fore reef to the beach, and is sensitive to the influence of flow acceleration. We find that a projected 12% decrease in annual wave energy by 2100 CE has negligible influence on reef flat hydrodynamics. However, 0.5–2.0 m of sea level rise increases wave heights, skewness, and shear stress on the reef flat and decreases wave skewness and shear stress on the fore reef. These hydrodynamic changes decrease potential sediment inputs onshore from the fore reef where coral production is greatest but increase potential cross‐reef sediment transport from the outer reef flat to the beach. Assuming sediment production on the fore reef remains constant or decreases due to increasing ocean temperatures and acidification, these processes have the potential to decrease net sediment delivery to atoll islands, causing erosion.

Journal of Geophysical Research – Earth Surface↗

Endless forams: >34,000 modern planktonic foraminiferal images for taxonomic training and automated species recognition using convolutional neural networks

Accurate planktonic foraminiferal species identification is central to many paleoceanographic studies, from selecting specific species for geochemical research to elucidating the biotic dynamics of microfossil communities relevant to physical oceanographic processes and interconnected phenomena such as climate change. However, species identification varies among taxonomic schools, few resources exist to train students in the difficult task of discerning amongst closely related species, and the number of taxonomic experts is limited. Here, we take the first steps towards removing these rate-limiting steps by generating the first extensive image library of modern planktonic foraminifera, providing taxonomic training tools and resources, and automating species-level taxonomic identification of planktonic foraminifera via machine learning using convolution neural networks. Taxonomic experts identified 34,640 images of modern planktonic foraminifera to the species level. These images are served as species exemplars through the online portal Endless Forams (endlessforams.org) and a taxonomic training portal hosted on the citizen science platform Zooniverse (zooniverse.org/projects/ahsiang/endless-forams/). A supervised machine learning classifier was then trained with more than 24,000 images of planktonic foraminifera and tested using the remaining ~10,000 images (i.e., the validation set). The best classifier provided the correct species name for an image in the validation set 87.4% of the time. Together, these resources provide a rigorous set of training tools in modern planktonic foraminiferal taxonomy and a means of rapidly generating assemblage data via machine learning in future studies.

Paleoceanography and Paleoclimatology↗

The imminent calving retreat of Taku Glacier

Along the rugged Southeast Alaska coast, 30 kilometers northeast of the state capital Juneau, a tidewater glacier has largely defied global trends by steadily advancing for most of the past century while most glaciers on Earth retreated. This 55-kilometer-long and nearly 1,500-meter-thick tidewater glacier, named Taku Glacier, or T'aaḵú Ḵwáan Sít'i in the language of the Indigenous Tlingit people, has been the focus of continuous scientific study for more than 70 years. Some records even extend back to the mid-18th century. With this long observation record and the glacier’s year-round accessibility and proximity to Juneau and adjacent research facilities, Taku provides an unparalleled locale to study tidewater glaciers and their response to Earth’s rapidly changing climate.

Eos, American Geophysical Union↗

Credit where credit is due

Credit is the currency of science. Scientists are evaluated and promoted in their jobs and professional communities on the basis of their recognized contributions to science. Unlike a financial contribution, a scientific contribution is difficult to measure. Traditionally, credit for scientific contributions has been given through authorship and citations in scientific literature as well as awards and the naming of geographic features, instruments, and methods and other honorifics. However, these practices do not capture the breadth and depth of the contributions by all actors in modern, open science.

Eos Science News↗

Development and application of an Infragravity Wave (InWave) driver to simulate nearshore processes

Infragravity waves are key components of the hydro-sedimentary processes in coastal areas, especially during extreme storms. Accurate modeling of coastal erosion and breaching requires consideration of the effects of infragravity waves. Here, we present InWave, a new infragravity wave driver of the Coupled Ocean-Atmopshere-Waves-Sediment Transport (COAWST) modeling system. InWave computes the spatial and temporal variation of wave energy at the wave group scale and the associated incoming bound infragravity wave. Wave group-varying forces drive free infragravity wave growth and propagation within the hydrodynamic model of the coupled modeling system, which is the Regional Ocean Modeling System (ROMS) in this work. Since ROMS is a three-dimensional model, this coupling allows for the combined formation of undertow currents and infragravity waves. We verified the coupled InWave-ROMS with one idealized test case, one laboratory experiment, and one field experiment. The coupled modeling system correctly reproduced the propagation of gravity wave energy with acceptable numerical dissipation. It also captured the transfer of energy from the gravity band to the infragravity band, and within the different infragravity bands in the surf zone, the measured three-dimensional flow structure, and dune morphological evolution satisfactorily. The idealized case demonstrated that the infragravity wave variance depends on the directional resolution and horizontal grid resolution, which are known challenges with the approach taken here. The addition of InWave to COAWST enables novel investigation of nearshore hydro-sedimentary dynamics driven by infragravity waves using the strengths of the other modeling components, namely the three-dimensional nature of ROMS and the sediment transport routines.

Journal of Advances in Modeling Earth Systems↗

Local energy flux estimates for unstable conditions using variance data in semiarid rangelands

A network of meteorological stations was installed during the Monsoon '90 field campaign in the Walnut Gulch experimental watershed. The study area has a fairly complex surface. The vegetation cover is heterogeneous and sparse, and the terrain is mildly hilly, but dissected by ephemeral channels. Besides measurement of some of the standard weather data such as wind speed, air temperature, and solar radiation, these sites also contained instruments for estimating the local surface energy balance. The approach utilized measurements of net radiation ( R n ), soil heat flux ( G ) and Monin-Obukhov similarity theory applied to first- and second-order turbulent statistics of wind speed and temperature for determining the sensible heat flux ( H ). The latent heat flux ( LE ) was solved as a residual in the surface energy balance equation, namely, LE = −( R n + G + H ). This procedure (VAR-RESID) for estimating the energy fluxes satisfied monetary constraints and the requirement for low maintenance and continued operation through the harsh environmental conditions experienced in semiarid regions. Comparison of energy fluxes using this approach with more traditional eddy correlation techniques showed differences were within 20% under unstable conditions. Similar variability in flux estimates over the study area was present in the eddy correlation data. Hence, estimates of H and LE using the VAR-RESID approach under unstable conditions were considered satisfactory. Also, with second-order statistics of vertical velocity collected at several sites, the local momentum roughness length was estimated. This is an important parameter used in modeling the turbulent transfer of momentum and sensible heat fluxes across the surface-atmosphere interface.

Water Resources Research↗

Quantifying Anderson's fault types

Anderson [1905] explained three basic types of faulting (normal, strike-slip, and reverse) in terms of the shape of the causative stress tensor and its orientation relative to the Earth's surface. Quantitative parameters can be defined which contain information about both shape and orientation [ Célérier , 1995], thereby offering a way to distinguish fault-type domains on plots of regional stress fields and to quantify, for example, the degree of normal-faulting tendencies within strike-slip domains. This paper offers a geometrically motivated generalization of Angelier's [1979, 1984, 1990] shape parameters ϕ and ψ to new quantities named A ϕ and A ψ . In their simple forms, A ϕ varies from 0 to 1 for normal, 1 to 2 for strike-slip, and 2 to 3 for reverse faulting, and Aψ/ranges from 0° to 60°, 60° to 120°, and 120° to 180°, respectively. After scaling, A ϕ and A ψ agree to within 2% (or 1°), a difference of little practical significance, although A ψ has smoother analytical properties. A formulation distinguishing horizontal axes as well as the vertical axis is also possible, yielding an A ϕ ranging from −3 to +3 and A ψ from −180° to +180°. The geometrically motivated derivation in three-dimensional stress space presented here may aid intuition and offers a natural link with traditional ways of plotting yield and failure criteria. Examples are given, based on models of Bird [1996] and Bird and Kong [1994], of the use of Anderson fault parameters A ϕ and A ψ for visualizing tectonic regimes defined by regional stress fields.

Journal of Geophysical Research B: Solid Earth↗

Observations at the Mars Pathfinder site: Do they provide "unequivocal" evidence of catastrophic flooding?

After Mars Pathfinder landed at the mouth of Ares Valus, a large channel that drains into the Chryse Planitia basin, the mission reports unanimously supported the interpretation that the lander site is the locus of catastrophic flooding by noting that all aspects of the scene are consistent with this interpretation. However, alternatives cannot be ruled out by any site observations, as all aspects of the scene are equally consistent with other interpretations of origin, namely, ice and mass-flow processes subsequently modified by wind erosion. The authors discuss alternative explanations for the geologic history of the channel based on a regional view of the circum-Chryse channels from Viking images (our best broad-scale information to date) and the local view from the recent Pathfinder landing site. Mega-indicators of channel origin, the regional geomorphology, geology, and planetary climatic conditions, taken together suggest some combination of flood, mass flow, glacial, and eolian processes. The macro-indicators of channel origin (sedimentologic) are also not indicative of one process of emplacement, either as single criteria or taken cumulatively. Finally, the micro-indicators of channel origin (geochemical and mineralogic composition) do not provide very tight constraints on the deposits' possible origins other than that water was in some way involved.

Journal of Geophysical Research E: Planets↗

Submarine landslides and volcanic features on Kohala and Mauna Kea volcanoes and the Hana Ridge, Hawaii

The deep submarine eastern flanks of Mauna Kea, Kohala, and Haleakala volcanoes were mapped for the first time with a multibeam bathymetric and sidescan sonar system during joint Japan-US cruises aboard the JAMSTEC vessel R/V Yokosuka in 1999. The Pololu slump off northeast Kohala is overlain by a carbonate platform in the shallow region and the deeper areas are incised by downslope oriented channels. It is cut by several faults and slump scars and appears to override an older slide located farther east, here named the Laupahoehoe slump. The structures characteristic of the Laupahoehoe slump are NW-SE oriented scarp-and-bench topographic features analogous to the Hilina slump on the mobile SE flank of Kilauea. Enclosed basins lie at 3000-5000 m, fronted by ridges on their seaward sides. The basins may result from local rotational slumps or from uplift above discontinuous thrust faults. The Laupahoehoe slump appears to be overlapped by shield margins of both Kohala and Mauna Kea and thus may have been derived from an elongate Kohala edifice. A large debris apron continues from the base of the two-slide complex, abutting the distal Hana Ridge (Haleakala east rift zone). The tip of the Hana Ridge displays a curious steepsided arcuate rift tip that resembles the classic amphitheater scarp of a landslide. Numerous flat-topped volcanic domes are distributed along the broad crest of the lower rift zone. Similar cones, though fewer in number, are also present on the Hilo Ridge, which may be the continuation of a Kohala rift zone.

Hawaii↗

Investigation of internal friction in fused quartz, steel, Plexiglass, and Westerly granite from 0.01 to 1.00 Hertz at 10-8 to 10-7 strain amplitude

A detailed evaluation on the method of internal friction measurement by the stress-strain hysteresis loop method from 0.01 to 1 Hz at 10 −8 to 10 −7 strain amplitude and 23.9°C is presented. Significant systematic errors in relative phase measurement can result from convex end surfaces of the sample and stress sensor and from end surface irregularities such as nicks and asperities. Preparation of concave end surfaces polished to optical smoothness having a radius of curvature >3.6×10 4 cm reduces the systematic error in relative phase measurements to <(5.5±2.2)×10 −4 radians. The values of Q E −1 (internal friction under uniaxial compression) determined from the relative phase measurements are | Q E −1 – Q s −1 |< 2.8×10 −3 for the tool steel sample and | Q E −1 – Q s −1 |< 2.2×10 −3 for the Westerly granite sample, where Q s −1 is the internal friction of the fused quartz stress sensor under uniaxial compression. These values are consistent with those inferred from the relative modulus dispersion data also presented in this paper. The polymethyl methacrylate (PMM, trade name Plexiglass) sample shows high values of internal friction ( Q E −1 ≅5×10 −2 ) with strong frequency dependence and with a maximum in Q E −1 at ≅0.4 Hz.

Journal of Geophysical Research Solid Earth↗

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗

Evidence for a new geomagnetic reversal from lava flows in Idaho: Discussion of short polarity reversals in the Brunhes and late Matuyama polarity chrons

K-Ar ages and paleomagnetic data for basalt samples from a new core hole (site E) at the Idaho National Engineering Laboratory (INEL) indicate that the age of the reversed polarity event recorded in Snake River Plain lavas is older than 465±50 ka (1000 years before present) reported previously by Champion et al. (1981). Nine basalt flows, eight with normal polarity and one with reversed polarity, were recognized in the site E core hole. The flows above and below the reversed flow have ages of 491±80 ka and 580±93 ka, respectively. The inclination of the paleomagnetic field direction of the reversed flow at site E agrees with the inclination of reversed flows elsewhere at INEL which have an age of 565±14 ka. These reversed flows were previously thought to be correlative with the Emperor event. We suggest that this polarity event is an older event which we name the Big Lost Reversed Polarity Subchronozone and Subchron. A review of data documenting short reversal records from volcanic and sedimentary rocks shows that there is evidence for eight polarity subchrons in the Brunhes and two besides the Jaramillo in the late Matuyama. These 10 short subchrons begin to indicate the many short events that Cox (1968) hypothesized must exist if polarity interval lengths have a Poisson distribution. These events are true subchrons, not excursions, and may or may not have low associated paleointensities, although low field strengths might explain why the reversal process aborts. The mean sustained polarity interval length since late Matuyama Chron time is 90,000 years. The similarity of this number with the 10 5 -year period of the Earth's orbital eccentricity suggests anew that linkage between geomagnetic, paleoclimatic, and possible underlying Earth orbital parameters should be evaluated.

Journal of Geophysical Research Solid Earth↗