Field Locations Research Articles

Addressing Missing Data in Slope Displacement Monitoring: Comparative Analysis of Advanced Imputation Methods

Slope displacement monitoring is essential for assessing slope stability and preventing catastrophic failures, particularly in geotechnically sensitive areas. However, continuous data collection is often disrupted by environmental factors, sensor malfunctions, and communication issues, leading to missing data that can compromise analysis accuracy and reliability. This study addresses these challenges by evaluating advanced machine learning models—SAITS, ImputeFormer, and BRITS (Bidirectional Recurrent Imputation for Time Series)—for missing data imputation in slope displacement datasets. Sensors installed at two field locations, Yangyang and Omi, provided high-resolution displacement data, with approximately 34,000 data points per sensor. We simulated missing data scenarios at rates of 1%, 3%, 5%, and 10%, reflecting both random and block missing patterns to mimic realistic conditions. The imputation performance of each model was evaluated using Mean Absolute Error, Mean Squared Error, and Root Mean Square Error to assess accuracy and robustness across varying levels of data loss. Results demonstrate that each model has distinct advantages under specific missingness patterns, with the ImputeFormer model showing strong performance in capturing long-term dependencies. These findings underscore the potential of machine learning-based imputation methods to maintain data integrity in slope displacement monitoring, supporting reliable slope stability assessments even in the presence of significant data gaps. This research offers insights into the optimal selection and application of imputation models for enhancing the quality and continuity of geotechnical monitoring data.

Wild bees show local spatial and temporal dynamics in southeastern US blueberry farmscapes.

Wild bee communities are the target of various conservation and ecological restoration programs. Strategic conservation can influence bee communities visiting fields and help mitigate pollinator limitations in fruit production. However, planning compatible conservation strategies and gauging their effectiveness requires understanding how local communities vary across space and time in crops and adjacent semi-natural areas. Here, we assessed the spatiotemporal changes in the composition of wild bee communities in blueberry fields and adjacent forests. In partnership with commercial farms in southeast Georgia, USA, we deployed blue vane traps at the interior and edge of blueberry fields and within adjacent forests of 8 fields, from March to October over 2 yr. We identified 72 wild bee species across 26 genera. The most common were Melissodes communis (Cresson, Hymenoptera: Apidae), Bombus bimaculatus (Cresson, Apidae), Melissodes bimaculatus (Lepeletier), Ceratina floridana (Mitchell, Apidae), Lasioglossum pectorale (Smith, Halictidae), and Lasioglossum nymphale (Smith), which accounted for nearly 60% of the wild bees caught. Bee diversity and richness fluctuated over time, with peaks in all 3 habitat types occurring after the blueberry blooming. Bee abundance in the adjacent forest was relatively constant throughout the season, while in the field interior and edge assemblages, abundance peaked between May and June. We observed dissimilarity in species composition related to month and field location. This difference was explained by fluctuation in the identity and relative abundance of the most abundant species. Together, our study advances foundational knowledge of wild bee community dynamics and species identity in blueberry fields which will help inform and prioritize conservation practices.

Description of Encyonemataceae fam. nov. and Witkowskiaceae fam. nov. (Bacillariophyceae, Cymbellales) based on molecular and morphological analyses

In this study, based on the results of molecular analysis of SSU rRNA and rbcL sequences, we propose the descriptions of two new families in the order Cymbellales. Molecular data demonstrates that diatoms of the genus Encyonema constitute an independent monophyletic clade, which represents a new family described herein – Encyonemataceae fam. nov. Another family introduced in this article, with regards to molecular data, is Witkowskiaceae fam. nov. In our research, it is comprised of three closely related genera – Geissleria, Paraplaconeis and Witkowskia. The results of performed molecular investigation are supported by morphological analysis. Morphological diagnoses of the new families are based on a combination of features: type of valve symmetry, raphe structure, presence and position of stigmata or stigmoids, number and location of apical pore fields and morphology of chloroplast. In addition, we discuss the phylogeny of selected genera of the order Cymbellales. Based on their morphology, 3 and 5 genera are assigned with provisionary position in Encyonemataceae fam. nov. and Witkowskiaceae fam. nov., respectively.

MATHEMATICAL MODEL OF THE VOLUMETRIC STRESS-STRAIN STATE OF GEARS WITH ARBITRARY TOOTH SHAPE

This paper presents an effective method for numerical analysis of the bulk stress-strain state of a gear of a real configuration. Namely, the determination of the bending stress fields in the gear teeth fillets with regard to geometric parameters. The methods of the classical theory of elasticity in a three-dimensional formulation with mixed boundary conditions for a region with a complex boundary surface are used. The solution is given in displacements. The components of the elastic displacement vector at an arbitrary point in the domain are determined, followed by the components of the stress tensor at the same point using Cauchy's differential dependences and the generalized Hooke's law. One of the main difficulties in solving this problem was the problem of representing the entire elastic section of the wheel together with the teeth in the form of an implicit continuous function of a continuous argument, which decisively determines the possibility of applying the methods of elasticity theory and constructing boundary conditions. These difficulties were overcome by using the theory of R-functions to describe the boundary surface of the gear as a whole. The force transmitted by the tooth is introduced in the problem in such a way that the configuration and size of the modeled meshing field are taken into account with the possibility of varying the law of its distribution over the meshing field. Different positions of the instantaneous contact patch in terms of the meshing phase over the entire period of meshing of a pair of teeth are also taken into account. The numerical solution to the problem is based on the Ritz method, where a linearly independent orthonormalized system of Legendre polynomials is used to develop coordinate sequences, and the geometry of the region and boundary conditions are taken into account. The proposed mathematical model of the volumetric stress-strain state can be used both in scientific research and engineering practice at the stage of designing or finishing gears with arbitrary tooth shapes for any real gear meshing system, rational selection of the size and location of the meshing field, and modeling the shape of the instantaneous contact area by the meshing phase for the entire period of meshing of a pair of teeth.

Bronze Age Fields in Suffolk: a Preliminary Survey

Farming developed in Britain during the Neolithic period but across much of England the earliest good archaeological evidence for fields and enclosures in which crops were grown and livestock kept dates from the Middle Bronze Age, c. 1600/1500 bc. While these Bronze Age sub-divided agricultural landscapes are widespread across southern and eastern England, Suffolk and Norfolk were, until recently, essentially a ‘blank’ in their distribution. Over the last 15 years an increasing number of such field systems have been excavated, particularly in Norfolk, and some have started to appear in print. This article adds to this developing picture by briefly describing parts of seven additional Bronze Age – and probable Bronze Age – field systems that have been investigated through recent development-led excavation in south-east Suffolk. Currently published and unpublished evidence from elsewhere in the county is also considered, with the aims of identifying how widespread such land divisions were and establishing the current state of knowledge regarding the location, date, development, layout, and agricultural function of Bronze Age fields in the county. Some of the implications are of wider interest for understanding Bronze Age landscape organisation and land use in lowland England.

Visual Field Asymmetries in Responses to ON and OFF Pathway Biasing Stimuli.

Recent reports suggest the ON and OFF pathways are differentially susceptible to selective vision loss in glaucoma. Thus, perimetric assessment of ON- and OFF-pathway function may serve as a useful diagnostic. However, this necessitates a developed understanding of normal ON/OFF pathway function around the visual field and as a function of input intensity. Here, using electroencephalography, we measured ON- and OFF-pathway biased contrast response functions in the upper and lower visual fields. Using the steady-state visually evoked potential paradigm, we flickered achromatic luminance probes according to a saw-tooth waveform, the fast phase of which biased responses towards the ON or OFF pathways. Neural responses from the upper and lower visual fields were simultaneously measured using frequency tagging - probes in the upper visual field modulated at 3.75Hz, while those in the lower visual field modulated at 3Hz. We find that responses to OFF/decrements are larger than ON/increments, especially in the lower visual field. In the lower visual field, both ON and OFF responses were well described by a sigmoidal non-linearity. In the upper visual field, the ON pathway function was very similar to that of the lower, but the OFF pathway function showed reduced saturation and more cross-subject variability. Overall, this demonstrates that the relationship between the ON and OFF pathways depends on the visual field location and contrast level, potentially reflective of natural scene statistics.

Heading performance outcomes at the FIFA Women's World Cup France 2019™

This study evaluated the performance outcomes of headers during the FIFA Women's World Cup France 2019™. Video analysis was used to code performance outcomes of headers (uncontested and contested) and their descriptors (e.g., playing position, match situation, field location, and the distance the ball travelled). Descriptive statistics, and odds ratios (ORs) (odds of a successful outcome) are reported for headers based on their descriptors. Less than half of all headers resulted in a successful outcome, with headers observed to result in a turnover of possession 53% of the time (uncontested: 51%, contested: 57%). Headed goal conversion rates ranged from 0–38% across countries/teams (mean: 13%), with variability in the frequency of headed shots (range n = 1–22). Headed shots were most efficient from free kicks with 24% of shots resulting in a goal. Odds of a successful heading outcome was lowest from long balls (>20 m) in all areas of the pitch. Uncontested headers had greater odds of a successful heading outcome than contested headers from corner kicks (OR: 2.33, p = 0.004) and free play (OR: 1.30, p = 0.001), but had lower odds of success from goal kicks (OR: 0.62, p = 0.017). Central defensive midfielders (OR: 1.45, p = 0.002) and centrebacks (OR: 1.25, p = <0.001) had significantly greater odds of successful heading outcomes, and strikers (uncontested) (OR: 0.82, p = 0.043) and wingers (contested) (OR: 0.72, p = 0.041) had the lowest. The findings of the current study suggest that heading commonly results in lost possession, particularly from long balls. These findings may help to guide future heading coaching frameworks.

STRUKTUR POPULASI, KOMPOSISI PROKSIMAT DAN MINERAL KEONG GONDANG (Pila scutata) PADA EKOSISTEM PERSAWAHAN DI WILAYAH KABUPATEN GIANYAR

The Pila scutata is a species of Mollusca Phylum that is also known as the kakul in Bali and the gondang snail in Indonesia. The environment of rice fields is one of the snail's habitats. Snails are utilized in religious ceremonies in Bali as a ceremonial component and as an alternate source of protein. The purpose of this study was to ascertain the snail population structure, including density, frequency of presence, and dispersion patterns in Gianyar Regency's rice field ecosystem. The snail sample's mineral and proximate compositions were also examined. Samples were collected from three plots of rice fields with rice plants that were less than a month old at each of the several rice field locations. Sampling using a 1m × 1m square positioned in every rice field corner. Environmental factor measurements are made for every sample. Samples of snails collected from every site were brought to the laboratory for morphological analysis and counting. According to research result, there were 0,83–2,0 ind/m2 of snails, 58,33–91,67% of the frequency of present, and their distribution pattern was regular. There was more than 50% protein. Iron 18,987–38,70; potassium 162,886-216,121, 407,082–567,960 mg/100g; copper and lead not found. Calcium content ranges 90,347–122,97 mg/100g. Keywords: Freshwater snail; native species; nutrition; rice field

Geoinformation monitoring of the condition of rice fields in Giang Province (Vietnam) according to multispectral ERS data and field spectroradiometering

The spatial distribution of areas of territory used for agricultural work is of great importance for the development of measures for managing territories and planning the rational use of land and water resources. As a result of the high development of land and its use for agricultural production, timely assessment of both the condition of soils and the growing season of crops in the fields plays an important role. Since in the conditions of the study area it is possible to choose the timing of sowing, growing and harvesting, spatial data on the location of fields for growing rice can be used to estimate the volume of water consumed for its cultivation and develop a crop rotation model for different volumes of available water based on the level of water reserves of the main crop. source. Geoinformation classification of Earth remote sensing data and the use of spectral indices can be used to monitor the dynamics of rice crop formation under existing conditions. Refinement of the results of geoinformation processing of satellite images is carried out using field standardization methods, including photo standardization, showing the real value of reflected energy during large-scale photography, and spectroradiometry, which makes it possible to determine the characteristics of the reflected energy by these objects in various spectral ranges. The results obtained made it possible to obtain statistical data on the values of the area of plots based on 4844 measurements; the average values of the area of plots were established – 0.447, standard deviation – 0.309, maximum area – 5.84 hectares, minimum 0.02 hectares. The obtained statistical results make it possible to determine on average the estimated number of plots in the territory of An Giang province, which is 442 thousand fields. Thus, when deciphering images of rice fields on satellite images, local geoinformation cartographic layers of such fields are developed, taking into account the actual characteristics and stage of plant vegetation, as well as the stage of harvesting and the stage of tillage.

Agronomic biofortification of genetically biofortified wheat genotypes with zinc, selenium, iodine, and iron under field conditions.

Inherently low concentrations of zinc (Zn), iron (Fe), iodine (I), and selenium (Se) in wheat (Triticum aestivum L.) grains represent a major cause of micronutrient malnutrition (hidden hunger) in human populations. Genetic biofortification represents a highly useful solution to this problem. However, genetic biofortification alone may not achieve desirable concentrations of micronutrients for human nutrition due to several soil- and plant-related factors. This study investigated the response of genetically biofortified high-Zn wheat genotypes to soil-applied Zn and foliarly applied Zn, I, and Se in India and Pakistan. The effect of soil-applied Zn (at the rate of 50 kg ha-1 as ZnSO4·7H2O) and foliar-applied Zn (0.5% ZnSO4·7H2O), I (0.04% KIO3), Se (0.001% Na2SeO4), and a foliar cocktail (F-CT: combination of the above foliar solutions) on the grain concentrations of Zn, I, Se, and Fe of high-Zn wheat genotypes was investigated in field experiments over 2 years. The predominantly grown local wheat cultivars in both countries were also included as check cultivars. Wheat grain yield was not influenced by the micronutrient treatments at all field locations, except one location in Pakistan where F-CT resulted in increased grain yield. Foliar-applied Zn, I, and Se each significantly enhanced the grain concentration of the respective micronutrients. Combined application of these micronutrients was almost equally effective in enhancing grain Zn, I, and Se, but with a slight reduction in grain yield. Foliar-applied Zn, Zn+I, and F-CT also enhanced grain Fe. In India, high-Zn genotypes exhibited a minor grain yield penalty as compared with the local cultivar, whereas in Pakistan, high-Zn wheat genotypes could not produce grain yield higher than the local cultivar. The study demonstrates that there is a synergism between genetic and agronomic biofortification in enrichment of grains with micronutrients. Foliar Zn spray to Zn-biofortified genotypes provided additional increments in grain Zn of more than 15 mg kg-1. Thus, combining agronomic and genetic strategies will raise grain Zn over 50 mg kg-1. A combination of fertilization practice with plant breeding is strongly recommended to maximize accumulation of micronutrients in food crops and to make significant progress toward resolving the hidden hunger problem in human populations.

Развитие методических подходов при моделировании свойств многокомпонентных гранулитов, применяемых в условиях низких температур и других осложняющих факторов

The article addresses the questions of modelling the properties of multicomponent granulites to be used in conditions of low temperatures and other complicating factors. A methodical approach is used which is based on identifying a set of basic and indirect features that characterize the modelled multicomponent simple explosives as a complex object with essential features, distinctive properties, qualitative and quantitative characteristics, and technological indicators. An algorithm to modify the properties of multicomponent granulites was developed in order to solve the problem. This algorithm structurally separates the varying indicators and parameters in the manufacturing and application processes. Specific features of the used raw material components are shown as well as limitations in low-temperature applications. The effects of physical, mechanical, structural and other characteristics of granulite components on the explosive properties are shown, as well as the technical and technological aspects of production. The material mix of recycling materials formed in the processes of mining, concentration and processing of mineral raw materials is provided and the comparative efficiency of granulites is shown with account of cost indicators of the materials used. A method of desulphurisation during explosion of multi-component granulites is proposed which is aimed at neutralization of sulphur dioxide in the blast products. Natural and climatic factors, economic and geographical characteristics of the field location, geocryological conditions and characteristics of the permafrost soils are considered, as well as mining, geological and technical conditions of the development that affect blasting operations in the permafrost soils.

Functional multiple-point simulation

We present a new paradigm, called functional multiple-point simulation, in which multiple-point geostatistical simulation can be performed when functions or curves are observed at each location of a random field. Multiple-point simulation is a non-parametric method used for conditional geostatistical simulation of complex spatial patterns by inferring multiple-point statistics from a training image, rather than from a two-point variogram or covariance model. When the observable at each spatial location is a functional random variable, such multiple-point simulation must take into account not only the spatial correlation among locations but also the similarity of functions or curves observed at each location. The data events to be compared in this case are now functional, in the sense that they consist of spatial arrangements of functions. Consequently, we propose four distances, inspired by the functional data analysis literature, for measuring similarities between functional data events and use these to extend the direct sampling method to perform multiple-function geostatistical simulation with functional fields. We coin the new method Functional Direct Sampling and carry out extensive qualitative and quantitative performance comparison between the four proposed distances using simulation techniques on two well-known applications of multiple-point simulation: simulating copies of a functional random field and gap-filling of locations in a functional random field. We apply the proposed method to a gap-filling task of simulated wind profiles spatial functions over the Arabian Peninsula.

Optimizing liquid‐applied iron sulfate rate and application interval for dollar spot suppression on golf course fairways

AbstractIron sulfate applications reduce dollar spot severity on golf course turfgrass, but further optimizing the use of iron sulfate may increase efficacy and encourage more widespread adoption by golf course superintendents. Three iron sulfate rates (9, 18, 36 kg ha−1), two application intervals (7 and 14 days), three water volumes (306, 612, and 1223 L ha−1), and combinations with fertilizer and fungicide products were evaluated at field locations in Madison, WI, Bloomington, MN, and Blacksburg, VA, in 2018 and 2019. Iron sulfate rate decreased dollar spot at all locations in a dose‐dependent manner, and 7‐day applications were more efficacious than 14‐day applications. Water volume had no impact on disease severity at any location and combining iron sulfate with fertilizer had minimal impact as well. These results demonstrate improved iron sulfate strategies for dollar spot management and can be implemented by golf course superintendents in areas where fungicide resistance and fungicide regulations limit the use of traditional fungicides.

Retinotopic coding organizes the interaction between internally and externally oriented brain networks.

The human brain seamlessly integrates internally generated thoughts with incoming sensory information, yet the networks supporting internal (default network, DN) and external (dorsal attention network, dATN) processing are traditionally viewed as antagonistic. This raises a crucial question: how does the brain integrate information between these seemingly opposed systems? Here, using precision neuroimaging methods, we show that these internal/external networks are not as dissociated as traditionally thought. Using densely-sampled 7T fMRI data, we defined individualized whole-brain networks from participants at rest and calculated the retinotopic preferences of individual voxels within these networks during an visual mapping task. We show that while the overall network activity between the DN and dATN is independent at rest, considering a latent retinotopic code reveals a complex, voxel-scale interaction stratified by visual responsiveness. Specifically, the interaction between the DN and dATN at rest is structured at the voxel-level by each voxel's retinotopic preferences, such that the spontaneous activity of voxels preferring similar visual field locations is more anti-correlated than that of voxels preferring different visual field locations. Further, this retinotopic scaffold integrates with the domain-specific preferences of subregions within these networks, enabling efficient, parallel processing of retinotopic and domain-specific information. Thus, DN and dATN are not independent at rest: voxel-scale interaction between these networks preserves and encodes information in both positive and negative BOLD responses, even in the absence of visual input or task demands. These findings suggest that retinotopic coding may serve as a fundamental organizing principle for brain-wide communication, providing a new framework for understanding how the brain balances and integrates internal cognition with external perception.

Response of n-Heptane and Dodecane Swirl-Stabilized Spray Flames to Transverse Forcing Characterized By Flame Describing Functions Based on Downstream Acoustic Pressure or Velocity Signals

Abstract Predicting thermoacoustic instabilities in annular combustors requires knowledge of the impact of acoustic oscillations on heat release rate oscillations. Flame Describing Functions (FDF) measured at the burner exit using acoustic forcing are key elements of thermoacoustic instability analyses. FDFs based on acoustic pressure measurements, FP' or on axial velocity measurements FU', are compared here. This study is done on the TACC-Spray bench, an original linear array of spray flames stabilized by a strong swirling flow, representing an unfolded sector of a self-unstable annular combustor. Acoustic forcing of a standing transverse chamber mode is applied downstream of the injectors. Experiments are conducted with liquid n-heptane or dodecane, with the flames placed at a pressure or an intensity antinode of the transverse mode. FP' does not depend on the measurement location for acoustically-compact flames provided that this location remains in the flame vicinity. FU' can lead to significant discrepancies, as swirling flows present strong velocity gradients, which can be minimized by carefully choosing the measurement location. The injector admittance linking the two FDFs is shown to be quasi-independent of the forcing amplitude here. Consequently, both FDFs show a similar dependence on the forcing amplitude. FP' indicates constructive combustion-acoustics interference whatever the flame location in the acoustic field and the fuel, consistent with self-sustained instabilities observed in the annular combustor. An analysis using the Rayleigh criterion corroborates the results derived from the FDFs. So, FP' appears as a powerful and practical tool for characterizing combustion dynamics.

Optical trapping using quasi-bound states in the continuum of photonic crystal slab.

This paper explores efficient and stable optical trapping using quasi-bound states in the continuum (quasi-BICs) in photonic crystal slabs. By breaking inversion symmetry by transforming square holes into isosceles trapezoidal holes, we create quasi-BICs with finite quality factors. Calculations show that optical forces at enhanced electric field locations effectively trap particles, with significant potential wells at these sites. Multipole expansion analysis indicates that particles at sharp corners suppress radiation and enhance trapping stability. Moreover, an increased number of trapped particles causes a redshift in resonance frequency and strengthens optical forces, especially at sharp corners. These insights are crucial for designing photonic crystal slabs for practical optical trapping applications, where positioning more particles at sharp corners improves quasi-BIC mode and trapping efficiency.

Positive and Negative Retinotopic Codes in the Human Hippocampus.

The hippocampus is thought to coordinate sensory-mnemonic information streams in the brain, representing both the apex of the visual processing hierarchy and the central hub of mnemonic processing. Yet, the mechanisms underlying sensory-mnemonic interactions in the hippocampus are poorly understood. Recent work in cortex suggests that a retinotopic code - typically thought to be exclusive to visual areas - may help organize internal and external information at the cortical apex via opponent interactions. Here, we leverage high-resolution 7T functional MRI to test whether a bivalent retinotopic code structures activity within the human hippocampus and mediates hippocampal-cortical interactions. In seven densely-sampled individuals, we defined the retinotopic preferences of individual voxels within the hippocampus and cortex during a visual mapping task, as well as their functional connectivity during independent runs of resting-state fixation. Our findings reveal a robust retinotopic code in the hippocampus, characterized by stable population receptive fields (pRFs) with consistent preferred visual field locations across experimental runs. Notably, this retinotopic code is comprised of roughly equal proportions of positive and negative pRFs, aligning with the hypothesized role of negative pRFs in mnemonic processing. Finally, the signed amplitude of hippocampal pRFs predicts functional connectivity between retinotopic hippocampal and cortical voxels. Taken together, these results suggest that retinotopic coding may scaffold internal mnemonic and external sensory information processing within the hippocampus, and across hippocampal-cortical interactions.

Player and match characteristics associated with head acceleration events in elite-level men’s and women’s rugby union matches

ObjectiveTo examine the likelihood of head acceleration events (HAEs) as a function of previously identified risk factors: match time, player status (starter or substitute) and pitch location in elite-level men’s...

Frequency-of-seeing curves (psychometric functions) for perimetric stimuli in age-related macular degeneration.

Frequency-of-seeing (FoS) curves (psychometric functions) for perimetric stimuli have been widely used in computer simulations of new visual field test procedures. FoS curves for age-related macular degeneration (AMD) are not available in the literature and are needed for the development of improved microperimetry test procedures, which are of particular interest for use as clinical trial endpoints. Data were refitted from a previous study to generate FoS curves for 20 participants with AMD, each tested at nine locations within the central 10°. Stimulus parameters, background luminance and dB scale were matched to the MAIA-2 microperimeter, and stimuli were presented in a method of constant stimuli to build up FoS curves over multiple runs. FoS curves were fitted with a modified cumulative Gaussian function. The relationship between sensitivity and slope of fitted FoS curves was modelled by robust linear regression, producing models both with and without an eccentricity parameter. FoS curves were satisfactorily fitted to data from 174 visual field locations in 20 participants (age 65-83 years, 11 female). Each curve was made up of a median of 243 (range 177-297) stimulus presentations over a median of 12 (range 9-32) levels. Median sensitivity was 25.5 dB (range 3.8-31.4 dB). The median slope (SD of fitted function) was 1.6 dB (range 0.5-8.5 dB). As in previous studies of other conditions, the slope of fitted FoS curves increased as sensitivity decreased (p < 0.001). FoS are provided for participants with AMD, as well as models of the relationship between sensitivity and slope. These fitted models and data may be useful for computer simulation studies of microperimetry procedures. Full details of the fitted curves are provided as supporting information.

Extreme thunderstorm ground enhancements registered on Aragats in 2023

In 2023, a series of intense thunderstorm ground enhancements (TGEs) were recorded on Mount Aragats in Armenia, with seven events exceeding the fair-weather cosmic ray flux by more than 75%. This study comprehensively analyzes these TGEs, investigating the atmospheric conditions and electric fields contributing to their occurrence. Key insights include discovering relationships between TGE electron content and the location of atmospheric electric fields, recovering electron and gamma-ray energy spectra, and the impact of the cloud charge structure. The findings offer a deeper understanding of TGEs’ role in atmospheric physics and its synergy with high-energy astrophysics. Published by the American Physical Society 2024