این تحقیق یک مدل عددی سه بعدی برای پیش بینی گردش آب های ناشی از وزش باد متغیر ارایه می دهد و نتایج حاصل از آن در مقایسه با مشاهدات، مورد بحث قرار می گیرد. مدل عددی مذکور براساس حل معادلات حاکم بر جریان و توزیع دما با روش تفاضل محدود (finite difference) با فرض DY, DX ثابت و DZ متغیر استوار است. معادلات تکانه، تاوایی، پیوستگی و روابط تنش در سطح، کف و لایه های میانی مورد استفاده قرار گرفته اند. لزجت افقی و عمودی و ضریب پخش تلاطمی براساس فرمول های استفاده شده بنت در دریاچه آنتریو (1997) برآورد شده است. با توجه به ابعاد دریای خزر در مدل عددی، از پارامتر کوریولیس (coriolis) متغیر و همچنین از تقریب هیدروستاتیک (hydrostatic)، صفحه بتا و تقریب بوسینسک (boussinesq) استفاده شده است.برای آزمایش مدل از چند مورد ساده با شرایط مرزی مشابه دریای خزر، استفاده شده است. برای مدل سازی دریای خزر، شبکه بندی حوزه در سطح افقی با ابعاد 10×10 کیلومتر و در جهت قائم در 10 لایه با ضخامت های متغیر از سطح به کف به ترتیب: 5، 10، 20، 30، 60، 100، 150، 250، 500 و بالاتر از 500 متر در نظر گرفته شده اند. داده های باد شامل سرعت، جهت و دمای آب که مربوط به روز پانزدهم شهریور سال 1374 در ساعت های 6، 12 و 18 بود که از طریق ایستگاه های همدیده ای حاشیه دریای خزر و گشت دریایی حاجی علی اف، گردآوری شده است. اطلاعات باد ساحلی اندازه گیری شده با تصحیحات اعمالی به روش (SPM, shore protection manual) به باد دور از ساحل تبدیل شد و به روش درون یابی با استفاده از عکس مربع فاصله، توزیع مکانی سرعت باد در سطح حوزه دریای خزر به دست آمد.مدل بر اساس مشاهدات و گزارش های موجود مورد ارزیابی قرار گرفته است. وضعیت جریان در چند لایه بررسی شده و سرعت ها در برش طولی و برش های عرضی در ناحیه شمالی، میانی و جنوبی مورد بحث قرار گرفته است. در کل حوزه، اندازه سرعت جریان متفاوت است. سرعت جریان محاسبه شده از 1.6cm.s-1 تا 15cm.s-1 است که بیشترین مقدار سرعت در ناحیه میانی به ازای سرعت باد حداکثر 7 متر بر ثانیه و تنش باد tx=ty=0.06N.m-2. همچنین نتایج به دست آمده، وجود سلول های گردشی را به طور پاد ساعتگرد در سه ناحیه فوق نشان می دهد و اینکه گردش با عمق به طور قابل ملاحظه ای کاهش می یابد. نتایج حاصل از حل عددی معادله دما با اطلاعات حاصل از مشاهدات گشت دریایی مقایسه شده است و تغییرات دما در لایه های متفاوت برش های عرضی و طولی، هم خوانی نسبتا قابل قبولی از مدل با مشاهدات را بیان می کند.