Search

XTE J1810-189 is a Low-Mass X-ray Binary transient system hosting a neutron star, which underwent a three-month-long outburst in 2020. In order to study its spectral evolution during this outburst, we analysed all the available observations performed by NICER, in the 1-10 keV energy band. Firstly, we fitted the spectra with a thermal Comptonisation model. Our analysis revealed the lack of a significant direct emission from a black-body-like component, therefore we calculated the optical depth of the Comptonising region, deriving an upper limit of 4.5, which suggests the presence of a moderately thick corona. We also attempted to fit the spectrum with an alternative model, i.e. a cold Comptonised emission from a disc and a direct thermal component from the neutron star, finding a similarly good fit. The source did not enter a full high luminosity/soft state throughout the outburst, with a photon index ranging from 1.7 to 2.2, and an average unabsorbed flux in the 1-10 keV band of 3.6x10^(-10) erg cm^(-2) s^(-1). We searched for the presence of Fe K-shell emission lines in the range 6.4-7 keV, significantly detecting a broad component only in a couple of observations. Finally, we conducted a time-resolved spectral analysis of the detected type-I X-ray burst, observed during the outburst, finding no evidence of a photospheric radius expansion. The type-I burst duration suggests a mix of H/He fuel., Comment: 11 pages, 8 figures, accepted for publication on MNRAS

The Accreting Millisecond X-ray Pulsar IGR J17591-2342 is a LMXB system that went in outburst on August 2018 and it was monitored by the NICER observatory and partially by other facilities. We aim to study how the spectral emission of this source evolved during the outburst, by exploiting the whole X-ray data repository of simultaneous observations. The continuum emission of the combined broad-band spectra is on average well described by an absorbed Comptonisation component scattering black-body-distributed photons peaking at (0.8+/-0.5) keV, by a moderately optically thick corona (tau=2.3+/-0.5) with temperature of (34+/-9) keV. A black-body component with temperature and radial size of (0.8+/-0.2) keV and (3.3+/-1.5) km respectively is required by some of the spectra and suggests that part of the central emission, possibly a fraction of the neutron star surface, is not efficiently scattered by the corona. The continuum at low energies is characterised by significant residuals suggesting the presence of an absorption edge of O VIII and of emission lines of Ne IX ions. Moreover, broad Fe I and Fe XXV K-alpha emission lines are detected at different times of the outburst, suggesting the presence of reflection in the system., Comment: 13 pages, 7 figures, accepted for publication on MNRAS

The neutron star low-mass X-ray binary SWIFT J1749.4-2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this manuscript we perform a spectral characterization of the system throughout its 2021, two-week-long outburst, analyzing 11 NICER observations and quasi-simultaneous XMM-Newton and NuSTAR single observations at the outburst peak. The broadband spectrum is well-modeled with a black body component with a temperature of $\sim$0.6 keV, most likely consistent with a hot spot on the neutron star surface, and a Comptonisation spectrum with power-law index $\Gamma \sim 1.9$, arising from a hot corona at $\sim$12 keV. No direct emission from the disc was found, possibly due to it being too cool. A high truncation radius for the disc, i.e., at $\sim$20--30 R$_{G}$ , was obtained from the analysis of the broadened profile of the Fe line in the reflection component. The significant detection of a blue-shifted Fe XXVI absorption line at $\sim$7 keV indicates weakly relativistic X-ray disc winds, which are typically absent in the hard state of X-ray binaries. By comparing the low flux observed during the outburst and the one expected in a conservative mass-transfer, we conclude that mass-transfer in the system is highly non-conservative, as also suggested by the wind detection. Finally, using the Nicer spectra alone, we followed the system while it was fading to quiescence. During the outburst decay, as the spectral shape hardened, the hot spot on the neutron star surface cooled down and shrank, a trend which could be consistent with the pure power-law spectrum observed during quiescence., Comment: 15 pages, 9 figure; accepted for publication on MNRAS

We present the pulsar timing analysis of the accreting millisecond X-ray pulsar SWIFT J1749.4-2807 monitored by NICER and XMM-Newton during its latest outburst after almost eleven years of quiescence. From the coherent timing analysis of the pulse profiles, we updated the orbital ephemerides of the system. Large phase jumps of the fundamental frequency phase of the signal are visible during the outburst, consistent with what was observed during the previous outburst. Moreover, we report on the marginally significant evidence for non-zero eccentricity ($e\simeq 4\times 10^{-5}$) obtained independently from the analysis of both the 2021 and 2010 outbursts and we discuss possible compatible scenarios. Long-term orbital evolution of SWIFT J1749.4-2807 suggests a fast expansion of both the NS projected semi-major axis $(x)$, and the orbital period $(P_{\rm orb})$, at a rate of $\dot{x}\simeq 2.6\times 10^{-13}\,\text{lt-s}\,\text{s}^{-1}$ and $\dot{P}_{\rm orb}\simeq 4 \times 10^{-10}\,\text{s}\,\text{s}^{-1}$, respectively. SWIFT J1749.4-2807 is the only accreting millisecond X-ray pulsar, so far, from which the orbital period derivative has been directly measured from appreciable changes on the observed orbital period. Finally, no significant secular deceleration of the spin frequency of the compact object is detected, which allowed us to set a constraint on the magnetic field strength at the polar caps of $B_{PC}<1.3\times 10^{8}~\text{G}$, in line with typical values reported for AMXPs., Comment: 13 pages, 3 figures, 1 table. Accepted for publication in MNRAS

The X-ray source 4U 1822-371 is an eclipsing low-mass X-ray binary and X-ray pulsar, hosting a NS that shows periodic pulsations in the X-ray band. The inclination angle of the system is so high that in principle, it should be hard to observe both the direct thermal emission of the central object and the reflection component of the spectrum because they are hidden by the outer edge of the accretion disc. Assuming that the source accretes at the Eddington limit, we analysed non-simultaneous XMM-Newton and NuSTAR observations and studied the average broadband spectrum, with the aim to investigate the presence of a reflection component. No such component has been observed before in a high-inclination source such as 4U 1822-371. We modelled the spectral emission of the source using two different reflection models, Diskline plus Pexriv and the self-consistent model RfxConv. In our analysis, we find significant evidence of a reflection component in the spectrum, in addition to two lines associated with neutral or mildly ionised iron. The continuum spectrum is well fitted by a saturated Comptonisation model and a thermal black-body component emitted by the accretion disc at a lower temperature. We updated the ephemeris, adding two new eclipse times to the most recent ephemeris reported in literature. In our proposed scenario, the source is accreting at the Eddington limit with an intrinsic luminosity of $10^{38}$ erg/s, while the observed luminosity is two orders of magnitude lower. Despite the high inclination, we find that a reflection component is required to fit residuals at the Fe line range and the hard excess observed in the spectrum. The best-fit value of the inner disc radius is still uncertain and model dependent. More observations are therefore needed to confirm these results, which can give important information on this enigmatic and peculiar source.

Low-mass X-ray binaries hosting a low-magnetised neutron star, which accretes matter via Roche-lobe overflow, are generally grouped in two classes, named Atoll and Z sources after the path described in their X-ray colour-colour diagrams. Scorpius X-1 is the brightest persistent low-mass X-ray binary known so far, and it is the prototype of the Z sources. We analysed the first NuSTAR observation of this source to study its spectral emission exploiting the high statistics data collected by this satellite. Examining the colour-colour diagram, the source was probably observed during the lower normal and flaring branches of its Z-track. We separated the data from the two branches in order to investigate the evolution of the source along the track. We fitted the 3-60 keV NuSTAR spectra using the same models for both the branches. We adopted two description for the continuum: in the first case we used a blackbody and a thermal Comptonisation with seed photons originating in the accretion disc; in the second one, we adopted a disc-blackbody and a Comptonisation with a blackbody-shaped spectrum of the incoming seed photons. A power-law fitting the high energy emission above 20 keV was also required in both cases. The two models provide the same physical scenario for the source in both the branches: a blackbody temperature between 0.8 and 1.5 keV, a disc-blackbody with temperature between 0.4 and 0.6 keV, and an optically thick Comptonising corona with optical depth between 6 and 10 and temperature about 3 keV. Furthermore, two lines related to the K$\alpha$ and K$\beta$ transitions of the He-like Fe XXV ions were detected at 6.6 keV and 7.8 keV, respectively. A hard tail modelled by a power law with a photon index between 2 and 3 was also required for both the models., Comment: 12 pages, 7 fugures, accepted for publication in A&A

The dipping source XB 1916-053 is a compact binary system with an orbital period of 50 min harboring a neutron star. Using ten new {\it Chandra} observations and one {\it Swift/XRT} observation, we are able to extend the baseline of the orbital ephemeris; this allows us to exclude some models that explain the dip arrival times. The Chandra observations provide a good plasma diagnostic of the ionized absorber and allow us to determine whether it is placed at the outer rim of the accretion disk or closer to the compact object. From the available observations we are able to obtain three new dip arrival times extending the baseline of the orbital ephemeris from 37 to 40 years. From the analysis of the dip arrival times we confirm an orbital period derivative of $\dot{P}=1.46(3) \times 10^{-11}$ s s$^{-1}$. We show that the $\dot{P}$ value and the luminosity values are compatible with a mass accretion rate lower than 10\% of the mass transfer rate. We show that the mass ratio $q=m_2/m_1$ of 0.048 explains the apsidal precession period and the nodal precession period. The observed absorption lines are associated with the presence of \ion{Ne}{x}, \ion{Mg}{xii}, \ion{Si}{xiv}, \ion{S}{xvi,} and \ion{Fe}{xxvi} ions. We observe a redshift in the absorption lines between $1.1 \times 10^{-3}$ and $1.3 \times 10^{-3}$. By interpreting it as gravitational redshift, as recently discussed in the literature, we find that the ionized absorber is placed at a distance of $10^8$ cm from the neutron star with a mass of 1.4 M$_{\odot}$ and has a hydrogen atom density greater than $10^{15}$ cm$^{-3}$. (Abstract abridged), Comment: 17 pages, 10 figures, accepted by Astronomy & Astrophysics

GX 9+9 (4U 1728-16) is a low mass X-ray binary (LMXB) source harboring a neutron star. Although it belongs to the subclass of the bright Atoll sources together with GX 9+1, GX 3+1, and GX 13+1, its broadband spectrum is poorly studied and apparently does not show reflection features in the spectrum. To constrain the continuum well and verify whether a relativistic smeared reflection component is present, we analyze the broadband spectrum of GX 9+9 using {\it BeppoSAX} and \textit{XMM-Newton} spectra covering the 0.3-40 keV energy band. We fit the spectrum adopting a model composed of a disk-blackbody plus a Comptonized component whose seed photons have a blackbody spectrum (Eastern Model). A statistically equivalent model is composed of a Comptonized component whose seed photons have a disk-blackbody distribution plus a blackbody that mimics a saturated Comptonization likely associated with a boundary layer (Western model). Other trials did not return a good fit. The spectrum of GX 9+9 was observed in a soft state and its luminosity is $2.3 \times 10^{37}$ erg s$^{-1}$ assuming a distance to the source of 5 kpc. In the Eastern Model scenario, we find the seed-photon temperature and electron temperature of the Comptonized component to be $1.14^{+0.10}_{-0.07}$ keV and $2.80^{+0.09}_{-0.04}$ keV, respectively, while the optical depth of the Comptonizing corona is $8.9\pm0.4$. The color temperature of the inner accretion disk is $0.86^{+0.08}_{-0.02}$ keV and $0.82 \pm 0.02$ keV for the {\it BeppoSAX} and \textit{XMM-Newton} spectrum, respectively. In the Western Model scenario, instead, we find that the seed-photon temperature is $0.87 \pm 0.07$ keV and $1.01 \pm 0.08$ keV for the {\it BeppoSAX} and \textit{XMM-Newton} spectrum, respectively. (Abridged abstract), Comment: 10 pages, 8 figures, accepted by Astronomy & Astrophysics

The X-ray transient eclipsing source MXB 1659-298 went in outburst in 1999 and 2015, respectively, during which it was observed by XMM-Newton, NuSTAR and Swift. Using these observations we studied the broadband spectrum of the source to constrain the continuum components and to verify the presence of a reflection component. We analysed the soft and hard state of the source, finding that the soft state can be modelled with a thermal component associated with the inner accretion disc plus a Comptonised component. A smeared reflection component and the presence of an ionised absorber are also requested in the best-fit model. On the other hand, the direct continuum emission in the hard state can be described by a Comptonised component with a temperature larger than 150 keV. Also in this case a reflection component and a ionised absorber are observed., Comment: Proceedings of the 12th INTEGRAL conference and 1st AHEAD Gamma-ray Workshop, Geneva (Switzerland), 11-15 February 2019, Ed. C. Ferrigno, E. Bozzo, P. von Balmoos; 4 pages, 1 figure, 2 tables, submitted to Mem. S.A.It

Context. Since the discovery of the first Accreting Millisecond X-ray Pulsar SAX J1808.4-3658 in 1998, the family of these sources kept growing on. Currently, it counts 22 members. All AMXPs are transients with usually very long quiescence periods, implying that mass accretion rate in these systems is quite low and not constant. Moreover, for at least three sources, a non-conservative evolution was also proposed. Aims. Our purpose is to study the long term averaged mass-accretion rates in all the Accreting Millisecond X-ray Pulsars discovered so far, to investigate a non-conservative mass-transfer scenario. Methods. We calculated the expected mass-transfer rate under the hypothesis of a conservative evolution based on their orbital periods and on the (minimum) mass of the secondary (as derived from the mass function), driven by gravitational radiation and/or magnetic braking. Using this theoretical mass-transfer, we determined the expected accretion luminosity of the systems. Thus, we achieved the lower limit to the distance of the sources by comparing the computed theoretical luminosity and the observed flux averaged over a time period of 20 years. Then, the lower limit to the distance of the sources has been compared to the value of the distance reported in literature to evaluate how reasonable is the hypothesis of a conservative mass-transfer. Results. Based on a sample of 18 sources, we found strong evidences of a non-conservative mass-transfer for five sources, for which the estimated distance lower limits are higher than their known distances. We also report hints for mass outflows in other six sources. The discrepancy can be fixed under the hypothesis of a non-conservative mass-transfer in which a fraction of the mass transferred onto the compact object is swept away from the system, likely due to the (rotating magnetic dipole) radiation pressure of the pulsar., Comment: 12 pages, 2 figures, accepted for publication in A&A

The source X 1822-371 is an eclipsing compact binary system with a period close to 5.57 hr and an orbital period derivative $\dot{P}_{\rm orb}$ of 1.51(7)$\times 10^{-10}$ s s$^{-1}$. The very large value of $\dot{P}_{\rm orb}$ is compatible with a super-Eddington mass transfer rate from the companion star, as suggested by X-ray and optical data. The XMM-Newton observation taken in 2017 allows us to update the orbital ephemeris and verify whether the orbital period derivative has been stable over the last 40 yr. We added to the X-ray eclipse arrival times from 1977 to 2008 two new values obtained from the RXTE and XMM-Newton observations performed in 2011 and 2017, respectively. We estimated the number of orbital cycles and the delays of our eclipse arrival times spanning 40 yr using as reference time the eclipse arrival time obtained from the Rossi-XTE observation taken in 1996. Fitting the delays with a quadratic model, we found an orbital period $P_{\rm orb}=5.57062957(20)$ hr and a $\dot{P}_{\rm orb}$ value of $1.475(54) \times 10^{-10}$ s s$^{-1}$. The addition of a cubic term to the model does not significantly improve the quality of the fit. We also determined a spin-period value of $P_{\rm spin}=0.5915669(4)$ s and its first derivative $\dot{P}_{\rm spin}= -2.595(11) \times 10^{-12}$ s s$^{-1}$. The obtained results confirm the scenario of a super-Eddington mass transfer rate; we also exclude a gravitational coupling between the orbit and the change in the oblateness of the companion star triggered by the nuclear luminosity of the companion star., Comment: 5 pages, 3 figures, accepted for publication in A&A Letters

Context: XB 1916-053 is a low mass X-ray binary system (LMXB) hosting a neutron star (NS) and showing periodic dips. The spectrum of the persistent emission was modeled with a blackbody component having a temperature between 1.31 and 1.67 keV and with a Comptonization component with an electron temperature of 9.4 keV and a photon index $\Gamma$ between 2.5 and 2.9. The presence of absorption features associated with highly ionized elements suggested the presence of partially ionized plasma in the system. Aims: In this work we performed a study of the spectrum of XB 1916-053, which aims to shed light on the nature of the seed photons that contribute to the Comptonization component. Methods: We analyzed three Suzaku observations of XB 1916-053: the first was performed in November 2006 and the others were carried out in October 2014. We extracted the persistent spectra from each observation and combined the spectra of the most recent observations, obtaining a single spectrum with a higher statistic. We also extracted and combined the spectra of the dips observed during the same observations. Results: On the basis of the available data statistics, we infer that the scenario in which the corona Comptonizes photons emitted both by the innermost region of the accretion disk and the NS surface is not statistically relevant with respect to the case in which only photons emitted by the NS surface are Comptonized. We find that the source is in a soft spectral state in all the analyzed observations. We detect the K$\alpha$ absorption lines of \ion{Fe}{xxv} and \ion{Fe}{xxvi}, which have already been reported in literature, and for the first time the K$\beta$ absorption lines of the same ions. We also detect an edge at 0.876 keV, which is consistent with a \ion{O}{viii} K absorption edge. (Abridged), Comment: 12 pages, 4 figures, accepted for publication in A&A

Most of the X-ray binary systems containing neutron stars classified as Atoll sources show two different spectral states, called soft and hard. Moreover, a large number of these systems show a reflection component relativistically smeared in their spectra, which gives information on the innermost region of the system. Our aim is to investigate the poorly studied broadband spectrum of the low mass X-ray binary system 4U 1702-429, which was recently analysed combining XMM-Newton and INTEGRAL data. The peculiar value of the reflection fraction brought us to analyse further broadband spectra of 4U 1702-429. We re-analysed the spectrum of the XMM-Newton/INTEGRAL observation of 4U 1702-429 in the 0.3-60 keV energy range and we extracted three 0.1-100 keV spectra of the source analysing three observations collected with the BeppoSAX satellite. We find that the XMM-Newton/INTEGRAL spectrum is well fitted using a model composed of a disc blackbody plus a Comptonised component and a smeared reflection component. We used the same spectral model for the BeppoSAX spectra, finding out that the addition of a smeared reflection component is statistically significant. The best-fit values of the parameters are compatible to each other for the BeppoSAX spectra. We find that the reflection fraction is $0.05^{+0.03}_{-0.01}$ for the XMM-Newton/INTEGRAL spectrum and between 0.15 and 0.4 for BeppoSAX ones. The relative reflection fraction and the ionisation parameter are incompatible between the XMM-Newton/INTEGRAL and the BeppoSAX observations and the characteristics of the Comptonising corona suggest that the source was in a soft state in the former observation and in a hard state in the latter., Comment: 10 pages, 14 figures, submitted to A&A

We discuss the spectral and timing properties of the accreting millisecond X-ray pulsar SWIFT J1756.9-2508 observed by XMM-Newton, NICER and NuSTAR during the X-ray outburst occurred in April 2018. The spectral properties of the source are consistent with a hard state dominated at high energies by a non-thermal power-law component with a cut-off at ~70 keV. No evidence of iron emission lines or reflection humps has been found. From the coherent timing analysis of the pulse profiles, we derived an updated set of orbital ephemerides. Combining the parameters measured from the three outbursts shown by the source in the last ~11 years, we investigated the secular evolution of the spin frequency and the orbital period. We estimated a neutron magnetic field of 3.1E+8 G < B_pc< 4.5E+8 G and measured an orbital period derivative of -4.1E-12 s/s < P_dot_orb < 7.1E-12 s/s. We also studied the energy dependence of the pulse profile by characterising the behaviour of the pulse fractional amplitude in the energy range 0.3-80 keV. These results are compared with those obtained from the previous outbursts of SWIFT J1756.9-2508 and other previously known accreting millisecond X-ray pulsars., Comment: 9 pages, 6 figures, 4 tables. Accepted for publication in MNRAS

The X-ray transient eclipsing source MXB 1659-298 went into outburst in 1999 and 2015. During these two outbursts the source was observed by XMM-Newton, nuSTAR, and Swift/XRT. Using these observations, we studied the broadband spectrum of the source to constrain the continuum components and to verify whether it had a reflection component, as is observed in other X-ray eclipsing transient sources. We combined the available spectra to study the soft and hard state of the source in the 0.45-55 keV energy range. We report a reflection component in the soft and hard state. The direct emission in the soft state can be modeled with a thermal component originating from the inner accretion disk plus a Comptonized component associated with an optically thick corona surrounding the neutron star. On the other hand, the direct emission in the hard state is described only by a Comptonized component with a temperature higher than 130 keV; this component is associated with an optically thin corona. We observed narrow absorption lines from highly ionized ions of oxygen, neon, and iron in the soft spectral state. We investigated where the narrow absorption lines form in the ionized absorber. The equivalent hydrogen column density associated with the absorber is close to $6 \times 10^{23}$ cm$^{-2}$ and $1.3 \times 10^{23}$ cm$^{-2}$ in the soft and hard state, respectively., Comment: 12 pages, 8 Figures. Accepted by Astronomy & Astrophysics

We report on the detection of X-ray pulsations at 2.1 ms from the known X-ray burster IGR J17379-3747 using XMM-Newton. The coherent signal shows a clear Doppler modulation from which we estimate an orbital period of ~1.9 hours and a projected semi-major axis of ~8 lt-ms. Taking into account the lack of eclipses (inclination angle of < 75 deg) and assuming a neutron star mass of 1.4 Msun, we estimated a minimum companion star of ~0.06 Msun. Considerations on the probability distribution of the binary inclination angle make less likely the hypothesis of a main-sequence companion star. On the other hand, the close correspondence with the orbital parameters of the accreting millisecond pulsar SAX J1808.4-3658 suggests the presence of a bloated brown dwarf. The energy spectrum of the source is well described by a soft disk black-body component (kT ~0.45 keV) plus a Comptonisation spectrum with photon index ~1.9. No sign of emission lines or reflection components is significantly detected. Finally, combining the source ephemerides estimated from the observed outbursts, we obtained a first constraint on the long-term orbital evolution of the order of dP_orb/dt = (-2.5 +/- 2.3)E-12 s/s., Comment: 5 pages, 2 figures, 1 table. Accepted for publication in A&A

During last years theoretical works shed new light and proposed new hypothesis on the mechanisms which regulate the time behaviour of biological populations in different natural systems. Despite of this, the role of environmental variables in ecological systems is still an open question. Filling this gap of knowledge is a crucial task for a deeper comprehension of the dynamics of biological populations in real ecosystems. In this work we study how the dynamics of food spoilage bacteria influences the sensory characteristics of fresh fish specimens. This topic is crucial for a better understanding of the role played by the bacterial growth on the organoleptic properties, and for the quality evaluation and risk assessment of food products. We therefore analyze the time behaviour, in fresh fish specimens, of sensory characteristics starting from the growth curves of two spoilage bacterial communities. The theoretical study, initially based on a deterministic model, exploits experimental temperature profiles. As a first step, a model of predictive microbiology is used to reproduce the experimental behaviour of the two bacterial populations. Afterwards, the theoretical bacterial growths are converted, through suitable differential equations, into "sensory" scores, based on the Quality Index Method (QIM), a scoring system for freshness and quality sensory estimation of fishery products. As a third step, the theoretical curves of QIM scores are compared with the experimental data obtained by sensory analysis. Finally, the differential equations for QIM scores are modified by adding terms of multiplicative white noise, which mimics the effects of uncertainty and variability in sensory analysis. A better agreement between experimental and theoretical QIM scores is observed, in some cases, in the presence of suitable values of noise intensity respect to the deterministic analysis., Comment: 25 pages, 12 figures, 1 table; to be published in Math. Model. Nat. Phenom

Noise through its interaction with the nonlinearity of the living systems can give rise to counter-intuitive phenomena. In this paper we shortly review noise induced effects in different ecosystems, in which two populations compete for the same resources. We also present new results on spatial patterns of two populations, while modeling real distributions of anchovies and sardines. The transient dynamics of these ecosystems are analyzed through generalized Lotka-Volterra equations in the presence of multiplicative noise, which models the interaction between the species and the environment. We find noise induced phenomena such as quasi-deterministic oscillations, stochastic resonance, noise delayed extinction, and noise induced pattern formation. In addition, our theoretical results are validated with experimental findings. Specifically the results, obtained by a coupled map lattice model, well reproduce the spatial distributions of anchovies and sardines, observed in a marine ecosystem. Moreover, the experimental dynamical behavior of two competing bacterial populations in a meat product and the probability distribution at long times of one of them are well reproduced by a stochastic microbial predictive model., Comment: 23 pages, 8 figures; to be published in Math. Model. Nat. Phenom. (2016)

In this paper, by using a stochastic reaction-diffusion-taxis model, we analyze the picophytoplankton dynamics in the basin of the Mediterranean Sea, characterized by poorly mixed waters. The model includes intraspecific competition of picophytoplankton for light and nutrients. The multiplicative noise sources present in the model account for random fluctuations of environmental variables. Phytoplankton distributions obtained from the model show a good agreement with experimental data sampled in two different sites of the Sicily Channel. The results could be extended to analyze data collected in different sites of the Mediterranean Sea and to devise predictive models for phytoplankton dynamics in oligotrophic waters., Comment: To be published in Ecological Complexity