Source author record

Cristóbal López

Cristóbal López appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2021arXiv

Sinking microplastics in the water column: simulations in the Mediterranean Sea

We study the vertical dispersion and distribution of negatively buoyant rigid microplastics within a realistic circulation model of the Mediterranean sea. We first propose an equation describing their idealized dynamics. In that framework, we evaluate the importance of some relevant physical effects: inertia, Coriolis force, small-scale turbulence and variable seawater density, and bound the relative error of simplifying the dynamics to a constant sinking velocity added to a large-scale velocity field. We then calculate the amount and vertical distribution of microplastic particles on the water column of the open ocean if their release from the sea surface is continuous at rates compatible with observations in the Mediterranean. The vertical distribution is found to be almost uniform with depth for the majority of our parameter range. Transient distributions from flash releases reveal a non-Gaussian character of the dispersion and various diffusion laws, both normal and anomalous. The origin of these behaviors is explored in terms of horizontal and vertical flow organization.

preprint2020arXiv

Accumulated densities of sedimenting particles in turbulent flows

We study the effect of turbulence on a sedimenting layer of particles by means of direct numerical simulations. A Lagrangian model in which particles are considered as tracers with an additional downward settling velocity is integrated together with an isotropic homogeneous turbulent flow. We study the spatial distribution of particles when they are collected on a plane at non-asymptotic times. We relate the resulting coarse-grained particle density to the history of the stretching rate along the particle trajectory and the projection of the density onto the accumulation plane, and analyse the deviation from homogeneity in terms of the Reynolds number and the settling velocity. We identify two regimes that arise during the early and during the well-mixed stage of advection. In the former regime, more inhomogeneity in the particle distribution is introduced for decreasing settling velocity or increasing Reynolds number, while the tendencies are opposite in the latter regime. A resonant-like crossover is found between these two regimes, where inhomogeneity is maximal.

preprint2020arXiv

Classical analogies for the force acting on an impurity in a Bose-Einstein condensate

We study the hydrodynamic forces acting on a small impurity moving in a two-dimensional Bose-Einstein condensate at non-zero temperature. The condensate is modelled by the damped-Gross Pitaevskii (dGPE) equation and the impurity by a Gaussian repulsive potential coupled to the condensate. For weak coupling, we obtain analytical expressions for the forces acting on the impurity, and compare them with those computed through direct numerical simulations of the dGPE and with the corresponding expressions for classical forces. For non-steady flows, there is a time-dependent force dominated by inertial effects and which has a correspondence in the Maxey-Riley theory for particles in classical fluids. In the steady-state regime, the force is dominated by a self-induced drag. Unlike at zero temperature, where the drag force vanishes below a critical velocity, at low temperatures the impurity experiences a net drag even at small velocities, as a consequence of the energy dissipation through interactions of the condensate with the thermal cloud. This dissipative force due to thermal drag is similar to the classical Stokes' drag. There is still a critical velocity above which steady-state drag is dominated by acoustic excitations and behaves non-monotonically with impurity's speed.

preprint2016arXiv

Correlation networks from flows. The case of forced and time-dependent advection-diffusion dynamics

Complex network theory provides an elegant and powerful framework to statistically investigate different types of systems such as society, brain or the structure of local and long-range dynamical interrelationships in the climate system. Network links in climate networks typically imply information, mass or energy exchange. However, the specific connection between oceanic or atmospheric flows and the climate network's structure is still unclear. We propose a theoretical approach for verifying relations between the correlation matrix and the climate network measures, generalizing previous studies and overcoming the restriction to stationary flows. Our methods are developed for correlations of a scalar quantity (temperature, for example) which satisfies an advection-diffusion dynamics in the presence of forcing and dissipation. Our approach reveals that correlation networks are not sensitive to steady sources and sinks and the profound impact of the signal decay rate on the network topology. We illustrate our results with calculations of degree and clustering for a meandering flow resembling a geophysical ocean jet.

preprint2016arXiv

Pattern formation with repulsive soft-core interactions: discrete particle dynamics and Dean-Kawasaki equation

Brownian particles interacting via repulsive soft-core potentials can spontaneously aggregate, despite repelling each other, and form periodic crystals of particle clusters. We study this phenomenon in low-dimensional situations (one and two dimensions) at two levels of description: performing numerical simulations of the discrete particle dynamics, and by linear and nonlinear analysis of the corresponding Dean-Kawasaki equation for the macroscopic particle density. Restricting to low dimensions and neglecting fluctuation effects we gain analytical insight into the mechanisms of the instability leading to clustering which turn out to be the interplay between diffusion, the intracluster forces and the forces between neighboring clusters. We show that the deterministic part of the Dean-Kawasaki equation provides a good description of the particle dynamics, including width and shape of the clusters, in a wide range of parameters, and analyze with weakly nonlinear techniques the nature of the pattern-forming bifurcation in one and two dimensions. Finally, we briefly discuss the case of attractive forces.

preprint2015arXiv

Boundaries of the Peruvian Oxygen Minimum Zone shaped by coherent mesoscale dynamics

Dissolved oxygen in sea water is a major factor affecting marine habitats and biogeochemical cycles. Oceanic zones with oxygen deficits represent significant portions of the area and volume of the oceans and are thought to be expanding. The Peruvian oxygen minimum zone is one of the most pronounced and lies in a region of strong mesoscale activity in the form of vortices and frontal regions, whose effect in the dynamics of the oxygen minimum zone is largely unknown. Here, we study this issue from a modeling approach and a Lagrangian point of view, using a coupled physical-biogeochemical simulation of the Peruvian oxygen minimum zone and finite-size Lyapunov exponent fields to understand the link between mesoscale dynamics and oxygen variations. Our results show that, at depths between 380 and 600 meters, mesoscale structures have a relevant dual role. First, their mean positions and paths delimit and maintain the oxygen minimum zone boundaries. Second, their high frequency fluctuations entrain oxygen across these boundaries as eddy fluxes that point towards the interior of the oxygen minimum zone and are one order of magnitude larger than mean fluxes. We conclude that these eddy fluxes contribute to the ventilation of the oxygen minimum zone.

preprint2015arXiv

Dominant transport pathways in an atmospheric blocking event

A Lagrangian flow network is constructed for the atmospheric blocking of eastern Europe and western Russia in summer 2010. We compute the most probable paths followed by fluid particles which reveal the {\it Omega}-block skeleton of the event. A hierarchy of sets of highly probable paths is introduced to describe transport pathways when the most probable path alone is not representative enough. These sets of paths have the shape of narrow coherent tubes flowing close to the most probable one. Thus, even when the most probable path is not very significant in terms of its probability, it still identifies the geometry of the transport pathways.

preprint2015arXiv

Pattern Formation in Populations with Density-Dependent Movement and Two Interaction Scales

We study the spatial patterns formed by a system of interacting particles where the mobility of any individual is determined by the population crowding at two different spatial scales. In this way we model the behavior of some biological organisms (like mussels) that tend to cluster at short ranges as a defensive strategy, and strongly disperse if there is a high population pressure at large ranges for optimizing foraging. We perform stochastic simulations of a particle-level model of the system, and derive and analyze a continuous density description (a nonlinear diffusion equation). In both cases we show that this interplay of scale-dependent-behaviors gives rise to a rich formation of spatial patterns ranging from labyrinths to periodic cluster arrangements. In most cases these clusters have the very peculiar appearance of ring-like structures, i.e., organisms arranging in the perimeter of the clusters, that we discuss in detail.

preprint2014arXiv

Hydrodynamic provinces and oceanic connectivity from a transport network help designing marine reserves

Oceanic dispersal and connectivity have been identified as crucial factors for structuring marine populations and designing Marine Protected Areas (MPAs). Focusing on larval dispersal by ocean currents, we propose an approach coupling Lagrangian transport and new tools from Network Theory to characterize marine connectivity in the Mediterranean basin. Larvae of different pelagic durations and seasons are modeled as passive tracers advected in a simulated oceanic surface flow from which a network of connected areas is constructed. Hydrodynamical provinces extracted from this network are delimited by frontiers which match multi-scale oceanographic features. By examining the repeated occurrence of such boundaries, we identify the spatial scales and geographic structures that would control larval dispersal across the entire seascape. Based on these hydrodynamical units, we study novel connectivity metrics for existing reserves. Our results are discussed in the context of ocean biogeography and MPAs design, having ecological and managerial implications.

preprint2013arXiv

Characterization of coherent structures in three-dimensional flows using the finite-size Lyapunov exponent

In this paper we use the finite size Lyapunov Exponent (FSLE) to characterize Lagrangian coherent structures in three-dimensional (3d) turbulent flows. Lagrangian coherent structures act as the organizers of transport in fluid flows and are crucial to understand their stirring and mixing properties. Generalized maxima (ridges) of the FSLE fields are used to locate these coherent structures. Three-dimensional FSLE fields are calculated in two phenomenologically distinct turbulent flows: a wall-bounded flow (channel flow) and a regional oceanic flow obtained by numerical solution of the primitive equations where two-dimensional turbulence dominates. In the channel flow, autocorrelations of the FSLE field show that the structure is substantially different from the near wall to the mid-channel region and relates well to the more widely studied Eulerian coherent structure of the turbulent channel flow. The ridges of the FSLE field have complex shapes due to the 3d character of the turbulent fluctuations. In the oceanic flow, strong horizontal stirring is present and the flow regime is similar to that of 2d turbulence where the domain is populated by coherent eddies that interact strongly. This in turn results in the presence of high FSLE lines throughout the domain leading to strong non-local mixing. The ridges of the FSLE field are quasi-vertical surfaces, indicating that the horizontal dynamics dominates the flow. Indeed, due to rotation and stratification, vertical motions in the ocean are much less intense than horizontal ones. This suppression is absent in the channel flow, as the 3d character of the FSLE ridges shows.

preprint2013arXiv

Clustering determines the survivor for competing Brownian and Lévy walkers

The competition between two ecologically similar species that use the same resources and differ from each other only in the type of spatial motion they undergo is studied. The latter is assumed to be described either by Brownian motion or Lévy flights. Competition is taken into account by assuming that individuals reproduce in a density-dependent fashion. It is observed that no influence of the type of motion occurs when the two species are in a well-mixed unstructured state. However, as soon as the species develop spatial clustering, the one forming more concentrated clusters gets a competitive advantage and eliminates the other. Similar competitive advantage would occur between walkers of the same type but with different diffusivities if this leads also to different clustering. Coexistence of both species is also possible under certain conditions.

preprint2013arXiv

Lagrangian transport in a microtidal coastal area: the Bay of Palma, island of Mallorca, Spain

Coastal transport in the Bay of Palma, a small region in the island of Mallorca, Spain, is characterized in terms of Lagrangian descriptors. The data sets used for this study are the output for two months (one in autumn and one in summer) of a high resolution numerical model, ROMS, forced atmospherically and with a spatial resolution of 300 m. The two months were selected because its different wind regime, which is the main driver of the sea dynamics in this area. Finite-size Lyapunov Exponents (FSLEs) were used to locate semi-persistent Lagrangian coherent structures (LCS) and to understand the different flow regimes in the Bay. The different wind directions and regularity in the two months have a clear impact on the surface Bay dynamics, whereas only topographic features appear clearly in the bottom structures. The fluid interchange between the Bay and the open ocean was tudied by computing particle trajectories and Residence Times (RT) maps. The escape rate of particles out of the Bay is qualitatively different, with a 32$%$ more of escape rate of particles to the ocean in October than in July, owing to the different geometric characteristics of the flow. We show that LCSs separate regions with different transport properties by displaying spatial distributions of residence times on synoptic Lagrangian maps together with the location of the LCSs. Correlations between the time-dependent behavior of FSLE and RT are also investigated, showing a negative dependence when the stirring characterized by FSLE values moves particles in the direction of escape.

preprint2013arXiv

Optimizing the search for resources by sharing information: Mongolian gazelles as a case study

We investigate the relationship between communication and search efficiency in a biological context by proposing a model of Brownian searchers with long-range pairwise interaction. After a general study of the properties of the model, we show an application to the particular case of acoustic communication among Mongolian gazelle, for which data are available, searching for good habitat areas. Using Monte Carlo simulations and density equations, our results point out that the search is optimal (i.e. the mean first hitting time among searchers is minimum) at intermediate scales of communication, showing that both an excess and a lack of information may worsen it.

preprint2013arXiv

The reduction of plankton biomass induced by mesoscale stirring: a modeling study in the Benguela upwelling

Recent studies, both based on remote sensed data and coupled models, showed a reduction of biological productivity due to vigorous horizontal stirring in upwelling areas. In order to better understand this phenomenon, we consider a system of oceanic flow from the Benguela area coupled with a simple biogeochemical model of Nutrient-Phyto-Zooplankton (NPZ) type. For the flow three different surface velocity fields are considered: one derived from satellite altimetry data, and the other two from a regional numerical model at two different spatial resolutions. We compute horizontal particle dispersion in terms of Lyapunov Exponents, and analyzed their correlations with phytoplankton concentrations. Our modelling approach confirms that in the south Benguela there is a reduction of biological activity when stirring is increased. Two-dimensional offshore advection and latitudinal difference in Primary Production, also mediated by the flow, seem to be the dominant processes involved. We estimate that mesoscale processes are responsible for 30 to 50% of the offshore fluxes of biological tracers. In the northern area, other factors not taken into account in our simulation are influencing the ecosystem. We suggest explanations for these results in the context of studies performed in other eastern boundary upwelling areas.

preprint2012arXiv

Oceanic three-dimensional Lagrangian Coherent Structures: A study of a mesoscale eddy in the Benguela ocean region

We study three dimensional oceanic Lagrangian Coherent Structures (LCSs) in the Benguela region, as obtained from an output of the ROMS model. To do that we first compute Finite-Size Lyapunov exponent (FSLE) fields in the region volume, characterizing mesoscale stirring and mixing. Average FSLE values show a general decreasing trend with depth, but there is a local maximum at about 100 m depth. LCSs are extracted as ridges of the calculated FSLE fields. They present a "curtain-like" geometry in which the strongest attracting and repelling structures appear as quasivertical surfaces. LCSs around a particular cyclonic eddy, pinched off from the upwelling front are also calculated. The LCSs are confirmed to provide pathways and barriers to transport in and out of the eddy.

preprint2010arXiv

How reliable are Finite-Size Lyapunov Exponents for the assessment of ocean dynamics?

Much of atmospheric and oceanic transport is associated with coherent structures. Lagrangian methods are emerging as optimal tools for their identification and analysis. An important Lagrangian technique which is starting to be widely used in oceanography is that of Finite-Size Lyapunov Exponents (FSLEs). Despite this growing relevance there are still many open questions concerning the reliability of the FSLEs in order to analyse the ocean dynamics. In particular, it is still unclear how robust they are when confronted with real data. In this paper we analyze the effect on this Lagrangian technique of the two most important effects when facing real data, namely noise and dynamics of unsolved scales. Our results, using as a benchmarch data from a primitive numerical model of the Mediterranean Sea, show that even when some dynamics is missed the FSLEs results still give an accurate picture of the oceanic transport properties.