Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (10,105)

Search Parameters:
Keywords = coupling mechanism

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
26 pages, 4319 KiB  
Article
POI Data–Driven Identification and Representation of Production–Living–Ecological Spaces at the Urban and Peri–Urban Scale: A Case Study of the Hohhot–Baotou–Ordos–Yulin Urban Agglomeration
by Shuai Zhang, Yixin Fang and Xiuqing Zhao
Sustainability 2025, 17(5), 2235; https://doi.org/10.3390/su17052235 - 4 Mar 2025
Abstract
The identification of the multifunctional combination of production–living–ecological spaces (PLES) in urban agglomerations, particularly in urban cores and peri–urban areas, is a critical issue in the urbanization process. This study, using the Hohhot–Baotou–Ordos–Yulin (HBOY) urban agglomeration, a key node in China’s “Two Horizontals [...] Read more.
The identification of the multifunctional combination of production–living–ecological spaces (PLES) in urban agglomerations, particularly in urban cores and peri–urban areas, is a critical issue in the urbanization process. This study, using the Hohhot–Baotou–Ordos–Yulin (HBOY) urban agglomeration, a key node in China’s “Two Horizontals and Three Verticals” urbanization strategy, proposes a hexagonal grid–based PLES quantification framework using POI data. A three–level POI classification system was developed, with functional element weights determined via the Analytic Hierarchy Process and public perception surveys. The framework quantifies PLES within hexagonal grids and analyzes its patterns and functional coupling mechanisms using spatial overlay, Average Nearest Neighbor Index (ANNI), kernel density analysis, and spatial autocorrelation analysis. The following results were obtained. (1) PLES classification accuracy reached 90.83%, confirming the reliability of the method. (2) The HBOY urban agglomeration exhibits a dominant production space (40.84%), balanced living and ecological spaces (29.37% and 29.36%, respectively), and a severe shortage of mixed spaces (0.43%). (3) Production and living spaces show significant clustering ( ≤ 0.581), mixed spaces follow ( = 0.660), and ecological spaces are relatively evenly distributed (= 0.870). (4) The spatial distribution patterns show that production and living spaces exhibit “core concentration with peripheral dispersion”, ecological spaces show “block concentration with point–like distribution”, and mixed spaces show “point–like dispersion”. (5) Production and living spaces exhibit strong spatial autocorrelation ( > 0.7) and the highest spatial correlation (= 0.692), while the spatial correlation with ecological spaces is weakest ( = 0.150). The proposed PLES identification framework, with its efficiency and dynamic updating potential, provides an innovative approach to urban spatial governance under the global Sustainable Development Goals. The findings offer integrated decision–making support for spatial diagnosis and functional regulation in the ecologically vulnerable areas of northwest China’s new urbanization. Full article
21 pages, 3900 KiB  
Article
A Study of Forced Vibrations with Nonlinear Springs and Dry Friction: Application to a Mechanical Oscillator with Very Large Vibrating Blades for Soil Cutting
by Dario Friso
Vibration 2025, 8(1), 10; https://doi.org/10.3390/vibration8010010 - 4 Mar 2025
Viewed by 65
Abstract
To cut a clod of soil containing the roots of trees in nurseries, a semi-circular vibrating blade digging machine with diameters up to 1.2 m is increasingly used. The heart of the machine is the mechanical oscillator that produces an excitation torque supplied [...] Read more.
To cut a clod of soil containing the roots of trees in nurseries, a semi-circular vibrating blade digging machine with diameters up to 1.2 m is increasingly used. The heart of the machine is the mechanical oscillator that produces an excitation torque supplied to the blade together with the cutting torque of the soil. The advantage of the vibrating blade is a reduction in the cutting torque up to 70%. This advantage led to the present study of the extension to blades of even 1.8 m for the digging of very large trees. To build an oscillator suitable for all blade sizes (from 0.6 to 1.8 m), it was necessary to equip it with nonlinear (quadratic) springs, since with traditional linear springs, it would not be versatile. The design and simulation of its operation required the development of a new mathematical model. Therefore, an approximate solution of the differential equation of the forced vibration with quadratic springs and dry friction between the blade and soil was developed, aimed at calculating the maximum blade displacement and the phase lag. These quantities, together with the cutting time, had to satisfy certain values to ensure the maximum reduction in the cutting torque (−70%). After the construction of the oscillator, it was coupled with all the blades (0.6, 0.9, 1.2, and 1.8 m) for experimental tests. For all diameters, the oscillator was able to optimally vibrate the blades, preventing the springs from reaching the end of the stroke when cutting the soil. Measuring the maximum blade displacement compared with the calculated one provided a good accuracy of the mathematical modeling, resulting in a mean error of 5.6% and a maximum error of 7.2%. Full article
Show Figures

Figure 1

Figure 1
<p>The semicircular blade and the mechanical oscillator (right) of the tree digger machine.</p>
Full article ">Figure 2
<p>The oscillator blade system: (1) gear wheels; (2) eccentric masses; (3) gear housing; (4) shaft; (5) horizontal butterfly bush; (6) semicircular blade; (7) springs; (8) vertical butterfly bush; (9) worm screw; (10) worm gear housing.</p>
Full article ">Figure 3
<p>Sketch of conical coil spring. The quantity <span class="html-italic">p</span> is the variable pitch of the coil spring.</p>
Full article ">Figure 4
<p>Diagram of spring force <span class="html-italic">F</span> versus shortening <span class="html-italic">x</span> of the conical coil spring with variable pitch of <a href="#vibration-08-00010-f003" class="html-fig">Figure 3</a>. It appears as a polygonal chain (black), next to which a regression parabola (red) is drawn.</p>
Full article ">Figure 5
<p>Single degree of freedom system to which the oscillating blade within the soil can be traced under the effect of the excitation torque <span class="html-italic">T<sub>em</sub></span>, the friction torque <span class="html-italic">T<sub>F</sub></span>, the springs torque <span class="html-italic">k·</span>(<span class="html-italic">α</span> + <span class="html-italic">α<sub>cut</sub></span>) + <span class="html-italic">K</span>·(<span class="html-italic">α</span> + <span class="html-italic">α<sub>cut</sub></span>)<sup>2</sup>, and the cutting torque <span class="html-italic">T<sub>cut</sub></span>.</p>
Full article ">Figure 6
<p>Calculated maximum angular displacement <span class="html-italic">α<sub>0</sub></span> vs. diameter of the blade <span class="html-italic">D</span>.</p>
Full article ">Figure 7
<p>Calculated maximum blade displacement <span class="html-italic">A</span><sub>0</sub> vs. diameter of the blade <span class="html-italic">D</span>.</p>
Full article ">Figure 8
<p>Calculated phase angle <span class="html-italic">ϕ</span> vs. diameter of the blade <span class="html-italic">D</span>.</p>
Full article ">Figure 9
<p>Calculated maximum spring shortening <math display="inline"><semantics> <mrow> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mi>c</mi> <mi>u</mi> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </mfenced> <mo>=</mo> <mfenced separators="|"> <mrow> <msub> <mrow> <mi>α</mi> </mrow> <mrow> <mi>c</mi> <mi>u</mi> <mi>t</mi> </mrow> </msub> <mo>+</mo> <msub> <mrow> <mi>α</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </mfenced> <mo>·</mo> <msub> <mrow> <mi>b</mi> </mrow> <mrow> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> during the soil cutting vs. total available space <span class="html-italic">S<sub>t</sub></span> and vs. diameter of the blade <span class="html-italic">D</span>.</p>
Full article ">Figure 10
<p>Calculated soil cutting time <span class="html-italic">t<sub>cut</sub></span> for the entire hemispherical soil clod vs. diameter of the blade <span class="html-italic">D</span>.</p>
Full article ">Figure 11
<p>Calculated oscillating power <span class="html-italic">P<sub>o</sub></span>, cutting power <span class="html-italic">P<sub>cut</sub></span>, and total power <span class="html-italic">P<sub>t</sub></span> vs. blade diameter <span class="html-italic">D</span>.</p>
Full article ">Figure 12
<p>Calculated total energy <span class="html-italic">W<sub>t</sub></span> during the cutting of the entire hemispherical soil clod vs. blade diameter <span class="html-italic">D</span>.</p>
Full article ">
16 pages, 3320 KiB  
Article
A Spike Train Production Mechanism Based on Intermittency Dynamics
by Stelios M. Potirakis, Fotios K. Diakonos and Yiannis F. Contoyiannis
Entropy 2025, 27(3), 267; https://doi.org/10.3390/e27030267 - 4 Mar 2025
Viewed by 19
Abstract
Spike structures appear in several phenomena, whereas spike trains (STs) are of particular importance, since they can carry temporal encoding of information. Regarding the STs of the biological neuron type, several models have already been proposed. While existing models effectively simulate spike generation, [...] Read more.
Spike structures appear in several phenomena, whereas spike trains (STs) are of particular importance, since they can carry temporal encoding of information. Regarding the STs of the biological neuron type, several models have already been proposed. While existing models effectively simulate spike generation, they fail to capture the dynamics of high-frequency spontaneous membrane potential fluctuations observed during relaxation intervals between consecutive spikes, dismissing them as random noise. This is eventually an important drawback because it has been shown that, in real data, these spontaneous fluctuations are not random noise. In this work, we suggest an ST production mechanism based on the appropriate coupling of two specific intermittent maps, which are nonlinear first-order difference equations. One of these maps presents small variation in low amplitude values and, at some point, bursts to high values, whereas the other presents the inverse behavior, i.e., from small variation in high values, bursts to low values. The suggested mechanism proves to be able to generate the above-mentioned spontaneous membrane fluctuations possessing the associated dynamical properties observed in real data. Moreover, it is shown to produce spikes that present spike threshold, sharp peak and the hyperpolarization phenomenon, which are key morphological characteristics of biological spikes. Furthermore, the inter-spike interval distribution is shown to be a power law, in agreement with published results for ST data produced by real biological neurons. The use of the suggested mechanism for the production of other types of STs, as well as possible applications, are discussed. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Biological membrane potential <math display="inline"><semantics> <mrow> <mi>V</mi> <mfenced separators="|"> <mrow> <mi>t</mi> </mrow> </mfenced> </mrow> </semantics></math> of neuron 089s1c1 from in vitro intracellular recordings of CA1 pyramidal neurons of Wistar male rats (adopted from [<a href="#B17-entropy-27-00267" class="html-bibr">17</a>]). (<b>b</b>) Zoom-in to a spike from <a href="#entropy-27-00267-f001" class="html-fig">Figure 1</a>a, along with pre- and after-spike high-frequency fluctuations. The horizontal dashed red line, denoting the firing threshold, highlights the fact that the after-spike mean level is lower than the pre-spike one (hyperpolarization phenomenon).</p>
Full article ">Figure 2
<p>An intermittent time series in which it has been considered that the region from 0 up to 0.2 (bounded upwards by the red horizontal line) is the laminar region and all values above this zone correspond to bursts.</p>
Full article ">Figure 3
<p>Typical examples of the return plots of the <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>1</mn> </mrow> </semantics></math>, given by Equation (1) (green), and the <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>2</mn> </mrow> </semantics></math>, given by Equation (3) (blue). The laminar region in both maps is the part of the trajectory that closely follows (is almost parallel to) the bisector (red line). A zoom-in to each laminar region, <math display="inline"><semantics> <mrow> <mi>l</mi> <mi>a</mi> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mi>a</mi> <mi>r</mi> <mn>1</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>l</mi> <mi>a</mi> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mi>a</mi> <mi>r</mi> <mn>2</mn> </mrow> </semantics></math>, is presented in the corresponding insets. As soon as the trajectory begins to move away from the bisector, the map has entered the bursts region. Such plots can be constructed for various combinations of the parameters’ values in Equations (1) and (3), moving the two laminar regions closer or further away from each other.</p>
Full article ">Figure 4
<p>Suggested spike train production mechanism based on intermittency dynamics.</p>
Full article ">Figure 5
<p>Return plots of <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>1</mn> </mrow> </semantics></math> (green) and <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>2</mn> </mrow> </semantics></math> (blue), uncoupled, for the maps parameters values used in our numerical experiment (see text). Red line denotes the bisector.</p>
Full article ">Figure 6
<p>Pseudocode for the production of the ST time series presented in <a href="#entropy-27-00267-f007" class="html-fig">Figure 7</a>.</p>
Full article ">Figure 7
<p>(<b>a</b>) The time series produced by the ST production mechanism suggested in <a href="#sec3dot1-entropy-27-00267" class="html-sec">Section 3.1</a>, using the parameters: <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>1</mn> </mrow> </semantics></math> {<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>z</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>4</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>u</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>0.011</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ε</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>0.0175</mn> </mrow> </semantics></math>}, <math display="inline"><semantics> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> <mn>2</mn> </mrow> </semantics></math> {<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>z</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>u</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>17</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ε</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>0.07</mn> </mrow> </semantics></math>}, switching thresholds {<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>T</mi> <mi>h</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>0.31</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>T</mi> <mi>h</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>}. Only the first 300,000 points of the produced 3,000,000-points-long time series are shown, so that the spike pattern is clearly visible. (<b>b</b>) A spike from the time series of <a href="#entropy-27-00267-f007" class="html-fig">Figure 7</a>a where the spike threshold and the hyperpolarization phenomenon are marked.</p>
Full article ">Figure 8
<p>(<b>a</b>) The ST time series produced by the numerical experiment of <a href="#sec3dot2-entropy-27-00267" class="html-sec">Section 3.2</a> (also depicted in <a href="#entropy-27-00267-f007" class="html-fig">Figure 7</a>a). The laminar region of the high-frequency fluctuations of the relaxation intervals was found to be bound between <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>r</mi> <mi>e</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mo>−</mo> <mn>0.64</mn> </mrow> </semantics></math> (red horizontal line) and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>ϕ</mi> </mrow> <mrow> <mi>b</mi> <mi>l</mi> <mi>u</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mn>0.07</mn> </mrow> </semantics></math> (blue horizontal line). (<b>b</b>) The distribution of laminar lengths resulting from the laminar region marked in <a href="#entropy-27-00267-f008" class="html-fig">Figure 8</a>a. The estimated exponent values by fitting Equation (7) are: <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>p</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mn>1.372</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>p</mi> </mrow> <mrow> <mn>3</mn> </mrow> </msub> <mo>=</mo> <mn>0.001</mn> </mrow> </semantics></math>, with goodness of fit <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>R</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> <mo>=</mo> <mn>1.000</mn> </mrow> </semantics></math>. Although only the first 300,000 points of the analyzed time series are shown in <a href="#entropy-27-00267-f008" class="html-fig">Figure 8</a>a, the distribution was calculated using the total length of the time series (3,000,000 points).</p>
Full article ">Figure 9
<p>The distribution of the inter-spike intervals of the 3,000,000-points-long ST time series produced in the numerical experiment of <a href="#sec3dot2-entropy-27-00267" class="html-sec">Section 3.2</a> is a power law of the form <math display="inline"><semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>s</mi> <mo>)</mo> <mo>~</mo> <msup> <mrow> <mi>s</mi> </mrow> <mrow> <mo>−</mo> <mn>1.372</mn> </mrow> </msup> </mrow> </semantics></math>.</p>
Full article ">
10 pages, 585 KiB  
Article
The Quantum Relative Entropy of the Schwarzschild Black Hole and the Area Law
by Ginestra Bianconi
Entropy 2025, 27(3), 266; https://doi.org/10.3390/e27030266 - 4 Mar 2025
Viewed by 29
Abstract
The area law obeyed by the thermodynamic entropy of black holes is one of the fundamental results relating gravity to statistical mechanics. In this work, we provide a derivation of the area law for the quantum relative entropy of the Schwarzschild black hole [...] Read more.
The area law obeyed by the thermodynamic entropy of black holes is one of the fundamental results relating gravity to statistical mechanics. In this work, we provide a derivation of the area law for the quantum relative entropy of the Schwarzschild black hole for an arbitrary Schwarzschild radius. The quantum relative entropy between the metric of the manifold and the metric induced by the geometry and the matter field has been proposed in G. Bianconi as the action for entropic quantum gravity leading to modified Einstein equations. The quantum relative entropy generalizes Araki’s entropy and treats the metrics between zero-forms, one-forms, and two-forms as quantum operators. Although the Schwarzschild metric is not an exact solution of the modified Einstein equations of the entropic quantum gravity, it is an approximate solution valid in the low-coupling, small-curvature limit. Here, we show that the quantum relative entropy associated to the Schwarzschild metric obeys the area law for a large Schwarzschild radius. We provide a full statistical mechanics interpretation of the results. Full article
Show Figures

Figure 1

Figure 1
<p>Diagrammatic sketch of the entropic quantum gravity approach. In this approach the action is given by the quantum relative entropy between the metric <math display="inline"><semantics> <mover accent="true"> <mi>g</mi> <mo stretchy="false">˜</mo> </mover> </semantics></math> and the metric <math display="inline"><semantics> <mover accent="true"> <mi mathvariant="bold">G</mi> <mo stretchy="false">˜</mo> </mover> </semantics></math> induced by the matter fields and the geometry of the manifold.</p>
Full article ">Figure 2
<p>The quantum relative entropy of the Schwarzschild metric <math display="inline"><semantics> <mi mathvariant="script">S</mi> </semantics></math> divided by its asymptotic expression <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="script">S</mi> <mi>A</mi> </msub> <mo>=</mo> <mi mathvariant="script">C</mi> <mi>A</mi> <mo>/</mo> <mrow> <mo>(</mo> <mn>4</mn> <mi>G</mi> <mo>)</mo> </mrow> </mrow> </semantics></math>, obeying the area law, is plotted as a function of the Schwarzschild radius <math display="inline"><semantics> <msub> <mi>R</mi> <mi>s</mi> </msub> </semantics></math> for <math display="inline"><semantics> <mrow> <mi>G</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">
20 pages, 47140 KiB  
Article
Analysis of the Dominant Factors and Interannual Variability Sensitivity of Extreme Changes in Water Use Efficiency in China from 2001 to 2020
by Shubing Hou, Wenli Lai, Jie Zhang, Yichen Zhang, Wenjie Liu, Feixiang Zhang and Shuqi Zhang
Forests 2025, 16(3), 454; https://doi.org/10.3390/f16030454 - 4 Mar 2025
Viewed by 79
Abstract
Ecosystem water use efficiency (WUE) is a key indicator of the coupling between carbon and water cycles. With the increasing frequency of extreme climate events, WUE may also show trends of extremization. Understanding the dominant drivers behind extreme WUE variations is crucial for [...] Read more.
Ecosystem water use efficiency (WUE) is a key indicator of the coupling between carbon and water cycles. With the increasing frequency of extreme climate events, WUE may also show trends of extremization. Understanding the dominant drivers behind extreme WUE variations is crucial for assessing the impact of climate variability on WUE. We investigate the main drivers and regional sensitivity of extreme WUE variations across seven geographical regions in China. The results reveal that extreme WUE variations are collectively influenced by gross primary productivity (GPP) and evapotranspiration (ET) (43.72%). GPP controls extreme WUE variations in 36.00% of the areas, while ET controls 20.17%. Furthermore, as the climate shifts from arid to humid regions, the area where GPP dominates extreme WUE variations increases, while the area dominated by ET decreases, suggesting a relationship with precipitation. Ridge regression analysis shows that vapor pressure deficit (VPD) is the primary driver of interannual WUE variation in China, with an average relative contribution of 38.64% and an absolute contribution of 0.025 gC·m−2·mm−1·a−1. We studied the changes in WUE and its driving mechanisms during extreme disaster events, providing a perspective focused on extreme conditions. In the future, these results may help regulate the carbon–water cycle in different regions, such as by guiding vegetation planting and land use planning based on the spatial characteristics of the dominant factors influencing extreme WUE variations to improve vegetation WUE. Full article
(This article belongs to the Section Forest Hydrology)
Show Figures

Figure 1

Figure 1
<p>Seven distinct geographical regions (divided by black lines) and different vegetation types (represented by different colors) are shown, with vegetation distribution illustrated using 2010 as an example (the midpoint year). WB: water body, ENF: evergreen needleleaf forest, EBF: evergreen broadleaf forest, DNF: deciduous needleleaf forest, DBF: deciduous broadleaf forest, MXF: mixed forest, WSN: woody savanna, SN: savanna, GL: grassland, CL: cropland, UB: urban and built-up land, BN: barren (The vegetation type data in the figure is the result of mode resampling of the land use data from <a href="#sec2dot2dot2-forests-16-00454" class="html-sec">Section 2.2.2</a>, with a resolution of 0.5°.).</p>
Full article ">Figure 2
<p>Spatial distribution of the dominant factor (either GPP, ET, or GPP and ET) controlling annual WUE extreme variation during 2001–2020. (<b>a</b>) WUE<sub>max</sub>; (<b>c</b>) WUE<sub>min</sub>. Area covered (in percentage) by contributions of the driving factors. (<b>b</b>) WUE<sub>max</sub>; (<b>d</b>) WUE<sub>min</sub>.</p>
Full article ">Figure 3
<p>The area of regions controlled solely by GPP and ET under extreme vegetation WUE conditions across seven geographical regions from 2001 to 2020, along with their fitted estimates against regional annual precipitation. The green line represents GPP, and the red line represents ET. (<b>a</b>) WUE<sub>max</sub>, GPP controlled; (<b>b</b>) WUE<sub>min</sub>, GPP controlled; (<b>c</b>) WUE<sub>max</sub>, ET controlled; (<b>d</b>) WUE<sub>min</sub>, ET controlled.</p>
Full article ">Figure 4
<p>The spatial pattern of the SPEI-12 corresponding to extreme WUE (<b>a</b>) WUE<sub>max</sub> and (<b>c</b>) WUE<sub>min</sub> at annual scale. The regional proportions of statistical distribution for drought years, normal years, and humid years. (<b>b</b>) WUE<sub>max</sub>, 2001–2020; (<b>d</b>) WUE<sub>min</sub>, 2001–2020.</p>
Full article ">Figure 5
<p>Regional statistics of the dominant factor of extreme variations in vegetation WUE in China. (<b>a</b>) WUE<sub>max</sub>, 2001–2010; (<b>c</b>) WUE<sub>max</sub>, 2011–2020; (<b>b</b>) WUE<sub>min</sub>, 2001–2010; (<b>d</b>) WUE<sub>min</sub>, 2011–2020.</p>
Full article ">Figure 6
<p>The spatial distribution of the relative contributions of (<b>a</b>) Prec, (<b>b</b>) Srad, (<b>c</b>) Temp, and (<b>d</b>) VPD to interannual variations in WUE.</p>
Full article ">Figure 7
<p>(<b>a</b>) Relative contributions rate and (<b>b</b>) absolute contributions of changes in Prec, Srad, Temp, and VPD to WUE for sub-regions.</p>
Full article ">Figure 8
<p>The spatial distribution of the absolute contributions of (<b>a</b>) Prec, (<b>b</b>) Srad, (<b>c</b>) Temp, and (<b>d</b>) VPD to interannual variations in WUE.</p>
Full article ">
14 pages, 2028 KiB  
Article
Metabolically Healthy Obesity Is Characterized by a Distinct Proteome Signature
by Fayaz Ahmad Mir, Houari B. Abdesselem, Farhan Cyprian, Ahmad Iskandarani, Asmma Doudin, Mutasem AbdelRahim Shraim, Bader M. Alkhalaf, Meis Alkasem, Ibrahem Abdalhakam, Ilham Bensmail, Hamza A. Al Halabi, Shahrad Taheri and Abdul-Badi Abou-Samra
Int. J. Mol. Sci. 2025, 26(5), 2262; https://doi.org/10.3390/ijms26052262 - 4 Mar 2025
Viewed by 99
Abstract
Obesity is commonly associated with metabolic diseases including type 2 diabetes, hypertension, and dyslipidemia. Moreover, individuals with obesity are at increased risk of cardiovascular disease. However, a subgroup of individuals within the obese population presents without concurrent metabolic disorders. Even though this group [...] Read more.
Obesity is commonly associated with metabolic diseases including type 2 diabetes, hypertension, and dyslipidemia. Moreover, individuals with obesity are at increased risk of cardiovascular disease. However, a subgroup of individuals within the obese population presents without concurrent metabolic disorders. Even though this group has a stable metabolic status and does not exhibit overt metabolic disease, this status may be transient; these individuals may have subclinical metabolic derangements. To investigate the latter hypothesis, an analysis of the proteome signature was conducted. Plasma samples from 27 subjects with obesity but without an associated metabolic disorder (obesity only (OBO)) and 15 lean healthy control (LHC) subjects were examined. Fasting samples were subjected to Olink proteomics analysis targeting 184 proteins enriched in cardiometabolic and inflammation pathways. Our results distinctly delineated two groups with distinct plasma protein expression profiles. Specifically, a total of 24 proteins were differentially expressed in individuals with obesity compared to LHC. Among these, 13 proteins were downregulated, whereas 11 proteins were upregulated. The pathways that were upregulated in the OBO group were related to chemoattractant activity, growth factor activity, G protein-coupled receptor binding, chemokine activity, and cytokine activity, whereas the pathways that were downregulated include regulation of T cell differentiation, leukocyte differentiation, reproductive system development, inflammatory response, neutrophil, lymphocyte, monocyte and leukocyte chemotaxis, and neutrophil migration. The study identifies several pathways that are altered in individuals with obesity compared to healthy control subjects. These findings provide valuable insights into the underlying mechanisms, potentially paving the way for the identification of therapeutic targets aimed at improving metabolic health in individuals with obesity. Full article
(This article belongs to the Special Issue Advance of Cell Metabolism in Endocrine Diseases)
Show Figures

Figure 1

Figure 1
<p>The experimental study designs. Abbreviations: BMI (body mass index), OBO (obesity only (no metabolic disease)), LHC (lean healthy controls).</p>
Full article ">Figure 2
<p>Differential protein expression analysis between OBO and LHC. Differentially expressed proteins (DEPs) were identified from combined Olink inflammation and cardiometabolic panels (184 unique proteins) defined as DEPs with more than a 1.25-fold change with FDR &lt; 0.05. (<b>A</b>) Hierarchical clustering based on all 184 proteins assayed using the two Olink panels showed a separation between obese patients compared to controls, with a cutoff Z score of 4. (<b>B</b>) A volcano plot summarizing DEPs based on log2 fold changes across the two groups. The red and green circles show significantly differential expressed proteins ≥ 0.32 log2 fold change (upregulated) or ≥−0.32 log2 fold change (downregulated) and FDR &lt; 0.05. Differential expression analysis accounted for obesity as the main effect, while correcting for age, sex, and race. (<b>C</b>) OBO signature protein list.</p>
Full article ">Figure 3
<p>Protein–protein interactions and pathway analysis of DEPs in OBO versus LHC. (<b>A</b>) The 24 DEPs were run in the STRING network to assess their protein–protein interactions. (<b>B</b>,<b>C</b>) Pathway analysis was performed using Gene Ontology (GO): (<b>A</b>) shows GO biological processes and (<b>B</b>) GO molecular functions. The size of the dots corresponds to adjusted <span class="html-italic">p</span>-values for statistically significant upregulation or downregulation represented by the red and green circles, respectively.</p>
Full article ">Figure 4
<p>A schematic drawing representing the potential outcome of the metabolic signature and associated signaling pathways altered in OBO individuals.</p>
Full article ">
40 pages, 656 KiB  
Article
The Impact of Digital–Green Synergy on Total Factor Productivity: Evidence from Chinese Listed Companies
by Dongfeng Chen, Junpeng Wang, Bin Li, Huihui Luo and Guangming Hou
Sustainability 2025, 17(5), 2200; https://doi.org/10.3390/su17052200 - 3 Mar 2025
Viewed by 206
Abstract
Driven by the dual imperatives of global economic green transformation and the advancement of digital technologies, achieving synergistic enhancement through digitalization and greenization to promote sustainable development has become a focal point for both academia and practical fields. This study, utilizing a sample [...] Read more.
Driven by the dual imperatives of global economic green transformation and the advancement of digital technologies, achieving synergistic enhancement through digitalization and greenization to promote sustainable development has become a focal point for both academia and practical fields. This study, utilizing a sample of Chinese A-share listed companies from 2010 to 2023, aims to explore the transformative potential of digital–green synergy (DGS) for enhancing enterprise sustainable development within the realm of production efficiency improvement. Employing a coupling coordination model based on the entropy-weighted TOPSIS method, the research measures the DGS levels of enterprises. Grounded in strategic synergy theory, the resource-based view, and dynamic capability theory, this study thoroughly investigates the direct impacts of DGS on corporate TFP, intermediary mechanisms, moderating effects, and heterogeneous roles. The research findings robustly demonstrate that DGS can significantly improve enterprise TFP through optimizing resource allocation, reducing cost stickiness, and enhancing operational efficiency, thereby facilitating the dynamic reorganization of production factors and the creation of sustainable value. Furthermore, external factors, such as financing constraints and environmental regulation, alongside internal organizational factors like executive characteristics, are shown to exert significant moderating effects on the effectiveness of DGS. In summary, this research not only highlights the crucial role of DGS in enhancing production efficiency as a driver for high-quality corporate development and the pursuit of sustainable goals but also provides important theoretical guidance for policymakers to incentivize digital and green transformation. It also offers practical insights for enterprise managers to strategically formulate synergistic development strategies, enhance economic benefits, and achieve long-term sustainable performance. Beyond these practical implications, this study further enriches the theoretical landscape by first extending strategic synergy theory to firm-level digital–green synergy in emerging markets; second by enhancing sustainability research by adopting a broader “environment-society” framework; methodologically innovating by developing a novel “goal-strategy-input-technology” synergy measurement framework; and finally, deepening the theoretical understanding of DGS-TFP relationships through mechanism and moderator exploration. Full article
(This article belongs to the Special Issue Sustainable Digital Transformation and Corporate Practices)
Show Figures

Figure 1

Figure 1
<p>PSM comparison: (<b>a</b>) Before matching (<b>b</b>) After matching.</p>
Full article ">
18 pages, 7897 KiB  
Article
The Dynamic Process of CO2 Leakage Along Wellbores Under Different Sequestration Conditions
by Baolei Zhu, Tianfu Xu, Xi Zhang, Chenglong Zhang and Guanhong Feng
Energies 2025, 18(5), 1237; https://doi.org/10.3390/en18051237 - 3 Mar 2025
Viewed by 158
Abstract
Abandoned production and monitoring wells in depleted oil and gas fields can readily serve as primary leakage pathways for stored CO2. The temperature, pressure conditions around the wellbore bottom, and CO2 concentration influence the phase behavior of CO2 during [...] Read more.
Abandoned production and monitoring wells in depleted oil and gas fields can readily serve as primary leakage pathways for stored CO2. The temperature, pressure conditions around the wellbore bottom, and CO2 concentration influence the phase behavior of CO2 during leakage. This study establishes a 3D wellbore–reservoir coupled model using CO2 injection data from 1 December 2009, in the DAS area, eastern Cranfield oilfield, Mississippi, USA, to analyze the dynamic evolution of CO2 leakage along wellbores. Simulations are conducted using the collaboration of ECO2M and ECO2N v2.0 modules. The study examines leakage regimes under varying distances from the injection well and different reservoir temperatures. The results indicate that CO2 phase changes occur primarily in wells near the injection point or under high-pressure and high CO2 saturation conditions, usually with a short leakage period due to ice formation at the wellhead. In areas with low CO2 saturation, prolonged leakage periods lead to significant pressure drops at the bottom, as well as the temperature as a result of the Joule–Thomson effect. Lower reservoir temperatures facilitate smoother and more gradual leakage. These findings provide a theoretical foundation for ensuring the safe implementation of CCUS projects and offer insights into the mechanical explanation of CO2 geyser phenomena. Full article
(This article belongs to the Section B3: Carbon Emission and Utilization)
Show Figures

Figure 1

Figure 1
<p>The reservoir schematic diagram of the DAS area from Jung et al. [<a href="#B24-energies-18-01237" class="html-bibr">24</a>].</p>
Full article ">Figure 2
<p>3D sketch of the wellbore–reservoir coupled model and the plane view.</p>
Full article ">Figure 3
<p>Porosity and permeability distribution in the reservoir of Cranfield; φ is porosity, and k is horizontal permeability (modified from Jung et al. (2020) [<a href="#B24-energies-18-01237" class="html-bibr">24</a>]).</p>
Full article ">Figure 4
<p>CO<sub>2</sub> injection rate from Pan et al. (2018) [<a href="#B30-energies-18-01237" class="html-bibr">30</a>] and simplified data for this simulation.</p>
Full article ">Figure 5
<p>Comparison between the simulated and measured pressure (<b>a</b>) and temperature (<b>b</b>) at the injection well bottom.</p>
Full article ">Figure 6
<p>The distribution of (<b>a</b>) pressure, (<b>b</b>) temperature, and (<b>c</b>) CO<sub>2</sub> saturation in the reservoir after 1552 days of injection (three black solid lines represent the monitoring wells).</p>
Full article ">Figure 7
<p>The profile evolution of (<b>a</b>) pressure, (<b>b</b>) temperature, (<b>c</b>) CO<sub>2</sub> saturation, and (<b>d</b>) mass fraction of CO<sub>2</sub> in the aqueous phase in the MW1 well.</p>
Full article ">Figure 8
<p>The profile evolution of (<b>a</b>) pressure, (<b>b</b>) temperature, (<b>c</b>) CO<sub>2</sub> saturation, and (<b>d</b>) mass fraction of CO<sub>2</sub> in the aqueous phase in the MW2 well.</p>
Full article ">Figure 9
<p>The profile evolution of (<b>a</b>) pressure, (<b>b</b>) temperature, (<b>c</b>) CO<sub>2</sub> saturation, and (<b>d</b>) mass fraction of CO<sub>2</sub> in the aqueous phase in the MW3 well.</p>
Full article ">Figure 10
<p>Pressure–temperature evolution at the well bottom of MW1 and MW2, with the background of Joule–Thomson coefficient distribution of CO<sub>2</sub>.</p>
Full article ">Figure 11
<p>Pressure–temperature profiles of different monitoring wells.</p>
Full article ">Figure 12
<p>The profile evolution of (<b>a</b>) pressure, (<b>b</b>) temperature, (<b>c</b>) CO<sub>2</sub> saturation, and (<b>d</b>) mass fraction of CO<sub>2</sub> in the aqueous phase in the MW1 well, under the condition that the ground temperature gradient is 16.2 °C/km.</p>
Full article ">Figure 13
<p>The profile evolution of Sg × Sl in the MW1 well; (<b>a</b>) base model and (<b>b</b>) geothermal gradient reduced by half to 16.2 °C/km.</p>
Full article ">Figure A1
<p>Primary variables in ECO2M and ECO2N v2.0 modules.</p>
Full article ">
27 pages, 4814 KiB  
Article
Evaluation of Urban Infrastructure Resilience Based on Risk–Resilience Coupling: A Case Study of Zhengzhou City
by Wenli Dong, Yunhan Zhou, Dongliang Guo, Zhehui Chen and Jiwu Wang
Land 2025, 14(3), 530; https://doi.org/10.3390/land14030530 - 3 Mar 2025
Viewed by 142
Abstract
The frequent occurrence of disasters has brought significant challenges to increasingly complex urban systems. Resilient city planning and construction has emerged as a new paradigm for dealing with the growing risks. Infrastructure systems like transportation, lifelines, flood control, and drainage are essential to [...] Read more.
The frequent occurrence of disasters has brought significant challenges to increasingly complex urban systems. Resilient city planning and construction has emerged as a new paradigm for dealing with the growing risks. Infrastructure systems like transportation, lifelines, flood control, and drainage are essential to the operation of a city during disasters. It is necessary to measure how risks affect these systems’ resilience at different spatial scales. This paper develops an infrastructure risk and resilience evaluation index system in city and urban areas based on resilience characteristics. Then, a comprehensive infrastructure resilience evaluation is established based on the risk–resilience coupling mechanism. The overall characteristics of comprehensive infrastructure resilience are then identified. The resilience transmission level and the causes of resilience effects are analyzed based on the principle of resilience scale. Additionally, infrastructure resilience enhancement strategies under different risk scenarios are proposed. In the empirical study of Zhengzhou City, comprehensive infrastructure resilience shows significant clustering in the city area. It is high in the central city and low in the periphery. Specifically, it is relatively high in the southern and northwestern parts of the airport economy zone (AEZ) and low in the center. The leading driving factors in urban areas are risk factors like flood and drought, hazardous materials, infectious diseases, and epidemics, while resilience factors include transportation networks, sponge city construction, municipal pipe networks, and fire protection. This study proposes a “risk-resilience” coupling framework to evaluate and analyze multi-hazard risks and the multi-system resilience of urban infrastructure across multi-level spatial scales. It provides an empirical resilience evaluation framework and enhancement strategies, complementing existing individual dimensional risk or resilience studies. The findings could offer visualized spatial results to support the decision-making in Zhengzhou’s resilient city planning outline and infrastructure special planning and provide references for resilience assessment and planning in similar cities. Full article
20 pages, 11785 KiB  
Article
IRFNet: Cognitive-Inspired Iterative Refinement Fusion Network for Camouflaged Object Detection
by Guohan Li, Jingxin Wang, Jianming Wei and Zhengyi Xu
Sensors 2025, 25(5), 1555; https://doi.org/10.3390/s25051555 - 3 Mar 2025
Viewed by 130
Abstract
Camouflaged Object Detection (COD) aims to identify objects that are intentionally concealed within their surroundings through appearance, texture, or pattern adaptations. Despite recent advances, extreme object–background similarity causes existing methods struggle with accurately capturing discriminative features and effectively modeling multiscale patterns while preserving [...] Read more.
Camouflaged Object Detection (COD) aims to identify objects that are intentionally concealed within their surroundings through appearance, texture, or pattern adaptations. Despite recent advances, extreme object–background similarity causes existing methods struggle with accurately capturing discriminative features and effectively modeling multiscale patterns while preserving fine details. To address these challenges, we propose Iterative Refinement Fusion Network (IRFNet), a novel framework that mimics human visual cognition through progressive feature enhancement and iterative optimization. Our approach incorporates the following: (1) a Hierarchical Feature Enhancement Module (HFEM) coupled with a dynamic channel-spatial attention mechanism, which enriches multiscale feature representations through bilateral and trilateral fusion pathways; and (2) a Context-guided Iterative Optimization Framework (CIOF) that combines transformer-based global context modeling with iterative refinement through dual-branch supervision. Extensive experiments on three challenging benchmark datasets (CAMO, COD10K, and NC4K) demonstrate that IRFNet consistently outperforms fourteen state-of-the-art methods, achieving improvements of 0.9–13.7% across key metrics. Comprehensive ablation studies validate the effectiveness of each proposed component and demonstrate how our iterative refinement strategy enables progressive improvement in detection accuracy. Full article
(This article belongs to the Special Issue Transformer Applications in Target Tracking)
Show Figures

Figure 1

Figure 1
<p>Overview of our proposed network architecture.</p>
Full article ">Figure 2
<p>Illustration of our proposed Hierarchical Feature Fusion (HFF) module.</p>
Full article ">Figure 3
<p>Illustration of our proposed Dynamic Channel-Spatial Attention (DCSA) module.</p>
Full article ">Figure 4
<p>Architecture of our Global Context Awareness Module (GCAM).</p>
Full article ">Figure 5
<p>Comparison of precision–recall curves on three benchmark datasets.</p>
Full article ">Figure 6
<p>Visual comparison with state-of-the-art methods. Our proposed IRFNet generates more accurate and complete segmentation maps across various challenging scenarios.</p>
Full article ">Figure 7
<p>Qualitative comparison of different model variants.</p>
Full article ">Figure 8
<p>Qualitative comparison of HFEM variants.</p>
Full article ">Figure 9
<p>Qualitative comparison of CIOF components.</p>
Full article ">Figure 10
<p>Evolution of performance metrics during the iterations in our optimal three-iteration model on the CAMO, COD10K, and NC4K datasets. The plots demonstrate consistent improvement across iterations, with MAE decreasing and other metrics increasing.</p>
Full article ">Figure 11
<p>Visualization of the iterative refinement process in our optimal three-iteration model. From left to right: input image, ground truth (GT), and prediction after first iteration (Iter 1), second iteration (Iter 2), and third iteration (Iter 3). The prediction quality progressively improves with each iteration.</p>
Full article ">Figure 12
<p>Visual comparison of polyp segmentation results.</p>
Full article ">
21 pages, 8306 KiB  
Article
Magmatic–Hydrothermal Processes of the Pulang Giant Porphyry Cu (–Mo–Au) Deposit, Western Yunnan: A Perspective from Different Generations of Titanite
by Mengmeng Li, Xue Gao, Guohui Gu and Sheng Guan
Minerals 2025, 15(3), 263; https://doi.org/10.3390/min15030263 - 3 Mar 2025
Viewed by 165
Abstract
The Yidun island arc was formed in response to the Late Triassic westward subduction of the Ganzi–Litang oceanic plate, a branch of the Paleo-Tethys Ocean. The Zhongdian arc, located in the south of the Yidun island arc, has relatively large number of porphyry [...] Read more.
The Yidun island arc was formed in response to the Late Triassic westward subduction of the Ganzi–Litang oceanic plate, a branch of the Paleo-Tethys Ocean. The Zhongdian arc, located in the south of the Yidun island arc, has relatively large number of porphyry (skarn) type Cu–Mo ± Au polymetallic deposits, the largest of which is the Pulang Cu (–Mo–Au) deposit with proven Cu reserves of 5.11 Mt, Au reserves of 113 t, and 0.17 Mt of molybdenum. However, the relationship between mineralization and the potassic alteration zone, phyllic zone, and propylitic zone of the Pulang porphyry deposit is still controversial and needs further study. Titanite (CaTiSiO5) is a common accessory mineral in acidic, intermediate, and alkaline igneous rocks. It is widely developed in various types of metamorphic rocks, hydrothermally altered rocks, and a few sedimentary rocks. It is a dominant Mo-bearing phase in igneous rocks and contains abundant rare earth elements and high-field-strength elements. As an effective geochronometer, thermobarometer, oxybarometer, and metallogenic potential indicator mineral, titanite is ideal to reveal the magmatic–hydrothermal evolution and the mechanism of metal enrichment and precipitation. In this paper, major and trace element contents of the titanite grains from different alteration zones were obtained using electron probe microanalysis (EPMA) and laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) to define the changes in physicochemical conditions and the behavior of these elements during the process of hydrothermal alteration at Pulang. Titanite in the potassic alteration zone is usually shaped like an envelope. It occurs discretely or is enclosed by feldspar, with lower contents of CaO, Al, Sr, Zr and Hf; a low Nb/Ta ratio; high ∑REE + Y, U, Th, Ta, Nb, and Ga content; and high FeO/Al2O3 and LREE/HREE ratios. This is consistent with the characteristics of magmatic titanite from fresh quartz monzonite porphyry in Pulang and other porphyry Cu deposits. Titanite in the potassium silicate alteration zone has more negative Eu anomaly and a higher U content and Th/U ratio, indicating that the oxygen fugacity decreased during the transformation to phyllic alteration and propylitic alteration in Pulang. High oxygen fugacity is favorable for the enrichment of copper, gold, and other metallogenic elements. Therefore, the enrichment of copper is more closely related to the potassium silicate alteration. The molybdenum content of titanite in the potassium silicate alteration zone is 102–104 times that of the phyllic alteration zone and propylitic alteration zone, while the copper content is indistinctive, indicating that molybdenum was dissolved into the fluid or deposited in the form of sulfide before the medium- to low-temperature hydrothermal alteration, which may lead to the further separation and deposition of copper and molybdenum. Full article
Show Figures

Figure 1

Figure 1
<p>Regional tectonic map of the Yidun Island arc collision orogenic belt. Modified from Deng et al., 2014 [<a href="#B4-minerals-15-00263" class="html-bibr">4</a>]. (<b>a</b>) The location map of the study area; (<b>b</b>) Geological schematic map of the Yidun Island arc; (<b>c</b>) Geological map of the southern Yidun arc (Zhongdian) region.</p>
Full article ">Figure 2
<p>Geological map (<b>a</b>) and AA’ profile (<b>b</b>) of the Pulang porphyry copper deposit, showing the orebodies and magmatic–hydrothermal alteration zones. Modified from Li et al., 2011 [<a href="#B8-minerals-15-00263" class="html-bibr">8</a>].</p>
Full article ">Figure 3
<p>Petrographical characteristics of the Pulang porphyry copper deposit showing the different alteration zones and mineral assemblages. (<b>a</b>) A hand specimen of quartz monzonite porphyry with potassium silicate alteration overprinted by quartz veins with obvious mineralization; (<b>b</b>) A hand specimen of quartz monzonite porphyry with phyllic alteration overlapping the potassium silicate alteration; (<b>c</b>) A hand specimen of quartz diorite porphyry with propylitic alteration and phyllic alteration; (<b>d</b>) A hand specimen of quartz monzonite porphyry with phyllic alteration; (<b>e</b>) A hand specimen of granite porphyry with propylitic alteration; (<b>f</b>) A hand specimen of quartz monzonite porphyry with potassium silicate alteration. Ccp-chalcopyrite; Mo-molybdenite; Po-pyrrhotite; Py-pyrite; Qtz-quartz.</p>
Full article ">Figure 4
<p>Photomicrograph of the three styles of titanites in different alteration zones. (<b>a</b>–<b>c</b>) are from the potassium silicate alteration zone; (<b>d</b>,<b>e</b>) are from the phyllic alteration zone; (<b>f</b>,<b>g</b>) are from the propylitic alteration zone; (<b>h</b>–<b>n</b>) are backscatter electron images correspond to a–g, respectively. The red spot is the electron probe point. Qz-quartz; Bt-biotite; Kfs-K-feldspar; Ttn-titanite; Scr-sericite; Rt-rutile; Ccp-chalcopyrite; Mt-magnetite; Py-pyrite; Po-pyrrhotite; llm-ilmenite.</p>
Full article ">Figure 5
<p>Variation range of major elements of titanite in different alteration zones.</p>
Full article ">Figure 6
<p>Variation range of trace elements and REEs of titanite in different alteration zones.</p>
Full article ">Figure 7
<p>Chondrite-normalized pattern of REEs. Normalization ratios from McDonough and Sun, 1995 [<a href="#B41-minerals-15-00263" class="html-bibr">41</a>].</p>
Full article ">Figure 8
<p>Al versus Fe (apfu) diagram of titanite in different alteration zones. Background map is from Aleinikoff et al., 2002 [<a href="#B43-minerals-15-00263" class="html-bibr">43</a>].</p>
Full article ">Figure 9
<p>The substitution mechanism of Fe<sub>2</sub>O<sub>3</sub> + Al<sub>2</sub>O<sub>3</sub> − TiO<sub>2</sub> (<b>a</b>), CaO − TiO<sub>2</sub> (<b>b</b>), F − TiO<sub>2</sub> (<b>c</b>), Na<sub>2</sub>O − CaO (<b>d</b>), and Ce<sub>2</sub>O<sub>3</sub> + Y<sub>2</sub>O<sub>3</sub> − CaO (<b>e</b>) in titanite from different alteration zones.</p>
Full article ">Figure 10
<p>Sr-(La/Yb)<sub>N</sub> (<b>a</b>), Sr-(Sm/Yb)<sub>N</sub> (<b>b</b>), and Sr-(La/Sm)<sub>N</sub> (<b>c</b>) diagrams of titanite in different alteration zones.</p>
Full article ">Figure 11
<p>CaO−Al<sub>2</sub>O<sub>3</sub> (<b>a</b>) and CaO−FeO/Al<sub>2</sub>O<sub>3</sub> (<b>b</b>) diagrams of titanite in different alteration zones.</p>
Full article ">Figure 12
<p>Cu-Mo separation mechanism of the Pulang copper deposit.</p>
Full article ">
13 pages, 1841 KiB  
Article
The N-Linked Glycosylation Asn191 and Asn199 Sites Are Controlled Differently Between PKA Signal Transduction and pEKR1/2 Activity in Equine Follicle-Stimulating Hormone Receptor
by Sung-Hoon Kim, Munkhzaya Byambaragchaa, Sei Hyen Park, Myung-Hum Park, Myung-Hwa Kang and Kwan-Sik Min
Curr. Issues Mol. Biol. 2025, 47(3), 168; https://doi.org/10.3390/cimb47030168 - 2 Mar 2025
Viewed by 77
Abstract
Equine follicle-stimulating hormone receptor (eFSHR) contains four extracellular N-linked glycosylation sites, which play important roles in agonist-induced signal transduction. Glycosylation regulates G protein-coupled receptor mechanisms by influencing folding, ligand binding, signaling, trafficking, and internalization. Here, we examined whether the glycosylated sites in eFSHR [...] Read more.
Equine follicle-stimulating hormone receptor (eFSHR) contains four extracellular N-linked glycosylation sites, which play important roles in agonist-induced signal transduction. Glycosylation regulates G protein-coupled receptor mechanisms by influencing folding, ligand binding, signaling, trafficking, and internalization. Here, we examined whether the glycosylated sites in eFSHR are necessary for cyclic adenosine monophosphate (cAMP) signal transduction and the phosphate extracellular signal-regulated kinase 1/2 (pERK1/2) response. We constructed mutants (N191Q, N199Q, N268Q, and N293Q) of the four N-linked glycosylation sites in eFSHR using site-directed mutagenesis. In wild-type (wt) eFSHR, the cAMP response gradually increased dose-dependently, displaying a strong response at the EC50 and Rmax. Two mutants (N191Q and N199Q) considerably decreased the cAMP response. Both EC50 values were approximately 0.46- and 0.44-fold compared to that of the eFSHR-wt, whereas Rmax levels were 0.29- and 0.45-fold compared to eFSHR-wt because of high-ligand treatment. Specifically, the EC50 and Rmax values in the N268Q mutant were increased 1.23- and 1.46-fold, respectively, by eFSHR-wt. pERK1/2 activity in eFSHR-wt cells was rapid, peaked within 5 min, consistently sustained until 15 min, and then sharply decreased. pERK1/2 activity in the N191Q mutant showed a pattern similar to that of the wild type, despite impaired cAMP responsiveness. The N199Q mutant showed low pERK1/2 activity at 5 and 15 min. Interestingly, pERK1/2 activity in the N268Q and N298Q mutants was similar to that of eFSHR-wt at 5 min, but neither mutant showed any signaling at 15 min, despite displaying high cAMP responsiveness. Overall, eFSHR N-linked glycosylation sites can signal to pERK1/2 via PKA and the other signals, dependent on G protein coupling and β-arrestin-dependent recruitment. Our results provide strong evidence for a new paradigm in which cAMP signaling is not activated, yet pERK1/2 cascade remains strongly induced. Full article
(This article belongs to the Special Issue Hormonal Regulation in Germ Cell Development)
Show Figures

Figure 1

Figure 1
<p>Schematic representation of the eFSHR structure. The N-linked glycosylation sites (N191, N199, N268, and N293) in the extracellular domain regions are indicated. The green circles indicate the putative four N-glycosylated sites. The extracellular region comprises 365 amino acids, indicating the longest extracellular region among G protein-coupled receptors. The three intracellular loops comprise 9, 25, and 25 amino acids. The intracellular regions have 65 amino acids, and the 10 potential phosphorylation sites (serine and threonine residues) are S641, T655, T657, S658, S659, T660, S664, T674, T684, and S689, indicated by yellow circles. EL, extracellular loop; IL, intracellular loop.</p>
Full article ">Figure 2
<p>Total cAMP levels induced by stimulation with the FSH in CHO-K1 cells transiently transfected with the N-linked glycosylation site mutants of eFSHR. Empty circles denote eFSHR-wt and black circles denote each mutant. The value of ΔF% was recalculated as cAMP concentration (nM). A representative dataset was obtained from three independent experiments. The figure depicts the results of a representative experiment performed with the indicated mutants.</p>
Full article ">Figure 3
<p>Rmax levels in the N-linked glycosylation mutants of eFSHR. The maximal cAMP responses presented in <a href="#cimb-47-00168-f002" class="html-fig">Figure 2</a> are displayed using a bar graph. * Statistically significant differences (<span class="html-italic">p</span> &lt; 0.05) compared to the Rmax level of the eFSHR wild type.</p>
Full article ">Figure 4
<p>pERK1/2 activation in HEK-293 cells transfected with eFSHR-wt and mutants after stimulation with FSH. HEK-293 cells transiently expressing eFSHR-wt or mutants were serum-starved for 4–6 h before stimulation with 500 ng/mL FSH for the indicated time periods. Whole-cell lysates were analyzed for pERK1/2 and total ERK levels. (<b>A</b>–<b>C</b>) Western blot results for phospho-ERK1/2 and total ERK. (<b>D</b>–<b>G</b>) Quantified phosphor-ERK1/2 levels, normalized to total ERK, are expressed as a percentage of the maximal response (100% for eFSHR-wt at 5 min). Densitometry was used to quantify the phospho-ERK1/2 band. Representative data are shown, and graphs represent the mean ± SEM of three independent experiments. Statistical significance was determined using one-way ANOVA, followed by Tukey’s comparison test. * <span class="html-italic">p</span> &lt; 0.05 compared with eFSHR-wt cells at the corresponding time point.</p>
Full article ">
29 pages, 8539 KiB  
Article
Three-Dimensional FEM Analysis of the Protective Effects of Isolation Piles on Tunnels Under Adjacent Excavations
by Libo Xu, Junneng Ye, Yanming Yao, Chi Liu and Xiaoli Liu
Appl. Sci. 2025, 15(5), 2673; https://doi.org/10.3390/app15052673 - 2 Mar 2025
Viewed by 261
Abstract
Isolation piles are critical for mitigating excavation-induced tunnel displacements, yet two unresolved challenges persist in tunnel engineering: (1) controversies regarding the influence of key parameters (e.g., pile head depth, pile length, and pile-to-pit distance) on their performance, and (2) insufficient understanding of the [...] Read more.
Isolation piles are critical for mitigating excavation-induced tunnel displacements, yet two unresolved challenges persist in tunnel engineering: (1) controversies regarding the influence of key parameters (e.g., pile head depth, pile length, and pile-to-pit distance) on their performance, and (2) insufficient understanding of the effects on both horizontal and vertical displacement control of tunnel. These challenges stem from the current research focus on isolated displacement components or simplified scenarios, which fails to address the complex interactions between key parameters and the deformation mechanisms. To address these gaps, this study proposes a hybrid validation framework integrating a three-dimensional finite element model (HS-Small constitutive model) with field monitoring data. A concept of “control efficiency” is introduced to quantify the effectiveness of isolation piles, complemented by a parametric sensitivity analysis framework. By synergizing the mirror image method and statistical theory, the research reveals a dual-path control mechanism involving displacement blocking and tunnel geometric reconfiguration. The findings advance the state of the art by resolving controversies over critical parameters and establishing a unified theoretical framework for coupled displacement control, providing actionable insights for optimizing isolation pile design in engineering practice. Full article
(This article belongs to the Special Issue New Challenges in Urban Underground Engineering)
Show Figures

Figure 1

Figure 1
<p>A workflow diagram of this study.</p>
Full article ">Figure 2
<p>Stress–strain relationship in a standard drained triaxial test.</p>
Full article ">Figure 3
<p>Schematic layout of the project.</p>
Full article ">Figure 4
<p>Profile of the relative positions of the foundation pit, isolation pile, and tunnel.</p>
Full article ">Figure 5
<p>Engineering geological cross-section.</p>
Full article ">Figure 6
<p>Three-dimensional finite element model: (<b>a</b>) model, (<b>b</b>) structure, and (<b>c</b>) mesh map.</p>
Full article ">Figure 7
<p>Evolutions of the horizontal displacement of the retaining wall (C1) and the depth.</p>
Full article ">Figure 8
<p>Evolutions of the horizontal displacement of the soil at a distance of 15 m away from the excavation (C2) and the depth.</p>
Full article ">Figure 9
<p>Evolutions of the horizontal displacement (<b>a</b>) and vertical displacement (<b>b</b>) of the tunnel with the monitoring point along the tunnel.</p>
Full article ">Figure 10
<p>Comparison between the predicted and measured values.</p>
Full article ">Figure 11
<p>Simplified finite-element model for parametric study.</p>
Full article ">Figure 12
<p>Three-dimensional finite-element model for parametric analysis.</p>
Full article ">Figure 13
<p>Relationships between control efficiency and pile position (x = 2 m) for varying pile lengths: 5 m (<b>a</b>), 10 m (<b>b</b>), 15 m (<b>c</b>), 20 m (<b>d</b>), 25 m (<b>e</b>), 30 m (<b>f</b>), and 35 m (<b>g</b>).</p>
Full article ">Figure 14
<p>Relationships between control efficiency and pile position (x = 5.5 m) for varying pile lengths: 5 m (<b>a</b>), 10 m (<b>b</b>), 15 m (<b>c</b>), 20 m (<b>d</b>), 25 m (<b>e</b>), 30 m (<b>f</b>), and 35 m (<b>g</b>).</p>
Full article ">Figure 15
<p>Relationships between control efficiency and pile position (x = 9 m) for varying pile lengths: 5 m (<b>a</b>), 10 m (<b>b</b>), 15 m (<b>c</b>), 20 m (<b>d</b>), 25 m (<b>e</b>), 30 m (<b>f</b>), and 35 m (<b>g</b>).</p>
Full article ">Figure 16
<p>Evolutions of the maximum horizontal displacement of soil at the isolation piles with the depth of the pile head for cases with x = 2 (<b>a</b>) and 9 m (<b>b</b>).</p>
Full article ">Figure 17
<p>Evolutions of vertical (<b>a</b>) and horizontal (<b>b</b>) displacement of the tunnel with the center angel (α).</p>
Full article ">Figure 18
<p>Schematic diagram of the tunnel deformation due to excavation with and without isolation piles.</p>
Full article ">Figure 19
<p>Relationships between control efficiency and pile head depth for cases with x = 2 (<b>a</b>), 5.5 (<b>b</b>), and 9 m (<b>c</b>).</p>
Full article ">Figure 20
<p>Relationships between control efficiency and pile length for cases with x = 2 (<b>a</b>), 5.5 (<b>b</b>), and 9 m (<b>c</b>).</p>
Full article ">Figure 21
<p>Relationships between the maximum horizontal displacement of soil at isolation piles and the position of the pile head and toe for cases with x = 2 m.</p>
Full article ">Figure 22
<p>Schematic diagram of the tunnel deformation due to excavation with various pile lengths.</p>
Full article ">Figure 23
<p>Relationships between vertical control efficiency and distance between the foundation pit and isolation piles for varying pile lengths: 10 m (<b>a</b>), 15 m (<b>b</b>), 20 m (<b>c</b>), 25 m (<b>d</b>), 30 m (<b>e</b>), 35 m (<b>f</b>), and 40 m (<b>g</b>).</p>
Full article ">Figure 24
<p>Relationships between horizontal control efficiency and distance between the foundation pit and isolation piles for varying pile lengths: 10 m (<b>a</b>), 15 m (<b>b</b>), 20 m (<b>c</b>), 25 m (<b>d</b>), 30 m (<b>e</b>), 35 m (<b>f</b>), and 40 m (<b>g</b>).</p>
Full article ">Figure 25
<p>Variation pattern of the control efficiency of isolation piles for a constant pile length (<b>a</b>) and pile toe depth (<b>b</b>).</p>
Full article ">
16 pages, 6479 KiB  
Article
Vat Photopolymerization of CeO2-Incorporated Hydrogel Scaffolds with Antimicrobial Efficacy
by Nelly Aimelyne Mpuhwe, Gyu-Nam Kim and Young-Hag Koh
Materials 2025, 18(5), 1125; https://doi.org/10.3390/ma18051125 - 2 Mar 2025
Viewed by 219
Abstract
We herein demonstrate the utility of gelatin methacryloyl (GelMA)/poly(ethylene glycol) diacrylate (PEGDA)–cerium oxide (CeO2) hydrogel inks for manufacturing hydrogel scaffolds with antimicrobial efficacy by vat photopolymerization. For uniform blending with GelMA/PEGDA hydrogels, CeO2 nanoparticles with a round shape were synthesized [...] Read more.
We herein demonstrate the utility of gelatin methacryloyl (GelMA)/poly(ethylene glycol) diacrylate (PEGDA)–cerium oxide (CeO2) hydrogel inks for manufacturing hydrogel scaffolds with antimicrobial efficacy by vat photopolymerization. For uniform blending with GelMA/PEGDA hydrogels, CeO2 nanoparticles with a round shape were synthesized by the precipitation method coupled with calculation at 600 °C. In addition, they had highly crystalline phases and the desired chemical structures (oxidation states of Ce3+ and Ce4+) required for outstanding antimicrobial efficacy. A range of GelMA/PEGDA-CeO2 hydrogel scaffolds with different CeO2 contents (0% w/v, 0.1% w/v, 0.5% w/v, 1% w/v, and 5% w/v with respect to distilled water content) were manufactured. The photopolymerization behavior, mechanical properties, and biological properties (swelling and biodegradation behaviors) of hydrogel scaffolds were characterized to optimize the CeO2 content. GelMA/PEGDA-CeO2 hydrogel scaffolds produced with the highest CeO2 content (5% w/v) showed reasonable mechanical properties (compressive strength = 0.56 ± 0.09 MPa and compressive modulus = 0.19 ± 0.03 MPa), a high swelling ratio (1063.3 ± 10.9%), and the desired biodegradation rate (remaining weight after 28 days = 39.6 ± 2.3%). Furthermore, they showed outstanding antimicrobial efficacy (the number of colony-forming units = 76 ± 44.6 (×103)). In addition, macroporous GelMA/PEGDA-CeO2 hydrogel scaffolds with tightly controlled porous structures could be manufactured by vat photopolymerization. Full article
(This article belongs to the Section Biomaterials)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram showing the manufacturing of CeO<sub>2</sub>-incorporated GelMA/PEGDA hydrogel scaffolds by VP-based 3D printing using photo-crosslinking process.</p>
Full article ">Figure 2
<p>Schematic diagram of the DLP 3D printing technique with GelMA-based hydrogel.</p>
Full article ">Figure 3
<p>Representative FE-SEM images of synthesized CeO<sub>2</sub> nanoparticles (<b>A</b>) at low and (<b>B</b>) at high magnifications.</p>
Full article ">Figure 4
<p>(<b>A</b>) XRD pattern and (<b>B</b>) XPS spectrum of synthesized CeO<sub>2</sub> nanoparticles.</p>
Full article ">Figure 5
<p>(<b>A</b>) Photo-DSC results and (<b>B</b>) cure depths obtained at 60 °C as a function of exposure time for photopolymerization from various GelMA/PEGDA-CeO<sub>2</sub> hydrogel inks (CeO<sub>2</sub> contents = 0% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 0.1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 0.5% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, and 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>).</p>
Full article ">Figure 6
<p>Representative FE-SEM images of GelMA/PEGDA-CeO<sub>2</sub> hydrogel inks with various CeO<sub>2</sub> contents of (<b>A</b>) 0% <span class="html-italic">w</span>/<span class="html-italic">v</span>, and (<b>B</b>) 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>, where white contrast represents CeO<sub>2</sub> nanoparticles. The yellow arrows indicate GelMA/PEGDA frameworks, and the inside marked by dashed line corresponds to pores.</p>
Full article ">Figure 7
<p>(<b>A</b>) Representative compressive stress–strain curves and (<b>B</b>) compressive strengths and moduli of various GelMA/PEGDA-CeO<sub>2</sub> hydrogel inks (CeO<sub>2</sub> contents = 0% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, and 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>). Different letters in each color represent statistical significance (<span class="html-italic">p</span>-value &lt; 0.05).</p>
Full article ">Figure 8
<p>Swelling ratios of various GelMA/PEGDA-CeO<sub>2</sub> hydrogel scaffolds (CeO<sub>2</sub> contents = 0.1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 0.5% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, and 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>). Same letters (a) represent no statistical significance (<span class="html-italic">p</span>-value &gt; 0.05).</p>
Full article ">Figure 9
<p>(<b>A</b>) Remaining weights of various GelMA/PEGDA-CeO<sub>2</sub> hydrogel scaffolds (CeO<sub>2</sub> contents = 0.1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 0.5% <span class="html-italic">w</span>/<span class="html-italic">v</span>, 1% <span class="html-italic">w</span>/<span class="html-italic">v</span>, and 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>) during biodegradation in PBS-collagenase solutions as a function of time and (<b>B</b>) remaining weights after 28 days. Different letters represent statistical significance (<span class="html-italic">p</span>-value &lt; 0.05).</p>
Full article ">Figure 10
<p>(<b>A</b>) Optical images of <span class="html-italic">S. mutans</span> adhered to GelMA/PEGDA and GelMA/PEGDA-CeO<sub>2</sub> hydrogel scaffolds (CeO<sub>2</sub> contents = 5% <span class="html-italic">w</span>/<span class="html-italic">v</span>) and (<b>B</b>) the numbers of colony-forming units.</p>
Full article ">Figure 11
<p>(<b>A</b>) A unit cell used to manufacture honeycomb-inspired scaffolds and (<b>B</b>) optical image of the GelMA/PEGDA-CeO<sub>2</sub> hydrogel scaffold manufactured by our VP printer. Letter ‘a’ represents the channel size, and letter ‘b’ represents wall thickness of honeycomb-inspired structure.</p>
Full article ">
16 pages, 16056 KiB  
Article
Spatial–Temporal Evolution and Coupling and Coordination of “Production–Life–Ecological” in Energy-Rich Area: A Perspective on Structure and Function
by Lin Zhang, Xingyue Ji, Yumeng Su and Zhaohua Lu
Land 2025, 14(3), 520; https://doi.org/10.3390/land14030520 - 1 Mar 2025
Viewed by 241
Abstract
The conflict between socio-economic development and ecological protection is prominent, as the practice framework for territorial spatial planning and the rational layout and function coordination of production–life–ecological (PLE) spaces are crucial for achieving regional sustainable development. However, the dynamic evolution of PLE structure [...] Read more.
The conflict between socio-economic development and ecological protection is prominent, as the practice framework for territorial spatial planning and the rational layout and function coordination of production–life–ecological (PLE) spaces are crucial for achieving regional sustainable development. However, the dynamic evolution of PLE structure and function, as well as the driving mechanisms for the sustainable development of PLE, are still understudied. Therefore, this study takes the Ji-shaped bend Energy-Rich Area (ERA) of the Yellow River basin as a case study, classifies the PLE spaces based on land use data, and develops a PLE function indicator system consistent with the regional characteristics of an ERA. This paper characterizes PLE from both structure and function perspectives and explores the coupling and coordinated degree (CCD) among PLE functions and their driving factors. The results show the following: (1) From 2000 to 2020, the area of living space increased by 35.86%, while areas of production and ecological space decreased by 2.10% and 0.08%, respectively. (2) The PLE function increased, with the production function performing better in the typical ERA and the ecological function performing well in the atypical ERA. (3) From 2000 to 2020, the CCD of the PLE function increased by 24.85%, with atypical ERA demonstrating a higher CCD than typical ERA. (4) Factors in production function had the most significant impact on the CCD of PLE function, followed by living drivers. These results provide valuable insights and guidance for regulating PLE and promoting sustainable development. Full article
(This article belongs to the Section Land Systems and Global Change)
Show Figures

Figure 1

Figure 1
<p>Spatial relationship diagram of production–living–ecological.</p>
Full article ">Figure 2
<p>Study area. (<b>a</b>) Study area location. (<b>b</b>) Digital Elevation Model (DEM). (<b>c</b>) Energy-rich areas.</p>
Full article ">Figure 3
<p>Spatial and temporal changes in production–living–ecological space.</p>
Full article ">Figure 4
<p>Spatial and temporal changes in production–living–ecological function.</p>
Full article ">Figure 5
<p>Characterization changes in CCD of production–living–ecological.</p>
Full article ">
Back to TopTop