SEARCH · Search NASA

Results for “Mixing layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Shallow Remineralization in the Sargasso Sea Estimated from Seasonal Variations in Oxygen and Dissolved Inorganic Carbon

A diagnostic model of the mean annual cycles of dissolved inorganic carbon (DIC) and oxygen below the mixed layer at the Bermuda Atlantic Time-series Study (BATS) site is presented and used to estimate organic carbon remineralization in the seasonal thermocline. The model includes lateral and vertical advection as well as vertical, diffusion. Very good agreement is found for the remineralization estimates based on oxygen and DIC. Net remineralization averaged from mid-spring to early fall is found to be a maximum between 120 and 140 in. Remineralization integrated between 100 (the compensation depth) and 250 m during this period is estimated to be about 1 mol C/sq m. This flux is consistent with independent estimates of the loss of particulate and dissolved organic carbon.

Ono, S.↗

NASA NTRS1999

A Comparison of Three Navier-Stokes Solvers for Exhaust Nozzle Flowfields

A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer.

Georgiadis, Nicholas J.↗

NASA NTRS1996

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗

NASA NTRS1999

Satellite Estimated Surface Heat Budgets in the Pacific Warm Pool During TOGA COARE

Chou (1998) produced daily mean surface net shortwave and longwave radiative fluxes for the westem Pacific warm pool (from the Japan GMS data) for the COARE IOP (intensive observing period). In addition, daily mean air-sea turbulent fluxes have been produced from the SSM/I data using the method of Chou (1997). In this paper, we focus on three areas. The first area is to further assess accuracy of the retrieved daily mean surface fluxes using those from the IMET buoy, RV Moana Wave, and RV IWecoma. The second area is to describe the temporal and spatial variability of the surface heat budget (as well as the related parameters) over the warm pool during the IOP. The third area is to discuss the relation of SST (sea surface temperature) change to surface fluxes, including the tranmission of solar flux through the ocean mix-layer.

Chou, Shu-Hsien↗

NASA NTRS1997

Eulerian Time-Domain Filtering for Spatial LES

Eulerian time-domain filtering seems to be appropriate for LES (large eddy simulation) of flows whose large coherent structures convect approximately at a common characteristic velocity; e.g., mixing layers, jets, and wakes. For these flows, we develop an approach to LES based on an explicit second-order digital Butterworth filter, which is applied in,the time domain in an Eulerian context. The approach is validated through a priori and a posteriori analyses of the simulated flow of a heated, subsonic, axisymmetric jet.

Pruett, C. David↗

NASA NTRS1999

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Canuto, V. M.↗

NASA NTRS1999

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗

NASA NTRS1999

Ocean Turbulence: One-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 2

We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models have adopted three adjustable "background diffusivities" for momentum, heat and salt. We propose a model that avoids such adjustable diffusivities. We assume that below the ML, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking.measured by Gargett et al. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Canuto, V. M.↗

NASA NTRS1999

Development and Validation of a Supersonic Helium-Air Coannular Jet Facility

Data are acquired in a simple coannular He/air supersonic jet suitable for validation of CFD (Computational Fluid Dynamics) codes for high speed propulsion. Helium is employed as a non-reacting hydrogen fuel simulant, constituting the core of the coannular flow while the coflow is composed of air. The mixing layer interface between the two flows in the near field and the plume region which develops further downstream constitute the primary regions of interest, similar to those present in all hypersonic air breathing propulsion systems. A computational code has been implemented from the experiment's inception, serving as a tool for model design during the development phase.

Carty, Atherton A.↗

NASA NTRS1997

Climate Variability in the Amundsen and Bellingshausen Seas

Satellite data reveal a 20% decline in sea ice extent in the Amundsen and Bellingshausen Seas in the two decades following 1973. This change is negatively correlated with surface air temperatures on the west side of the Antarctic Peninsula, which have increased -0.5 C / decade- since the mid-1940s, The recession was strongest during summer, when monthly average minima in 1991-92 removed much of the incipient multiyear ice over the continental shelf. This would have lowered the regional-mean ice thickness, impacting snow ice formation, brine production, and vertical heat flux. The northern ice edge contracted by 1 deg of latitude in all seasons from 1973-79 to 1987-93, returning toward mean conditions in 1993-95. The decline included multiyear cycles of several years in length, superimposed on high interannual variability. A review of atmospheric forcing shows winds consistent with mean and extreme ice extents, and suggests links to larger-scale circulation changes in the South Pacific. Historical ocean measurements are sparse in this sector, but mixed-layer depths and upper pycnoclines beneath the sea ice resemble those in the Weddell Sea. Weaker surface currents or changes in the upwelling of Circumpolar Deep Water on the continental shelf could have contributed to the anomaly persistence.

Jacobs, Stanley S.↗

NASA NTRS1997

El Nino Southern Oscillation and Tuna in the Western Pacific

Nearly 70% of the world's annual tuna harvest, currently 3.2 million tonnes, comes from the Pacific Ocean. Skipjack tuna (Katsuwonus pelamis) dominate the catch. Although skipjack are distributed in the surface mixed layer throughout the equatorial and subtropical Pacific, catches are highest in the western equatorial Pacific warm pool, a region characterized by low primary productivity rates that has the warmest surface waters of the world's oceans. Assessments of tuna stocks indicate that recent western Pacific skipjack catches approaching one million tonnes annually are sustainable. The warm pool, which is fundamental to the El Nino Southern Oscillation (ENSO) and the Earth's climate in general, must therefore also provide a habitat capable of supporting this highly productive tuna population. Here we show that apparent spatial shifts in the skipjack population are linked to large zonal displacements of the warm pool that occur during ENSO events. This relationship can be used to predict (several months in advance) the region of highest skipjack abundance, within a fishing ground extending over 6,000 km along the Equator.

Lehodey, P.↗

NASA NTRS1997

Solitary Waves in the Western Equatorial Pacific Ocean

During the spring tides of early January and February 1993, groups of solitary internal waves were observed propagating through the Intensive Flux Array of the TOGA COARE experiment. The waves appear to originate near the islands of Nugarba (3 deg S 30 deg S - 154 deg 30'E). They travel north-eastward at 2.5-3 m/s, closely coupled with the semi-diurnal baroclinic tide. Peak amplitudes exceed 60 m. Velocities are in excess of .8 m/s. Sea-surface vertical displacements of order.3 m can be inferred directly from the lateral acceleration of surface waters. The Equatorial Undercurrent is displaced by soliton passage but apparently is unaffected otherwise. The intrinsic shear of the solitary crests is small compared to ambient equatorial shears. The crests, while not themselves unstable, are effective at triggering instabilities on the background flow. The motions potentially contribute 10-15 Watts/sq m to the flux of heat into the mixed layer.

Pinkel, R.↗

NASA NTRS1997

Can TOPEX/Poseidon Data be Used for Monitoring the Indonesian Throughflow?

The effects of the Indonesian throughflow(ITF) are studied in a reduced gravity, primitive equation, sigma coordinate model. The surface heat fluxes are provided by coupling the ocean GCAF to an advective atmospheric mixed layer model. There is an ENSO related signal in the ITF but the correlation with the Southern Oscillation Index (SOI) is only -0.31. When the winds over the Indian Ocean are held to climatology, this correlation jumps to -0.65 indicating that the non-EASO signal in the ITF is caused by the downstream winds. On interannual time-scales the ITF can be explained in terms of sea level differences between the western Pacific and eastern Indian Oceans when appropriate representative locations are chosen as demonstrated in both model and TOPEX data. It is shown that the main climatological effect of the ITF is to warm the Indian Ocean and cool the Pacific. Mile the main ENSO indices, NIN03 and NIN04, are almost identical with and without the ITF, the total SSTs show significant ENSO dependence.

Murtugudde, Ragu↗

NASA NTRS1997

Evolution of Nonlinear Internal Waves in China Seas

Synthetic Aperture Radar (SAR) images from ERS-I have been used to study the characteristics of internal waves of Taiwan in the East China Sea, and east of Hainan Island in the South China Sea. Rank-ordered packets of internal solitons propagating shoreward from the edge of the continental shelf were observed in the SAR images. Based on the assumption of a semidiurnal tidal origin, the wave speed can be estimated and is consistent with the internal wave theory. By using the SAR images and hydrographic data, internal waves of elevation have been identified in shallow water due to a thicker mixed layer as compared with the bottom layer on the continental shelf. The generation mechanism includes the influences of the tide and the Kuroshio intrusion across the continental shelf for the formations of elevation internal waves. The effects of water depth on the evolution of solitons and wave packets are modeled by nonlinear Kortweg-deVries (KdV) type equation and linked to satellite image observations. The numerical calculations of internal wave evolution on the continental shelf have been performed and compared with the SAR observations. For a case of depression waves in deep water, the solitons first disintegrate into dispersive wave trains and then evolve to a packet of elevation waves in the shallow water area after they pass through a turning point of approximately equal layer depths has been observed in the SAR image and simulated by numerical model.

Liu, Antony K.↗

NASA NTRS1998

A System Concept for the Advanced Post-TRMM Rainfall Profiling Radars

Atmospheric latent heating field is fundamental to all modes of atmospheric circulation and upper mixed layer circulations of the ocean. The key to understanding the atmospheric heating process is understanding how and where precipitation occurs. The principal atmospheric processes which link precipitation to atmospheric circulation include: (1) convective mass fluxes in the form of updrafts and downdrafts; (2) microphysical. nucleation and growth of hydrometeors; and (3) latent heating through dynamical controls on the gravitation-driven vertical mass flux of precipitation. It is well-known that surface and near-surface rainfall are two of the key forcing functions on a number of geophysical parameters at the surface-air interface. Over ocean, rainfall variation contributes to the redistribution of water salinity, sea surface temperature, fresh water supply, and marine biology and eco-system. Over land, rainfall plays a significant role in rainforest ecology and chemistry, land hydrology and surface runoff. Precipitation has also been closely linked to a number of atmospheric anomalies and natural hazards that occur at various time scales, including hurricanes, cyclones, tropical depressions, flash floods, droughts, and most noticeable of all, the El Ninos. From this point of view, the significance of global atmospheric precipitation has gone far beyond the science arena - it has a far-reaching impact on human's socio-economic well-being and sustenance. These and many other science applications require the knowledge of, in a global basis, the vertical rain structures, including vertical motion, rain intensity, differentiation of the precipitating hydrometeors' phase state, and the classification of mesoscale physical structure of the rain systems. The only direct means to obtain such information is the use of a spaceborne profiling radar. It is important to mention that the Tropical Rainfall Measuring Mission (TRMM) have made a great stride forward towards this ultimate goal. The Precipitation Radar (PR) aboard the TRMM satellite is the first ever spaceborne radar dedicated to three-dimensional, global precipitation measurements over the tropics and the subtropics, as well as the detailed synopsis of a wide range of tropical rain storm systems. In only twelve months since launch, the PR, together with other science instruments abroad the satellite have already provided unprecedented insights into the rainfall systems. It is anticipated the a lot more exciting and important rain observations would be made by TRMM throughout its mission duration. While TRMM has provided invaluable data to the user community, it is only the first step towards advancing our knowledge on rain processes and its contributions to climate variability. It is envisioned that a TRMM follow-on mission is needed in such a way to capitalize on the pioneering information provided by TRMM, and its instrument capability must be extended beyond TRMM in such a way to fully address the key science questions from microphysical to climatic time scale. In fact, a number of new and innovative mission concepts have recently put forth for this purpose. Almost all of these new concepts have suggested the utility of a more advanced, high-resolution, Doppler-enabled, vertical profiling radar that can provide multi-parameter observations of precipitation. In this paper, a system concept for a second- gene ration precipitation radar (PR-2) which addresses the above requirements will be described.

Im, Eastwood↗

NASA NTRS1999

Ocean-Atmosphere Interaction in Climate Changes

The diagram, which attests the El Nino teleconnection observed by the NASA Scatterometer (NSCAT) in 1997, is an example of the results of our research in air-sea interaction - the core component of our three-part contribution to the Climate Variability Program. We have established an interplay among scientific research, which turns spacebased data into knowledge, a push in instrument technology, which improves observations of climate variability, and an information system, which produces and disseminates new data to support our scientific research. Timothy Liu led the proposal for advanced technology, in response to the NASA Post-2002 Request for Information. The sensor was identified as a possible mission for continuous ocean surface wind measurement at higher spatial resolution, and with the unique capability to measure ocean surface salinity. He is participating in the Instrument Incubator Program to improve the antenna technology, and is initiating a study to integrate the concept on Japanese missions. He and his collaborators have set up a system to produce and disseminate high level (gridded) ocean surface wind/stress data from NSCAT and European missions. The data system is being expanded to produce real-time gridded ocean surface winds from Quikscat, and precipitation and evaporation from the Tropical Rain Measuring Mission. It will form the basis for a spacebased data analysis system which will include momentum, heat and water fluxes. The study on 1997 El Nino teleconnection illustrates our interdisciplinary and multisensor approach to study climate variability. The diagram shows that the collapse of trade wind and the westerly wind anomalies in the central equatorial Pacific led to the equatorial ocean warming. The equatorial wind anomalies are connected to the anomalous cyclonic wind pattern in the northeast Pacific. The anomalous warming along the west coast of the United States is the result of the movement of the pre-existing warm sea surface temperature anomalies with the cyclonic wind anomalies toward the coast. The results led to a new study which identifies decadal ocean variations in the Northeast Pacific. Three studies of oceanic responses to wind forcing caused by the seasonal change of monsoons, the passage of a typhoon, and the 1997 El Nino, were successfully conducted. Besides wind forcing, we continue to examine new techniques for estimating thermal and hydrologic fluxes, through the inverse ocean mixed-layer model, through divergence of atmospheric water transport, and by direct retrieval from radiances observed by microwave radiometers. Greenhouse warming has been linked to water vapor measured by two spaceborne sensors in two studies. In the first study, strong baroclinicity and deep convection were found to transport water vapor to the upper atmosphere and increase greenhouse trapping over the storm tracks of the North Pacific and Atlantic. In another study, the annual cycle of greenhouse warming were related to sea surface temperature (SST) and integrated water vapor, and the latitudinal dependence of the magnitudes and phases of the annual cycles were compared.

Liu, W. Timothy↗

NASA NTRS2000

Significant Findings: Tracking the SeaWiFS Record with a Coupled Physical/Biogeochemical/Radiative Model of the Global Oceans

The Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) has observed 2.5 years of routine global chlorophyll observations from space. The mission was launched into a record El Nino event, which eventually gave way to one of the most intensive and longest-lasting La Nina events ever recorded. The SeaWiFS chlorophyll record captured the response of ocean phytoplankton to these significant events in the tropical Indo-Pacific basins, but also indicated significant interannual variability unrelated to the El Nino/La Nina events. This included large variability in the North Atlantic and Pacific basins, in the North Central and equatorial Atlantic, and milder patterns in the North Central Pacific. This SeaWiFS record was tracked with a coupled physical/biogeochemical/radiative model of the global oceans using near-real-time forcing data such as wind stresses, sea surface temperatures, and sea ice. This provided an opportunity to offer physically and biogeochemically meaningful explanations of the variability observed in the SeaWiFS data set, since the causal mechanisms and interrelationships of the model are completely understood. The coupled model was able to represent the seasonal distributions of chlorophyll during the SeaWiFS era, and was capable of differentiating among the widely different processes and dynamics occurring in the global oceans. The model was also reasonably successful in representing the interannual signal, especially when it was large, such as, the El Nino and La Nina events in the tropical Pacific and Indian Oceans. The model provided different phytoplankton group responses for the different events in these regions: diatoms were predominant in the tropical Pacific during the La Nina but other groups were predominant during El Nino. The opposite condition occurred in the tropical Indian Ocean. Both situations were due to the different responses of the basins to El Nino. The interannual variability in the North Atlantic, which was exhibited in SeaWiFS data as a decline in the spring/summer bloom in 1999 relative to 1998, resulted in the model from a more slowly shoaling mixed layer, allowing herbivore populations to keep pace with increasing phytoplankton populations. However, several aspects of the interannual cycle were not well-represented by the model. Explanations ranged from inherent model deficiencies, to monthly averaging of forcing fields, to biases in SeaWiFS atmospheric correction procedures.

Watson, Gregg W.↗

NASA NTRS1999

Comparison of AGE and Spectral Methods for the Simulation of Far-Wakes

Turbulent flow simulation methods based on finite differences are attractive for their simplicity, flexibility and efficiency, but not always for accuracy or stability. This report demonstrates that a good compromise is possible with the Advected Grid Explicit (AGE) method. AGE has proven to be both efficient and accurate for simulating turbulent free-shear flows, including planar mixing layers and planar jets. Its efficiency results from its localized fully explicit finite difference formulation (Bisset 1998a,b) that is very straightforward to compute, outweighing the need for a fairly small timestep. Also, most of the successful simulations were slightly under-resolved, and therefore they were, in effect, large-eddy simulations (LES) without a sub-grid-scale (SGS) model, rather than direct numerical simulations (DNS). The principle is that the role of the smallest scales of turbulent motion (when the Reynolds number is not too low) is to dissipate turbulent energy, and therefore they do not have to be simulated when the numerical method is inherently dissipative at its resolution limits. Such simulations are termed 'auto-LES' (LES with automatic SGS modeling) in this report.

Bisset, D. K.↗