High-resolution simulations by the Regional Ocean Modeling System (ROMS) were used to investigate the dispersal of the San Francisco Bay (SFB) plume over the northern-central California continental shelf during the period of 2011 to 2012. The modeled bulk dynamics of surface currents and state variables showed many similarities to corresponding observations. After entering the Pacific Ocean through the Golden Gate, the SFB plume is dispersed across the shelf via three pathways: (i) along the southern coast towards Monterey Bay, (ii) along the northern coast towards Point Arena, and (iii) an offshore pathway restricted within the shelf break. On the two-year mean timescale, the along-shore zone of impact of the northward-dispersed plume is about 1.5 times longer than that of the southern branch. Due to the opposite surface Ekman transports induced by the northerly or southerly winds, the southern plume branch occupies a broader cross-shore extent, roughly twice as wide as the northern branch which extends roughly two times deeper due to coastal downwelling. Besides these mean characteristics, the SFB plume dispersal also shows considerable temporal variability in response to various forcings, with wind and surface-current forcing most strongly related to the dispersing direction. Applying constituent-oriented age theory, we determine that it can be as long as 50 days since the SFB plume was last in contact with SFB before being flushed away from the Gulf of the Farallones. This study sheds light on the transport and fate of SFB plume and its impact zone with implications for California’s marine ecosystems.

The Great Slave Lake shear zone (GSLsz) is a type example for deeply eroded continental transform boundaries located in the Northwest Territories, Canada. Formed during the oblique convergence of the Archean Rae and Slave cratons, the GSLsz has accommodated up to 700 km of dextral shear. Here we present the results of in situ U-Pb apatite and titanite geochronology from 11 samples that were collected across the strike of the shear zone. Both geochronometers record a near-continuous history of ductile shear during crustal cooling and exhumation that spans ca. 1920–1740 Ma. By integrating the geochronological data with structural and metamorphic observations across the structure, we propose a tectonic model for the shear zone that consists of three stages. The first stage (ca. 1920–1880 Ma) is characterized by strain accommodation along two coeval fault strands. During the second stage (ca. 1880–1800 Ma), ductile shear ceases along the northernmost fault strand and the locus of strain migrates southwards towards the hinterland of the Rae cratonic margin. In the third stage (ca. 1800–1740 Ma), ductile strain localizes back along the southern of the two original fault strands, after which the present-day surface level of the shear zone transitions to brittle shear. Our results highlight both the significance of the lateral migration of the zone of active deformation in major crustal shear zones as well as the localization of strain along existing crustal structures.

Basin-scale quasi-geostrophic gyres are common features of large lakes subject to Coriolis force. Cyclonic gyres are often characterized by dome-shaped thermoclines that form due to pelagic upwelling which takes place in their center. At present, dynamics of pelagic upwelling in the Surface Mixed Layer (SML) of oceans and lakes are poorly documented. A unique combination of high-resolution 3D numerical modeling, satellite imagery and field observations allowed confirming for the first time in a lake, the existence of intense pelagic upwelling in the center of cyclonic gyres under strong shallow (summer) and weak deep (winter) stratified conditions/thermocline. Field observations in Lake Geneva revealed that surprisingly intense upwelling from the thermocline to the SML and even to the lake surface occurred as chimney-like structures of cold water within the SML, as confirmed by Advanced Very High-Resolution Radiometer data. Results of a calibrated 3D numerical model suggest that the classical Ekman pumping mechanism cannot explain such pelagic upwelling. Analysis of the contribution of various terms in the vertically-averaged momentum equation showed that the nonlinear (advective) term dominates, resulting in heterogeneous divergent flows within cyclonic gyres. The combination of nonlinear heterogeneous divergent flow and 3D ageostrophic strain caused by gyre distortion is responsible for the chimney-like upwelling in the SML. The potential impact of such pelagic upwelling on long-term observations at a measurement station in the center of Lake Geneva suggests that caution should be exercised when relying on limited (in space and/or time) profile measurements for monitoring and quantifying processes in large lakes.

The solar radio flux at 10.7 cm, known as F10.7, is a critical operational space weather index. However, without a clear backup, any interruption to the service can result in substantial errors in model outputs. In this paper we show the impact of one such outage in March 2022 and present a number of alternative solutions for any future outages. The approach resulting in the smallest reconstruction error of F10.7 uses the solar radio flux observations at alternative wavelengths (the best giving a percentage error of 3.1%). Alternatively, use of Sunspot Number, a regular, robust alternative observation, results in a mean percentage error of 8.2% and is also a reliable fallback solution. Additionally, analysis of the error on the use of the conversion between the 12-month rolling sunspot number (R12) and its conversion to F10.7 as used by the IRI is included.

This paper reports on the standing whistler waves upstream of Mercury’s quasi-perpendicular bow shock. Using MESSENGER’s magnetometer data, 36 wave events were identified during interplanetary coronal mass ejections (ICMEs). These elliptic or circular polarized waves were characterized by: (1) a constant phase with respect to the shock, (2) propagation along the normal direction to the shock surface, and (3) rapid damping over a few wave periods. We inferred the speed of Mercury’s bow shock as ~31 km/s and a shock width of 1.76 ion inertial length. These events were observed in 20% of the MESSENGER orbits during ICMEs. We conclude that standing whistler wave generations at Mercury are generic to ICME impacts and the low Alfvén Mach number (MA) collisionless shock, and are not affected by the absolute dimensions of its bow shock. Our results further support the theory that these waves are generated by the current in the shock.

The use of heavy machinery is increasing in agricultural industries in particular cotton farming systems, which induces an increased risk of soil compaction and yield reduction. Hence, there is a need for a technical solution to use available tools to measure projected soil compaction due to farm machinery traffics. The aim of this work was to compare the effects of static and dynamic loads on soil compaction. In this study, three vertosols (common soil for cotton production) were selected to examine soil compaction under a range of static and dynamic loads using uniaxial compression equipment and a modified proctor test, respectively. In general, soils behaved similarly under static and dynamic loads with no significant difference between bulk density values for all moisture contents with a high index of agreement (d=0.96, RMSE= 0.056). The results further indicate better agreement between soil compaction for static and dynamic loads Uniaxial compression test (static loads) produced higher compaction compared to the modified proctor test (dynamic loads) in particular at moisture contents lower than the plastic limit condition. The variation in soil compaction for static and loads was often evident for loads ≥600 kPa, with the highest soil compaction induced under loads ≥1200 kPa. The findings of this study confirm the suitability of a modified proctor method to assess soil compaction as an alternative tool under a range of moisture contents and machinery loads for vertosols.

Climate change has altered the frequency, intensity, and timing of mean and extreme precipitation events. Extreme precipitation has caused tremendous socio-economic losses and displays strong regional variability. Although many previous studies have addressed daily extreme precipitation, hourly extreme rainfall still needs to be thoroughly investigated. In this study, we investigated the trends, spatio-temporal variability, and long-term variations in mean and extreme precipitation over South Korea using daily and hourly observational data. During the past 50 years (1973–2022), there has been a notable escalation in maximum hourly precipitation, although the boreal summer mean precipitation has increased only marginally. Regionally, an increase in mean and extreme rainfall occurred in the northern part of the central region. Moreover, increased intensity and frequency of extreme precipitation have contributed more to the total summer precipitation in recent years. Our findings provide scientific insights into the progression of extreme summer precipitation events in South Korea.

We investigate the effect of extremely rough bathymetry on energy dissipation and mixing in a coastal region characterized by small-scale seafloor features penetrating a strongly-stratified density interface of comparable vertical scale. Our data from the non-tidal Baltic Sea include shear microstructure measurements and observations from a broadband echosounder, here used to resolve the extreme variability and intermittency of stratified turbulence in the vicinity of obstacles. Scale analysis and acoustic imaging of small-scale turbulent motions suggest that the underlying mixing mechanisms are related to topographic wake eddies and, to a smaller extent, to breaking internal waves near the bathymetric features. Vertical diffusivities exceed those at a nearby reference station with smooth bathymetry by up to two orders of magnitude. Our study emphasizes the importance of rough small-scale (< 1 km) bathymetric features for energy dissipation and vertical turbulent transport in coastal areas shaped by e.g., glacial, tectonic, or volcanic processes.

This research targets disentangling shallow causes of anthropogenically-induced subsidence in a reclaimed and urbanized coastal plain. The study area is around the city of Almere, in the South Flevoland polder, the Netherlands, which is among the countries’ fastest subsiding areas. The procedure consists of integrating synthetic Interferometric Synthetic Aperture Radar (InSAR) data with high-resolution phreatic groundwater and lithoclass models, and a database containing construction details. The two main parts of the workflow are isolation of the InSAR points of structures without a pile foundation and a data assimilation procedure by Ensemble Smoothing with Multiple Data Assimilation. The shrinkage of surficial clay beds by phreatic groundwater level lowering is identified to be the main cause of shallow subsidence in the area, with an average contribution of 6 mm/year. The history-matched physics-based model predicts that one meter drop in phreatic groundwater level now translates into 10 millimeter of subsidence in the next five years. Also, this study showed that a groundwater deficiency due to severe dry periods should be considered as an accelerator of subsidence in both the short- and long-term planning. To ensure a robust network to estimate future subsidence, we advise on a consistent monitoring strategy of the phreatic groundwater level.

The observed retreat and anticipated further decline in Arctic sea ice hold strong climate, environmental, and societal implications. In predicting climate evolution, ensembles of coupled climate models have demonstrated appreciable accuracy in simulating sea ice area and volume trends throughout the historical period. However, individual climate models still show significant differences in simulating the sea ice thickness distribution. To better understand individual model performance in sea ice simulation, nine climate models previously identified to provide plausible sea ice decline and global temperature change were evaluated in comparison with Arctic satellite and reanalysis derived sea ice thickness data, sea ice extent records, and atmospheric reanalysis data of surface wind and air temperature. Assessment found that the simulated spatial distribution of historical sea ice thickness varies greatly between models and that several key limitations persist among models. Primarily, most models do not capture the thickest regimes of multi-year ice present in the Wandel and Lincoln Seas; those that do, often possess erroneous positive bias in other regions such as the Laptev Sea or along the Eurasian Arctic Shelf. From analysis, no model could be identified as performing best overall in simulating historic sea ice, as model bias varies regionally and seasonally. Nonetheless, the bias maps and statistical measures derived from this analysis should enhance understanding of the limitations of each climate model. This research is motivated in-part to inform future usage of coupled climate model projection for regional modeling efforts and enhance climate change preparedness and resilience in the Arctic.

Parameterizations for bottom shear stress are required to predict sediment resuspension from field observations and within numerical models that do not resolve flow within the viscous sublayer. This study assessed three observation-based bottom shear stress (τb) parameterizations, including (1) the sum of surface wave stress and mean current (quadratic) stress (τb= τw +τc); (2) the log-law (τb= τL); and (3) the turbulent kinetic energy (τb= τTKE); using two years of observations from a large shallow lake. For this system, the parameterization τb= τw +τc was sufficient to qualitatively predict resuspension, since bottom currents and surface wave orbitals were the two major processes found to resuspend bottom sediments. However, the τL and τTKE parameterizations also captured the development of a nepheloid layer within the hypolimnion associated with high-frequency internal waves. Reynolds-averaged Navier-Stokes (RANS) equation models parameterize τb as the summation of modeled current-induced bottom stress (τc,m) and modelled surface wave-induced bottom stress (τw,m). The performance of different parameterizations for τc,m and τw,m in RANS models was assessed against the observations. The optimal parameterizations yielded root-mean-square errors of 0.031 and 0.025 Pa, respectively, when τc,m, and τw,m were set using a constant canonical drag coefficient. A RANS-based τL parameterization was developed; however, the grid-averaged modelled dissipation did not always match local observations, leading to O(10) errors in prediction of bottom stress. Turbulence-based parameterizations should be further developed for application to flows with mean shear-free boundary turbulence.

The Indian plate underthrusting the Himalaya is considered to be segmented along the collision belt arc and seismic images of the Indian mantle lithosphere (IML) suggest along-arc variations in the angle of underthrusting and its northern limit beneath Tibet. The pre-existing transverse tectonic structures of the Indian plate mapped in the Ganga foreland basin have been related to these segmentation boundaries. These segmentations imply changes in mechanical properties of adjoining blocks which should manifest in the form of spatial variations in topography build-up. We have analysed a geomorphic index, normalized channel steepness (ksn), along the Himalayan arc using the ALOS elevation dataset to test whether there is any correlation between the and these segmentation boundaries. Our results bring out spatial variability in the along the arc. Based on these results, the arc can be segmented into five blocks, similar to the ones delineated based on correlation between the width of the Ganga foreland basin and the disposition of major Himalayan thrusts from the foothills. Thus, the ksn can be used as a proxy to demarcate different tectonic blocks along the Himalayan arc. Further, we have found a good correlation between the basin width and the northern limit of the IML for all block except the Uttarakhand block. We infer that transverse crustal heterogeneities in this block due to the continuation of different litho-units of the Aravalli-Delhi Fold Belt could be a plausible cause for this anti-correlation.

This study summarizes paleomagnetic secular variation (PSV) in five published Holocene records from Eastern North America. We have developed 100-year increment time series for the declinations and inclinations for all sites and compared their directional variability. We see evidence of ten correlatable features in both inclination and declination. We focus on the clockwise or counter-clockwise motion of paleomagnetic directions (termed circularity) in these PSV records. We have first calculated the incremental rate and direction of motion (clockwise or counter-clockwise) for each record over the last 4000-8000 years. We have separately looked for discernable looping in individual records. We estimate the loop sizes, durations, and circularity direction. We see the same pattern of circularity in both measurement techniques. There are seven intervals of oscillating circularity and looping in all five sites. Both techniques suggest a distinctive oscillating, teeter-totter like, behavior to PSV circularity that must be due to the pattern of fluid flow in the outer core. This teeter-totter behavior is unbalanced with more time spent in clockwise motion than in counter-clockwise motion. We think the teeter-totter oscillation may be due to torsional oscillation in the outer core fluid flow. The loops have a distribution of sizes and durations with smaller loops being shorter in duration and bigger loops having longer durations. All five PSV records show 5 intervals of ~102 yr significant acceleration in circularity rate and PSV rate combined with change in circularity direction. These features are broadly analogous to historic geomagnetic jerks.

Agricultural management strategies are crucial in regulating the soil-atmosphere interaction. The crop landscape is influenced by farmers through different field practices, and further impacts the variations of soil temperature, soil moisture, and field microclimate. To examine how different management strategies affect the soil properties and the aforementioned interaction, two observation systems were installed in an organic-certified (ORG) tea field and a conventional (CONV) tea field in northern Taiwan. The results show that the variation of canopy temperature was more significant in CONV while the difference in soil diurnal temperature range was minor. However, the daily loss rate of soil water content in ORG was two times faster than that in CONV (0.93% d−1 vs. 0.46% d−1). These findings suggest that the appropriate management strategies could assist farmers in adapting to environmental fluctuations and provide quantitative references for assessing soil characteristics under different agricultural applications and climatic conditions.

Laterite is a red weathering crust developed with various rocks and Quaternary loose sediments as its parent material in the tropics and subtropics regions of the world, and it is also the most widely distributed Quaternary earthy accumulation in China. Since the 1930s, most researchers have believed that the fluvial reticulated laterite in southern China was influenced by the warm and humid climate of the Middle Pleistocene. In recent years, the remains of Paleolithic human activities are often found in the reticulated laterite of southern China. However, the study of laterite chronology is sporadic or there is no critical chronological analysis, which causes uncertainty in the identification and discussion of the ages of reticulated laterite and Paleolithic sites in South China. In this study, a paleolithic site found in fluvial reticulated laterite in South China was systematically tested by quartz optical luminescence dating and geomorphic process analysis. The results show that, (1) The T3 terrace, an archive of hominin activity in the study area, primarily formed between 56 and 11 ka. (2) Reticulated laterite cannot be used simply to determine the ages of the Paleolithic sites found in this stratum, and typical reticulated laterite cannot be used as a marker for climatic stratigraphy and chronostratigraphy. The fluvial reticulated laterite in the southern tropics, under suitable hydrothermal conditions, can form within tens of thousands of years or even within 10 ka. (3) Human activities can also lead to an inversion in the age of reticulated red soil.

Assessing space weather modeling capability is a key element in improving existing models and developing new ones. In order to track improvement of the models and investigate impacts of forcing, from the lower atmosphere below and from the magnetosphere above, on the performance of ionosphere-thermosphere models, we expand our previous assessment for 2013 March storm event [Shim et al., 2018]. In this study, we evaluate new simulations from upgraded models (Coupled Thermosphere Ionosphere Plasmasphere Electrodynamics (CTIPe) model version 4.1 and Global Ionosphere Thermosphere Model (GITM) version 21.11) and from NCAR Whole Atmosphere Community Climate Model with thermosphere and ionosphere extension (WACCM-X) version 2.2 including 8 simulations in the previous study. A simulation of NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model version 2 (TIE-GCM 2) is also included for comparison with WACCM-X. TEC and foF2 changes from quiet-time background are considered to evaluate the model performance on the storm impacts. For evaluation, we employ 4 skill scores: Correlation coefficient (CC), root-mean square error (RMSE), ratio of the modeled to observed maximum percentage changes (Yield), and timing error(TE). It is found that the models tend to underestimate the storm-time enhancements of foF2 (F2-layer critical frequency) and TEC (Total Electron Content) and to predict foF2 and/or TEC better in the North America but worse in the Southern Hemisphere. The ensemble simulation for TEC is comparable to results from a data assimilation model (Utah State University-Global Assimilation of Ionospheric Measurement (USU-GAIM)) with differences in skill score less than 3% and 6% for CC and RMSE, respectively.

A document by Srishti Singh. Click on the document to view its contents.

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture and area, stiffness proxy, under dynamic stressing. Experiments are conducted on rough, tensile-fractured Westerly granite specimen subject to triaxial stresses. Fracture permeability is measured from steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1~MPa at 1~Hz frequency. During dynamic stressing we transmit acoustic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor the evolution of fracture interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10, 12.5, 15, 17.5, and 20~MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and they map contact area at stresses from 9 to 21~MPa. The resulting ‘true’ area of contact measurements made for the entire fracture surface and within the calculated PZT sensor footprints, numerical modeling of Fresnel zone. We compare the elastodynamic response of the the fracture using the stress-induced changes ultrasonic wave velocities for a range of transmitter-receiver pairs to image spatial variations in contact properties, which is informed by fracture contact area measurements. These measurements of the nonlinear elasticity are related to the fluid-flow, permeability, in response to dynamic stressing and similar comparisons are made for the slow-dynamics, recovery, of the fracture interface following the stress perturbations.