Point Source Research Articles

Revealing asymmetry on mid-plane of protoplanetary disc through modelling of axisymmetric emission: methodology

ABSTRACT This study proposes an analytical framework for deriving the surface brightness profile and geometry of a geometrically thin axisymmetric disc from interferometric observation of continuum emission. Such precise modelling facilitates the exploration of faint non-axisymmetric structures, such as spirals and circumplanetary discs. As a demonstration, we simulate interferometric observations of geometrically thin axisymmetric discs. The proposed method can reasonably recover the injected axisymmetric structures, whereas Gaussian fitting of the same data yielded larger errors in disc orientation estimation. To further test the applicability of the method, it was applied to the mock data for $m=1,2$ spirals and a point source, which are embedded in a bright axisymmetric structure. The injected non-axisymmetric structures were reasonably recovered except for the innermost parts, and the disc geometric parameter estimations were better than Gasussian fitting. The method was then applied to the real data of Elias 20 and AS 209, and it adequately subtracted the axisymmetric component, notably in Elias 20, where substantial residuals remained without our method. We also applied our method to continuum data of PDS 70 to demonstrate the effectiveness of the method. We successfully recovered emission from PDS 70 c consistently with previous studies, and also tentatively discovered new substructures. The current formulation can be applied to any data for disc continuum emission, and aids in the search of spirals and circumplanetary discs, whose detection is still limited.

T-ReX: The Tarantula—Revealed by X-Rays

The Tarantula Nebula (30 Doradus, 30 Dor) is the most important star-forming complex in the Local Group, offering a microscope on starburst astrophysics. At its heart lies the exceptionally rich young stellar cluster R136, containing the most massive stars known. Stellar winds and supernovae have carved 30 Dor into an amazing display of arcs, pillars, and bubbles. We present first results and advanced data-processing products from the 2 Ms Chandra X-ray Visionary Project, “The Tarantula—Revealed by X-rays” (T-ReX). The 3615 point sources in the T-ReX catalog include massive stars, compact objects, binaries, bright pre-main-sequence stars, and compact young stellar (sub)clusters in 30 Dor. After removing point sources and excluding the exceptionally bright supernova remnant N157B (30 Dor B), the global diffuse X-ray maps reveal hot plasma structures resolved at 1–10 pc scales, with an absorption-corrected total-band (0.5–7 keV) X-ray luminosity of 2.110 × 1037 erg s−1. Spatially resolved spectral modeling provides evidence for emission lines enhanced by charge-exchange processes at the interfaces. We identify a candidate for the oldest X-ray pulsar detected to date in 30 Dor, PSR J0538-6902, inside a newly resolved arcuate X-ray wind nebula, the Manta Ray. The long temporal baseline of T-ReX allowed monitoring of dozens of massive stars, several showing periodic variability tied to binary orbital periods, and captured strong flares from at least three low-mass Galactic foreground stars.

Understanding Limitations in Electrochemical Conversion to CO at Low CO2 Concentrations.

Low-temperature electrochemical CO2 reduction has demonstrated high selectivity for CO when devices are operated with pure CO2 streams. However, there is currently a dearth of knowledge for systems operating below 30% CO2, a regime interesting for coupling electrochemical devices with CO2 point sources. Here we examine the influence of ionomer chemistry and cell operating conditions on the CO selectivity at low CO2 concentrations. Utilizing advanced electrochemical diagnostics, values for cathode catalyst layer ionic resistance and electrocatalyst capacitance as a function of relative humidity (RH) were extracted and correlated with selectivity and catalyst utilization. Staying above 20% CO2 concentration with at least a 50% cathode RH resulted in >95% CO/H2 selectivity regardless of the ionomer chemistry. At 10% CO2, however, >95% CO/H2 selectivity was only obtained at 95% RH under scenarios where the resulting electrode morphology enabled high catalyst utilization.

Optimal multi-period source-sink matching in carbon capture, utilization, and storage (CCUS) systems under neutrosophic environment

Carbon capture, utilization, and storage (CCUS) is an important technology for mitigating the effect of climate change. It involves the capture of CO2 from a point source, transport it through pipelines and store it in geological reservoirs for permanent storage. CCUS enables the use of fossil fuel-based power plant sources that produce a steady amount of energy and the reduction of CO2 emissions from these sources simultaneously. While CO2 sources are easy to characterize, the storage parameters of geological reservoirs are more difficult to measure, thus, estimates are provided and used for planning large-scale CCUS systems. CO2 utilization can provide a means to delay CO2 emissions by converting it highly valuable products, however, its sequestration potential is highly uncertain. A systematic mathematical approach is needed to effectively plan CCUS systems considering the risks associated with the estimates in the characteristics of the geological reservoir and the sequestration benefits from CO2 utilization. In this study, an optimization model is developed for planning multi-period CCUS systems under neutrosophic environment. The model considers managing storage risk by maximizing satisfaction from minimizing risks, minimizing dissatisfaction from consequences brought by these risks, and minimizing indeterminacy from inaccuracy of storage characteristics estimates. For utilization, the model considers the satisfaction for CO2 utilization benefits, dissatisfaction for opportunity losses, and indeterminacy in sequestration potential. Two case studies are used to illustrate the model where a sensitivity analysis is performed to illustrate the effect of the different risk appetite on the optimal solution. These cases studies reveal important insights on scheduling operations between CO2 sources and sinks and the allocation of captured CO2 for storage and utilization. The model can be used to develop policies for considering the risks associated with these uncertainties and how it will be managed.

Double wavefront reversal during six-wave interaction on Kerr nonlinearity in a waveguide with infinitely conducting surfaces

Background. Generation of a wave with a double reversed wavefront in multimode waveguides increases the efficiency of six-wave radiation converters and expands the possibilities of its use in adaptive optics problems and the conversion of complex spatially inhomogeneous waves. Aim. The quality of double wavefront reversal during six-wave interaction in a waveguide with infinitely conducting surfaces with Kerr nonlinearity is analyzed for the ratio of the wave numbers of the pump waves equal to 2 and 0,5, and the condition that one of the pump waves excites the zero mode of the waveguide, and the amplitude distribution of the other pump wave excites the edges of the waveguide are described by a Gaussian function. Methods. The influence of pump wave parameters on the half-width and contrast of the amplitude modulus of the object wave was studied using numerical methods. A wave from a point source located on the front face of the waveguide was used as a signal wave. Results. The dependences of the half-width and contrast of the amplitude modulus of the object wave on the ratio between the width of the waveguide and the width of the Gaussian pump wave are obtained. Conclusion. It is shown that the maximum change in the characteristics of a wave with a double reversed wavefront is observed when the width of the Gaussian pump waves changes in the range from 0,3 to 2 half-widths of the waveguide.

In-situ Measurement of Unsaturated Hydraulic Conductivity Function by a New Point Source Field Dripper Method (NPSFDM)

In-situ Measurement of Unsaturated Hydraulic Conductivity Function by a New Point Source Field Dripper Method (NPSFDM)

Comparison of patch-based and synoptic connectivity algorithms with graph theory metrics. A case study in Reggio Calabria metropolitan area (South Italy)

ContextPredicting and mapping connectivity between habitats and populations is critical to addressing habitat loss and biodiversity issues. Several strategies in the literature exist to understand, restore, and preserve ecological connectivity. The main issue of the current research is to identify which connectivity modeling strategies are the most reliable for planning purposes. ObjectivesOur goals in this paper were to compare connectivity predictions using a wide variety of commonly used approaches to improve the understanding of the similarities and differences in the predictions of these methods. Specifically, we investigated the differences in connectivity predictions related to the connectivity algorithm, the number and distribution of source points, and the threshold distance at which connectivity is allowed between locations. First, we separately applied different strategies and methods commonly used in the literature to model connectivity in the same study area. Then, going through a series of hypotheses, we compared the different models to confirm or disprove the initial hypotheses. MethodsWe proposed 4 main hypotheses and 14 combinations of them, hypothesizing that what most influences the results of connectivity models are different dispersal distance thresholds; differences in connectivity algorithms, especially kernel, path, and graph theory-based approaches; differences in predictions produced by two different software tools, UNICOR and Graphab; use of source points derived from a synoptic or patch-based perspective. ResultsWe found that the dominant pattern of differences in the predictions of different connectivity analyses was related to the method of analysis, with clear differences between kernel, path, and graph-theory approaches and relatively little effect due to the density and distribution of source points or the distance threshold used to define dispersal capability. ConclusionsThis work provides one of the first comparisons of spatial predictions of different methods, frameworks, and parameterizations of connectivity models. Our results support environmental planning by clarifying what most influences predictions of movement patterns and how the predicted connectivity networks differ between different analytical frameworks.

Ultra‐Inclined Nanocolumnar ZnO Films Sputtered Using a Novel Masking Configuration Providing Controlled and Restricted Oblique Angle Deposition for Enhanced Sensing Platforms

AbstractOblique angle deposition (OAD) of inclined thin films is mainly performed using electron beam evaporation due to its accurate point source control over the incoming evaporated flux angle α, leading to thin films with a nanocolumnar inclination angle β. However, the utilization of magnetron sputtering (MS) with an extended source for OAD is not extensively studied and reported. This work presents a thorough analysis of ZnO inclined thin films deposited by a novel restricted DC‐reactive MS‐OAD technique. OAD‐inclined films are deposited at α ranged 60°‐88°, where incoming flux is restricted using a patented masking configuration enabling tunable control of deposited nanocolumn angular range. The described technique provides accurate control over the resulting β (99.5% reproducibility), allowing demonstrated βmax of 47.3°, close to theoretical limits predicted for ZnO. The approach discussed here probes enhanced control of β comparable to that observed in evaporation, however using an extended source, resulting in high‐quality reproducible nanocolumnar‐inclined films. The mentioned improvements result from the exploration of operational parameters such as magnetron power, working pressure, and chamber temperature, as well as the design of the restricting configuration and substrate holders and their influence on the resulting inclined thin film crystallinity, and morphology.

Monitoring microplastics in the Seine River in the Greater Paris area

This data paper presents microplastic (MP) concentration and composition data from four campaigns in a one-year period in the Seine River within the Greater Paris area. Detailed data on MP occurrence, physical water parameters, hydrological conditions, and potential sewer point sources are presented. Initial analysis of the data set shows that MP concentrations were significant (median: 600 particles m−3 and 100 mg m−3) and MP fluxes showed variations across sampling sites and campaigns. Three polymer types (polypropylene - PP, polyethylene - PE, polystyrene - PS) account for 96% of the total number of MP. This dataset provides valuable insights for understanding river MP pollution under anthropogenic pressure.

Low mercury risks in paddy soils across the Pakistan

Mercury (Hg) is a globally distributed heavy metal. Here, we study Hg concentration and isotopic composition to understand the status of Hg pollution and its sources in Pakistan's paddy soil. The collected paddy soils (n = 500) across the country have an average THg concentration of 22.30 ± 21.74 ng/g. This low mean concentration suggests Hg pollution in Pakistan was not as severe as previously thought. Meanwhile, samples collected near brick kilns and industrial areas were significantly higher in THg than others, suggesting the influence of Hg emitted from point sources in certain areas. Soil physicochemical properties showed typical characteristic of mineral soils due to the study area's arid to semi-arid climate. Hg stable isotopes analysis, depicted mean Δ199Hg of −0.05 ± 0.12‰ and mean δ202Hg −0.45 ± 0.35‰, respectively, for contaminated sites, depicting Hg was primarily sourced from coal combustion by local anthropogenic sources. While uncontaminated sites show mean Δ199Hg of 0.15 ± 0.08‰, mean Δ200Hg of 0.06 ± 0.07‰ and mean δ202Hg of −0.32 ± 0.28‰, implying long-range transboundry Hg transport through wet Hg(II) deposition as a dominant Hg source. This study fills a significant knowledge gap regarding the Hg pollution status in Pakistan and suggests that the Hg risk in Pakistan paddies is generally low.

Estimating internal phosphorus loading for a water quality model using chemical characterisation of sediment phosphorus and contrasting oxygen conditions

The Finnish Archipelago Sea (AS) has long been subject to intensive anthropogenic phosphorus (P) loading. The area suffers from seasonal hypoxia and cyanobacterial blooms despite reductions in nutrient discharge from the catchment and point sources. Internal loading may even dominate the P budget. Previous estimates of internal P loading have limitations (e.g., in spatial coverage and infrequent measurements). We present the first area-wide estimates of the magnitude of internal P loading based on the long-term release of P stored in the sediments.Modelling the internal P loading in the AS is challenging due to the complexity of biogeochemical processes in the sediment-water interface, as well as the heterogenic topography of the seafloor. Instead, we calculated estimates of internal P loading based on data from previous studies on sequential chemical extraction of sediment P, sediment physical characteristics (e.g., organic content, location of muddy seabed substrates), and near-bottom oxygen (O2) conditions. The estimates in three scenarios of contrasting O2 conditions were based on potentially mobile P pools in the sediments, recycled from sediment to water (i.e., loosely-bound or exchangeable P, P bound to reducible iron oxy(hydr)oxides, and labile organic P). The potentially mobile P pools were determined by chemical extraction methods (modified from Psenner et al., 1984 and Ruttenberg, 1992). The internal P loading under presumable O2 conditions was estimated to be fivefold that of waterborne P input to the AS; comparable to previous estimates for hypoxic areas in the Baltic Sea.Our estimates revealed wide spatial variability in the internal P loading, depending on O2 conditions and seabed sediment substrate. The site-specific P release estimates are included in a water quality model used by regional authorities, which increases the model's reliability for estimating the impact of human activities on the water quality across the AS.

Controlled Discovery and Localization of Signals via Bayesian Linear Programming

Scientists often must simultaneously localize and discover signals. For instance, in genetic fine-mapping, high correlations between nearby genetic variants make it hard to identify the exact locations of causal variants. So the statistical task is to output as many disjoint regions containing a signal as possible, each as small as possible, while controlling false positives. Similar problems arise, for example, when locating stars in astronomical surveys and in changepoint detection. Common Bayesian approaches to these problems involve computing a posterior distribution over signal locations. However, existing procedures to translate these posteriors into credible regions for the signals fail to capture all the information in the posterior, leading to lower power and (sometimes) inflated false discoveries. We introduce Bayesian Linear Programming (BLiP), which can efficiently convert any posterior distribution over signals into credible regions for signals. BLiP overcomes an extremely high-dimensional and nonconvex problem to verifiably nearly maximize expected power while controlling false positives. Applying BLiP to existing state-of-the-art analyses of UK Biobank data (for genetic fine-mapping) and the Sloan Digital Sky Survey (for astronomical point source detection) increased power by 30%–120% in just a few minutes of additional computation. BLiP is implemented in pyblip (Python) and blipr (R). Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work.

A Note on Exact Closed-Form Solutions for Lamb’s Problem with an Arbitrary Poisson’s Ratio

Abstract Lamb’s Problem, a cornerstone in theoretical seismology, investigates the displacement of an elastic half-space due to surficial or buried point sources. Although the exact closed-form solution for surface displacements from surficial sources is well established, the solutions for buried sources in the 3D case have remained limited, either confined to integral forms or applicable only for certain Poisson’s ratios. This study addresses this gap by generalizing the existing closed-form solutions to accommodate Poisson’s ratio for buried sources. Our solutions are expressed in terms of elementary algebraic expressions, elliptic integrals, and Gauss hypergeometric functions. The derived solution is validated by numerical tests. With the newly derived solution, we also explore the properties of the leaking mode, which stems from the complex conjugate roots of the Rayleigh equation.

Probabilistic and progressive deblended far-infrared and sub-millimetre point source catalogues. I. Methodology and first application in the COSMOS field

Probabilistic and progressive deblended far-infrared and sub-millimetre point source catalogues. I. Methodology and first application in the COSMOS field

A study on decompression wave propagation characteristics during CO2 pipeline leakage with consideration of gas‐liquid transition

AbstractPipelines stand as the most cost‐effective method for large‐scale transportation of CO2 from a source point to the storage site, especially over extensive distances. The potential for crack propagation following a pipeline rupture highlights the need for precise analysis of decompression wave propagation. To accurately model this, understanding the decompression wave's propagation laws becomes imperative. Although previous studies have predominantly focused on pipeline leaks within the dense phase or supercritical state, the transition from liquid to gas during leakage significantly affects the decompression wave propagation. When a gaseous CO2 pipeline ruptures, the high Joule‐–Thomson coefficient causes a swift temperature plunge, potentially leading to a gas–liquid transition. However, research on how this phase transition impacts the decompression wave characteristics is limited. To address this gap, this study proposes a transition computational fluid dynamics model to predict the decompression wave behavior. The model is validated with an industrial‐scale full‐bore rupture experiment. The results reveal that the gaseous CO2 leakage induces a pressure plateau at a certain distance from the leakage due to the gas‐liquid phase transition. The influences of initial conditions on this pressure plateau and decompression wave are also explored. This study provides valuable insights into understanding the decompression wave behaviors of gaseous CO2 pipelines, which are essential for ensuring the safety and reliability of CO2 transportation within the carbon capture and storage technology chain. © 2024 Society of Chemical Industry and John Wiley & Sons, Ltd.

ON THE PROPERTIES OF A FAST ALGORITHM FOR SOLVING A THREE-DIMENSIONAL INVERSE PROBLEM OF SCALAR ACOUSTICS

t Theinverse single-frequency problem of scalar acoustics in three-dimensional space is considered. It consists in determining the characteristics of acoustic inhomogeneities lying in a flat layer from the distribution of the complex amplitude of the acoustic field in a flat recording layer. To effectively solve it, a special numerical algorithm is used, previously proposed and justified by the authors. As its component part, the algorithm uses solutions to special one-dimensional Fredholm equations of the 1st kind, and in the case of ambiguity in the solution of the latter, their normal solutions (solutions with a minimum norm) are sought. To regularize this ill-posed inverse problem, Tikhonov regularization and the TSVD method are used. A systematic numerical study is carried out of the influence of various parameters in the data recording scheme on the accuracy of the approximate solutions obtained using the algorithm. In particular, we study the dependence of this accuracy on the position of point sources and composite sources that cause sound vibrations, its dependence on the position of the data recording layer, and on the number of planes in which the recording sensors lie. It is shown that acceptable accuracy of approximate solutions can be obtained even with two such layers. In addition, an approach to the numerical study of the ambiguity of the solution to the inverse problem under consideration is proposed. It is based on special theoretical statements given in the article, and on the same numerical algorithm for solving the inverse problem. The approach is focused on comparing several normal solutions of one-dimensional integral equations from the problem under consideration. These normal solutions are calculated with respect to various elements. An example of using this theory to numerically estimate the non-uniqueness of a solution is given.

Extended Aquifer System Pressure Behavior under Carbon Storage

Summary Most reported carbon storage projects have involved inexpensive carbon dioxide (CO2) capture from gas processing plants or ethanol refineries. However, widespread carbon capture and storage (CCS) application must avoid any risk that high capital investment cost for carbon capture from stationary point sources leads to unanticipated issues related to the aquifer storage. This paper reviews successful and unsuccessful carbon storage projects and explains simple extended aquifer system fundamentals that must be considered in selecting a storage aquifer. This study begins by evaluating reported carbon storage projects in the context of an extended aquifer system with specific attention to initial formation pore pressure and potential or known hydraulic vertical or lateral communication with hydrocarbon accumulations and/or fresh water. Further study focuses on how the contrast between injection well and aquifer pressure evolution enables understanding of the overall aquifer material balance. Finally, we consider implications of brine migration during and after long-term CO2 injection in unconfined aquifers. Experience in the petroleum industry with aquifer behavior includes presence or lack of water influx and production from hydrocarbon reservoirs that share a common aquifer. Of particular importance is the observation that hydrostatic initial formation pressure indicates the possibility that a petroleum system, or an extended aquifer system without hydrocarbon accumulation(s), connects to atmospheric pressure through an unconfined aquifer. In such cases, indefinite injection will never increase the regional aquifer pressure. Furthermore, initial formation pressure that exceeds hydrostatic pressure implies a petroleum system or an extended aquifer system that is volumetrically limited. In such cases, injection will increase the system pressure, and pressure monitoring can detect leakage from the system. Finally, CO2 injection into an aquifer will displace brine in the direction of lower pressure that could relate to distant production from the same aquifer or from hydrocarbon reservoirs with which it communicates. Reasons for known carbon storage project interruptions have included unexpected lateral plume migration or aquifer pressure increase during CO2 injection that might have been anticipated with attention to straightforward consideration of aquifer-enabled hydraulic communication. Such extended aquifer dynamics must be included in long-term models for permanent CO2 storage during and after injection.

Periodicity from X-ray sources within the inner Galactic disk

Aims. For many years it had been claimed that the Galactic ridge X-ray emission at the Galactic Center (GC) is truly diffuse in nature. However, with the advancement of modern X-ray satellites, it has been found that most of the diffuse emission actually comprises thousands of previously unresolved X-ray point sources. Furthermore, many studies suggest that a vast majority of these X-ray point sources are magnetic cataclysmic variables (CVs) and active binaries. One unambiguous way to identify these magnetic CVs and other sources is by detecting their X-ray periodicity. Therefore, we systematically searched for periodic X-ray sources in the inner Galactic disk, including the GC region. Methods. We used data from our ongoing XMM-Newton Heritage Survey of the inner Galactic disk (350° ≲l ≲ +7° and −1° ≲b ≲ +1°) plus archival XMM-Newton observations of the GC. We computed the Lomb-Scargle periodogram for the soft (0.2–2 keV), hard (2–10 keV), and total (0.2–10 keV) band light curves to search for periodicities. Furthermore, we modeled the power spectrum using a power-law model to simulate 1000 artificial light curves and estimate the detection significance of the periodicity. We fitted the energy spectra of the sources using a simple power-law model plus three Gaussians, at 6.4, 6.7, and 6.9 keV, for the iron K emission complex. Results. We detected periodicity in 26 sources. For 14 of them, this is the first discovery of periodicity. For the other 12 sources, we found periods similar to those already known, indicating no significant period evolution. The intermediate polar (IP) type sources display relatively hard spectra compared to polars. We also searched for the Gaia counterparts of the periodic sources to estimate their distances using the Gaia parallax. We found a likely Gaia counterpart for seven sources. Conclusions. Based on the periodicity, hardness ratio, and the equivalent width of Fe K line emission, we have classified the sources into four categories: IPs, polars, neutron star X-ray binaries, and unknown. Of the 14 sources for which we detect the periodicity for the first time, four are likely IPs, five are likely polars, two are neutron star X-ray binaries, and three are of an unknown nature.

Quantifying water effluent violations and enforcement impacts using causal AI

AbstractIn the landscape of environmental governance, controlling water pollution through the regulation of point sources is vital as it preserves ecosystems, protects human health, ensures legal compliance, and fulfills global environmental responsibilities. Under the Clean Water Act, the integrated compliance information system monitors the compliance and enforcement status of facilities regulated by the National Pollutant Discharge Elimination System (NPDES) permit program. This study assesses temporal and geographic trends for effluent violations within the United States and introduces a novel metric for quantifying violation trends at the facility level. Furthermore, we utilize a linear parametric approach for Conditional Average Treatment Effect (CATE) causal analysis to quantify the heterogeneous effects of EPA and state enforcement actions on effluent violation trends at facilities with NPDES permits. Our research reveals insights into national pollutant discharge trends, regional clustering of all pollutant violation types in Ohio ( Z‐score of 2.15), and priority pollutants in West Virginia ( Z‐score of 3.07). The trend metric identifies regulated facilities that struggle with severe and recurring violations. The causal model highlights variations in state compliance and enforcement effectiveness, underscoring the successful moderation of violation trends by states such as Montana and Maryland, among others.

Perspectives on total phosphorus response in rivers: Examining the influence of rainfall extremes and post-dry rainfall

Eutrophication is a significant environmental problem caused by nutrient loads from both point and non-point sources. Weather variables, particularly precipitation, affect the concentration of nutrients in water bodies, particularly those from non-point sources, in two contrasting ways. Heavy precipitation causes surface runoff which transports pollutants to rivers and increases nutrient concentration. Conversely, increased river flow can dilute the concentration, lowering it. This study investigates the impact of extreme precipitation, prolonged precipitation, and precipitation after a dry period on the total phosphorus concentration in the Moehne and Erft rivers in Germany, given the projected increase in frequency of extreme precipitation events and long drought periods due to climate change. The study comprises two parts: selecting extreme weather days from 2001 to 2021 and comparing observed Total Phosphorus concentrations with estimated concentrations derived from Generalized Additive Models and linear regression based on the discharge-concentration relationship. Changes in river TP concentration in response to continuous precipitation and precipitation after a dry period were also studied. Our results showed that during wet extreme and post-dry period rainfall events, TP concentration consistently surpassed expected values, underscoring the profound influence of intense rainfall on nutrient mobilization. However, we observed the impact of continuous rainfall to be non-unidirectional. Our work is distinguished by three key innovations: 1) addressing limitations in studying the effects of extreme weather on water quality due to limited temporal resolution, 2) incorporating both linear and non-linear modeling approaches for discharge-concentration relationships, and 3) performing a comprehensive analysis of temporal and spatial patterns of Total Phosphorus concentrations in response to varying rainfall patterns.