**Acknowledgements**

This work was supported by the German Federal Ministry of Education and Research (BMBF).

[9] ChenJ.Y., TangC.Y., ShenY.J., SakuraY., KondonA.,ShimadaJ. Use of water balance calculationand tritium to examine the dropdown of groundwatertable in the pied‐

Model Based Sustainable Management of Regional Water Supply Systems

http://dx.doi.org/10.5772/51973

47

mont of the north china plain (ncp). Environmental Geology, 44:564-571, 2003.

mental Geology,44:106-111, 2003.

wasser, 4:195-207, 2008.

SIAM Press, 2005.

gramming, 2006: 106(1):25–57.

[10] Z. Han. Groundwater resources protection andaquifer recovery in china. Environ‐

[11] Zippel M., Hannappel S. Ermittlung des Grundwasserdargebotes der Berliner Was‐ serwerke mittels regionaler numerischer Grundwasserströmungsmodelle. Grund‐

[12] FosterS., GardunoH., EvansR., OlsonD., TianY.,ZhangW., HanZ. Quarternary aquifer

[13] Antoulas A. C. Approximation of Large-Scale Dynamical Systems. Philadelphia:

[14] Wolfrum P., Vargas A., Gallivan M., Allgöwer F.: Complexity reduction of a thin film deposition model using a trajectory based nonlinear model reduction technique.

[15] Wächter A., Biegler L. T. On the implementation of a primal-dual interior point filter line search algorithm for large-scale nonlinear programming. Mathematical Pro‐

[16] Franke R., Arnold E. The solver Omuses/HQP for structured large-scale constrained optimization: algorithm, implementation and example application, Sixth SIAM Con‐

[17] Griewank A., Juedes D., Utke J. ADOL-C: A package for the automatic differentiation of algorithms written in C/C++. ACM Trans. Math. Software, 1996: 22(2): 131-167.

of thenorth china plain. Hydrogeology Journal, 12:81-93,2004.

American Control Conference, Portland, USA, 2005.

ference on Optimization, Atlanta, USA, 1999.
