Marshall John C.

Last Name

Marshall

First Name

John C.

ORCID

0000-0001-9230-3591

Full item page

Search Results

Now showing 1 - 18 of 18

Hydrothermal plume dynamics on Europa : implications for chaos formation

(American Geophysical Union, 2004-03-20) Goodman, Jason C. ; Collins, Geoffrey C. ; Marshall, John C. ; Pierrehumbert, Raymond T.

Hydrothermal plumes may be responsible for transmitting radiogenic or tidally generated heat from Europa's rocky interior through a liquid ocean to the base of its ice shell. This process has been implicated in the formation of chaos regions and lenticulae by melting or exciting convection in the ice layer. In contrast to earlier work, we argue that Europa's ocean should be treated as an unstratified fluid. We have adapted and expanded upon existing work describing buoyant plumes in a rotating, unstratified environment. We discuss the scaling laws governing the flow and geometry of plumes on Europa and perform a laboratory experiment to obtain scaling constants and to visualize plume behavior in a Europa-like parameter regime. We predict that hydrothermal plumes on Europa are of a lateral scale (at least 25–50 km) comparable to large chaos regions; they are too broad to be responsible for the formation of individual lenticulae. Plume heat fluxes (0.1–10 W/m2) are too weak to allow complete melt-through of the ice layer. Current speeds in the plume (3–8 mm/s) are much slower than indicated by previous studies. The observed movement of ice blocks in the Conamara Chaos region is unlikely to be driven by such weak flow.
Suppressed pCO(2) in the Southern Ocean due to the interaction between current and wind

(American Geophysical Union, 2021-11-15) Kwak, Kyungmin ; Song, Hajoon ; Marshall, John C. ; Seo, Hyodae ; McGillicuddy, Dennis J.

The Southern Ocean, an important region for the uptake of anthropogenic carbon dioxide (CO2), features strong surface currents due to substantial mesoscale meanders and eddies. These features interact with the wind and modify the momentum transfer from the atmosphere to the ocean. Although such interactions are known to reduce momentum transfer, their impact on air-sea carbon exchange remains unclear. Using a 1/20° physical-biogeochemical coupled ocean model, we examined the impact of the current-wind interaction on the surface carbon concentration and the air-sea carbon exchange in the Southern Ocean. The current-wind interaction decreased winter partial pressure of CO2 (pCO2) at the ocean surface mainly south of the northern subantarctic front. It also reduced pCO2 in summer, indicating enhanced uptake, but not to the same extent as the winter loss. Consequently, the net outgassing of CO2 was found to be reduced by approximately 17% when including current-wind interaction. These changes stem from the combined effect of vertical mixing and Ekman divergence. A budget analysis of dissolved inorganic carbon (DIC) revealed that a weakening of vertical mixing by current-wind interaction reduces the carbon supply from below, and particularly so in winter. The weaker wind stress additionally lowers the subsurface DIC concentration in summer, which can affect the vertical diffusive flux of carbon in winter. Our study suggests that ignoring current-wind interactions in the Southern Ocean can overestimate winter CO2 outgassing.
Direct estimate of lateral eddy diffusivity upstream of Drake Passage

(American Meteorological Society, 2014-10) Tulloch, Ross ; Ferrari, Raffaele ; Jahn, Oliver ; Klocker, Andreas ; LaCasce, Joseph H. ; Ledwell, James R. ; Marshall, John C. ; Messias, Marie-Jose ; Speer, Kevin G. ; Watson, Andrew J.

The first direct estimate of the rate at which geostrophic turbulence mixes tracers across the Antarctic Circumpolar Current is presented. The estimate is computed from the spreading of a tracer released upstream of Drake Passage as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The meridional eddy diffusivity, a measure of the rate at which the area of the tracer spreads along an isopycnal across the Antarctic Circumpolar Current, is 710 ± 260 m2 s−1 at 1500-m depth. The estimate is based on an extrapolation of the tracer-based diffusivity using output from numerical tracers released in a one-twentieth of a degree model simulation of the circulation and turbulence in the Drake Passage region. The model is shown to reproduce the observed spreading rate of the DIMES tracer and suggests that the meridional eddy diffusivity is weak in the upper kilometer of the water column with values below 500 m2 s−1 and peaks at the steering level, near 2 km, where the eddy phase speed is equal to the mean flow speed. These vertical variations are not captured by ocean models presently used for climate studies, but they significantly affect the ventilation of different water masses.
The CLIMODE field campaign : observing the cycle of convection and restratification over the Gulf Stream

(American Meteorological Society, 2009-09) Marshall, John C. ; Ferrari, Raffaele ; Forget, Gael ; Andersson, A. ; Bates, Nicholas R. ; Dewar, William K. ; Doney, Scott C. ; Fratantoni, David M. ; Joyce, Terrence M. ; Straneo, Fiamma ; Toole, John M. ; Weller, Robert A. ; Edson, James B. ; Gregg, M. C. ; Kelly, Kathryn A. ; Lozier, M. Susan ; Palter, Jaime B. ; Lumpkin, Rick ; Samelson, Roger M. ; Skyllingstad, Eric D. ; Silverthorne, Katherine E. ; Talley, Lynne D. ; Thomas, Leif N.

A major oceanographic field experiment is described, which is designed to observe, quantify, and understand the creation and dispersal of weakly stratified fluid known as “mode water” in the region of the Gulf Stream. Formed in the wintertime by convection driven by the most intense air–sea fluxes observed anywhere over the globe, the role of mode waters in the general circulation of the subtropical gyre and its biogeo-chemical cycles is also addressed. The experiment is known as the CLIVAR Mode Water Dynamic Experiment (CLIMODE). Here we review the scientific objectives of the experiment and present some preliminary results.
Impact of near-inertial waves on vertical mixing and air-sea CO2 fluxes in the Southern Ocean

(American Geophysical Union, 2019-06-17) Song, Hajoon ; Marshall, John C. ; Campin, Jean-Michel ; McGillicuddy, Dennis J.

We report the significant impact of near‐inertial waves (NIWs) on vertical mixing and air‐sea carbon dioxide (CO2) fluxes in the Southern Ocean using a biogeochemical model coupled to an eddy‐rich ocean circulation model. The effects of high‐frequency processes are quantified by comparing the fully coupled solution (ONLINE) to two offline simulations based on 5‐day‐averaged output of the ONLINE simulation: one that uses vertical mixing archived from the ONLINE model (CTRL) and another in which vertical mixing is recomputed from the 5‐day average hydrodynamic fields (5dAVG). In this latter simulation, processes with temporal variabilities of a few days including NIWs are excluded in the biogeochemical simulation. Suppression of these processes reduces vertical shear and vertical mixing in the upper ocean, leading to decreased supply of carbon‐rich water from below, less CO2 outgassing in austral winter, and more uptake in summer. The net change amounts up to one third of the seasonal variability in Southern Ocean CO2 flux. Our results clearly demonstrate the importance of resolving high‐frequency processes such as NIWs to better estimate the carbon cycle in numerical model simulations.
Evaluation of ocean carbon cycle models with data-based metrics

(American Geophysical Union, 2004-04-02) Matsumoto, K. ; Sarmiento, Jorge L. ; Key, Robert M. ; Aumont, Olivier ; Bullister, John L. ; Caldeira, Ken ; Campin, J.-M. ; Doney, Scott C. ; Drange, Helge ; Dutay, J.-C. ; Follows, Michael J. ; Gao, Y. ; Gnanadesikan, Anand ; Gruber, Nicolas ; Ishida, Akio ; Joos, Fortunat ; Lindsay, Keith ; Maier-Reimer, Ernst ; Marshall, John C. ; Matear, Richard J. ; Monfray, Patrick ; Mouchet, Anne ; Najjar, Raymond G. ; Plattner, Gian-Kasper ; Schlitzer, Reiner ; Slater, Richard D. ; Swathi, P. S. ; Totterdell, Ian J. ; Weirig, Marie-France ; Yamanaka, Yasuhiro ; Yool, Andrew ; Orr, James C.

New radiocarbon and chlorofluorocarbon-11 data from the World Ocean Circulation Experiment are used to assess a suite of 19 ocean carbon cycle models. We use the distributions and inventories of these tracers as quantitative metrics of model skill and find that only about a quarter of the suite is consistent with the new data-based metrics. This should serve as a warning bell to the larger community that not all is well with current generation of ocean carbon cycle models. At the same time, this highlights the danger in simply using the available models to represent the state-of-the-art modeling without considering the credibility of each model.
Evaluating global ocean carbon models : the importance of realistic physics

(American Geophysical Union, 2004-09-15) Doney, Scott C. ; Lindsay, Keith ; Caldeira, Ken ; Campin, J.-M. ; Drange, Helge ; Dutay, J.-C. ; Follows, Michael J. ; Gao, Y. ; Gnanadesikan, Anand ; Gruber, Nicolas ; Ishida, Akio ; Joos, Fortunat ; Madec, G. ; Maier-Reimer, Ernst ; Marshall, John C. ; Matear, Richard J. ; Monfray, Patrick ; Mouchet, Anne ; Najjar, Raymond G. ; Orr, James C. ; Plattner, Gian-Kasper ; Sarmiento, Jorge L. ; Schlitzer, Reiner ; Slater, Richard D. ; Totterdell, Ian J. ; Weirig, Marie-France ; Yamanaka, Yasuhiro ; Yool, Andrew

A suite of standard ocean hydrographic and circulation metrics are applied to the equilibrium physical solutions from 13 global carbon models participating in phase 2 of the Ocean Carbon-cycle Model Intercomparison Project (OCMIP-2). Model-data comparisons are presented for sea surface temperature and salinity, seasonal mixed layer depth, meridional heat and freshwater transport, 3-D hydrographic fields, and meridional overturning. Considerable variation exists among the OCMIP-2 simulations, with some of the solutions falling noticeably outside available observational constraints. For some cases, model-model and model-data differences can be related to variations in surface forcing, subgrid-scale parameterizations, and model architecture. These errors in the physical metrics point to significant problems in the underlying model representations of ocean transport and dynamics, problems that directly affect the OCMIP predicted ocean tracer and carbon cycle variables (e.g., air-sea CO2 flux, chlorofluorocarbon and anthropogenic CO2 uptake, and export production). A substantial fraction of the large model-model ranges in OCMIP-2 biogeochemical fields (±25–40%) represents the propagation of known errors in model physics. Therefore the model-model spread likely overstates the uncertainty in our current understanding of the ocean carbon system, particularly for transport-dominated fields such as the historical uptake of anthropogenic CO2. A full error assessment, however, would need to account for additional sources of uncertainty such as more complex biological-chemical-physical interactions, biases arising from poorly resolved or neglected physical processes, and climate change.
Anomalous chlorofluorocarbon uptake by mesoscale eddies in the Drake Passage region

(John Wiley & Sons, 2015-02-23) Song, Hajoon ; Marshall, John C. ; Gaube, Peter ; McGillicuddy, Dennis J.

The role of mesoscale eddies in the uptake of anthropogenic chlorofluorocarbon-11 (CFC-11) gas is investigated with a 1/20° eddy-resolving numerical ocean model of a region of the Southern Ocean. With a relatively fast air-sea equilibrium time scale (about a month), the air-sea CFC-11 flux quickly responds to the changes in the mixed layer CFC-11 partial pressure (pCFC-11). At the mesoscale, significant correlations are observed between pCFC-11 anomaly, anomalies in sea surface temperature (SST), net heat flux, and mixed layer depth. An eddy-centric analysis of the simulated CFC-11 field suggests that anticyclonic warm-core eddies generate negative pCFC-11 anomalies and cyclonic cold-core eddies generate positive anomalies of pCFC-11. Surface pCFC-11 is modulated by mixed layer dynamics in addition to CFC-11 air-sea fluxes. A negative cross correlation between mixed layer depth and surface pCFC-11 anomalies is linked to higher CFC-11 uptake in anticyclones and lower CFC-11 uptake in cyclones, especially in winter. An almost exact asymmetry in the air-sea CFC-11 flux between cyclones and anticyclones is found.
Ocean system science to inform the exploration of ocean worlds

(Oceanography Society, 2022-05-23) German, Christopher R. ; Blackman, Donna K. ; Fisher, Andrew T. ; Girguis, Peter R. ; Hand, Kevin P. ; Hoehler, Tori M. ; Huber, Julie A. ; Marshall, John C. ; Pietro, Kathryn R. ; Seewald, Jeffrey S. ; Shock, Everett ; Sotin, Christophe ; Thurnherr, Andreas M. ; Toner, Brandy M.

Ocean worlds provide fascinating opportunities for future ocean research. They allow us to test our understanding of processes we consider fundamental to Earth’s ocean and simultaneously provide motivation to explore our ocean further and develop new technologies to do so. In parallel, ocean worlds research offers opportunities for ocean scientists to provide meaningful contributions to novel investigations in the coming decades that will search for life beyond Earth. Key to the contributions that oceanographers can make to this field is that studies of all other ocean worlds remain extremely data limited. Here, we describe an approach based on ocean systems science in which theoretical modeling can be used, in concert with targeted laboratory experimentation and direct observations in Earth’s ocean, to predict what processes (including those essential to support life) might be occurring on other ocean worlds. In turn, such an approach would help identify new technologies that might be required for future space missions as well as appropriate analog studies that could be conducted on Earth to develop and validate such technologies. Our approach is both integrative and interdisciplinary and considers multiple domains, from processes active in the subseafloor to those associated with ocean-ice feedbacks.
Divergent behavior of hydrothermal plumes in fresh versus salty icy ocean worlds

(American Geophysical Union, 2023-11-21) Bire, Suyash ; Mittal, Tushar ; Kang, Wanying ; Ramadhan, Ali ; Tuckman, Philip J. ; German, Christopher R. ; Thurnherr, Andreas M. ; Marshall, John C.

Water parcels close to their freezing point contract and become heavy on warming if they are sufficiently fresh (salinity less than 22g kg−1 for earth's ocean), but expand and become buoyant when salty (salinity greater than 22g kg−1). We explore the resulting divergent behavior of hydrothermal plumes in fresh versus salty icy ocean worlds, with particular emphasis on Europa and Enceladus. Large, salty, putative Europa-like oceans, develop buoyant plumes which rise upwards in the water column when energized by localized hydrothermal vents. Instead, small, fresher, putative Enceladus-like oceans, can develop bottom-hugging gravity currents when heated near the freezing point, due to the anomalous contraction of fluid parcels on warming. Such a bottom-filling regime would most likely be a transient stage in the evolution of an icy moon over geological time. The contrasting dynamics are highlighted and rationalized in terms of key non-dimensional numbers with a focus on the ability of ocean to carry bio-markers from the hydrothermal activity at the bottom to the ice shell at the top. Finally, the implications of our study for prioritizing future missions to icy moons are discussed. An advantage of a mission to a large icy moon (e.g., Europa), rather than a smaller target (e.g., Enceladus), is that a larger moon's ocean would likely support buoyant convection, which could bring signatures of seafloor venting to the outer ice-shell regardless of that ocean's salinity. For smaller icy moons, the nature of convection would hinge on its assumed salinity.
“Climate response functions” for the Arctic Ocean : a proposed coordinated modelling experiment

(Copernicus Publications on behalf of the European Geosciences Union, 2017-07-21) Marshall, John C. ; Scott, Jeffery ; Proshutinsky, Andrey

A coordinated set of Arctic modelling experiments, which explore how the Arctic responds to changes in external forcing, is proposed. Our goal is to compute and compare "climate response functions" (CRFs) – the transient response of key observable indicators such as sea-ice extent, freshwater content of the Beaufort Gyre, etc. – to abrupt "step" changes in forcing fields across a number of Arctic models. Changes in wind, freshwater sources, and inflows to the Arctic basin are considered. Convolutions of known or postulated time series of these forcing fields with their respective CRFs then yield the (linear) response of these observables. This allows the project to inform, and interface directly with, Arctic observations and observers and the climate change community. Here we outline the rationale behind such experiments and illustrate our approach in the context of a coarse-resolution model of the Arctic based on the MITgcm. We conclude by summarizing the expected benefits of such an activity and encourage other modelling groups to compute CRFs with their own models so that we might begin to document their robustness to model formulation, resolution, and parameterization.
Observed mesoscale eddy signatures in Southern Ocean surface mixed-layer depth

(John Wiley & Sons, 2017-01-27) Hausmann, Ute ; McGillicuddy, Dennis J. ; Marshall, John C.

Combining satellite altimetry with Argo profile data a systematic observational estimate of mesoscale eddy signatures in surface mixed-layer depth (MLD) is provided across the Southern Ocean (SO). Eddy composite MLD anomalies are shallow in cyclones, deep in anticyclones, and increase in magnitude with eddy amplitude. Their magnitudes show a pronounced seasonal modulation roughly following the depth of the climatological mixed layer. Weak eddies of the relatively quiescent SO subtropics feature peak late-winter perturbations of ±10 m. Much larger MLD perturbations occur over the vigorous eddies originating along the Antarctic Circumpolar Current (ACC) and SO western boundary current systems, with late-winter peaks of −30 m and +60 m in the average over cyclonic and anticyclonic eddy cores (a difference of ≈ 100 m). The asymmetry between modest shallow cyclonic and pronounced deep anticyclonic anomalies is systematic and not accompanied by corresponding asymmetries in eddy amplitude. Nonetheless, the net deepening of the climatological SO mixed layer by this asymmetry in eddy MLD perturbations is estimated to be small (few meters). Eddies are shown to enhance SO MLD variability with peaks in late winter and eddy-intense regions. Anomalously deep late-winter mixed layers occur disproportionately within the cores of anticyclonic eddies, suggesting the mesoscale heightens the frequency of deep winter surface-mixing events along the eddy-intense regions of the SO. The eddy modulation in MLD reported here provides a pathway via which the oceanic mesoscale can impact air-sea fluxes of heat and carbon, the ventilation of water masses, and biological productivity across the SO.
Seasonal variation in the correlation between anomalies of sea level and chlorophyll in the Antarctic Circumpolar Current

(John Wiley & Sons, 2018-05-30) Song, Hajoon ; Long, Matthew C. ; Gaube, Peter ; Frenger, Ivy ; Marshall, John C. ; McGillicuddy, Dennis J.

The Antarctic Circumpolar Current has highly energetic mesoscale phenomena, but their impacts on phytoplankton biomass, productivity, and biogeochemical cycling are not understood well. We analyze satellite observations and an eddy‐rich ocean model to show that they drive chlorophyll anomalies of opposite sign in winter versus summer. In winter, deeper mixed layers in positive sea surface height (SSH) anomalies reduce light availability, leading to anomalously low chlorophyll concentrations. In summer with abundant light, however, positive SSH anomalies show elevated chlorophyll concentration due to higher iron level, and an iron budget analysis reveals that anomalously strong vertical mixing enhances iron supply to the mixed layer. Features with negative SSH anomalies exhibit the opposite tendencies: higher chlorophyll concentration in winter and lower in summer. Our results suggest that mesoscale modulation of iron supply, light availability, and vertical mixing plays an important role in causing systematic variations in primary productivity over the seasonal cycle.
Mesoscale modulation of air-sea CO2 flux in Drake Passage

(John Wiley & Sons, 2016-09-10) Song, Hajoon ; Marshall, John C. ; Munro, David R. ; Dutkiewicz, Stephanie ; Sweeney, Colm ; McGillicuddy, Dennis J. ; Hausmann, Ute

We investigate the role of mesoscale eddies in modulating air-sea CO2 flux and associated biogeochemical fields in Drake Passage using in situ observations and an eddy-resolving numerical model. Both observations and model show a negative correlation between temperature and partial pressure of CO2 (pCO2) anomalies at the sea surface in austral summer, indicating that warm/cold anticyclonic/cyclonic eddies take up more/less CO2. In austral winter, in contrast, relationships are reversed: warm/cold anticyclonic/cyclonic eddies are characterized by a positive/negative pCO2 anomaly and more/less CO2 outgassing. It is argued that DIC-driven effects on pCO2 are greater than temperature effects in austral summer, leading to a negative correlation. In austral winter, however, the reverse is true. An eddy-centric analysis of the model solution reveals that nitrate and iron respond differently to the same vertical mixing: vertical mixing has a greater impact on iron because its normalized vertical gradient at the base of the surface mixed layer is an order of magnitude greater than that of nitrate.
Impact of current-wind interaction on vertical processes in the Southern Ocean

(American Geophysical Union, 2020-03-17) Song, Hajoon ; Marshall, John C. ; McGillicuddy, Dennis J. ; Seo, Hyodae

Momentum input from westerly winds blowing over the Southern Ocean can be modulated by mesoscale surface currents and result in changes in large‐scale ocean circulation. Here, using an eddy‐resolving 1/20 degree ocean model configured near Drake Passage, we evaluate the impact of current‐wind interaction on vertical processes. We find a reduction in momentum input from the wind, reduced eddy kinetic energy, and a modification of Ekman pumping rates. Wind stress curl resulting from current‐wind interaction leads to net upward motion, while the nonlinear Ekman pumping term associated with horizontal gradients of relative vorticity induces net downward motion. The spatially averaged mixed layer depth estimated using a density criteria is shoaled slightly by current‐wind interaction. Current‐wind interaction, on the other hand, enhances the stratification in the thermocline below the mixed layer. Such changes have the potential to alter biogeochemical processes including nutrient supply, biological productivity, and air‐sea carbon dioxide exchange.
Mechanisms controlling the SST air-sea heat flux feedback and its dependence on spatial scale

( 2016-04-05) Hausmann, Ute ; Czaja, Arnaud ; Marshall, John C.

The turbulent air-sea heat flux feedback (α, in W m-2 K-1) is a major contributor to setting the damping timescale of sea surface temperature (SST) anomalies. In this study we compare the spatial distribution and magnitude of α in the North Atlantic and the Southern Ocean, as estimated from the ERA-Interim reanalysis dataset. The comparison is rationalized in terms of an upper bound on the heat flux feedback, associated with \fast" atmospheric export of temperature and moisture anomalies away from the marine boundary layer, and a lower bound associated with "slow" export. It is found that regions of cold surface waters (≤10°C) are best described as approaching the slow export limit. This conclusion is not only valid at the synoptic scale resolved by the reanalysis data, but also on basin scales. In particular, it applies to the heat flux feedback acting as circumpolar SST anomaly scales are approached in the Southern Ocean, with feedbacks of ≤10 W m-2 K-1. In contrast, the magnitude of the heat flux feed-back is close to that expected from the fast export limit over the Gulf Stream and its recirculation with values on the order of ≈40 W m-2 K-1. Further analysis suggests that this high value reflects a compensation between a moderate thermo-dynamic adjustment of the boundary layer, which tends to weaken the heat flux feedback, and an enhancement of the surface winds over warm SST anomalies, which tend to enhance the feedback.
Observational inferences of lateral eddy diffusivity in the halocline of the Beaufort Gyre

(John Wiley & Sons, 2017-12-26) Meneghello, Gianluca ; Marshall, John C. ; Cole, Sylvia T. ; Timmermans, Mary-Louise

Using Ekman pumping rates mediated by sea ice in the Arctic Ocean's Beaufort Gyre (BG), the magnitude of lateral eddy diffusivities required to balance downward pumping is inferred. In this limit—that of vanishing residual-mean circulation—eddy-induced upwelling exactly balances downward pumping. The implied eddy diffusivity varies spatially and decays with depth, with values of 50–400 m2/s. Eddy diffusivity estimated using mixing length theory applied to BG mooring data exhibits a similar decay with depth and range of values from 100 m2/s to more than 600 m2/s. We conclude that eddy diffusivities in the BG are likely large enough to balance downward Ekman pumping, arresting the deepening of the gyre and suggesting that eddies play a zero-order role in buoyancy and freshwater budgets of the BG.
Seasonally derived components of the Canada Basin halocline

(John Wiley & Sons, 2017-05-26) Timmermans, Mary-Louise ; Marshall, John C. ; Proshutinsky, Andrey ; Scott, Jeffery

The Arctic halocline stratification is an important barrier to the transport of deep ocean heat to the underside of sea ice. Surface water in the Chukchi Sea, warmed in summer by solar radiation, ventilates the Canada Basin halocline to create a warm layer below the mixed-layer base. The year-round persistence of this layer is shown to be consistent with the seasonal cycle of halocline ventilation. We present hydrographic observations and model results to show how Chukchi Sea density outcrops migrate seasonally as surface fluxes modify salinity and temperature. This migration is such that in winter, isopycnals bounding the warm halocline are blocked from ventilation, while the cool, relatively salty and deeper halocline layers are ventilated. In this way, the warm halocline is isolated by stratification (both vertically and laterally) each winter. Results shed light on the fate and impact to sea ice of the warm halocline under future freshening and warming of the surface Arctic Ocean.

Marshall John C.

Last Name

First Name

ORCID

Filters

Author

Subject

Date

Type

Has files

Settings

Sort By

Results per page

Search Results