Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals

Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals are investigated by considering the idealized problem of a magnetic multipole of order m and degree n located at the centre of a spherical cavity surrounded by a boundless perfect diamagnetic medium. In this illustrative idealization, the fixed spherical (magnetopause) boundary layer behaves as a perfectly conducting surface that shields the external diamagnetic medium from the compressed multipole magnetic field, which is therefore confined within the spherical cavity. For a general magnetic multipole of degree n, the non-radial components of magnetic induction just inside the magnetopause are increased by the factor 1 + [(n + 1)/n] relative to their corresponding values in the absence of the perfectly conducting spherical magnetopause. An exact equation is derived for the magnetic field lines of an individual zonal (m = 0), or axisymmetric, magnetic multipole of arbitrary degree n located at the centre of the magnetospheric cavity. For such a zonal magnetic multipole, there are always two neutral points and n – 1 neutral rings on the spherical magnetopause surface. The two neutral points are located at the poles of the spherical magnetopause. If n is even, one of the neutral rings is coincident with the equator; otherwise, the neutral rings are located symmetrically with respect to the equator. The actual existence of idealized higher-degree (n > 1) axisymmetric magnetospheres would necessarily imply multiple (n + 1) magnetospheric cusps and multiple (n) ring currents. Exact equations are also derived for the magnetic field lines of an individual non-axisymmetric magnetic multipole, confined by a perfectly conducting spherical magnetopause, in two special cases; namely, a symmetric sectorial multipole (m = n) and an antisymmetric sectorial multipole (m = n – 1). For both these non-axisymmetric magnetic multipoles, there exists on the spherical magnetopause surface a set of neutral points linked by a network of magnetic field lines. Novel magnetospheric processes are likely to arise from the existence of magnetic neutral lines that extend from the magnetopause to the surface of the Earth. Finally, magnetic field lines that are confined to, or perpendicular to, either special meridional planes or the equatorial plane, when the multipole is in free space, continue to be confined to, or perpendicular to, these same planes when the perfectly conducting magnetopause is present.Also Honorary Research Associate, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK and Visiting Reader in Physics. University of Sussex, Falmer, Brighton BN1 9QH, UK     

Response analysis of hysteretic multi-storey frames under earthquake excitation

Local incremental stiffness relations are formulated for a class of elasto-plastic beam elements. The earthquake acceleration is modelled as a filtered white noise process. The Itǒ differential equations of the integrated system made up of the structural system and the excitation process are then formulated. Instead of the original system an equivalent nonlinear system is considered, in which the drift vector is given by a series expansion of order n ≧ 1, where n = 1 represents the well-known case of equivalent linearization. Only components of the drift vector representing the non-analytieal constitutive equations are replaced by a polynomial expansion. The coefficients of this expansion are determined from a least mean square criterion, and are sequentially updated. Especially an equivalent system with a cubic expansion to the drift vector is investigated. The hierarchy of statistical moment equations is closed by a cumulant neglect closure scheme. The method has been applied to a two-storey frame. The results are compared to those of numerical simulation, and provide substantial improvements compared to equivalent linearization.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Analysis of gravity source moment inversion. Part II: Usinga priori information about the source

The method of moment inversion, based on the approximation of the gravity anomaly by thetruncated series obtained from its multipole expansion, uses, implicitly,a priori information about the anomalous body. The series truncation imposes a regularizing condition on the equipotential surfaces (produced by the anomalous body), allowing the unique determination of some moments and linear combinations of moments that are the coefficients of the basis functions in the multipole expansion series. These moments define a class of equivalent distributions of mass. The equivalence criterion is based on the misfit between the observations and the field produced by the series truncated at a prefixed maximum order for the moments. The estimates of the moments of the equivalent distribution are shown to compose the stationary solution of a system of first-order linear differential equations for which uniqueness and asymptotic stability are guaranteed. Specifically for the series retaining moments up to second order, the implicita priori information introduced requires that the source have finite volume, be sufficiently distant from the measurement plane and that its spatial distribution of mass present three orthogonal planes of symmetry intersecting at the center of mass. Subject to these hypotheses, it is possible to estimate uniquely and simultaneously the total excess of mass, the position of the center of mass and the directions of the three principal axes of the anomalous body.     

On the relations between the hydrological dynamical systems of water budget,travel time,response time and tracer concentrations

Using previous results on extended Petri Nets (EPN), we present the relations between various hydrological dynamical systems (HDSys) derived from the water budget. Once the water budget has been implemented, there is a consistent way of getting the equations for backward travel time distributions, for forward response time distributions and for the concentration of a solute or tracer. We show that the water budget has a correspondence of one to many with the backward travel time distributions. In fact, to any one of the water budget equations there correspond as many equations as there are input precipitation events. The backward travel time distributions are related to the response time distributions by Niemi's relationship and, if there are n outputs, by the definition of n − 1 partition functions. These determine what fraction of the water volume injected into the control volume at a specific time t_in goes asymptotically to a specific output. Given the backward travel time distributions, the output concentration of a solute or tracer also depends on the input concentration. The paper clarifies the complicated relations described above by taking [Hydrology and Earth System Sciences, 20, 299–328] as an example from literature. Once the appropriate information is presented, it is shown how these HDSys can be solved simultaneously without duplicating calculations. Then, it is also shown that, under the hypothesis of uniform mixing of water ages within each reservoir, these systems can be solved exactly.     

各向异性双相介质中多极源声波测井理论研究

将各向异性地层模拟为横向各向同性双相介质 ,使其对称主轴与井孔中心轴线相垂直 .本文针对井孔模型下的多极源声波测井问题提出了一种解析的摄动方法 .将各向异性双相介质看作是由对应的各向同性双相介质加上摄动部分而组成 ,并引入 5个关于弹性模量偏离各向同性时的摄动量 ,对对称主轴与井轴垂直的各向异性双相介质的多极源声波测井问题进行了系统的理论研究和分析 ,严格求解了井孔内外的零级和一级近似声波场 .发现在n级多极源激发下 ,井孔中不仅出现n级多极声场 ,而且还能激发出与各向同性地层时不相同的高于和低于n级的多极声波场     

各向异性双相介质中多极源声波测井理论研究

将各向异性地层模拟为横向各向同性双相介质 ,使其对称主轴与井孔中心轴线相垂直 .本文针对井孔模型下的多极源声波测井问题提出了一种解析的摄动方法 .将各向异性双相介质看作是由对应的各向同性双相介质加上摄动部分而组成 ,并引入 5个关于弹性模量偏离各向同性时的摄动量 ,对对称主轴与井轴垂直的各向异性双相介质的多极源声波测井问题进行了系统的理论研究和分析 ,严格求解了井孔内外的零级和一级近似声波场 .发现在n级多极源激发下 ,井孔中不仅出现n级多极声场 ,而且还能激发出与各向同性地层时不相同的高于和低于n级的多极声波场     

Coefficient estimation for a generalized equation of earthquake magnitude by means of the linear programming

Summary Beginning with a definition to the equation of the earthquake magnitude a generalization of this equation is obtained by approximating the magnitude and depth — distance functions with the help ofmth andnth degree polynomial functions. The determination of the polynomial function coefficients is proposed to be done either by the method of the least squares or by the linear programming. Owing to the solving proportions the use of electronic computers is considered as being necessary. By means of this generalized equation, magnitude equations formerly studied by different authors, are obtained through particularisation.     

Planar charged-particle trajectories in multipole magnetic fields

This paper provides a complete generalization of the classic result that the radius of curvature () of a charged-particle trajectory confined to the equatorial plane of a magnetic dipole is directly proportional to the cube of the particles equatorial distance () from the dipole (i.e. ³). Comparable results are derived for the radii of curvature of all possible planar chargedparticle trajectories in an individual static magnetic multipole of arbitrary order m and degree n. Such trajectories arise wherever there exists a plane (or planes) such that the multipole magnetic field is locally perpendicular to this plane (or planes), everywhere apart from possibly at a set of magnetic neutral lines. Therefore planar trajectories exist in the equatorial plane of an axisymmetric (m = 0), or zonal, magnetic multipole, provided n is odd: the radius of curvature varies directly as ⁿ⁼². This result reduces to the classic one in the case of a zonal magnetic dipole (n = 1). Planar trajectories exist in 2m meridional planes in the case of the general tesseral (0 < m < n) magnetic multipole. These meridional planes are defined by the 2m roots of the equation cos[m()–_n^m)] = 0, where _n^m = (1/m) arctan (h_n^m/g_n^m); g_n^m and h_n^m denote the spherical harmonic coefficients. Equatorial planar trajectories also exist if (n – m) is odd. The polar axis ( = O,) of a tesseral magnetic multipole is a magnetic neutral line if m > I. A further 2_m(n – m) neutral lines exist at the intersections of the 2_m meridional planes with the (n – m) cones defined by the (n – m) roots of the equation P_n^m(cos ) = 0 in the range 0 < 9 < , where P_n^m(cos ) denotes the associated Legendre function. If (n – m) is odd, one of these cones coincides with the equator and the magnetic field is then perpendicular to the equator everywhere apart from the 2m equatorial neutral lines. The radius of curvature of an equatorial trajectory is directly proportional to ⁿ⁼² and inversely proportional to cos[m(–)]. Since this last expression vanishes at the 2m equatorial neutral ines, the radius of curvature becomes infinitely large as the particle approaches any one of these neutral lines. The radius of curvature of a meridional trajectory is directly proportional to rⁿ⁺², where r denotes radial distance from the multiple, and inversely proportional to P_n^m(cos )/sin . Hence the radius of curvature becomes infinitely large if the particle approaches the polar magnetic neutral ine (m > 1) or any one of the 2m(n – m) neutral ines located at the intersections of the 2m meridional planes with the (n – m) cones. Illustrative particle trajectories, derived by stepwise numerical integration of the exact equations of particle motion, are pressented for low-degree (n 3) magnetic multipoles. These computed particle trajectories clearly demonstrate the non-adiabatic scattering of charged particles at magnetic neutral lines. Brief comments are made on the different regions of phase space defined by regular and irregular trajectories.Also Visiting Reader in Physics, University of Sussex, Falmer, Brighton, BN1 9QH, UK     

A mathematical structure for a theory of electromagnetic induction in the Earth in low latitude

Summary In low latitude the spatial distribution functions of the source field over the surface and the dimensions of the source, are important in any theory of electromagnetic induction developed for studying the conductivity structure of the Earth. The author has built up a mathematical structure for a theory of electromagnetic induction in anyn-layered earth model in low latitude. No simple solution is assumed for the horizontal distribution function of the source field and no assumption is made about the horizontal gradients of the source. The mathematical structure involves the concept of downward continuity of the field equations inside then-layered earth model. The resulting mathematical functions derived for anyn-layered earth model are complex. Hence a new matrix algebra of complex numbers is introduced by the author and this is built into the theory. From the upward continuity of the field equations, an inequality equation is derived in order to determine the heighth ₀ at which the induction field of the earth becomes negligible compared with the source field. The comparison of such heights at two or more stations under the same influence of the source field can be used for the resolution of the lateral distribution of the earch conductivity structure at these stations. The application of the theory will follow in a subsequent paper.     

Explorations into the formal structure of drainage basins

Every basin of higher than first order is drained by a channel network composed of two subnetworks. Their basins are separated by a drainage divide line, called the basin divider, which is the primary organizing feature of the main basin. Each basin of magnitude n contains n – 1 subnetworks of higher order, and is therefore organized by a set of n – 1 dividers. The dividers and the basin boundary are interconnected in a graph called the divider network of the basin; in graph-theoretic terms this network forms a tree and has the same magnitude and link numbers as the channel network draining the basin. While the subbasins and subnetworks of a drainage basin form a nesting hierarchy, the corresponding dividers do not; indeed, any two dividers share at most one node in common, and whether they do so is independent of whether the corresponding subbasins are nesting or disjoint. However, the dividers of nesting basins are linked by recursive relationships which permit the derivation of a set of algebraic equations; these equations relate the dividers of a basin to other basin components; for example, their combined length is equal to half the length of all first-order basin boundaries minus the length of the main basin boundary. The second part of the paper explores the dependence of the divider length on other basin parameters. The expected length, as predicted by the assumption of topological randomness, is clearly rejected by the data. An alternative approach (regression) is based on the observed magnitudes of the subbasins separated by each divider, and is reasonably successful in estimating divider length. The last section introduces the concept of the standardized basin defined by a boundary length of unity; the estimated lengths of the basin divider and the basin boundary permit an approximate reconstruction of the idealized basin shape and the location of the divider in it.     

Pressure wind adjustment relationships during a multilevel primitive equation prediction process using tropical atmospheric data

Summary In numerical weather forecasting process, with primitive equations, the wind and pressure fields mutually adjust to each other until some form of balance is achieved. The type of balance so achieved by the mass and wind fields during the numerical integration of the primitive equations governing atmospheric motions is not knowna priori. This is particularly so in the case of tropical regions where the pressure wind adjustment laws prevailing in a tropical atmosphere are not well understood.In this study we perform a systematic investigation of the pressure wind adjustment relations during a numerical integration of the primitive equations governing atmospheric motions in a tropical atmosphere. Therefore, a two-day prediction experiment is carried out using the Florida State University Tropical Prediction (FSU) Model (Krishnamurti, 1969;Krishnamurti,et al. 1973;Kanamitsu, 1975). The 200 mb predicted motion (u, v) and height (z) fields are then extracted at 0, 12, 24, 36 and 48 hours of forecast time. Using these motion (u, v) fields three other 200 mb height (z) fields were computed from the inverse nonlinear, linear and quasigeostrophic balance equations. Each of these three diagnostic heights for the 200 mb pressure surface were compared with the respective 200 mb heights obtained from the Florida State University Tropical Preciction Model. The comparison is done by computing the root-mean-square differences between the predicted 200 mb height fields and each of the three 200 mb heights obtained from the inverse non-linear, linear and quasigeostrophic balance equations. The results show that the root-meansquare differences between thez fields from the FSU model and those obtained from the non-linear and linear balance equations lie within the ranges 23 to 44 and 25 to 50 metres respectively. The root-mean-square differences between the predicted heights and the heights computed from the quasigeostrophic balance equation lie in the range 54 to 62 metres. These root-mean-square differences are of significant magnitude since large-scale disturbances in the tropical atmosphere are associated with rather small pressure changes.The variations of these root-mean-square differences as one moves from one forecast time to another exhibit no clear increasing or decreasing trend. In fact the variations appear somewhat random. This rather unsystematic time variation of the root-mean-square differences is a manifestation of the constant changes of the physics in the model as different weather systems evolve in the course of the forecasting process. It seems therefore that the pressure-wind adjustments that take place during a numerical integration of the model equations are of complex nature and cannot simply be approximated by simple diagnostic relations like the ones used in this study.Most of this work was done while the author was at the Florida State University, tallahassee, USA.     

COMPUTATION OF INTERVAL VELOCITIES FROM COMMON REFLECTION POINT MOVEOUT TIMES FOR n LAYERS WITH ARBITRARY DIPS AND CURVATURES IN THREE DIMENSIONS WHEN ASSUMING SMALL SHOT-GEOPHONE DISTANCES*

It is well known that interval velocities can be determined from common-reflection-point moveout times. However, the mathematics becomes complicated in the general case of n homogeneous layers with curved interfaces dipping in three dimensions. In this paper the problem is solved by mathematical induction using the second power terms only of the Taylor series which represents the moveout time as a function of the coordinate differences between shot and geophone points. Moreover, the zero-offset reflection times of the nth interface in a certain area surrounding the point of interest have to be known. The n—I upper interfaces and interval velocities are known too on account of the mathematical induction method applied. Thus, the zero-offset reflection raypath of the nth interface can be supposed to be known down to the intersection with the (n—1)th interface. The method applied consists mainly in transforming the second power terms of the moveout time from one interface to the next one. This is accomplished by matrix algebra. Some special cases are discussed as e.g. uniform strike and small curvatures.     

ANALYSIS OF THE COMMON MIDPOINT GATHER BY DECOMPOSITION INTO ELEMENTARY WAVEFRONTS1

The imaginary wavefront of a common midpoint gather of seismic data is analysed in terms of elementary wavefronts. The wavefronts emitted by the sources are split into a sequence of elementary circular wavefronts. Each such wavefront is independently transmitted through the subsurface and received by a group of receivers associated with the sources. The theory is first demonstrated in the case of a one-layer medium and is then extended to the case of a multilayer medium. We show that pairs of imaginary stationary points F and F' occur in the first medium through which all the rays, corresponding to the group of sources and receivers associated with the elementary wavefront pass. The location of the points is determined by 2n equations with 2n unknowns; (2n - 1) of these equations correspond to n layers of the medium, and an additional equation is given by the geometry of the CMP gather. Applying Fermat's principle between the points F and F', we obtain an algorithm for CMP ray tracing. The wavefront sampling problem is considered but not theoretically solved. Results obtained in the general case (a multilayer medium with lateral velocity variation) are presented and compared with those obtained by other modelling techniques.     

Improved dynamic correction method in seismic analysis of both classically and non‐classically damped structures

A step‐by‐step approximate procedure taking into consideration high‐frequency modes, usually neglected in the modal analysis of both classically and non‐classically damped structures, is presented. This procedure can be considered as an extension of traditional modal correction methods, like the mode‐acceleration method and the dynamic correction method, which are very effective for structural systems subjected to forcing functions described by analytical laws. The proposed procedure, herein called improved dynamic correction method, requires two steps. In the first step, the number of differential equations of motion are reduced and consequently solved by using the first few undamped mode‐shapes. In the second step, the errors due to modal truncation are reduced by correcting the dynamic response and solving a new set of differential equations, formally similar to the original differential equations of motion. The difference between the two groups of differential equations lies in the forcing vector, which is evaluated in such a way as to correct the effects of modal truncation on applied loads. Copyright © 2001 John Wiley & Sons, Ltd.     

Earthquake analysis of a class of torsionally-coupled buildings

An effective procedure is summarized without derivation for estimating the maximum response of a class of torsionally-coupled, multi-storey buildings due to earthquake ground motion characterized by response spectra. The analysis procedure is shown to be efficient and to provide insight into torsional coupling effects because the maximum response of such a torsionally-coupled, N-storey building in its njth mode of vibration is determined exactly by analysing (i) the response in the jth vibration mode (j = 1,2, N) of the corresponding torsionally-uncoupled, N-storey system; and (ii) the response in the nth vibration mode (n = 1, 2 for a one-way symmetric plan) of an associated torsionally-coupled, one-storey system. It is also demonstrated that an earlier, simpler approach in which the total (considering all vibration modes) response of the building is determined from the total response of the two simpler systems is strictly valid only under restrictive conditions, which are identified, but the simpler approach leads to results that are accurate to a useful degree.     

Discrete direction method in nonlinear problem of incoherent scattering spectrum fitting

A new fitting technique based on functional partition by simplexes of the nth order is developed to determine ionospheric plasma parameters with spectral measurement data from the Irkutsk incoherent scatter radar. In this approach, the search for a minimum of a multidimensional functional is reduced to a minimized calculation of model functions in each iteration. The algorithm also allows simplification of the search-stop criterion.     

Hydraulic model calibration for extreme floods in bedrock‐confined channels: case study from northern Thailand

Palaeoflood reconstructions based on stage evidence are typically conducted in data‐poor field settings. Few opportunities exist to calibrate the hydraulic models used to estimate discharge from this evidence. Consequently, an important hydraulic model parameter, the roughness coefficient (e.g. Manning's n), is typically estimated by a range of approximate techniques, such as ‘visual estimation’ and semi‐empirical equations. These techniques contribute uncertainty to resulting discharge estimates, especially where the study reach exhibits sensitivity in the discharge–Manning's n relation. We study this uncertainty within a hydraulic model for a large flood of known discharge on the Mae Chaem River, northern Thailand. Comparison of the ‘calibrated’ Manning's n with that obtained from semi‐empirical equations indicates that these underestimate roughness. Substantial roughness elements in the extra‐channel zone, inundated during large events, contribute significant additional sources of flow resistance that are captured neither by the semi‐empirical equations, nor by existing models predicting stage–roughness variations. This bedrock channel exhibits a complex discharge–Manning's n relation, and reliable estimates of the former are dependent upon realistic assignment of the latter. Our study demonstrates that a large recent flood can provide a valuable opportunity to constrain this parameter, and this is illustrated when we model a palaeoflood event in the same reach, and subsequently examine the magnitude–return period consequences of discharge uncertainty within a flood frequency analysis, which contributes its own source of uncertainty. Copyright © 2005 John Wiley & Sons, Ltd.     

Mass balance calculations with end member compositional variability: applications to petrologic problems

Although end member compositional variability is quite common, most quantitative mass balance procedures cannot accommodate this variability in a systematic manner. By rearranging the traditional mass balance relations, a series of equations can be derived to account for such heterogeneity. While this approach can be applied ton-component systems, the results are difficult to represent graphically. Accordingly, the procedure is most useful for two- and three-component systems. For two components, the concentration of an element in a mixture is: and the isotopic signature is: where pm, hm and 1 refer to parental magma, hybrid magma and component 1, respectively, X is the proportion of parental magma, C is concentration, is isotopic ratio, and the concentration of the denominator isotope in the particular isotopic ratio of interest. These equations describe straight lines, termed isoproportional or IP lines, of fixed mixing proportion in the C_pm−C₁ concentration or planes. Only those IP lines intersecting the rectangle defined by observed end member compositions represent viable mixing proportions. For three-component systems, the mass balance equations are: and: where Y is now the proportion of component 1 and the proportion of component 2, Z, is simply 1 − X − Y. Using major, trace and rare earth element as well as isotopic data, a region of the X-Y plane representing possible mixing combinations can be defined. Due to the compositional variability of most magmatic end members, this new mass balance procedure should be applicable to a diverse range of petrologic problems. This procedure has been applied to three different petrologic processes: Aleutian parental magma genesis, assimilation/contamination, and crust-mantle differentiation.     

On the evolution and stability of interfacial ripples on and off the critical frequency

Abstract

Under consideration are interfaces between two media of different densities and which arise from the interaction between the Mth and Nth harmonics of the motion where 1 ≤ N < M. By means of the method of multiple scales in both space and time a pair of nonlinear coupled partial differential equations is derived which model the progression of the interface. The equations contain a detuning parameter [sgrave] which allow imperfections in the resonance to be taken into account. Stokes-type sinusoidal solutions to the equations were sought. It was found that solutions exist for all values of the interaction ratio M/N. In some situations interfaces exist at both exact and near resonance; while in others they are destroyed by amplifications in the detuning. In yet others, a quantity of detuning is actually necessary for the profiles to exist. In all cases, even when the parameters are fixed, a very large class of interface profiles is possible. Finally, the stability of the profiles is studied. It is found that some are quite stable, even to perturbations with wavenumbers close to the main flow.