Feed aggregator

Variance estimations in the presence of intermittent interference and their applications to incoherent scatter radar signal processing
Qihou Zhou, Yanlin Li, and Yun Gong
Atmos. Meas. Tech., 17, 4197–4209, https://doi.org/10.5194/amt-17-4197-2024, 2024
We discuss several robust estimators to compute the variance of a normally distributed random variable to deal with interference. Compared to rank-based estimators, the methods based on the geometric mean are more accurate and are computationally more efficient. We apply three robust estimators to incoherent scatter power and velocity processing, along with the traditional sample mean estimator. The best estimator is a hybrid estimator that combines the sample mean and a robust estimator.

Surface equilibrium vapor pressure of organic nanoparticles measured from the dynamic-aerosol-size electrical mobility spectrometer
Ella Häkkinen, Huan Yang, Runlong Cai, and Juha Kangasluoma
Atmos. Meas. Tech., 17, 4211–4225, https://doi.org/10.5194/amt-17-4211-2024, 2024
We report measurements of evaporation kinetics and surface equilibrium vapor pressures for various laboratory-generated organic nanoparticles using the dynamic-aerosol-size electrical mobility spectrometer (DEMS), a recent advancement in aerosol process characterization. Our findings align well with literature values, demonstrating DEMS's effectiveness. We suggest future improvements to DEMS and anticipate its potential for probing aerosol-related kinetic processes with unknown mechanisms.

Abstract

The accurate assessment of hydraulic transmissivity in rock fractures filled with particles is not only a scientific challenge but also a critical need for various industrial applications. However, the intricate dynamics of particle erosion and pore clogging that govern transmissivity evolution remain largely unexplored. In this study, we experimentally examine the fluid-driven particle migration behavior in filled fractures and its consequent impact on fracture transmissivity under various hydraulic gradients, normal stresses, and fracture apertures. We find that escalating hydraulic gradients not only intensify particle erosion through amplified fluid drag forces and hydro-mechanical coupling effects but also lead to an increase in the size of migrating particles, thereby augmenting pore clogging. The dynamics of erosion and clogging define four distinct migration phases within the filled fractures. Variations in normal stress and initial fracture aperture significantly alter the particle arrangement and the soil structure stability within the fractures, thereby modulating the progress of particle migration in response to hydraulic gradients. The pattern of particle migration in filled fractures dictates the development of the internal pore structure and normal deformation, ultimately affecting fracture transmissivity. We propose an empirical expression to encapsulate the comprehensive evolution of fracture transmissivity across different particle migration patterns. Our research advances the understanding of fluid-driven particle migration within filled fractures and provides a practical tool for the precise determination of hydraulic properties of fractured rocks amidst complex geological settings.

Abstract

Methane clumped isotope signatures of abiogenesis may be diagnostic of the origin of methane on Earth and other planetary bodies. We performed synthesis of abiogenic methane in hydrothermal conditions between 130 and 300°C and determined δ13C, δD, Δ13CH3D, and Δ12CH2D2. The experiments were performed by heating water in the presence of Fe0 powder and CO. The reduction of water on metallic iron led to the formation of H2. CO was reacted with both H2 and H2O, generating both CH4 and CO2. Methane δ13C values are isotopically depleted by ∼25‰ relative to the CO starting material. This is consistent with carbon isotopic equilibrium between methane, carbon monoxide and carbon dioxide in our experiments. In contrast, D/H ratios are inconsistent with equilibrium isotopic fractionation, as illustrated by δD values of methane fractionated by ∼500‰ relative to starting H2O. This suggests that under our experimental conditions, hydrogen additions to carbon may be governed by kinetics. Δ13CH3D values track experimental temperature, with values between +1.5‰ and +5.0‰ for most samples. In contrast, Δ12CH2D2 values are displaced from equilibrium. We find exclusively negative Δ12CH2D2 values, showing deficits down to 40‰ relative to thermodynamic equilibrium. We interpret the data as evidence for distinct, kinetically induced D/H pools contributing to methane assembly, that is, a combinatorial effect. The cumulative D/H fractionations associated with CO hydrogenation explain the direction and magnitude of Δ12CH2D2 values during abiotic methane formation. We suggest that near equilibrium Δ13CH3D with negative Δ12CH2D2 signatures will help identify methane formed abiotically in nature.

Abstract

In this study we apply electron tomography to characterize 3D dislocation microstructures in two quartz mylonite specimens from the Moine and Main Central Thrusts, both of which accommodated displacements by dislocation creep in the presence of water. Both specimens show dislocation activity with dislocation densities of the order of 3–4 × 1012 m−2 and evidence of recovery from the presence of subgrain boundaries. 〈a〉 slip occurs predominantly on pyramidal and prismatic planes (〈a〉 basal glide is not active). [c] Glide is not significant. On the other hand, we observe a very high level of activation of 〈c + a〉 glide on the 101‾0 $\left\{10\overline{1}0\right\}$, 101‾1 $\left\{10\overline{1}1\right\}$, 112‾n $\left\{11\overline{2}n\right\}$ (n = 1,2) and even 213‾1 $\left\{21\overline{3}1\right\}$ planes. Approximately 60% of all dislocations show evidence of climb with a predominance of mixed climb, a deformation mechanism characterized by dislocations moving in a plane intermediate between the glide and the climb planes. This atypical mode of deformation demonstrates comparable glide and climb efficiency under natural deformation conditions. It promotes dislocation glide in planes not expected for the quartz structure, probably by inhibiting lattice friction. Our quantitative characterization of the microstructure enables us to assess the strain that dislocations can generate. We show that glide systems indicated by the observed dislocations are sufficient to accommodate any arbitrary 3D strain by themselves. Although historically dislocation glide has been regarded as being primarily responsible for producing strain, activation of climb can also directly contribute to the finite strain. On the basis of this characterization, we propose a numerical modeling approach for attempting to characterize the local stress state that gave rise to the observed microstructure.

Abstract

Sea ice cools Earth by reducing its absorbed solar energy. We combine radiative transfer modeling with satellite-derived surface albedo, sea ice, and cloud distributions to quantify the top-of-atmosphere sea ice radiative effect (SIRE). Averaged over 1980–2023, Arctic and Antarctic SIREs range from −0.64 to −0.86 W m−2 and −0.85 to −0.98 W m−2, respectively, with different cloud data sets and assumptions of climatological versus annually-varying clouds. SIRE trends, however, are relatively insensitive to these assumptions. Arctic SIRE has weakened quasi-linearly at a rate of 0.04–0.05 W m−2 decade−1, implying a 21%–27% reduction in the reflective power of Arctic sea ice since 1980. Antarctic sea ice exhibited a regime change in 2016, resulting in 2016–2023 Antarctic and global SIRE being 0.08–0.12 and 0.22–0.27 W m−2 weaker, respectively, relative to 1980–1988. Global sea ice has therefore lost 13%–15% of its planetary cooling effect since the early/mid 1980s, and the implied global sea ice albedo feedback is 0.24–0.38 W m−2 K−1.

Author(s): Baihui Chen (陈百慧), Chengjie Zhan (湛承杰), Zhenhua Chai (柴振华), and Baochang Shi (施保昌)

In this work, we develop a phase-field-based lattice Boltzmann (LB) method for a two-scalar model of the two-phase flows with interfacial mass or heat transfer. Through the Chapman-Enskog analysis, we show that the present LB method can correctly recover the governing equations for phase field, flow…

[Phys. Rev. E 110, 015307] Published Wed Jul 17, 2024

Abstract

Updated summaries of the August 1972 and March 1989 space weather events have been constructed. The features of these two events are compared to the Carrington 1859 event and a few other major space weather events. It is concluded that solar active regions release energy in a variety of forms (X-rays, EUV photons, visible light, coronal mass ejection (CME) plasmas and fields) and they in turn can produce other energetic effects (solar energetic particles (SEPs), magnetic storms) in a variety of ways. It is clear that there is no strong one-to-one relationship between these various energy sinks. The energy is often distributed differently from one space weather event to the next. Concerning SEPs accelerated at interplanetary CME (ICME) shocks, it is concluded that the Fermi mechanism associated with quasi-parallel shocks is relatively weak and that the gradient drift mechanism (electric fields) at quasi-perpendicular shocks will produce harder spectra and higher fluxes. If the 4 August 1972 intrinsic magnetic cloud condition (southward interplanetary magnetic field instead of northward) and the interplanetary Sun to 1 au conditions were different, a 4 August 1972 magnetic storm and magnetospheric dawn-to-dusk electric fields substantially larger than the Carrington event would have occurred. Under these special interplanetary conditions, a Miyake et al. (2012), https://doi.org/10.1038/nature11123-like extreme SEP event may have been formed. The long duration complex 1989 storm was probably greater than the Carrington storm in the sense that the total ring current particle energy was larger.

Abstract

Electron density enhancements in the ionospheric D-region due to the precipitation of high-energy electrons (>30 keV) have been measured as increases in cosmic radio noise absorption (CNA) using ground-based riometers. CNA has been studied since the 1960s. However, there have been few studies of the spatiotemporal development of CNA at multi-point ground stations distributed in longitude at subauroral latitudes, where plasma particles with a wide energy range are intermingled. In this study, we analyzed the longitudinal development of CNA steep increases using simultaneous riometer observations at six stations at subauroral latitudes in Canada, Alaska, Russia, and Iceland over 3 years from 2017 to 2020. The results revealed that the occurrence rate of steep increases in CNA was highest at midnight at 22-08 magnetic local time (MLT), and lowest near dusk at 17–21 MLT. We also showed statistically that the CNA steep increases expanded eastward on the dawn side and westward on the dusk side. The CNA expansion velocity was slightly faster than the results of previous studies in the auroral zone. Correlation and superposed epoch analyses of CNA with solar wind and geomagnetic parameters revealed that CNA intensity was dependent on the Interplanetary Magnetic Field Bz, Interplanetary Electric Field Ey, SYM-H index, and SME index. These results indicate that the CNA at subauroral latitudes is closely related to solar wind and geomagnetic activities, and its propagation characteristics correspond to the dynamics of high energy electrons in the inner magnetosphere.

Abstract

Electromagnetic Ion Cyclotron (EMIC) wave scattering has been proved to be responsible for the fast loss of both radiation belt (RB) electrons and ring current (RC) protons. However, its role in the concurrent dropout of these two co-located populations remains to be quantified. In this work, we study the effect of EMIC wave scattering on both populations during the 27 February 2014 storm by employing the global physics-based RAM-SCB model. Throughout this storm event, MeV RB electrons and 100s keV RC protons experienced simultaneous dropout following the occurrence of intense EMIC waves. By implementing data-driven initial and boundary conditions, we perform simulations for both populations through the interplay with EMIC waves and compare them against Van Allen Probes observations. The results indicate that by including EMIC wave scattering loss, especially by the He-band EMIC waves, the model aligns closely with data for both populations. Additionally, we investigate the simulated pitch angle distributions (PADs) for both populations. Including EMIC wave scattering in our model predicts a 90° peaked PAD for electrons with stronger losses at lower pitch angles, while protons exhibit an isotropic PAD with enhanced losses at pitch angles above 40°. Furthermore, our model predicts considerable precipitation of both particle populations, predominantly confined to the afternoon to midnight sector (12 hr < MLT < 24 hr) during the storm's main phase, corresponding closely with the presence of EMIC waves.

Abstract

We comprehensively analyzed geomagnetic perturbations using ground magnetic records from over 400 stations spanning four solar cycles, from 1976 to 2023. We assess the perturbations in the three magnetic components separately. Our study covers low, middle, and high magnetic latitudes in the northern magnetic hemisphere, with the primary objective of quantifying extreme values and evaluating their variability on magnetic latitude, local time, and solar cycle phases “minimum, ascending, maximum, and declining.” Our findings reveal spatial patterns to be less discernible as perturbations intensify, with distinct responses at middle and high latitudes. The extreme values, defined as percentiles 0 and 100, were observed to be localized and randomly distributed in local time, especially in the east magnetic component. Additionally, we observed dusk-dawn asymmetries in the magnitude of perturbations related to the auroral electrojets, indicating complex interactions between the magnetosphere and ionosphere. Furthermore, the results reveal a preference for the most significant extreme values to occur in the declining phase of the solar cycle. These insights deepen our understanding of geomagnetic perturbations and their variability, contributing to space weather forecasting and mitigation strategies.

Abstract

Instrumental data set have revealed several summer precipitation patterns in eastern China, being summarized as “tripole,” “dipole” and “coast” modes. The former two have been found to persist at different time scales, leaving the latter unclear in geological records. Here we present 1300-year hydroclimate records in a tropical maar lake in the southern coast of China using archeal lipid GDGTs, which can reflect lower water redox conditions largely regulated by lake water depth. The down-core GDGTs reveal a relatively drier condition during the medieval climate anomaly compared to the Little Ice Age, in-phase with other records in southeast coast of China but opposite to the inland region, hence demonstrating a persistent “coast” mode in eastern China. The thermal state of equatorial Pacific is suggested to play an important role in shaping the “coast” mode by modulating the location and strength of the western Pacific subtropical high and tropical typhoons.

Abstract

We investigate similarities and differences between natural and simulated slow earthquakes using nonlinear dynamical system tools. We use spatio-temporal slip potency rate data derived from Global Navigation Satellite System (GNSS) position time series in the Cascadia subduction zone and numerical simulations intended to reproduce their pulse-like behavior and scaling laws. We provide metrics to evaluate the accuracy of simulations in mimicking slow earthquake dynamics. We investigate the influence of spatio-temporal coarsening as well as observational noise. Despite the use of many degrees of freedom, numerical simulations display a surprisingly low average dimension, akin to natural slow earthquakes. Instantaneous dynamical indices can reach large values (>10) instead, and differences persist between numerical simulations and natural observations. We propose to use the suggested metrics as an additional tool to narrow the divergence between slow earthquake observations and dynamical simulations.

Abstract

Paleomagnetic analyses have suggested that the lunar magnetic field underwent a significant change from 4.25 to 3.19 Ga, indicating the rapid transition of the lunar dynamo mechanism. We used the van der Pauw (vdP) method to measure the electrical resistivity of Fe-S-P alloys under conditions relevant to the lunar core and estimated the thermal conductivity of the Fe-S-P lunar core. These values were incorporated into thermal and dynamo models to investigate the evolution of the lunar core. Our model indicates that the inner core began to grow as early as 4.35 Ga, the solidification regime switched at 3.50 Ga, and the thermal dynamo ceased between 3.78 and 3.51 Ga. The cessation of the dynamo could be due to a low buoyancy flux and insufficient entropy dissipation. Thermal and compositional dynamos cannot sustain the ancient strength of the Moon's magnetic field, and require other energy sources.

Abstract

A comprehensive statistical study is conducted on O+ and H+ outflows obtained from the TEAMS/FAST data during the 23rd solar cycle (1996–2007). The study investigates interhemispheric asymmetry in ionospheric outflows during local summer, winter, and equinox seasons. Data are classified into two distinct periods: the pre-storm and geomagnetic storm phases. Numerous statistical asymmetries are identified. The findings indicate that the dayside cusp consistently demonstrates more outflow rates of O+ and H+ in the northern hemisphere than southern hemisphere during geomagnetic storms in all seasons as well as during the pre-storm period in the summer season with the exception of H+ during summer storms. Conversely, the nightside O+ and H+ outflow rates are higher in the southern hemisphere during pre-storm and storm periods in the summer season. Additionally, the dawnside and duskside outflow rates of O+ and H+ are predominantly stronger in the southern hemisphere.

Abstract

The application of global navigation satellite system in vehicle navigation imposes high requirements on its safety, which makes it extremely important to detect deceptive interference (DI). The integration of prior road position information is usually necessary for the vehicle navigation, but it requires extensive learning and data collection, resulting in high implementation costs. We focus on efficient and practical research on DI detection by utilizing the most traditional wheel speed sensor auxiliary information. The classical pseudorange consistency detection method avoids the use of absolute position information, and overcomes the limitation of clustering algorithms caused by the number of spoofing interferences. Hypothesis test statistics are constructed by using the differencing method in the east and north directions, and the probability distributions of detection values under normal and spoofing conditions are theoretically derived. A comprehensive deception detection scheme is introduced to enhance the detection probability, catering to both single and multiple satellite deception situations, with a central focus on the contribution change factor of positioning satellites. Experimental results demonstrate that within a time scale of several tens of seconds, effective identification and correction of DI involving more than one code chip bias can be achieved.

SummaryWe consider application of full waveform inversion (FWI) to radio-frequency electromagnetic (EM) data. Radio-frequency imaging (RIM) is a cross-borehole technique to image electromagnetic subsurface properties from measurements of transmitted radio-frequency waves. It is used in coal seam imaging, ore exploration and various engineering and civil engineering applications. RIM operates at frequencies from 50 kHz to several tens of MHz. It differs from other geophysical EM methods, because the frequency band includes the transition between the wave propagation and diffusion regimes. RIM data are acquired in two-dimensional cross-hole sections in a reciprocal manner. Traditionally, radio-frequency data are inverted by straight ray tomography because it is inexpensive and easy to implement. It is argued that due to attenuation, the sensitivity of the transmitted electric field is the strongest within the first Fresnel zone of the ray connecting the transmitter and receiver. While straight ray tomography is a simple to implement and fast method, the non-linearity in the relationship between model parameters and data is often strong enough to warrant non-linear inversion techniques. FWI is an iterative high resolution technique, in which the physical properties are updated to minimize the misfit between the measured and modelled wavefields. Full waveform techniques have been used and extensively studied for the inversion of seismic data, and more recently, they have been applied to the inversion of GPR data. Non-linear inversion methods for RIM data are less advanced. Their use has been hindered by the high cost of full wave modelling and the high conductivity contrasts of many RIM targets, and, to some extent, by the limitations of the measuring instruments. We present the first application of this methodology to simultaneous conductivity and permittivity inversion of RIM data. We implement the inversion in the frequency domain in two dimensions using L-BFGS optimization. We analyze the sensitivity of the data to the model parameters and the parameter trade-off and validate the proposed methodology on a synthetic example with moderate conductivity variations and localized highly conductive targets. We then apply the FWI methodology to a field data set from Sudbury, Canada. For the field data set, we determine the most appropriate preprocessing steps that take into account specific peculiarities of RIM: the insufficient prior information about the subsurface and the limitations of the measuring equipment. We show that FWI is applicable under the conditions of RIM and is robust to imperfect prior knowledge: we obtain satisfactory model recoveries starting from homogeneous initial models in all of our examples. Just as other methods, FWI underestimates large conductivity contrasts due to the loss of sensitivity of the transmitted electric field to the conductivity variations as the conductivity increases above a certain level. The permittivity inside high conductors can not be recovered, however, recovering permittivity variations in the resistive zones helps obtain better focused conductivity images with fewer artifacts. Overall, FWI produces cleaner, less noisy and higher resolution reconstructions than the methods currently used in practice.

Abstract

Saturn’s mid-sized icy moons have complex relationships with Saturn’s interior, the rings, and with each other, which can be expressed in their shapes, interiors, and geology. Observations of their physical states can, thus, provide important constraints on the ages and formation mechanism(s) of the moons, which in turn informs our understanding of the formation and evolution of Saturn and its rings. Here, we describe the cratering records of the mid-sized moons and the value and limitations of their use for constraining the histories of the moons. We also discuss observational constraints on the interior structures of the moons and geologically-derived inferences on their thermal budgets through time. Overall, the geologic records of the moons (with the exception of Mimas) include evidence of epochs of high heat flows, short- and long-lived subsurface oceans, extensional tectonics, and considerable cratering. Curiously, Mimas presents no clear evidence of an ocean within its surface geology, but its rotation and orbit indicate a present-day ocean. While the moons need not be primordial to produce the observed levels of interior evolution and geologic activity, there is likely a minimum age associated with their development that has yet to be determined. Uncertainties in the populations impacting the moons makes it challenging to further constrain their formation timeframes using craters, whereas the characteristics of their cores and other geologic inferences of their thermal evolutions may help narrow down their potential histories. Disruptive collisions may have also played an important role in the formation and evolution of Saturn’s mid-sized moons, and even the rings of Saturn, although more sophisticated modeling is needed to determine the collision conditions that produce rings and moons that fit the observational constraints. Overall, the existence and physical characteristics of Saturn’s mid-sized moons provide critical benchmarks for the development of formation theories.

Abstract

Passive microwave radiometers onboard satellites rely on the received upwelling radiation to retrieve precipitation, which is a mixed signal from the surface, atmosphere and precipitation hydrometeors. Liquid precipitation droplets increase the upwelling radiation from the surface at lower frequencies, while ice particles cause a decrease in upwelling radiation at higher frequencies. The task of the retrieval algorithm is to identify the precipitation phase and to isolate the signal of precipitation from that of the surface. This study develops a machine learning method to retrieve rainfall and snowfall rates based on observations from the Microwave Hydrometer Sounder and Microwave Temperature Sounder onboard FY-3E. Self-organized mapping (SOM) is selected to classify the precipitation and underlying surface types, and an artificial neural network (ANN) is subsequently used to relate the brightness temperature to the precipitation rate for the clusters derived from the SOM. The half-hour product of the Integrated Multi-Satellite Retrieval for Global Precipitation Measurement (IMERG) is used to train the ANN. To address the issue that number of heavy precipitation samples are not enough for training, the simulation of radiative transfer for TOVS is used as a supplement to heavy rain samples. The SOM-ANN algorithm outperforms the IMERG and Goddard profiling algorithm (GPROF) retrieval products in both rainfall and snowfall retrieval. Compared with the hourly observations at ∼4,400 stations during a 2-year period, the root mean square errors of SOM-ANN proposed here are 1.06 and 0.34 mm/hr for the rainfall and snowfall rates, which are better than those of IMERG (1.23 and 0.42 mm/hr) and GPROF (1.22 and 0.44 mm/hr).

From cloud formation to sea temperatures, sea foam plays many roles in the dynamic interactions that occur at the surface level of the world's oceans.