The Eastern United States has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader eastern United States. To produce reliable images of the subsurface, we simultaneously invert smoothed P-wave receiver functions, Rayleigh-wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear wave speed. Using surface-wave observations (3-250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear-wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment. A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western Mississippi Embayment. A comparison of crust thickness derived from our model against four recent published models shows first-order consistency. A relatively small upper mantle low-speed region correlates with a published P-waves analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

The agriculture sector consumes more than two-thirds of world’s limited freshwater resources. However, only a small part of the water (less than 5%) that is taken up by roots is used for plant growth, while the rest (above 95%) is lost due to transpiration through the stomatal apertures. Therefore, reducing the transpiration of agricultural plants will contribute to the preservation of precious water resources. However, reducing the transpiration rate artificially is difficult because most plants react delicately and negatively, resulting in water-stressed conditions that often cause different physiological disorders. The present study investigated the transpiration light response in tomato plants (Solanum lycopersicum) grown under LED lights and assessed different irradiation techniques’ ability to reduce transpiration and maintained proper plant growth in a controlled environment. Tomato plants were grown in three enclosed hydroponic units under blue (460 nm) and red (630 nm) LEDs inside an air-conditioned glasshouse. The test plants and multiple replicates were grown five consecutive times, and the irradiation intensity (photosynthetic photon flux density (PPFD)), irradiation pattern (simultaneous/alternate irradiation for red/blue LEDs) and LED combination (number/ratio of red/blue LEDs) were changed each time. The plants’ physiological parameters (transpiration, stomatal conductance, stem-diameter, stem height, and number of leaves) and daily transpiration rates were recorded periodically and analyzed. The results show that a typical photoperiod of 12 hours with simultaneous irradiation of red/blue LEDs produced balanced physiological growth for plants in general. However, when normalized against water use efficiency (transpiration), an alternate irradiation pattern (6 hours: blue LED on/off repeatedly for 15-minute intervals + 6 hours: red LED on/off repeatedly for 15-minute intervals) was the most suitable for tomato cultivation in controlled environments.

Extreme weather conditions are associated with a variety of water quality issues that can pose harm to humans and aquatic ecosystems. Under dry extremes, contaminants become more concentrated in streams with a greater potential for harmful algal blooms, while wet extremes can cause flooding and broadcast pollution. Developing appropriate interventions to improve water quality in a changing climate requires a better understanding of how extremes affect watershed processes, and which places are most vulnerable. We developed a Soil and Water Assessment Tool model of the Cape Fear River Basin (CFRB) in North Carolina, USA, representing contemporary land use, point and non-point sources, and weather conditions from 1979 to 2019. The CFRB is a large and complex river basin undergoing urbanization and agricultural intensification, with a history of extreme droughts and floods, making it an excellent case study. To identify intervention priorities, we developed a Water Quality Risk Index (WQRI) using the load average and load variability across normal conditions, dry extremes, and wet extremes. We found that the landscape generated the majority of contaminants, including 90.1% of sediment, 85.4% of total nitrogen, and 52.6% of total phosphorus at the City of Wilmington’s drinking water intake. Approximately 16% of the watershed contributed most of the pollutants across conditions—these represent high priority locations for interventions. The WQRI approach considering risks to water quality across different weather conditions can help identify locations where interventions are more likely to improve water quality under climate change.

Geomagnetically induced currents (GICs) are driven by the geoelectric field induced by fluctuations of Earth’s magnetic field. Drivers of intense GICs are often associated with large impulsive events such as coronal mass ejections. To a lesser extent fluctuations from regular oscillations of the geomagnetic field, or geomagnetic pulsations, have also been identified as possible drivers of GICs. In this work we show that these low-frequency pulsations are directly observed in measured GIC data from power networks. Due to the low-pass nature of GICs, Pc5 and lower frequency pulsations drive significant GICs for an extended duration at mid-latitudes. Longer period Ps6-type disturbances apparently not typical of mid-latitudes are seen with GIC amplitudes comparable to the peak GIC at storm sudden commencement. The quasi-ac nature of the sustained pulsation driving affects the power system response and cannot be properly modelled using only dc models. A further consideration is that the often used dB/dt GIC proxy is biased to the sampling rate of the geomagnetic field measurements used. The dB/dt metric does not adequately characterise GIC activity at frequencies in the low ULF range and a frequency weighted proxy akin to geoelectric field should be used instead.

Large earthquakes rapidly denude hillslopes by triggering thousands of coseismic landslides. The sediment produced by these landslides is initially quickly mobilised from the landscape by an interconnected cascade of processes. This cascade can dramatically but briefly enhance local erosion rates. Hillslope and channel processes, such as landsliding and debris flows, interact to influence the total mass, calibre, and rate of sediment transport through catchments. Calculating the sediment budget of an earthquake lends insight into the nature of these interactions. Using satellite imagery derived landslide inventories, channel surveys and a literature review combined with a Monte Carlo simulation approach we present a constrained sediment budget of the first decade after the 2008 Mw7.9 Wenchuan earthquake. With this sediment budget we demonstrate that debris flows are dominant process for delivering sediment into channels and that large volumes of sediment remain in the landscape. In our study area over 88% (469 Mega tonnes) of the coseismically generated sediment remains on the hillslopes in 2018. Of the 12% of the sediment that was mobilised, 69% (40.7 14 Mt) was mobilised by debris flows. Despite the large proportion of sediment remaining on the hillslope, the frequency of debris flows declined significantly over our observation period. The reduction in debris-flow frequency is not correlated to reductions in the frequency of triggering storms, suggesting changes in the mechanical properties of hillslope sediment may drive this observation. The stabilisation of coseismically generated sediment greatly extends its residence time and may influence catchment sediment yields for centuries or millennia.

The upper-level jet stream impacts surface-level trace gas variability, yet the cause of this relationship remains unclear. We investigate the mechanism(s) responsible for the relationship using idealized tracers with different source regions within a chemical transport model. All tracers’ daily variabilities are correlated with the meridional position of the jet in the mid-latitudes, but tracers emitted south (north) of the jet increase (decrease) in the mid-latitudes when the jet is shifted poleward. The jet stream regulates the near-surface meridional wind, and this coupling together with the meridional tracer gradient robustly predicts where the jet stream and tracers are in and out of phase. Our study elucidates a major driver of trace gas variability and links it to the location of the jet stream and emissions. These results are useful for understanding changes in trace gas variability if the jet stream’s position or major emission source regions change in the future.

Unless humanity achieves United Nations Sustainable Development Goals (SDGs) by 2030 and restores the relatively stable climate of pre-industrial CO2 levels (as early as 2110), species extinctions, starvation, drought/floods, and violence will exacerbate mass migrations. This paper presents conceptual designs and techno-economic analyses to calculate sustainable limits for growing high-protein seafood and macroalgae-for-biofuel. We review the availability of wet solid waste and outline the mass balance of carbon and plant nutrients passing through a hydrothermal liquefaction process. The paper reviews the availability of dry solid waste and dry biomass for bioenergy with CO2 capture and storage (BECCS) while generating Allam Cycle electricity. Sufficient wet-waste biomass supports quickly building hydrothermal liquefaction facilities. Macroalgae-for-biofuel technology can be developed and straightforwardly implemented on SDG-achieving high protein seafood infrastructure.). The analyses indicate a potential for (1) 0.5 billion tonnes/yr of seafood; (2) 20 million barrels/day of biofuel from solid waste; (3) more biocrude oil from macroalgae than current fossil oil; and (4) sequestration of 28 to 38 billion tonnes/yr of bio-CO2. Carbon dioxide removal (CDR) costs are between 25−33% of those for BECCS with pre-2019 technology or the projected cost of air-capture CDR.

A flowslide overriding liquefied substrate can vastly enhance its disaster after failure initiation, due to rapid velocity and long-runout distance during landslides mobilized into flows. It is crucial to provide improved understanding to the mechanism of these catastrophic flowslides for hazard mitigation and risk assessment. This study focuses on the Saleshan landslide of Gansu in China, which is a typically catastrophic flowslide overrode a liquefied sand substrate. Geomorphologic and topographic maps along with analysis of seismic signals confirm its dynamic features and mobilized behaviors. ERT surveying detected abundant groundwater in the landslide, which is fundamental to its rapid long-runout distance. Particle size distributions and triaxial shear behaviors affirmed more readily liquefied behavior of superficial loess and underlying alluvial sand than red soil sandwiched them. We also examined the liquefaction susceptibility of the alluvial sand under loading impact at undrained and drained conditions. The alluvial sand is readily liquefied in the undrained condition while it is difficult at drained condition due to rapid water pore pressure dissipation. The results showed that the landslide experienced a sudden transformation from slide on the steep slope where it originated to flow on a nearly flat terrace with abundant groundwater that it overrode. This transformation can be attributed to the liquefied alluvial sand substrate enhancing the whole landslide body mobility. Along with recent, similar findings from landslides worldwide, substrate liquefaction may present a widespread, significant increase in landslide hazard and consequent mobility and our study reveals conditions necessary for this phenomenon to occur.

The theory for recording of thermally blocked remanences predicts a quasi-linear relationship between low fields like the Earth’s in which rocks cool and acquire a magnetization. This serves as the foundation for estimating ancient magnetic field strengths. Addressing long-standing questions concerning Earth’s magnetic field require a global paleointensity dataset, but recovering the ancient field strength is complicated because the theory only pertains to uniformly magnetized particles. A key requirement of a paleointensity experiment is that a magnetization blocked at a given temperature should be unblocked by zero-field reheating to the same temperature. However, failure of this requirement occurs frequently and the causes and consequences of failure are poorly understood. Recent experiments demonstrate that the remanence in many samples typical of those used in paleointensity experiments is unstable, and exhibits an ”aging’ effect in which the unblocking temperature spectrum changes over only a few years resulting in non-ideal experimental behavior. While a fresh remanenence may conform to the requirement of equality of blocking and unblocking temperatures, aged remanences may not. Blocking temperature spectra can be unstable (fragile), which precludes reproduction of the conditions under which the original magnetization was acquired. This limits our ability to acquire accurate and precise ancient magnetic field strength estimates because differences between known and estimated fields can be significant (up to 10 μT) for individual specimens, with a low field bias. Fragility of unblocking temperature spectra appears to be related to grain size and may be related to features observed in first-order reversal curves.

The seismic activity of the middle strand of the North Anatolian Fault (MNAF), Northwestern Turkey, is debated because of its quiescence during the instrumental period, in contrast to a significant historical activity documented by several chronicles over the last two millennia. Here, we focus on Lake Iznik, bordered by the MNAF, to get a new insight into its long-term seismicity and its tectonic setting. The study of lacustrine sediment cores reveals fourteen earthquake-induced turbidites since their ages correspond to seismic events during the past two millennia. Bathymetry and high-resolution seismic reflection data allow describing two hitherto unknown subaquatic active fault structures (the South Boyalica and Iznik faults) that belong to the MNAF system. Sediment cores sampled on both sides of the Iznik Fault document an event deposit and a sedimentary unit vertically offset of ~ 40 cm interpreted as the last rupture during the 1065 CE destructive earthquake. Older events are supposed on this fault more than a thousand years ago. Further studies will help to estimate the horizontal coseismic offset of this oblique-slip fault and the calendar of older ruptures. The current seismic gap of thousand years on this segment greatly increases the seismic hazard in this region and must be considered in the seismic risk assessment of the NAF system.

Changes in tropical cyclone intensity prior to landfall represent a significant risk to human life and coastal infrastructure. Such changes can be influenced by shelf water temperatures through their role in mediating heat exchange between the ocean and atmosphere. However, the evolution of shelf sea surface temperature during a storm is dependent on the initial thermal conditions of the water column, information that is often unavailable. Here, observational data from multiple monitoring stations and satellite sensors were used to identify the sequence of events that led to the development of storm-favorable thermal conditions in the Mississippi Bight prior to the transit of Hurricane Sally (2020), a storm that rapidly intensified over the shelf. The annual peak in depth-average temperature of >29°C that occurred prior to the arrival of Hurricane Sally was the result of two distinct warming periods caused by a cascade of weather events. The event sequence transitioned the system from below average to above average thermal conditions over a 25-d period. The transition was initiated with the passage of Hurricane Marco (2020), which mixed the upper water column, transferring heat downward and minimizing the cold bottom water reserved over the shelf. The subsequent reheating of the upper ocean by a positive surface heat flux, followed by a period of downwelling winds, effectively elevated shelf-wide thermal conditions for the subsequent storm. The climatological coupling of warm sea surface temperature and downwelling winds suggest regions with such characteristics are at an elevated risk for storm intensification over the shelf.

On September 5-7, 2018, a series of tremors were reported in Nigeria's capital city, Abuja. These events followed a growing list of tremors felt in the stable intra-plate region, where earthquakes are not expected. Here, we review available seismological, geological, and geodetic data that may shed light on the origin of these tremors. First, we investigate the seismic records for parent location of the orphan tremors using a technique suitable when a single-seismic station is available such as the Western Africa region, which has a sparse seismic network. We find no evidence of the reported tremors within the seismic record of Western Africa. Next, we consider the possibility of a local amplification of earthquakes from regional tectonics, reactivation of local basement fractures by far-field tectonic stresses, landward continuation of oceanic fracture zones, or induced earthquakes triggered by groundwater extraction. Our assessments pose important implications for understanding Western Africa's intraplate seismicity and its potential connection to tectonic inheritance, active regional tectonics, and anthropogenic stress perturbation.

In closed-canopy forests, the availability of photosynthetically active light has been a focal point of research, emphasizing the role of light as a resource in limiting carbon assimilation and individual tree growth. However, light shapes the functioning of forest ecosystems through multiple mechanisms. Here, using a series of studies from a network of tree diversity experiments, we explore the multifaceted ways in which light---in terms of both quantity and quality---shapes productivity in mixed-species forests. Spectral reflectance from remote sensing of forest canopies is being increasingly used to detect how tree diversity influences productivity. We demonstrate that airborne imaging spectroscopy captures functionally important differences among canopies related to their structure, chemistry, and underlying biological interactions. Ground-based analyses can show in detail how photosynthetically active light is partitioned among species in mixed-species communities. We show that greater interception of light and greater efficiency of light use, generated by inter- and intra-specific differences, combine to enhance productivity in mixed-species forests. Light may shape forest function not only as a resource but also as a stressor and cue. Plants can perceive light at various wavelengths, use this information to assess their neighborhoods, and subsequently adjust their physiology and allocation. We characterize how light quality---from the ultraviolet to shortwave infrared---varies among and within canopies of differing diversity. We explore how these diversity-light quality relationships arise and connect across levels of biological organization from leaf-level trait expression to forest function. Together these studies lend insight into light-mediated mechanisms that drive relationships between biodiversity and productivity in forest ecosystems---insights that are crucial to predict how biodiversity change will affect future forest function.

Mixing parameters can be inaccurate in ocean data assimilation systems, even if there is close agreement between observations and mixing parameters in the same modeling system when data are not assimilated. To address this, we investigate whether there are additional observations that can be assimilated by ocean modeling systems to improve their representation of mixing parameters and thereby gain knowledge of the global ocean’s mixing parameters. Observationally-derived diapycnal diffusivities–using a strain-based parameterization of finescale hydrographic structure–are included in the Estimating the Circulation & Climate of the Ocean (ECCO) framework and the GEOS-5 coupled Earth system model to test if adding observational diffusivities can reduce model biases. We find that adjusting ECCO-estimated and GEOS-5-calculated diapycnal diffusivity profiles toward profiles derived from Argo floats using the finescale parameterization improves agreement with independent diapycnal diffusivity profiles inferred from microstructure data. Additionally, for the GEOS-5 hindcast, agreement with observed mixed layer depths and temperature/salinity/stratification (i.e., hydrographic) fields improves. Dynamic adjustments arise when we make this substitution in GEOS-5, causing the model’s hydrographic changes. Adjoint model-based sensitivity analyses suggest that the assimilation of dissolved oxygen concentrations in future ECCO assimilation efforts would improve estimates of the diapycnal diffusivity field. Observationally-derived products for horizontal mixing need to be validated before conclusions can be drawn about them through similar analyses.

Western Jezero crater of Mars exposes smooth- and rough-surfaced plains that bound a steep-sided fan-shaped plateau (i.e., western Jezero delta) in the west. The two plain terrains have a gradational contact and their end-member occurrences can be defined by the absence (smooth-surfaced) and presence (rough-surfaced) of littered surficial boulders and depressions. Early researchers interpret the plain terrains as morphological expression of volcanic flow and aeolian, fluvial, or lacustrine deposition, but these hypotheses have not been tested rigorously via detailed mapping. Here we show results of HiRISE-based mapping that reveals four landform units hosted by both smooth- and rough-surfaced terrains: (1) polygonal networks of double-ridged grooves that are 2-4 m wide, up to ~1700 m long, and spaced > ~100 m, (2) reticulate networks of grooves spaced mostly < ~10 m, (3) undulating surfaces hosting variously shaped depressions, and (4) NE-trending boulder-bearing ridges (~10 m wide and 70-300 m long) locally displaying NE-pointing streamlined shapes. Our mapping shows that the NE-trending boulder ridges formed first followed sequentially by the formation of depressions, double-ridged grooves, and reticulate grooves. None of the above observations can be explained as a whole by volcanic emplacement and/or earlier proposed aeolian, fluvial, and lacustrine depositional processes. Guided by Earth analogues, we interpret the streamlined boulder ridges as subglacial flutes, depressions including most circular pits with raised rims as kettle holes and/or thermokarsts, polygonal networks of double-ridged grooves as crevasse-filled ice-pressed moraine ridges and the polygonal pattern was inherited from the fractured glacier, and reticulate grooves as desiccation cracks during the final drying of glacial deposits at the end of a martian ice age.

This paper describes an implementation of the Combined Hybrid-Parallel Prediction (CHyPP) approach of Wikner et al. (2020) on a low-resolution atmospheric global circulation model (AGCM). The CHyPP approach combines a physics-based numerical model of a dynamical system (e.g., the atmosphere) with a computationally efficient type of machine learning (ML) called reservoir computing (RC) to construct a hybrid model. This hybrid atmospheric model produces more accurate forecasts of most atmospheric state variables than the host AGCM for the first 7-8 forecast days, and for even longer times for the temperature and humidity near the earth’s surface. It also produces more accurate forecasts than a model based only on ML, or a model that combines linear regression, rather than ML, with the AGCM. The potential of the approach for climate research is demonstrated by a 10-year long hybrid model simulation of the atmospheric general circulation, which shows that the hybrid model can simulate the general circulation with substantially smaller systematic errors and more realistic variability than the host AGCM.

Gyres are prominent surface structures in the global ocean circulation that often interact with the sea floor in a complex manner. Diagnostic methods, such as the depth-integrated vorticity budget, are needed to assess exactly how such model circulations interact with the bathymetry. Terms in the vorticity budget can be integrated over the area enclosed by streamlines to identify forces that spin gyres up and down. In this article we diagnose the depth-integrated vorticity budgets of both idealized gyres and the Weddell Gyre in a realistic global model. It is shown that spurious forces play a significant role in the dynamics of all gyres presented and that they are a direct consequence of the Arakawa C-grid discretization and the z-coordinate representation of the sea floor. The spurious forces include a numerical beta effect and interactions with the sea floor which originate from the discrete Coriolis force when calculated with the following schemes: the energy conserving scheme (ENE); the enstrophy conserving scheme (ENS); and the energy and enstrophy conserving scheme (EEN). Previous studies have shown that bottom pressure torques provide the main interaction between the depth-integrated flow and the sea floor. Bottom pressure torques are significant, but spurious interactions with bottom topography are similar in size. Possible methods for reducing the identified spurious topographic forces are discussed. Spurious topographic forces can be alleviated by using either a B-grid in the horizontal plane or a terrain-following vertical coordinate.

The ‘eddying’ ocean, recognized for several decades, has been the focus of much observational and theoretical research. We here describe a generalization for the analysis of eddy energy, based on the use of ensembles, that addresses two key related issues: the definition of an ‘eddy’ and the general computation of energy spectra. An ensemble identifies eddies as the unpredictable component of the flow, and permits the scale decomposition of their energy in inhomogeneous and non-stationary settings. We present two distinct, but equally valid, spectral estimates: one is similar to classical Fourier spectra, the other reminiscent of classical EOF analysis. Both satisfy Parseval’s equality and thus can be interpreted as length-scale dependent energy decompositions. The issue of ‘tapering’ or ‘windowing’ of the data, used in traditional approaches, is also discussed. We apply the analyses to a mesoscale ‘resolving’ (1/12$^\circ$) ensemble of the separated North Atlantic Gulf Stream. Our results reveal highly anisotropic spectra in the Gulf Stream and zones of both agreement and disagreement with theoretically expected spectral shapes. In general, we find spectral slopes that fall off faster than the steepest slope expected from quasi-geostrophic theory.