QP Models Research Articles

Inversion Modeling of Chlorophyll Fluorescence Parameters in Cotton Canopy via Moisture Data and Spectral Analysis

The study of chlorophyll fluorescence parameters is very important for understanding plant photosynthesis. Monitoring cotton chlorophyll fluorescence parameters via spectral technology can aid in understanding the photosynthesis, growth, and stress of cotton fields in real time and provide support for cotton growth regulation and planting management. In this study, cotton plot experiments with different water treatments were set up to obtain the spectral reflectance of the cotton canopy, the maximum photochemical quantum yield (Fv/Fm), and the photochemical quenching coefficient (qP) of leaves at different growth stages. Support vector machine regression (SVR), random forest regression (RFR), and artificial neural network regression (ANNR) were used to establish a fluorescence parameter inversion model of the cotton canopy leaves. The results show that the original spectrum was transformed by multivariate scattering correction (MSC), the standard normal variable (SNV), and continuous wavelet transform (CWT), and the model constructed with Fv/Fm passed accuracy verification. The SNV-SVR model at the budding stage, the MSC-SVR model at the early flowering stage, the SNV-SVR model at the full flowering stage, the MSC-SVR model at the flowering stage, and the CWT-SVR model at the full boll stage had the highest estimation accuracy. The accuracies of the three spectral preprocessing and qP models were verified, and the MSC-SVR model at the budding stage, SNV-SVR model at the early flowering stage, MSC-SVR model at the full flowering stage, SNV-SVR model at the flowering stage, and CWT-SVR model at the full boll stage presented the highest estimation accuracies.

High-Resolution 3D QP and QS Models of the Middle Eastern Boundary of the Sichuan–Yunnan Rhombic Block: New Insight into Implication for Seismogenesis

Abstract The eastern boundary of the Sichuan–Yunnan rhombic block (EB-SYRB) has complex structures and strong seismicity. Although multiple 3D high-resolution velocity models have been constructed for this region, its seismogenic environment has been controversial. Seismic wave attenuation (inversely proportional to Q) describes the anelastic properties of the Earth’s medium, and is more sensitive to changes in subsurface fluid and temperature than seismic wave velocity. Based on the data of a long-term dense array in downstream of the Jinsha River, this article uses local earthquake tomography to obtain 3D QP and QS models of the middle EB-SYRB with the highest resolution to date, improving the lateral resolution of the Q model from 100 km to 5–10 km and the depth resolution from 10 to 2 km. Combined with the existing high-resolution velocity and resistivity models and geochemical observation results, we can comprehensively understand the medium structure and the seismogenesis in the study area. The results show that the high-attenuation characteristics in the shallow layer of the Xiaojiang fault zone and the Zemuhe fault zone (within a depth of ∼5 km) are consistent with the topographic relief and the distribution of hot springs, which reveals the Quaternary sedimentary characteristics of the basins and the presence of shallow fluids in the fault zone. The columnar high-attenuation anomaly beneath Huize reveals the fluid channel created by deep melting. The Ludian earthquake sequence occurred in a prominent low-attenuation area, which is favorable for stress accumulation and has a seismogenic environment for strong earthquakes. The high attenuation near the southwest end of the Huize fault is closely related to the Huize earthquake cluster, which is driven by fluids in the upper crust.

Enhancing Short-Term Prediction of BDS-3 Satellite Clock Bias Based with BSO Optimized BP Neural Network

The satellite clock bias (SCB) prediction plays an important role in high-accuracy and real-time navigation and positioning. When predicting the SCB, the performance of the BP neural network is affected by the local optimum due to inaccurate initial parameters. Therefore, we propose an improved BP neural network based on the beetle swarm optimization (BSO-BP) algorithm to improve the performance of SCB prediction in third-generation Beidou satellite navigation system (BDS-3). The proposed model takes advantage of group learning strategy to optimize the initialization parameters of the BP neural network and obtains globally optimized parameters. In order to verify the proposed BSO-BP model, 15 BDS satellites are analyzed in terms of prediction accuracy and stability of SCB. The experimental results show that when predicting 1 hour SCB based on a 12 hours SCB data, the prediction accuracy of the BSO-BP model is the best, with an average accuracy of 0.064 ns. As compared with the LP, QP, and GM models, the average prediction accuracy of the proposed BSO-BP model increases by about 72.6%, 43.4%, and 86%, respectively. As the prediction time increases, the influence of the inaccurate initial parameters on SCB prediction gradually decreases, and the prediction accuracy improves. The proposed BSO-BP model has the best accuracy and stability when predicting the 1 h SCB based on the same data. The prediction stability of the proposed BSO-BP model improves by more than 36% as compared with LP, QP, and GM models. In addition, the prediction accuracies of PHM clock and Rb-II clock improved by more than 47%, as compared with that of the Rb clock. Therefore, the overall performance of the atomic clock based on BDS-3 is better than BDS-2. The positioning accuracy of the BSO-BP model can reach the centimeter level in east, north, and up directions.

Double-difference seismic attenuation tomography method and its application to The Geysers geothermal field, California

SUMMARY Knowledge of attenuation structure is important for understanding subsurface material properties. We have developed a double-difference seismic attenuation (DDQ) tomography method for high-resolution imaging of 3-D attenuation structure. Our method includes two main elements, the inversion of event-pair differential ${t^*}$ ($d{t^*}$) data and 3-D attenuation tomography with the $d{t^*}$ data. We developed a new spectral ratio method that jointly inverts spectral ratio data from pairs of events observed at a common set of stations to determine the $d{t^*}$ data. The spectral ratio method cancels out instrument and site response terms, resulting in more accurate $d{t^*}$ data compared to absolute ${t^*}$ from traditional methods using individual spectra. Synthetic tests show that the inversion of $d{t^*}$ data using our spectral ratio method is robust to the choice of source model and a moderate degree of noise. We modified an existing velocity tomography code so that it can invert $d{t^*}$ data for 3-D attenuation structure. We applied the new method to The Geyser geothermal field, California, which has vapour-dominated reservoirs and a long history of water injection. A new Qp model at The Geysers is determined using P-wave data of earthquakes in 2011, using our updated earthquake locations and Vp model. By taking advantage of more accurate $d{t^*}$ data and the cancellation of model uncertainties along the common paths outside of the source region, the DDQ tomography method achieves higher resolution, especially in the earthquake source regions, compared to the standard tomography method using ${t^*}$ data. This is validated by both the real and synthetic data tests. Our Qp and Vp models show consistent variations in a normal temperature reservoir that can be explained by variations in fracturing, permeability and fluid saturation and/or steam pressure. A prominent low-Qp and Vp zone associated with very active seismicity is imaged within a high temperature reservoir at depths below 2 km. This anomalous zone is likely partially saturated with injected fluids.

Distributed Energy Management of the Hybrid AC/DC Microgrid with High Penetration of Distributed Energy Resources Based on ADMM

This paper aims to investigate energy management of the hybrid AC/DC microgrid with the high penetration of distributed energy resources (DERs), such as electrical vehicles, heat pumps, and photovoltaics. In the previous studies, energy management of the hybrid microgrid is usually carried out by the system operator in a centralized manner, which suffers from the compromise of privacy information protection and the risk of single‐point failure. Therefore, this paper proposes a distributed energy management scheme of the hybrid microgrid using the projection function‐based alternating direction method of multipliers (P‐ADMM), which allows each subgrid, i.e., AC subgrid and DC subgrid, to make day‐ahead schedules independently with information exchanges while obtaining the optimal energy management solution. The energy management problem of the hybrid microgrid is formulated as a mixed‐integer quadratic programming (MIQP) model, considering DER and energy storage system operation constraints, system operation constraints, and converter operation constraints. Then, the MIQP model is decomposed and distributed into smaller‐scale QP models between subgrids using the P‐ADMM algorithm, which can handle binary variables through projection functions. The numerical results conducted on the hybrid microgrid demonstrate that the proposed distributed scheme can effectively achieve optimal energy management for the hybrid AC/DC microgrid in a distributed manner.

Joint local earthquake and teleseismic inversion for 3-D velocity and Q in New Zealand

We have developed a joint inversion of local earthquake and teleseismic data and applied it in New Zealand to obtain improved 3-D seismic velocity, Qp and Qs models. To accommodate sharper detail where there are dense local earthquake data, varied spatial resolution is enabled with linking of nodes into larger inversion elements where data are sparse. The teleseismic t∗ measurements are obtained by stacking and removing the average for the common source and path. The combined inversions used data from 126 teleseisms and 2365 local earthquakes that extend to 350 km depth. Vp results are similar to previous studies, but have enhanced detail near the slab provided by local seismicity. Sub-slab low velocity and low Q is imaged along the eastern North Island, which may represent asthenosphere with entrained water and melt. Along the western South Island, there is a high Vp, high Q feature distinct from the current slab, and interpreted to be a slab remnant from the period of Subduction-Transform Edge Propagator fault development and trench rollback. The axial South Island high Vp, high Q feature of the Hikurangi Plateau is imaged, and is bounded on the southeast by thinner lithosphere that may be related to the period of Great South Basin rifting.

The effect of crude protein reduction on performance and nitrogen metabolism in piglets (four to nine weeks of age) fed two dietary lysine levels1.

Lowering the CP level in piglet diets reduces the risk of postweaning diarrhea and N excretion to the environment. The question remains at what point CP becomes limiting. An experiment was designed with 2 standardized ileal digestible (SID) Lys levels (10 and 11 g) and 6 CP levels (140, 150, 160, 170, 180, 190 g/kg) in a 2 × 6 factorial design (with 6 pens of 6 animals each per treatment). Linear and quadratic (QP) mixed models of performance in function of CP were fitted to study the effect of SID Lys and CP and their interaction. To determine optima, QP models and broken line models with linear (BLL) or quadratic (BLQ) ascending portions were fitted through the data. It was hypothesized 1) that the response to a decreasing digestible CP level could be described with broken line models and 2) that the break point of these models is dependent on the dietary SID Lys level. Decreasing the CP level decreased ADG (P < 0.001). For G:F, the effect of decreasing CP level depended on the SID Lys level (P of the interaction = 0.028 in the linear model and P = 0.002 in the QP model). According to the BLL model, with 11 g SID Lys in the diet, G:F started to decline with CP levels < 176 g CP [SID Lys:CP = 0.062, SID Lys:apparent total tract digestible (ATTD) CP = 0.077], and with 10 g SID Lys, CP levels < 165 g/kg (SID Lys:CP = 0.061, SID Lys:ATTD CP = 0.075) depressed performance. Serum creatinine levels showed a linear decrease with increasing SID Lys:CP levels (P < 0.001). Across both SID Lys levels, when fitting a BLL model, minimal serum urea levels were reached at an SID Lys:CP ratio of 0.064. This seems to be the point where CP and not Lys limits muscle deposition. The small difference in break point between serum urea level and performance suggests that the composition of nonessential AA may also be at stake. The effect of decreasing CP level depends on SID Lys, and using a maximal SID Lys:CP ratio may be useful for optimizing the AA profile of dietary CP. When the SID Lys:CP ratio exceeds 0.064 (SID Lys:ATTD CP > 0.079), protein and not individual AA limits growth in most piglets between 4 and 9 wk of age.

Pseudo 2D elastic waveform inversion for attenuation in the near surface

Seismic waveform propagation could be significantly affected by heterogeneities in the near surface zone (0m–500m depth). As a result, it is important to obtain as much near surface information as possible. Seismic attenuation, characterized by QP and QS factors, may affect seismic waveform in both phase and amplitude; however, it is rarely estimated and applied to the near surface zone for seismic data processing. Applying a 1D elastic full waveform modelling program, we demonstrate that such effects cannot be overlooked in the waveform computation if the value of the Q factor is lower than approximately 100. Further, we develop a pseudo 2D elastic waveform inversion method in the common midpoint (CMP) domain that jointly inverts early arrivals for QP and surface waves for QS. In this method, although the forward problem is in 1D, by applying 2D model regularization, we obtain 2D QP and QS models through simultaneous inversion. A cross-gradient constraint between the QP and Qs models is applied to ensure structural consistency of the 2D inversion results. We present synthetic examples and a real case study from an oil field in China.

Entropy-inspired Lyapunov Functions and Linear First Integrals for Positive Polynomial Systems

Two partially overlapping classes of positive polynomial systems, chemical reaction networks with mass action law (MAL-CRNs) and quasi-polynomial systems (QP systems) are considered. Both of them have an entropy-like Lyapunov function associated to them which are similar but not the same. Inspired by the work of Prof. Gorban [12] on the entropy-functionals for Markov chains, and using results on MAL-CRN and QP-systems theory we characterize MAL-CRNs and QP systems that enable both types of entropy-like Lyapunov functions. The starting point of the analysis is the class of linear weakly reversible MAL-CRNs that are mathematically equivalent to Markov chains with an equilibrium point where various entropy level set equivalent Lyapunov functions are available. We show that non-degenerate linear kinetic systems with a linear first integral (that corresponds to conservation) can be transformed to linear weakly reversible MAL-CRNs using linear diagonal transformation, and the coefficient matrix of this system is diagonally stable. This implies the existence of the weighted version of the various entropy level set equivalent Lyapunov functions for non-degenerate linear kinetic systems with a linear first integral. Using translated X-factorable phase space transformations and nonlinear variable transformations a dynamically similar linear ODE model is associated to the QP system models with a positive equilibrium point. The non-degenerate kinetic property together with the existence of positive equilibrium point form a sufficient condition of the existence of the weighted version of the various entropy level set equivalent Lyapunov functions in this case. Further extension has been obtained by using the time re-parametrization transformation defined for QP models.

Seismic attenuation tomography of the Northeast Japan arc: Insight into the 2011 Tohoku earthquake (Mw9.0) and subduction dynamics

Detailed three-dimensional (3-D) P and S wave attenuation (Qp and Qs) models of the crust and upper mantle under the entire Northeast Japan (Tohoku) arc from the Japan Trench to the Japan Sea coast are determined, for the first time, using a large number of high-quality t* data measured precisely from P and S wave spectra of local earthquakes. The suboceanic earthquakes used in this work are relocated precisely using sP depth phases. Our results reveal a prominent landward dipping high-Q zone representing the subducting Pacific slab, a landward dipping intermediate- to high-Q zone in the mantle wedge between the Pacific coast and the volcanic front, and significant low-Q anomalies in the crust and mantle wedge between the volcanic front and the Japan Sea coast. Prominent high-Q patches surrounded by low-Q anomalies are revealed in the interplate megathrust zone under the Tohoku fore arc where the great 2011 Tohoku-oki earthquake (Mw 9.0) occurred. The high-Q patches in the megathrust zone generally exhibit large coseismic slips of megathrust earthquakes and large slip deficit on the plate interface. We think that these high-Q patches represent asperities in the megathrust zone, whereas the low-Q anomalies reflect weakly coupled areas. We also find that the hypocenters of the 2011 Tohoku-oki interplate earthquakes (Mw > 7.0) are located in areas where Qp, Qs, and Qp/Qs change abruptly. These results suggest that structural heterogeneities in the megathrust zone control the interplate seismic coupling and the nucleation of megathrust earthquakes.

Inferring water infiltration in the Longtan reservoir area by three-dimensional attenuation tomography

SUMMARY High-resolution 3-D QP and QS models at depths of 0∼14 km in the Longtan reservoir area, which were used to further infer the water infiltration around the reservoir, were constructed by the local earthquake tomography technique. The inversion used t* values measured from velocity spectral amplitudes of 16 134 P-wave and 17 027 S-wave observations of 2828 events recorded by 23 digital seismic stations in the Longtan reservoir area. 3-D QP and QS tomographic images suggest a close relation between the attenuation models and water infiltration. Low QP and QS exist down to 4∼7 km depths beneath the main rivers in the reservoir area, whereas high QP and QS are found both at deep layers and at shallow layers in regions other than the main river area. We inferred that the rocks beneath the main rivers are infiltrated by substantial water and water infiltration reaches 4∼7 km deep, and the ranges of water infiltration vary by region. Specifically, the zone of high water content is estimated to extend up to 10 km away from the Buliu river in the reservoir dam area where water depth in the reservoir reaches maximum, and the extent of water infiltration in rocks can only be about 5 km around from the Hongshui river and the Youla river. Low QP and high QS indicate that the carbonate rocks in the east of the Bala fault area may be saturated with water. Earthquake clusters are mainly concentrated on the low-QP and low-QS zones, indicating that water infiltration beneath the main rivers possibly plays an important role in inducing earthquakes in the Longtan reservoir area.

Three‐dimensional P wave attenuation and velocity upper mantle tomography of the southern Apennines–Calabrian Arc subduction zone

We propose a 3‐D crust–upper mantle seismic attenuation (QP) model of the southern Apennines–Calabrian Arc subduction zone together with a 3‐D velocity (VP) model. The QP model is calculated from relative t* using the spectral ratio method and the VP from traveltime data. The final data set used for the inversion of the VP model consists of 2400 traveltime arrivals recorded by 34 short‐period stations that are part of the Italian National Seismic Network, and for the QP model, 2178 Pn phases recorded by a subset of 32 stations. Traveltimes and waveforms come from 272 intermediate‐depth Calabrian slab events. This 3‐D model of attenuation, together with the 3‐D velocity model, improves our knowledge of the slab/mantle wedge structure and can be a starting point in determining the physical state of the asthenosphere (i.e., its temperature, the presence of melt and/or fluids) and its relation to volcanism found in the study area. Main features of the QP and VP models show that the mantle wedge/slab, in particular, the area of highest attenuation, is located in a volume underlying the Marsili Basin. The existence and shape of this main low‐QP (and low‐VP) anomaly points to slab dehydration and fluid/material flow, a process that may explain the strong geochemical affinities between the subduction‐related magmas from Stromboli and Vesuvius. Other interesting features in the models are strong lateral variations in QP and VP that are put in relation with known important tectonic structures and volcanic centers in the area.

P- andS-velocity images of the lithosphere-asthenosphere system in the Central Andes from local-source tomographic inversion

SUMMARY About 50 000 P and S arrival times and 25 000 values of t* recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P, S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp, Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting.

Attenuation models (QP and QS) in three dimensions of the southern California crust: Inferred fluid saturation at seismogenic depths

We analyze high dynamic range waveform spectra to determine t* values for both P and S waves from earthquakes in southern California. We invert the t* values for three‐dimensional (3‐D) frequency‐independent QP and QS regional models of the crust. The models have 15 km horizontal grid spacing and an average vertical grid spacing of 4 km, down to 22 km depth, and extend from the U.S.‐Mexico border to the Coast Ranges in the south and Sierra Nevada in the north. In general, QP and QS increase rapidly with depth, consistent with crustal densities and velocities. The 3‐D QP and QS models image the major tectonic structures and to a much lesser extent the thermal structure of the southern California crust. The near‐surface low QP and QS zones coincide with major sedimentary basins such as the San Bernardino, Chino, San Gabriel Valley, Los Angeles, Ventura, and Santa Maria basins and the Salton Trough. In contrast, at shallow depths beneath the Peninsular Ranges, southern Mojave Desert, and southern Sierras, we image high QP and QS zones, which correspond to the dense and high‐velocity rocks of the mountain ranges. Several clear transition zones of rapidly varying QP and QS coincide with major late Quaternary faults and connect regions of high and low QP and QS. At midcrustal depths, the QP and QS exhibit modest variation in slightly higher and lower QP or QS zones, which is consistent with reported crustal reflectivity. In general, for the southern California crust, QS/QP is greater than 1.0, suggesting partially fluid‐saturated crust. A few limited regions of QS/QP less than 1.0 correspond to areas mostly outside the major sedimentary basins, including areas around the San Jacinto fault, suggesting a larger reduction in the shear modulus compared to the bulk modulus or almost complete fluid saturation.

A parameterized multifrequency-polarization surface emission model

This study develops a parameterized bare surface emission model for the applications in analyses of the passive microwave satellite measurements from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). We first evaluated the capability of the advanced integral equation model (AIEM) in simulating wide-band and high-incidence surface emission signals in comparison with INRA's field experimental data obtained in 1993. The evaluation results showed a very good agreement. With the confirmed confidence, we generated a bare surface emission database for a wide range of surface dielectric and roughness properties under AMSR-E sensor configurations using the AIEM model. Through the evaluations of the commonly used semiempirical models with both the AIEM simulated and the field experimental data, we developed a parameterized multifrequency-polarization surface emission model-the Qp model. This model relates the effects of the surface roughness on the emission signals through the roughness variable Qp at the polarization p. The Qp can be simply described as a single-surface roughness property-the ratio of the surface rms height and the correlation length. The comparison of the emissivity simulations by the Qp and AIEM models indicated that the absolute error is extremely small at the magnitude of 10/sup -3/. The newly developed surface emission model should be very useful in modeling, improving our understanding, analyses, and predictions of the AMSR-E measurements.

Staggered-grid pseudospectral viscoacoustic wave field simulation in two-dimensional media

A viscoacoustic numerical modeling algorithm is developed, tested, and implemented for two-dimensional media. The loss mechanism is simulated according to a standard linear solid model. With anelastic effects defined through superposition of relaxation mechanisms, and incorporating them into the wave equation in terms of memory variables, the time-domain computation of responses of viscoacoustic media with arbitrary spatial variations is feasible. A numerical solution of viscoacoustic responses is obtained by solving either two-dimensional homogeneous or heterogeneous viscoacoustic wave equations. The model is parametrized through stress, strain relaxation times, density, relaxed and unrelaxed modulus, and velocity. On constructing P and S velocity, Qp, and Qs models separately, attenuated and decoupled propagating P and S wave fields are synthesized. According to this scheme one can model widely varying Q in highly heterogeneous media. Simulation of a pure acoustic wave field is regarded as a special case when temporally invariant stress and strain variables are considered. Through systematic analysis, a scheme with a staggered grid is highly desirable for stable computation and reduction of numerical artifacts such as grid dispersion and grid anisotropy. Considerable Nyquist error may arise when a regular grid for a first-derivative operator is used; a staggered grid enables significant improvement. Implementation of a pseudospectral technique with staggered grid provides a more general and more accurate description of wave propagation phenomena. Results from viscoacoustic simulation show that considerable variation of dynamic properties such as a decreased amplitude due to intrinsic Q and temporal shift due to physical wave dispersion phenomena can be analyzed quantitatively. Combined effects of source directivity, complex velocity and attenuation structures may alter significantly in the computed synthetic seismograms.