GMS:Temporal Discretization

Many of the parameters associated with feature objects can be specified as either constant or transient values. Transient values are defined as a simple list of time/data pairs using the XY Series editor. The time series represents a piece-wise linear curve indicating how the parameter varies with time. When the Map -> MODFLOW command is selected, these curves must undergo temporal discretization. Temporal discretization is a process of converting general time series into discrete values that apply over specific time ranges (stress periods). Transient parameters associated with feature objects are stored in an xy series. An xy series is a general-purpose object used in GMS to represent curves of data (in this case a time series). An xy series is manipulated by GMS with regards to feature objects in three different ways: extrapolation, interpolation, and integration.

Extrapolation

Because the user is free to enter any time values for the x parameter of an xy series, it is possible that the xy series as entered does not cover the same time range as the stress periods. In this case it may be necessary to extrapolate a value for the xy series at a time before or after the first or last entered value. In GMS the simplest approach has been used. If a value is required for a time previous to the times defined by the xy series, the first value is used. Likewise for a time that is later than the all of the times in the xy series, the last value is used. Since this behavior might hide an error in the input parameters, GMS will warn the user if any xy series does not cover the time range defined by the stress periods.

Interpolation

It is also sometimes necessary to create an xy series that is a composite of two other xy series. This is the case when obtaining transient values for an intermediate point along an arc segment that has differing transient parameters at both nodes at the ends of the arc. To perform this type of interpolation, a new xy series is constructed that is the union of the x times from the two original series. The y values that correspond to each of these new x values are obtained by evaluating the original series at the x value to get two y values and then interpolating these two y values using the following equation:

Where x_n and y_n are the points along the new xy series, F₀ and F₁ are the two original xy series and t is the interpolation weighting parameter.

An example of xy series interpolation. When one node of an arc has a constant parameter and the other has a transient parameter, the constant parameter is converted into an xy series with only one point. By using the extrapolation assumption above, it is then possible to perform a transient interpolation using two transient series.

Integration

MODFLOW and MT3DMS both use the concept of Stress Periods to define the times that stresses may be applied. A stress period is a time interval during which all external stresses are constant. Because an xy series is not constrained to be constant over a time interval, it is necessary to obtain a representative value from the xy series that will approximate this condition.

GMS uses integration of the curve defined by the xy series to obtain the average value over the stress period. Recall that integrating a function from a to b yields the area beneath a curve. Then, by dividing this area by the difference between a and b, the average value of the function can be determined. By interpreting an xy series as a piece-wise linear curve, F(x), the average value for a time series over a time range from a to b is:

This average value is then assigned to the stress period.

Time Units

When entering a time series in the Map module using the XY Series Editor, you can use relative times (i.e. 0.0, 3.2, 5.4 etc.), or dates/times (1/1/2004 12:00:00 AM, 2/13/2004 2:01:00 PM etc.).

See also

• MODF/MT3D Local Source/Sinks Coverage